Next Article in Journal
Melatonin Ameliorates Thermotolerance in Soybean Seedling through Balancing Redox Homeostasis and Modulating Antioxidant Defense, Phytohormones and Polyamines Biosynthesis
Next Article in Special Issue
1H NMR Study of the HCa2Nb3O10 Photocatalyst with Different Hydration Levels
Previous Article in Journal
Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein, Mpro and PLpro of SARS-CoV-2: Protein–Protein Interactions, Dynamics Simulations and Free Energy Calculations
Previous Article in Special Issue
Glycosylation States on Intact Proteins Determined by NMR Spectroscopy
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 17
10.3390/molecules26175115
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Gadolinium-Based Paramagnetic Relaxation Enhancement Agent Enhances Sensitivity for NUS Multidimensional NMR-Based Metabolomics

by
Chandrashekhar Honrao
1,
Nathalie Teissier
1,
Bo Zhang
1,
Robert Powers
2,3,* and
Elizabeth M. O’Day
1,*
1
Olaris, Inc., Waltham, MA 02451, USA
2
Department of Chemistry, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
3
Nebraska Center for Integrated Biomolecular Communication, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(17), 5115; https://doi.org/10.3390/molecules26175115
Submission received: 22 July 2021 / Revised: 16 August 2021 / Accepted: 20 August 2021 / Published: 24 August 2021
(This article belongs to the Special Issue Practical Applications of NMR to Solve Real-World Problems)
Browse Figures
Versions Notes

Abstract

:
Gadolinium is a paramagnetic relaxation enhancement (PRE) agent that accelerates the relaxation of metabolite nuclei. In this study, we noted the ability of gadolinium to improve the sensitivity of two-dimensional, non-uniform sampled NMR spectral data collected from metabolomics samples. In time-equivalent experiments, the addition of gadolinium increased the mean signal intensity measurement and the signal-to-noise ratio for metabolite resonances in both standard and plasma samples. Gadolinium led to highly linear intensity measurements that correlated with metabolite concentrations. In the presence of gadolinium, we were able to detect a broad array of metabolites with a lower limit of detection and quantification in the low micromolar range. We also observed an increase in the repeatability of intensity measurements upon the addition of gadolinium. The results of this study suggest that the addition of a gadolinium-based PRE agent to metabolite samples can improve NMR-based metabolomics.

1. Introduction

Metabolomics is a rapidly expanding field that relies on the detection and quantification of small molecular-weight (MW < 1500 Daltons) compounds present in a biological sample. Metabolite levels are often correlated with different disease states or phenotypic outcomes, which can lead to the development of highly valuable biomarkers and provide novel insights into human health and disease [1,2,3,4,5,6,7]. Nuclear magnetic resonance (NMR) spectroscopy has proven to be a powerful tool for metabolomics that meets the analytical requirements needed to achieve a robust and accurate characterization of the metabolome [8,9,10,11]. Conventional NMR-based approaches rely on one-dimensional (1D) 1H NMR experiments, which can facilitate the absolute quantification of metabolites. However, chemical shift overlap may limit the number of metabolites that can be accurately measured, which often relies on the application of peak-fitting algorithms. The size and completeness of the reference database used by these peak fitting algorithms will also limit the number of metabolites that can be quantified. Multi-dimensional techniques such as two-dimensional (2D) 1H-13C Heteronuclear Single Quantum Correlation (HSQC) spectroscopy can increase resolution by dispersing the chemical shifts along the carbon dimension, but necessitates long acquisition times due to the low natural abundance of 13C (1.1%), thus limiting the real-world practicality of this approach [12].
Expanding upon the work of Rai [13] and Von Schlippenbach [14], we recently demonstrated that non-uniform sampling (NUS) can be used to reduce the acquisition time of a 2D 1H-13C HSQC experiment to empower semi-quantitative metabolomics [15]. Indeed, a one-hour experiment using a 25% NUS 1H-13C HSQC led to 4-fold improvement in sensitivity, which also yielded highly linear and repeatable data. Further, we established guidelines based on a signal-to-noise ratio (S/N) to enable the reliable detection of a broad range of metabolites in the low micromolar range with a coefficient of variation (CV) of less than 20%. Using these results as our baseline, we sought herein to systematically evaluate the effects of relaxation delays in combination with paramagnetic relaxation enhancement (PRE) agents to further improve the sensitivity of 2D NMR experiments for metabolomics. First, we improved the mean signal intensity and S/N of a 25% NUS 1H-13C HSQC experiment by optimizing the relaxation delay and the number of scans. Then, we observed that the addition of a gadolinium-based PRE agent further improved the S/N of the 25% NUS 1H-13C HSQC spectra for both a model mixture and plasma samples. A lower limit of detection and quantification was achieved for most metabolites, but the most dramatic improvement in signal intensity was seen for the weakest peaks. We also observed that the addition of the PRE agent maintained linearity for all metabolites over a concentration range from 50 μM to 2 mM. These intensity measurements were highly repeatable, leading to smaller CVs. Overall, our results demonstrate that PRE agents can improve the sensitivity of 2D NUS NMR spectra routinely used in metabolomic studies.

2. Results

2.1. Optimizing the Relaxation Delay for Semi-Quantitative Metabolomics

A fundamental principle of NMR spectroscopy is that increasing the number of experimental scans (N) increases the S/N ratio by a factor of √N [16]. For pulsed NMR experiments, the relaxation delay, commonly known as d1, is the time required between scans to allow spins to return to equilibrium. The optimal d1 time depends on the longitudinal relaxation (T1) rate—the time required for full restoration of the nuclear spin to equilibrium along the direction of the polarizing magnetic field [16]. Each nuclei in a molecule has a different T1 value, and for small molecules like metabolites T1 values can be several seconds long. For example, formate has a T1 > 9 s at 600 MHz [17]. For quantitative NMR, it is advised to set d1 to 5 × T1 of the slowest relaxing nuclei in a sample [17]. This would require a d1 of upwards of a minute in length, leading to impractically long acquisition times that are not feasible for high-throughput NMR metabolomics. In practice, d1 is commonly set to a pre-determined value that allows for a relative quantitative comparison between spectra collected under identical conditions. It is important to note that only a comparison between the same metabolite can be made in this manner across the spectral dataset. A comparison between two or more different metabolites would be meaningless because of the d1-dependent variation in peak integrals that distorts the relationship between peak integral and metabolite concentration.
A model mixture (“Reference 1”) was composed of 29 commonly observed human metabolites, which included amino acids, organic acids, biogenic amines, sugars, etc., from the literature [18,19,20,21] as well as metabolites commonly observed in our own clinical studies. To find the optimal d1 for a model mixture of 29 metabolites (Reference 1), we recorded time equivalent experiments (4 min ± 8 s) with varying d1 values of 1.5 s, 1.2 s, 0.8 s and 0.6 s and observed the changes in both the 1D 1H NMR spectra and 2D 25% NUS 1H-13C HSQC spectra (Figure 1). At first, as the d1 decreased, the signal intensity for the majority of the metabolites increased, which is expected due to the increased number of scans (N = 64 to 92 for the 1D- and N = 36 to 84 for the 2D-experiments). For the 25% NUS 1H-13C HSQC spectra, we observed a steady increase in the mean intensity of metabolites from 2.9 × 107 to 4.9 × 107 as d1 decreased from 1.5 s to 0.8 s. Similarly, the mean S/N increased from 98.48 to 115.34. However, the mean S/N and intensity reached a maximum at a d1 of 0.8 s. As evident by the expanded regions of the 1D 1H NMR spectra (Figure 1a), peak intensities began to decrease at a d1 of 0.6 s despite the larger number of scans. This is consistent with the 25% NUS 1H-13C HSQC spectra at a d1 of 0.6 s, where the mean S/N and intensity decreased to 106.18 and 4.7 × 107, respectively. Furthermore, significant solvent artifacts were observed in the HSQC spectra relative to longer d1 values. Presumably, at a d1 of 0.6 s, factors related to T1 dominate spectral sensitivity, which could not be negated by the allowed increase in the number of scans. This led us to select 0.8 s as the optimal d1 value for improved S/N.

2.2. Gadolinium Provides Enhanced Sensitivity

The addition of PRE agents has been previously used to accelerate NMR data acquisition [13,17,22,23,24,25]. PRE agents contain unpaired electrons and decrease T1 relaxation times for all nuclei in a sample due to dipolar interactions between nuclear and electron spin states. The PRE effect is very large, owing to the large magnetic moment of an unpaired electron, and can be tunable by adjusting the concentration of the PRE agent [17]. By combining NUS with the relaxation enhancing agent, Cu(EDTA), Rai and colleagues demonstrated a 22-fold reduction in the 2D 1H-13C HSQC data collection time to quantify a handful of urine metabolites [13]. Gadolinium-based contrast agents have been widely used in MRI diagnostic imaging, for studying soluble proteins, for characterizing protein-protein, protein-oligosaccharides, and protein-nucleic acid complexes, and for investigating membrane proteins using NMR spectroscopy [26,27]. Sakol et al. have also shown the utility of the Gd-based contrast agent, Gd-DOTA, for cellular localization studies using NMR spectroscopy [28]. Similarly, Mulder and colleagues utilized gadolinium-based PRE agents and achieved a 3- to 4-fold improvement in acquisition time for quantifying several plasma metabolites [17]. We sought to expand upon these findings by focusing on parameters to increase spectral sensitivity for a fixed-time experiment (1 h ± 4 min) instead of accelerating acquisition times.
We first assessed the 1D 1H spectral changes for Reference 1 (Supplementary Material Table S1) with a d1 of 0.8 s with increasing concentrations (0.25 mM to 1 mM) of Cu (EDTA) and Gadobutrol (Gd) (Figure S1), a gadolinium-containing macrocyclic that has previously been shown to enhance the relaxation rates of urine metabolites [17,25] (Figure S2). In general, contrast agents containing Gd shorten T1 and T2 relaxation rates through a dipole–dipole interaction between the unpaired electron of Gd and nuclei in the compound. The decrease in T1 and T2 rates depends on the contrast agent used and its concentration, the charge state of the compound, the viscosity of the solution, and the protein affinity of either the compound or contrast agent, among other issues. The typical range of T1 values for nuclei of common metabolites such as glucose, lactate, citrate, acetate, glutamine, and alanine are between 0.9 and 4 s. Similarly, T2 values range from 100 to 600 ms [17,28,29]. In the presence of Gd, T1 values can decrease from 2- to 10-fold depending on the concentration of Gd. A similar reduction is observed for T2, but is more pronounced at higher Gd concentrations. Accordingly, NMR resonances will significantly broaden into the baseline with the increase in Gd concentration [28]. Experimentally, we observed that a concentration of Gd at 0.25 mM allowed us to decrease our recycle delay to 0.8 s and achieve an overall increase in sensitivity while avoiding substantial line broadening. As the concentration of the Gd agent increased, the decrease in T2 and the associated peak broadening eventually eclipsed the reduction in T1 and negated any intensity gains from a larger number of scans [17,28,30,31]. In agreement with these observations, at 0.25 mM Gd, we noted an increase in intensities for the majority of metabolite resonances. As the concentration of Gd increased to 0.5 mM, a handful of metabolite resonances continued to show an increase in intensity, while others began to broaden. At 1 mM Gd, the majority of resonances were diminished compared to the control that lacked Gd. Interestingly, our results are in line with the theoretical optimal recycle delay predictions of Rovnyak et al. [32]. To perform the comparison, we identified NMR relaxation times reported in the literature for metabolites included in our study. For example, the work by Mulder et al. [17] demonstrated that the addition of Gd at a concentration of 0.5 mM to a mixture of small molecules (glucose, creatinine, citrate, glutamine, acetate, alanine, etc.) greatly reduced the T1 relaxation times by 2- to 10-fold, resulting in an average T1 relaxation time of ~0.6 s. Using the equation derived from Rovnyak et al., in the presence of Gd the theoretical optimal recycle delay would be ~0.8 s (1.26 × 0.6 s), which is in perfect agreement with our experimental findings of an optimal d1 of 0.8 s. In the presence of Cu (EDTA), we observed a decrease in NMR resonance intensities and significant line broadening at all concentrations tested. These results suggest that the addition of 0.25 mM of Gd may offer an optimal improvement in S/N. Indeed, when we recorded a 25% NUS 1H-13C HSQC spectrum with a d1 of 0.8 s in the presence of 0.25 mM Gd, we observed an increase in both the mean intensity and mean S/N for Reference 1 (Figure 2). While the average fold-change increase in peak intensity due to the addition of Gd was modest (1.25-fold), we observed large fold-change increases (>2-fold) for the lowest intensity resonances (Figure 2c). Thus, the addition of Gd could improve the ability to detect low abundant metabolites. Of note, significant differences were observed in the intensity for individual metabolites, suggesting that Gd affects each metabolite to a different extent. Previous studies have suggested that a charge distribution, especially anionic metabolites, may be more affected by Gd [22,24]. We also verified that, for 0.25 mM Gd, the optimal d1 remained at 0.8 s as measured by both an increase in mean peak intensity and mean S/N (Figure S3). Taken together, our results suggest that the addition of Gd can improve both S/N and peak intensities, which will result in an overall sensitivity improvement, leading to a higher accuracy and precision in the measurement of metabolite concentrations.

2.3. Gadolinium Maintains Linearity

Metabolomics requires quantification across a broad range of concentrations and the ability to accurately detect changes in metabolite levels [33,34]. We previously demonstrated that NUS 1H-13C HSQC metabolite profiling is highly linear in the 0.05 μM to 2 mM range [15]. Rai and colleagues also observed that the addition of a PRE agent, Cu(EDTA), maintained linearity for an NUS 1H-13C HSQC experiment that measured four amino acids (glycine, alanine, valine and methionine) over a concentration range of 24 to 78 mM [13]. We first sought to confirm that the addition of Gd maintained linearity over a broad concentration range. A series of six NUS 1H-13C HSQC spectra were recorded for a mixture containing 29 metabolites (Reference 2) with concentrations ranging from 50 μM to 2 mM (Table S2). For each NMR resonance, the peak intensity was plotted as a function of concentration and the data were fit to a linear regression model (Figure S4). Example plots of the four NMR resonance peaks for leucine and the single resonance peak for pyruvic acid are shown in Figure 3. More than 98% of the metabolite resonances displayed a correlation coefficient of R2 > 0.9, indicating excellent linearity (Table 1). Interestingly, glucose resonances, which can be affected by isomers and conformational changes, had an R2 > 0.99 that was an improvement from our previous findings without Gd, where we observed an R2 of ~0.8 [15]. Overall, this analysis demonstrated that peak intensities are highly linear as a function of metabolite concentration for NUS 1H-13C HSQC spectra in the presence of Gd.

2.4. Gadolinium Improves the Lower Limit of Detection and Quantification

We next sought to determine the lower limit of detection (LOD) and lower limit of quantification (LOQ) for our NUS measurements in the presence of Gd. LOD and LOQ are defined as follows:
LOD = 3 × σ
LOQ = 10 × σ
where the variance of the noise (σ) was estimated by the median absolute deviation (MAD). MAD was calculated from the COLMAR database [35], where the positive values for all non-peak data (Xi) were used in the following equations:
MAD = mediani(|Xi − mediani(Xj)|)
σ = 1.4826 × MAD
Table 2 and Table 3 list the LOD and LOQ for each of the resonances detected in Reference 2. Metabolites with multiple resonances have an LOD/LOQ for each observed peak, and thus metabolites with multiple peaks will have a range of LOD/LOQ values. The average LOD and LOQ in the presence of Gd was 7.8 ± 0.3 μM and 26 ± 1 μM, respectively. This is a dramatic improvement over our previous findings that yielded an average LOD and LOQ of 19.1 μM and 65.6 μM, respectively [15]. These prior NMR experiments lacked the addition of Gd and used a longer d1 of 1.5 s. Thus, it is possible to detect lower abundant metabolites by adding Gd and decreasing d1. We also compared the effects of different NMR probes on LOD/LOQ. For the same d1 of 1.5 s, a TCI helium-cooled probe had a lower LOD/LOQ compared to a TXI nitrogen-cooled probe (Table S3).

2.5. Gadolinium Maintains Reproducibility

Highly reproducible measurements are required to detect changes in the large number of samples associated with metabolomics studies. We previously demonstrated that intensity measurements from NUS 1H-13C HSQC experiments with a d1 of 1.5 s were highly reproducible as evident by a percent coefficient of variation (%CV) of 14 ± 9% for a model mixture containing 15 metabolites at a concentration of 500 μM [15]. By decreasing the d1 to 0.8 s, we observed a decrease in the %CV to 8 ± 8% (Figure 4) for three replicates of Reference 1. This was expected, given that the increased number of scans would lead to an increase in peak intensities. We only observed a modest decrease in %CV to 7 ± 7% (Figure 4) by adding Gd to the samples while maintaining a d1 of 0.8 s. This suggests that the addition of Gd does not negatively impact the reproducibility of NUS 1H-13C HSQC experiments and may increase the reliability of these measurements.

2.6. Gadolinium Effect on Plasma Metabolites

We next assessed the effects of Gd on our ability to detect and quantify metabolites using a commercially available standard pooled human plasma sample. We recorded a 25% NUS 1H-13C HSQC with or without the addition of Gd, and with a relaxation delay of 0.8 s, a constant scan number of 72, and an acquisition time of ~1 h. The addition of Gd led to a 1.12-fold increase in overall mean peak intensities. This increase was slightly less pronounced than the fold change of 1.25 observed with the model mixture and could be due to the presence of additional anions and salts, which are known to influence the impact of PRE agents [22,24]. Nonetheless, as noted for the model mixtures, we observed that the largest increase in fold change was associated with low-intensity resonances. These results further suggest that the addition of Gd could improve our ability to detect low abundant metabolites (Figure 5c). Furthermore, the % CV was lowered from 15% to 10% for the pooled human plasma sample in the presence of Gd (Figure 5d). Collectively, these results suggest that the addition of Gd to plasma samples increases the S/N for metabolite NMR resonances, especially for low abundant metabolites, and increases the reproducibility of intensity measurements. Overall, the addition of Gd to a metabolomics sample could facilitate an increase in the confidence and reliability in the detection and quantification of metabolite NMR resonances.

3. Discussion

Metabolites are influenced both by the genome and the environment, and thus provide the most comprehensive readout for the state of an individual [36,37,38]. By monitoring changes in metabolites, it is possible to develop novel biomarkers that reveal important health information. Indeed, altered metabolite levels have been observed in many diseases, including diabetes [39], neurodegeneration [40], cancer [4], cardiovascular disease [6], and even aging [3]. Furthermore, in a series of separate studies, we have identified metabolite biomarkers of response (BoRs) that correlate with drug responsiveness for metastatic breast cancer patients treated with CDK4/6 inhibitors as well as the anti-HER2 therapy trastuzumab; and for gastrointestinal stromal tumor (GIST) patients treated with tyrosine kinase inhibitors [41,42,43]. While additional validation studies are required, these preliminary results suggest the exciting possibility that metabolite-based biomarkers have for designing optimal treatment strategies for individual patients, which is a major goal of precision medicine.
To uncover metabolite BoRs, it is first necessary to accurately measure metabolite levels in biospecimens collected from a large number of patients so that the relative metabolite concentration can be correlated with disease outcomes and/or a drug response. NMR and mass spectrometry (MS) are the two most commonly used analytical platforms for measuring metabolites. Traditionally, MS has been favored due to its high sensitivity, dynamic range, and potential for high throughput. There are numerous sensitive LC-MS methods reported in the literature for the identification of endogenous metabolites in human plasma [19,20,21,44,45,46]. For example, amino acids are routinely detected at submicromolar concentrations (0.01 to 0.04 µM) by these targeted LC-MS methods. In contrast, NMR-based approaches typically detect plasma concentrations in the micromolar (3–10 µM) range [44,45,46]. However, MS can suffer from reproducibility issues, requires chromatography because of the narrow molecular-weight distribution of metabolites, and still faces challenges in metabolite identification [9]. Conversely, NMR is highly reproducible and can reveal structural information to facilitate metabolite identification. However, NMR is limited by sensitivity and spectral overlap [47]. Multidimensional NMR can overcome some of these challenges but requires extremely long experimental times that are not practical for the large number of samples needed for BoR discovery. Efforts to increase the throughput of NMR are actively being explored. We and others have demonstrated that NUS can accelerate NMR acquisition times to meet the high-throughput demands of metabolomics [13,14,15]. In an approximate one-hour experiment, we verified that intensity measurements from NUS 1H-13C HSQC spectra are highly reproducible and can facilitate the detection of a wide variety of metabolites in the low micromolar range.
In this study, we sought to explore additional techniques to extend the limit of metabolite detection by multidimensional NMR. As a first step, we assessed the effect of the d1 relaxation delay on S/N. The relaxation delay is the experimental time between scans in an NMR experiment to allow the nuclear spins to return to equilibrium, which is influenced by T1 longitudinal relaxation rates of each nuclei in the sample. For quantitative NMR, it is suggested to set d1 to at least 5 times the slowest T1 [17]. For metabolomics, this is not practical as T1s can be several seconds in length or longer. Instead, d1 is commonly set to a shorter, predetermined value for semi-quantification. Herein, we demonstrated that a decrease in d1 from 1.5 s to 0.8 s enabled an increase in the number of scans from 36 to 72, which led to an overall improvement in S/N and an increase in the mean signal intensity for metabolite resonances. Notably, this was accomplished without increasing the total time to acquire the NMR spectrum. Any further reduction in d1 was observed to result in severe signal artifacts from the solvent.
With the optimal d1 selected, we next sought to manipulate the T1s of metabolite nuclei through PRE. PRE accelerates spin relaxation due to induced magnetic dipolar interactions with unpaired electrons. PRE-based applications have been used for macromolecular structure determination, characterizing long-range interactions and identifying transiently populated states of proteins and complexes [28,48,49]. PRE agents also provide the foundation for contrast agents in magnetic resonance imaging (MRI) [50]. Previous metabolomics studies have suggested that the addition of PRE-agents, Cu(EDTA) or Gd, can decrease T1 relaxation times for metabolites [13,17]. We observed similar trends using our standard 25% NUS 1H-13C HSQC experimental parameters and noted that the addition of Gd led to an overall improvement in S/N and mean signal intensity for metabolites in a model mixture and from plasma samples. Although the average increased fold change in intensity with Gd was relatively modest, we did observe a significant improvement (>2-fold increase) for metabolite resonances with the lowest signal intensity. For metabolomics studies, the ability to accurately detect and quantify a broad range of metabolites that span different chemical classes and concentration ranges is paramount. Thus, an increase in the NMR signal intensity for low abundant metabolites suggests that the addition of Gd could improve the coverage of the metabolome. Indeed, both the lower limit of detection and quantification (LoD/LoQ) were significantly improved in the presence of Gd. In our previous results, the LoD and LoQ for a model mixture of metabolites was 19.1 μM and 65.6 μM, respectively. For the same model mixture, the LoD and LoQ decreased by more than 2-fold to 7.8 μM and 26 μM, respectively, by decreasing the d1 and by the addition of Gd.

4. Materials and Methods

Commercially available analytical standards were used to prepare model mixtures of metabolites, Reference 1 and Reference 2 (Tables S1 and S2): acetylcholine chloride (C7H15NO2·HCl, >99%), L-arginine (C6H14N4O2, >98%), L-glutamine (C5H10N2O3, >99%), D-alpha-hydroxyglutaric acid disodium salt (C5H6Na2O5, >98%), α-ketoglutaric acid disodium salt dihydrate (C5H4Na2O5·2H2O, >98%), adenosine 5-monophosphate disodium (C10H12N5Na2O7P, >99%), D-(-)-fructose (C6H12O6, >99%), guanosine 5-triphosphate sodium salt (C10H16N5O14P3·xNa + yH2O, >95%), lithium potassium acetyl phosphate (C2H3KLiO5P, >97%), L-ornithine hydrochloride (C5H12N2O2·HCl, >98%), β-nicotinamide adenine dinucleotide hydrate (C21H27N7O14P2·xH2O, >98%), DL-malic acid (C4H6O5, >99%), D-ribose 5-phosphate disodium salt dihydrate (C5H9Na2O8P·2H2O, >99%), sodium succinate dibasic hexahydrate (C4H4Na2O4·6H2O, >99%), sodium acetate (C2H3NaO2, >99%), sodium L-lactate (C3H5NaO3, >99), sodium citrate tribasic dihydrate (C6H5O7Na3·2H2O, >99%), sodium fumarate dibasic (C4H2Na2O4, >98%), sodium pyruvate (C3H3NaO3, >99%), uridine 5-diphosphate (C9H12N2Na2O12P2·xH2O, >96%), L-alanine (C3H7NO2, >98%), L-cysteine (C3H7NO2S, >98%), D-(+)-glucosamine hydrochloride (C6H13NO5·HCl, >99%), D-(+)-glucose (C6H12O6, >99.5%), choline chloride (C5H13NO·HCl, >99%), cytidine (C9H13N3O5, >99%), L-leucine (C6H13NO2, >98.5%), L-glutamic acid monosodium salt monohydrate (C5H8NNaO4·H2O, >99%), L-histidine (C6H11N3O3·HCl, >98.5%), L-lysine, monohydrochloride (C6H14N2O2·HCl, >98.5%). All the compounds were obtained from Sigma-Aldrich. Deuterium oxide (D2O, 99.0%) was purchased from Cambridge Isotope Laboratory, Inc., Andover, MA. Pooled human plasma (apheresisderived, K2EDTA) was purchased from innovative research, Novi, MI. Paramagnetic relaxation agents gadobutrol (Gd) (C18H31GdN4O9, >99.9%) and copper (II) disodium ethylenediaminetetraacetate tetrahydrate (Cu- EDTA) (C10H12CuN2Na2O8·4H2O) were procured from MedChemExpress, Monmouth Junction, NJ and TCI America, Portland, OR respectively. The NMR reference standard, deuterated 3-(trimethylsilyl)-1-propanesulfonic acid sodium salt (DSS-d6, 98%) was purchased from Cambridge Isotope Laboratory, Andover, MA, USA.

4.1. NMR Sample Preparation

Reference 1 and Reference 2 were prepared as previously described [15]. Human plasma extraction: metabolites were extracted from 1 mL of human plasma via a methanol and chloroform liquid–liquid extraction. The aqueous phase was transferred to a 15 mL Falcon tube and freeze-dried. The powder was reconstituted in 180 µL of 50 mM phosphate buffer at pH 7.4 in D2O, and then immediately transferred to a 3 mm NMR tube for NMR data collection. The NMR standard, DSS-d6, was added to each sample for chemical shift referencing.

4.2. NMR Experiments and Data Processing

All NMR spectra were acquired on a Bruker AVANCE III solution-state NMR spectrometer equipped with a liquid helium-cooled TCI (H/F, C, N), deuterium lock, and a cryoprobe operating at a frequency of 599.773010 MHz for proton and 150.822998 MHz for carbon. NUS schedules were generated using a Poisson gap distribution with a sinusoidal weight of two and random seed generator [51]. The same 25% NUS schedule and seed were used for all experiments. All NMR data were collected at 298 K.
The spectral widths along the direct and the indirect dimensions were set at 9578.544 and 24,132.982 Hz, respectively. The number of complex points in the direct dimension was set at 512, and in the indirect dimension set at 32 with a 25% NUS sampling schedule. The number of scans for the 1D 1H experiments was set to 64 (d1 = 1.5 s), 72 (d1 = 1.2 s), 84 (d1 = 0.8 s), and 92 (d1 = 0.6 s). The number of scans for the 2D 1H-13C HSQC experiments was set to 36 (d1 = 1.5 s), 48 (d1 = 1.2 s), 72 (d1 = 0.8 s) and 84 (d1 = 0.6 s), respectively. The scan numbers were selected such that the total acquisition time for each 1D- and 2D- experiment was on average 246 s, and 69 min., respectively. The transmitter frequency offset was set to 75 ppm in the 13C dimension and 4.7 ppm in the 1H dimension.
The spectral data were processed using the NMRPipe software package, as previously described [52]. The NUS data were reconstructed using iterative soft thresholding according to the hmsIST algorithm [51] to generate the same number of direct dimension data points and twice the number of indirect dimension data points, 512 (N2) × 256 (N1). Both the NUS and US NMR data were zero-filled, Fourier-transformed and manually phase-corrected to yield a final digital resolution of 2048 (N2) × 2048 (N1) points. Chemical shift queries, metabolite identifications and quantifications were performed using the COLMARm NMR webserver (http://spin.ccic.ohio-state.edu/index.php/colmar (accessed on 04/02/21) [35]. The metabolite list is presented in the Supplementary Materials section for Reference 1 (Table S1) and Reference 2 (Table S2). The resonance assignments were used as previously reported [15].

5. Conclusions

In this study, we demonstrated the ability of Gd to improve the sensitivity of 2D NUS NMR spectra for the analysis of metabolomics samples. The addition of Gd led to an overall improvement in S/N and mean signal intensity for metabolites in both a model mixture and plasma samples. In the model mixture, the addition of Gd led to a 1.25-fold improvement in NMR signal intensities, which resulted in 1.7- and 1.6-fold improvements in LOD and LOQ, respectively. Interestingly, a significant improvement (>2-fold increase) was observed for metabolites with the lowest peak intensities, which suggests that the combination of Gd with NUS may improve the coverage of the plasma metabolome. The addition of Gd also maintained the highly linear intensity measurements that were correlated with a wide range of metabolite concentrations (50 μM to 2 mM). The reproducibility of intensity measurements, as noted by a decrease in %CV for both the model mixture (8% to 7%) and the plasma samples (15% to 10%), was similarly improved with the addition of Gd. Collectively, our results suggest that supplementing metabolomics samples with 0.25 mM Gd can improve the sensitivity of 2D NUS 1H-13C HSQC spectra and enhance the overall quality of the resulting data analysis. The routine adoption of PRE by the metabolomics community may expand the utility of multidimensional NMR to empower future biomarker discoveries.

Supplementary Materials

The following are available online. Figure S1. Chemical structure of Gadobutrol (Gd), Figure S2. Effect of addition of relaxation agents on 1H-13C HSQC 1D 1H spectral intensity for a model mixture (Reference 1) of metabolites in time equivalent experiments, Figure S3. Plots showing (a) mean peak intensity and (b) mean S/N for Reference 1 with and without addition of Gadobutrol, Figure S4. Linear regression curve of Reference 2 metabolites, Table S1. The list of metabolites in Reference 1 model mixture, Table S2. The list of metabolites in Reference 2 model mixture, Table S3. Comparison of LOD and LOQ at 0.8 D1 and at 1.5 D1 with and without Gd.

Author Contributions

Conceptualization, R.P. and E.M.O.; investigation, C.H., N.T. and B.Z.; data curation, C.H.; writing—original draft preparation, C.H. and E.M.O.; writing—review and editing, R.P. and E.M.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Olaris, Inc., and, in part, by funding from the National Science Foundation under Grant Number (1660921) to R.P. and the National Institutes of Health Grant Number (P20 GM113126, NIGMS, Nebraska Center for Integrated Biomolecular Communication) to R.P. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Data available on request due to proprietary restrictions. The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank Leonardo O. Rodrigues, Chen Dong, Srihari Raghavendra Rao, Tan Chui-Mae, Josephine J. Wolf, Kanchan Sonkar, Monty Mahantesh Kothiwale, Charles Sang, Gerhard Wagner, Amber Balazs and David Emmons, for their many valuable discussions that motivated the writing of this communication.

Conflicts of Interest

The authors report that this work (design of the study, collection, analysis, and interpretation of data) was supported by Olaris, Inc. Chandrashekhar Honrao, Nathalie Teissier and Elizabeth O’Day are employees at Olaris, Inc. Bo Zhang is a former employee. Robert Powers is a member of the Scientific Advisory Board at Olaris, Inc.

Sample Availability

Samples are or not available from the authors.

References

  1. Riekeberg, E.; Powers, R. New frontiers in metabolomics: From measurement to insight. F1000Research 2017, 6, 1148. [Google Scholar] [CrossRef] [Green Version]
  2. Li, L.; Krznar, P.; Erban, A.; Agazzi, A.; Martin-Levilain, J.; Supale, S.; Kopka, J.; Zamboni, N.; Maechler, P. Metabolomics Identifies a Biomarker Revealing In Vivo Loss of Functional β-Cell Mass Before Diabetes Onset. Diabetes 2019, 68, 2272–2286. [Google Scholar] [CrossRef] [Green Version]
  3. Shao, Y.; Le, W. Recent advances and perspectives of metabolomics-based investigations in Parkinson’s disease. Mol. Neurodegener. 2019, 14, 1–12. [Google Scholar] [CrossRef] [Green Version]
  4. Puchades-Carrasco, L. Metabolomics Applications in Precision Medicine: An Oncological Perspective. Curr. Top. Med. Chem. 2017, 17, 2740–2751. [Google Scholar] [CrossRef] [Green Version]
  5. Cirulli, E.T.; Guo, L.; Swisher, C.L.; Shah, N.; Huang, L.; Napier, L.A.; Kirkness, E.F.; Spector, T.D.; Caskey, C.T.; Thorens, B.; et al. Profound Perturbation of the Metabolome in Obesity Is Associated with Health Risk. Cell Metab. 2019, 29, 488–500. [Google Scholar] [CrossRef] [Green Version]
  6. Ussher, J.R.; Elmariah, S.; Gerszten, R.E.; Dyck, J.R. The Emerging Role of Metabolomics in the Diagnosis and Prognosis of Cardiovascular Disease. J. Am. Coll. Cardiol. 2016, 68, 2850–2870. [Google Scholar] [CrossRef] [PubMed]
  7. Darst, B.F.; Koscik, R.L.; Hogan, K.J.; Johnson, S.C.; Engelman, C.D. Longitudinal plasma metabolomics of aging and sex. Aging 2019, 11, 1262–1282. [Google Scholar] [CrossRef] [PubMed]
  8. Bhinderwala, F.; Powers, R. NMR Metabolomics Protocols for Drug Discovery. In Cardiovascular Development; Springer Science and Business Media LLC: Berlin, Germany, 2019; Volume 2037, pp. 265–311. [Google Scholar]
  9. Markley, J.L.; Brüschweiler, R.; Edison, A.; Eghbalnia, H.R.; Powers, R.; Raftery, D.; Wishart, D.S. The future of NMR-based metabolomics. Curr. Opin. Biotechnol. 2017, 43, 34–40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. McAlpine, J.B.; Chen, S.-N.; Kutateladze, A.; MacMillan, J.B.; Appendino, G.; Barison, A.; Beniddir, M.A.; Biavatti, M.W.; Bluml, S.; Boufridi, A.; et al. The value of universally available raw NMR data for transparency, reproducibility, and integrity in natural product research. Nat. Prod. Rep. 2019, 36, 35–107. [Google Scholar] [CrossRef] [Green Version]
  11. Dona, A.; Kyriakides, M.; Scott, F.; Shephard, E.; 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] [PubMed] [Green Version]
  12. Bingol, K.; Bruschweiler-Li, L.; Li, D.-W.; Zhang, B.; Xie, M.; Brüschweiler, R. Emerging new strategies for successful metabolite identification in metabolomics. Bioanalysis 2016, 8, 557–573. [Google Scholar] [CrossRef] [Green Version]
  13. Rai, R.; Sinha, N. Fast and Accurate Quantitative Metabolic Profiling of Body Fluids by Nonlinear Sampling of1H–13C Two-Dimensional Nuclear Magnetic Resonance Spectroscopy. Anal. Chem. 2012, 84, 10005–10011. [Google Scholar] [CrossRef]
  14. Von Schlippenbach, T.; Oefner, P.J.; Gronwald, W. Systematic Evaluation of Non-Uniform Sampling Parameters in the Targeted Analysis of Urine Metabolites by 1H,1H 2D NMR Spectroscopy. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef]
  15. Zhang, B.; Powers, R.; O’Day, E.M. Evaluation of Non-Uniform Sampling 2D 1H–13C HSQC Spectra for Semi-Quantitative Metabolomics. Metabolites 2020, 10, 203. [Google Scholar] [CrossRef] [PubMed]
  16. Sørensen, O.W. James Keeler. Understanding NMR Spectroscopy. Magn. Reson. Chem. 2006, 44, 820. [Google Scholar] [CrossRef]
  17. Mulder, F.A.A.; Tenori, L.; Luchinat, C. Fast and Quantitative NMR Metabolite Analysis Afforded by a Paramagnetic Co-Solute. Angew. Chem. Int. Ed. 2019, 58, 15283–15286. [Google Scholar] [CrossRef]
  18. Cai, X.; Li, R. Concurrent profiling of polar metabolites and lipids in human plasma using HILIC-FTMS. Sci. Rep. 2016, 6, 36490. [Google Scholar] [CrossRef]
  19. Simón-Manso, Y.; Lowenthal, M.S.; Kilpatrick, L.E.; Sampson, M.L.; Telu, K.H.; Rudnick, P.A.; Mallard, W.G.; Bearden, D.W.; Schock, T.; Tchekhovskoi, D.V.; et al. Metabolite Profiling of a NIST Standard Reference Material for Human Plasma (SRM 1950): GC-MS, LC-MS, NMR, and Clinical Laboratory Analyses, Libraries, and Web-Based Resources. Anal. Chem. 2013, 85, 11725–11731. [Google Scholar] [CrossRef]
  20. Dunn, W.B.; Broadhurst, D.; Begley, P.; Zelena, E.; Francis-McIntyre, S.; Anderson, N.; Brown, M.; Knowles, J.D.; Halsall, A.; Haselden, J.N.; et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protoc. 2011, 6, 1060–1083. [Google Scholar] [CrossRef] [PubMed]
  21. Telu, K.H.; Yan, X.; Wallace, W.E.; Stein, S.E.; Simón-Manso, Y. Analysis of human plasma metabolites across different liquid chromatography/mass spectrometry platforms: Cross-platform transferable chemical signatures. Rapid Commun. Mass Spectrom. 2016, 30, 581–593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Murphy, P.S.; Leach, M.; Rowland, I.J. Signal modulation in1H magnetic resonance spectroscopy using contrast agents: Proton relaxivities of choline, creatine, andN-acetylaspartate. Magn. Reson. Med. 1999, 42, 1155–1158. [Google Scholar] [CrossRef]
  23. Valkovič, L.; Lau, J.Y.C.; Abdesselam, I.; Rider, O.J.; Frollo, I.; Tyler, D.J.; Rodgers, C.T.; Miller, J.J.J. Effects of contrast agents on relaxation properties of 31 P metabolites. Magn. Reson. Med. 2021, 85, 1805–1813. [Google Scholar] [CrossRef]
  24. Botta, M.; Aime, S.; Barge, A.; Bobba, G.; Dickins, R.S.; Parker, D.; Terreno, E. Ternary Complexes between Cationic GdIII Chelates and Anionic Metabolites in Aqueous Solution: An NMR Relaxometric Study. Chem. A Eur. J. 2003, 9, 2102–2109. [Google Scholar] [CrossRef] [PubMed]
  25. Scott, L.J. Gadobutrol: A Review in Contrast-Enhanced MRI and MRA. Clin. Drug Investig. 2018, 38, 773–784. [Google Scholar] [CrossRef] [Green Version]
  26. Clore, G.M.; Iwahara, J. Theory, Practice, and Applications of Paramagnetic Relaxation Enhancement for the Characterization of Transient Low-Population States of Biological Macromolecules and Their Complexes. Chem. Rev. 2009, 109, 4108–4139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Giannoulis, A.; Ben-Ishay, Y.; Goldfarb, D. Characteristics of Gd(III) Spin Labels for the Study of Protein Conformations; Elsevier BV: Amsterdam, The Netherlands, 2021; Volume 651, pp. 235–290. [Google Scholar]
  28. Sakol, N.; Egawa, A.; Fujiwara, T. Gadolinium Complexes as Contrast Agent for Cellular NMR Spectroscopy. Int. J. Mol. Sci. 2020, 21, 4042. [Google Scholar] [CrossRef]
  29. Tang, H.; Wang, Y.; Nicholson, J.; Lindon, J. Use of relaxation-edited one-dimensional and two dimensional nuclear magnetic resonance spectroscopy to improve detection of small metabolites in blood plasma. Anal. Biochem. 2004, 325, 260–272. [Google Scholar] [CrossRef] [PubMed]
  30. Shahbazi-Gahrouei, D.; Williams, M.; Allen, B.J. In vitro study of relationship between signal intensity and gadolinium-DTPA concentration at high magnetic field strength. Australas. Radiol. 2001, 45, 298–304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Lee, M.-J.; Kim, M.-J.; Yoon, C.-S.; Song, S.Y.; Park, K.; Kim, W.S. The T2-Shortening Effect of Gadolinium and the Optimal Conditions for Maximizing the CNR for Evaluating the Biliary System: A Phantom Study. Korean J. Radiol. 2011, 12, 358–364. [Google Scholar] [CrossRef]
  32. Rovnyak, D.; Hoch, J.; Stern, A.; Wagner, G. Resolution and sensitivity of high field nuclear magnetic resonance spectroscopy. J. Biomol. NMR 2004, 30, 1–10. [Google Scholar] [CrossRef] [PubMed]
  33. Patterson, R.E.; Ducrocq, A.J.; McDougall, D.J.; Garrett, T.J.; Yost, R.A. Comparison of blood plasma sample preparation methods for combined LC–MS lipidomics and metabolomics. J. Chromatogr. B 2015, 1002, 260–266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Krug, S.; Kastenmüller, G.; Stückler, F.; Rist, M.J.; Skurk, T.; Sailer, M.; Raffler, J.; Römisch-Margl, W.; Adamski, J.; Prehn, C.; et al. The dynamic range of the human metabolome revealed by challenges. FASEB J. 2012, 26, 2607–2619. [Google Scholar] [CrossRef] [PubMed]
  35. 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]
  36. Johnson, C.; Ivanisevic, J.; Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 2016, 17, 451–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Clish, C.B. Metabolomics: An emerging but powerful tool for precision medicine. Mol. Case Stud. 2015, 1, a000588. [Google Scholar] [CrossRef] [Green Version]
  38. Fiehn, O. Metabolomics—The link between genotypes and phenotypes. Plant Mol. Biol. 2002, 48, 155–171. [Google Scholar] [CrossRef] [PubMed]
  39. Hameed, I.; Masoodi, S.R.; Mir, S.A.; Nabi, M.; Ghazanfar, K.; A Ganai, B. Type 2 diabetes mellitus: From a metabolic disorder to an inflammatory condition. World J. Diabetes 2015, 6, 598–612. [Google Scholar] [CrossRef]
  40. Procaccini, C.; Santopaolo, M.; Faicchia, D.; Colamatteo, A.; Formisano, L.; de Candia, P.; Galgani, M.; De Rosa, V.; Matarese, G. Role of metabolism in neurodegenerative disorders. Metab. Clin. Exp. 2016, 65, 1376–1390. [Google Scholar] [CrossRef]
  41. Zhang, B.; Warner, J.; Pinto, C.; Juric, D.; ODay, E. NMR-metabolite-resonance signature to predict HR+ breast cancer patient response to CDK4/6 inhibitors. J. Clin. Oncol. 2019, 37. [Google Scholar] [CrossRef]
  42. O’Day, E.; Leitzel, K.; Ali, S.M.; Zhang, B.; Dong, C.; Gu, H.; Shi, X.; Drabick, J.J.; Cream, L.; Vasekar, M.; et al. Abstract P4-10-25: Pretreatment serum metabolome predicts PFS in first-line trastuzumab-treated metastatic breast cancer. Poster Sess. Abstr. 2020, 80, 4–10. [Google Scholar]
  43. Honrao, C.; Rao, S.R.; Teissier, N.; Call, S.G.; Oday, E.M.; Janku, F. Abstract LB031: Plasma based metabolic profiling in metastatic gastrointestinal stromal tumors (GIST). Clin. Res. 2021, 81, LB031. [Google Scholar] [CrossRef]
  44. Gong, L.-L.; Yang, S.; Zhang, W.; Han, F.-F.; Xuan, L.-L.; Lv, Y.-L.; Liu, H.; Liu, L.-H. Targeted Metabolomics for Plasma Amino Acids and Carnitines in Patients with Metabolic Syndrome Using HPLC-MS/MS. Dis. Markers 2020, 2020, 1–8. [Google Scholar] [CrossRef] [PubMed]
  45. Tomita, R.; Todoroki, K.; Maruoka, H.; Yoshida, H.; Fujioka, T.; Nakashima, M.; Yamaguchi, M.; Nohta, H. Amino Acid Metabolomics Using LC-MS/MS: Assessment of Cancer-Cell Resistance in a Simulated Tumor Microenvironment. Anal. Sci. 2016, 32, 893–900. [Google Scholar] [CrossRef] [Green Version]
  46. Ni, J.; Xu, L.; Li, W.; Wu, L. Simultaneous determination of thirteen kinds of amino acid and eight kinds of acylcarnitine in human serum by LC-MS/MS and its application to measure the serum concentration of lung cancer patients. Biomed. Chromatogr. 2016, 30, 1796–1806. [Google Scholar] [CrossRef] [PubMed]
  47. Schätzlein, M.P.; Becker, J.; Schulze-Sünninghausen, D.; Pineda-Lucena, A.; Herance, J.R.; Luy, B. Rapid two-dimensional ALSOFAST-HSQC experiment for metabolomics and fluxomics studies: Application to a 13C-enriched cancer cell model treated with gold nanoparticles. Anal. Bioanal. Chem. 2018, 410, 2793–2804. [Google Scholar] [CrossRef]
  48. Clore, G.M.; Tang, C.; Iwahara, J. Elucidating transient macromolecular interactions using paramagnetic relaxation enhancement. Curr. Opin. Struct. Biol. 2007, 17, 603–616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Clore, G.M. Practical Aspects of Paramagnetic Relaxation Enhancement in Biological Macromolecules. In Biofilms; Elsevier BV: Amsterdam, The Netherlands, 2015; Volume 564, pp. 485–497. [Google Scholar]
  50. Xiao, Y.-D.; Paudel, R.; Liu, J.; Ma, C.; Zhang, Z.-S.; Zhou, S.-K. MRI contrast agents: Classification and application (Review). Int. J. Mol. Med. 2016, 38, 1319–1326. [Google Scholar] [CrossRef] [Green Version]
  51. Hyberts, S.G.; Milbradt, A.G.; Wagner, A.B.; Arthanari, H.; Wagner, G. Application of iterative soft thresholding for fast reconstruction of NMR data non-uniformly sampled with multidimensional Poisson Gap scheduling. J. Biomol. NMR 2012, 52, 315–327. [Google Scholar] [CrossRef] [Green Version]
  52. Delaglio, F.; Grzesiek, S.; Vuister, G.W.; Zhu, G.; Pfeifer, J.; Bax, A. NMRPipe: A multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 1995, 6, 277–293. [Google Scholar] [CrossRef]
Molecules 26 05115 g001 550
Figure 1. An optimized combination of relaxation delay (d1) and number of scans (N) improved the S/N for (a) 1D 1H and (b–e) 2D 1H-13C HSQC spectra of a model mixture of metabolites in time equivalent experiments. Four 1D 1H NMR spectra are overlaid and color-coded according to the d1 value: 1.5 s (blue), 1.2 s (red), 0.8 s (green), 0.6 s (purple). The boxed regions in the 1D 1H NMR spectra are shown as expanded inserts above each corresponding arrow.
Figure 1. An optimized combination of relaxation delay (d1) and number of scans (N) improved the S/N for (a) 1D 1H and (b–e) 2D 1H-13C HSQC spectra of a model mixture of metabolites in time equivalent experiments. Four 1D 1H NMR spectra are overlaid and color-coded according to the d1 value: 1.5 s (blue), 1.2 s (red), 0.8 s (green), 0.6 s (purple). The boxed regions in the 1D 1H NMR spectra are shown as expanded inserts above each corresponding arrow.
Molecules 26 05115 g001
Molecules 26 05115 g002 550
Figure 2. Gadolinium improves S/N and mean intensity of metabolite resonances. 2D 1H-13C HSQC spectra (a) without and (b) with the addition of Gd. (c) Fold change of the median normalized peak intensity in the presence of Gd.
Figure 2. Gadolinium improves S/N and mean intensity of metabolite resonances. 2D 1H-13C HSQC spectra (a) without and (b) with the addition of Gd. (c) Fold change of the median normalized peak intensity in the presence of Gd.
Molecules 26 05115 g002
Molecules 26 05115 g003 550
Figure 3. Linear regression analysis of NMR resonance intensities for (a) leucine and (b) pyruvic acid peaks.
Figure 3. Linear regression analysis of NMR resonance intensities for (a) leucine and (b) pyruvic acid peaks.
Molecules 26 05115 g003
Molecules 26 05115 g004 550
Figure 4. Gadolinium maintains the reproducibility of NMR intensity measurements. The percent coefficient of variation (%CV) in the peak intensities measured from NUS 1H-13C HSQC spectra. The %CV decreased for metabolite resonances in the presence of Gd (orange) compared to no Gd (blue).
Figure 4. Gadolinium maintains the reproducibility of NMR intensity measurements. The percent coefficient of variation (%CV) in the peak intensities measured from NUS 1H-13C HSQC spectra. The %CV decreased for metabolite resonances in the presence of Gd (orange) compared to no Gd (blue).
Molecules 26 05115 g004
Molecules 26 05115 g005a 550Molecules 26 05115 g005b 550
Figure 5. Gadolinium improves S/N and the mean intensity of plasma metabolite resonances. 1H-13C HSQC spectra (a) without and (b) with the addition of gadolinium. (c) Fold change of median normalized intensity in the presence of gadolinium. (d) Percent coefficient of variation (% CV) measured intensity decreases for metabolite resonances in the presence of Gd (orange) compared to no Gd (blue).
Figure 5. Gadolinium improves S/N and the mean intensity of plasma metabolite resonances. 1H-13C HSQC spectra (a) without and (b) with the addition of gadolinium. (c) Fold change of median normalized intensity in the presence of gadolinium. (d) Percent coefficient of variation (% CV) measured intensity decreases for metabolite resonances in the presence of Gd (orange) compared to no Gd (blue).
Molecules 26 05115 g005aMolecules 26 05115 g005b
Table
Table 1. Reference 2 metabolite resonances and their correlation coefficient (R2) for each resonance.
Table 1. Reference 2 metabolite resonances and their correlation coefficient (R2) for each resonance.
Metabolites/R212345678910111213
NAD0.9960.9990.9990.9971.0000.9990.9980.9991.0000.9991.0000.9990.999
NADPH1.0001.0000.8860.9290.9950.9660.9670.9660.9720.9680.950
Cytidine0.9120.9990.9970.9980.9970.9980.998
UDP0.9990.9990.9990.9990.9990.9991.000
Fructose1.0001.0001.0001.0000.9840.996
AMP0.9510.9730.9700.9920.9010.994
Lysine0.9980.9990.9990.9980.942
Histidine0.9981.0000.9980.9990.998
Glucose0.9950.9970.9930.9960.974
Ribose 5P0.9830.9900.9890.9960.994
Glucosamine0.9970.9590.9920.951
2-HG0.9940.9980.9990.984
Leucine0.9980.9970.9990.997
Nicotinic acid0.9980.9960.9940.997
Acetylcholine0.9990.9890.9980.998
Glutamic acid0.9980.9980.952
Malic acid0.9930.9890.996
Arginine0.9910.9970.999
Ornithine0.9980.9961.000
Choline0.9990.9970.996
Glutamine0.9980.9980.998
GTP0.9980.9980.967
Citrate0.9900.998
Alanine0.9980.995
Lactic acid0.9990.997
Pyruvic acid0.998
Acetic acid0.990
Fumaric acid0.989
Succinic acid0.998
Table
Table 2. Limit of detection (LOD) for measured metabolite resonances.
Table 2. Limit of detection (LOD) for measured metabolite resonances.
LOD (µM) per NMR Resonance
Metabolites12345678910111213Minimal Conc. (µM)
NAD12.7 ± 0.4912.58 ± 0.2626.05 ± 0.3212.2 ± 0.0410.61 ± 0.1915.54 ± 0.977.9 ± 0.0813.48 ± 0.1516.5 ± 0.0414.41 ± 0.1112.98 ± 0.169.2 ± 0.0912.9 ± 0.027.90 ± 0.08
NADPH26.85 ± 0.2121.76 ± 0.997.46 ± 0.1719.01 ± 0.3813.4 ± 0.4810.66 ± 0.4214.4 ± 0.1425.37 ± 0.519.46 ± 0.9612.98 ± 0.2724.4 ± 0.62 7.46 ± 0.17
Cytidine18.1 ± 0.668.44 ± 0.029.04 ± 0.087.4 ± 0.006.62 ± 0.057.86 ± 0.067.89 ± 0.05 6.62 ± 0.05
UDP6.78 ± 0.099.47 ± 0.159.82 ± 0.29.85 ± 0.018.69 ± 0.1110.65 ± 0.1510.69 ± 0.05 6.78 ± 0.09
Fructose19.58 ± 0.9313.82 ± 0.4513.85 ± 0.4418.89 ± 1.1718.71 ± 0.339.49 ± 0.09 9.49 ± 0.09
AMP8.35 ± 1.6812.14 ± 1.543.94 ± 0.1613.17 ± 2.9710.75 ± 1.714.96 ± 2.96 3.94 ± 0.16
Lysine23.24 ± 0.1910.07 ± 0.249.28 ± 0.136.59 ± 0.033.6 ± 0.1 3.60 ± 0.10
Histidine19.15 ± 2.1918.79 ± 1.9110.85 ± 1.667.9 ± 0.3811.3 ± 2.09 7.90 ± 0.38
Glucose21.62 ± 3.1818.61 ± 2.0116.35 ± 1.114.32 ± 0.1115.79 ± 2.56 14.32 ± 0.11
Ribose 5P41.21 ± 2.5830.31 ± 1.8221.22 ± 0.598.98 ± 0.2619.94 ± 0.71 8.98 ± 0.26
Glucosamine30.31 ± 1.1119.2 ± 0.570.5733.83 ± 1.5122.32 ± 1.71 19.20 ± 0.57
2-HG33.83 ± 1.7633.3 ± 1.8925.95 ± 0.7512.7 ± 0.3 12.70 ± 0.3
Leucine5.64 ± 0.985.58 ± 1.0720.77 ± 1.5411.49 ± 1.85 5.58 ± 1.07
Nicotinic acid10.78 ± 0.039.71 ± 0.198.04 ± 0.167.26 ± 0.01 7.26 ± 0.01
Acetylcholine8.56 ± 0.120.71 ± 0.028.42 ± 0.18.06 ± 0.07 0.71 ± 0.02
Glutamic acid22.1± 1.18.87 ± 0.134.09 ± 0.08 4.09 ± 0.08
Malic acid30.52 ± 2.0530.12 ± 0.418.98 ± 0.26 8.98 ± 0.26
Arginine25.49 ± 2.759.06 ± 0.096.09 ± 0.06 6.09 ± 0.06
Ornithine28.27 ± 1.774.59 ± 0.056.71 ± 0.05 4.59 ± 0.05
Choline1.41 ± 0.058.09 ± 0.37.76 ± 0.29 1.41 ± 0.05
Glutamine8.32 ± 0.139.46 ± 0.225.06 ± 0.07 5.06 ± 0.07
GTP11.66 ± 0.1815.58 ± 0.517.54 ± 1.44 7.54 ± 1.44
Citrate8.2 ± 1.118.15 ± 1.21 5.14 ± 0.01
Alanine5.14 ± 0.0113.4 ± 0.08 8.15 ± 1.21
Lactic acid5.4 ± 0.0213.89 ± 0.22 5.4 ± 0.02
Pyruvic acid30.79± 2.83 30.79 ± 2.83
Acetic acid7.21 ± 0.08 7.21 ± 0.08
Fumaric acid5.6 ± 0.14 5.67 ± 0.12
Succinic acid3.21 ± 0.03 3.21 ± 0.03
Table
Table 3. Limit of quantification (LOQ) for measured metabolite resonances.
Table 3. Limit of quantification (LOQ) for measured metabolite resonances.
LOQ (μM) per NMR Resonance
Metabolites12345678910111213Minimal Conc. (µM)
NAD42.34 ± 1.6241.92 ± 0.8886.83 ± 1.0640.66 ± 0.1335.37 ± 0.6251.78 ± 3.2226.35 ± 0.2744.92 ± 0.5155.01 ± 0.1448.03 ± 0.3843.26 ± 0.5530.66 ± 0.343.01 ± 0.0626.35 ± 0.27
NADPH89.49 ± 0.772.53 ± 3.2924.87 ± 0.5563.38 ± 1.2544.66 ± 1.6135.52 ± 1.3948.01 ± 0.4784.56 ± 1.6764.88 ± 3.1943.27 ± 0.9181.33 ± 2.07 35.52 ± 1.39
Cytidine60.34 ± 2.228.12 ± 0.0830.15 ± 0.2824.68 ± 0.0122.05 ± 0.1826.2 ± 0.1926.29 ± 0.15 22.05 ± 0.18
UDP22.59 ± 0.2931.55 ±0.5132.74 ± 0.6832.85 ± 0.0328.95 ± 0.3635.49 ± 0.4935.64 ± 0.15 22.59 ± 0.29
Fructose65.27 ± 3.146.06 ±1.546.18 ± 1.4762.98 ± 3.8962.38 ± 1.0931.63 ± 0.29 31.63 ± 0.29
AMP27.83 ± 5.6140.46 ± 5.1213.13 ± 0.5243.92 ± 9.8935.84 ± 5.6849.86 ± 9.88 13.13 ± 0.52
Lysine77.47 ± 0.6433.57 ± 0.830.92 ± 0.4421.95 ± 0.0812.01 ± 0.34 12.01 ± 0.34
Histidine63.83 ± 7.2962.63 ± 6.3836.18 ± 5.5326.32 ± 1.2537.67 ± 6.97 26.32 ± 1.25
Glucose72.07 ±10.5962.02 ± 6.6954.5 ± 3.6647.74 ± 0.3652.63 ± 8.54 47.74 ± 0.36
Ribose 5P137.37 ± 8.59101.02 ± 6.0870.73 ± 1.9730.11 ± 0.8966.45 ± 2.38 30.11 ± 0.89
Glucosamine101.02 ± 3.7142.87 ± 2.83112.76 ± 5.0274.4 ± 5.69 42.87 ± 2.83
2-HG112.76 ± 5.86110.99 ± 6.386.49 ± 2.5142.34 ±1.00 42.34 ± 1.00
Leucine18.78 ± 3.2718.61 ± 3.5869.24 ± 5.1338.28 ± 6.16 18.61 ± 3.58
Nicotinic acid35.93 ± 0.132.36 ± 0.6226.8 ± 0.5424.19 ± 0.04 24.19 ± 0.04
Acetylcholine28.54 ± 0.42.38 ± 0.0828.08 + BL86 ± 0.3326.88 ± 0.24 2.38 ± 0.08
Glutamic acid73.67 ± 3.6529.55 ± 0.4413.64 ± 0.28 13.64 ± 0.28
Malic acid101.74 ± 6.84100.41 ± 1.3629.92 ± 0.88 29.92 ± 0.88
Arginine84.96 ± 9.1630.2 ± 0.320.32 ± 0.19 20.32 ± 0.19
Ornithine94.24 ± 5.9115.3 ± 0.1622.35 ± 0.18 15.30 ± 0.16
Choline4.69 ± 0.1526.98 ± 0.9925.85 ± 0.97 4.69 ± 0.15
Glutamine27.72 ± 0.4431.54 ± 0.7316.85 ± 0.24 16.85 ± 0.24
GTP38.87 ± 0.651.94 ± 1.725.13 ± 4.8 25.13 ± 4.80
Citrate27.33 ± 3.7127.18 ± 4.03 27.18 ± 4.03
Alanine17.15 ± 0.0244.67 ± 0.26 17.15 ± 0.02
Lactic acid18.08 ± 0.1946.29 ± 0.74 18.08 ± 0.19
Pyruvic acid102.65 ± 9.42 102.65 ± 9.42
Acetic acid24.61 ± 0.35 24.61 ± 0.35
Fumaric acid18.89 ± 0.39 18.89 ± 0.39
Succinic acid10.69 ± 0.11 10.69 ± 0.11
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Honrao, C.; Teissier, N.; Zhang, B.; Powers, R.; O’Day, E.M. Gadolinium-Based Paramagnetic Relaxation Enhancement Agent Enhances Sensitivity for NUS Multidimensional NMR-Based Metabolomics. Molecules 2021, 26, 5115. https://doi.org/10.3390/molecules26175115

AMA Style

Honrao C, Teissier N, Zhang B, Powers R, O’Day EM. Gadolinium-Based Paramagnetic Relaxation Enhancement Agent Enhances Sensitivity for NUS Multidimensional NMR-Based Metabolomics. Molecules. 2021; 26(17):5115. https://doi.org/10.3390/molecules26175115

Chicago/Turabian Style

Honrao, Chandrashekhar, Nathalie Teissier, Bo Zhang, Robert Powers, and Elizabeth M. O’Day. 2021. "Gadolinium-Based Paramagnetic Relaxation Enhancement Agent Enhances Sensitivity for NUS Multidimensional NMR-Based Metabolomics" Molecules 26, no. 17: 5115. https://doi.org/10.3390/molecules26175115

APA Style

Honrao, C., Teissier, N., Zhang, B., Powers, R., & O’Day, E. M. (2021). Gadolinium-Based Paramagnetic Relaxation Enhancement Agent Enhances Sensitivity for NUS Multidimensional NMR-Based Metabolomics. Molecules, 26(17), 5115. https://doi.org/10.3390/molecules26175115

Article Metrics

Back to TopTop