Application of novel solid phase extraction-NMR protocols for metabolic profiling of human urine

Metabolite identification and annotation procedures are necessary for the discovery of biomarkers indicative of phenotypes or disease states, but these processes can be bottlenecked by the sheer complexity of biofluids containing thousands of different compounds. Here we describe low-cost novel SPE-NMR protocols utilising different cartridges and conditions, on both natural and artificial urine mixtures, which produce unique retention profiles useful for metabolic profiling. We find that different SPE methods applied to biofluids such as urine can be used to selectively retain metabolites based on compound taxonomy or other key functional groups, reducing peak overlap through concentration and fractionation of unknowns and hence promising greater control over the metabolite annotation/identification process

The association of fish consumption and its urinary metabolites with cardiovascular risk factors: The International Study of Macro-/Micronutrients and Blood Pressure (INTERMAP)

Background Results from observational studies regarding associations between fish (including shellfish) intake and cardiovascular disease risk factors, including blood pressure (BP) and BMI, are inconsistent. Objective To investigate associations of fish consumption and associated urinary metabolites with BP and BMI in free-living populations. Methods We used cross-sectional data from the International Study of Macro-/Micronutrients and Blood Pressure (INTERMAP), including 4680 men and women (40–59 y) from Japan, China, the United Kingdom, and United States. Dietary intakes were assessed by four 24-h dietary recalls and BP from 8 measurements. Urinary metabolites (2 timed 24-h urinary samples) associated with fish intake acquired from NMR spectroscopy were identified. Linear models were used to estimate BP and BMI differences across categories of intake and per 2 SD higher intake of fish and its biomarkers. Results No significant associations were observed between fish intake and BP. There was a direct association with fish intake and BMI in the Japanese population sample (P trend = 0.03; fully adjusted model). In Japan, trimethylamine-N-oxide (TMAO) and taurine, respectively, demonstrated area under the receiver operating characteristic curve (AUC) values of 0.81 and 0.78 in discriminating high against low fish intake, whereas homarine (a metabolite found in shellfish muscle) demonstrated an AUC of 0.80 for high/nonshellfish intake. Direct associations were observed between urinary TMAO and BMI for all regions except Japan (P < 0.0001) and in Western populations between TMAO and BP (diastolic blood pressure: mean difference 1.28; 95% CI: 0.55, 2.02 mmHg; P = 0.0006, systolic blood pressure: mean difference 1.67; 95% CI: 0.60, 2.73 mmHg; P = 0.002). Conclusions Urinary TMAO showed a stronger association with fish intake in the Japanese compared with the Western population sample. Urinary TMAO was directly associated with BP in the Western but not the Japanese population sample. Associations between fish intake and its biomarkers and downstream associations with BP/BMI appear to be context specific. INTERMAP is registered at www.clinicaltrials.gov as NCT00005271

Finding Correspondence between Metabolomic Features in Untargeted Liquid Chromatography-Mass Spectrometry Metabolomics Datasets

Integration of multiple datasets can greatly enhance bioanalytical studies, for example, by increasing power to discover and validate biomarkers. In liquid chromatography-mass spectrometry (LC-MS) metabolomics, it is especially hard to combine untargeted datasets since the majority of metabolomic features are not annotated and thus cannot be matched by chemical identity. Typically, the information available for each feature is retention time (RT), mass-to-charge ratio (m/z), and feature intensity (FI). Pairs of features from the same metabolite in separate datasets can exhibit small but significant differences, making matching very challenging. Current methods to address this issue are too simple or rely on assumptions that cannot be met in all cases. We present a method to find feature correspondence between two similar LC-MS metabolomics experiments or batches using only the features' RT, m/z, and FI. We demonstrate the method on both real and synthetic datasets, using six orthogonal validation strategies to gauge the matching quality. In our main example, 4953 features were uniquely matched, of which 585 (96.8%) of 604 manually annotated features were correct. In a second example, 2324 features could be uniquely matched, with 79 (90.8%) out of 87 annotated features correctly matched. Most of the missed annotated matches are between features that behave very differently from modeled inter-dataset shifts of RT, MZ, and FI. In a third example with simulated data with 4755 features per dataset, 99.6% of the matches were correct. Finally, the results of matching three other dataset pairs using our method are compared with a published alternative method, metabCombiner, showing the advantages of our approach. The method can be applied using M2S (Match 2 Sets), a free, open-source MATLAB toolbox, available at https://github.com/rjdossan/M2S

Finding correspondence between metabolomic features in untargeted liquid chromatography-mass spectrometry metabolomics datasets.

Integration of multiple datasets can greatly enhance bioanalytical studies, for example, by increasing power to discover and validate biomarkers. In liquid chromatography-mass spectrometry (LC-MS) metabolomics, it is especially hard to combine untargeted datasets since the majority of metabolomic features are not annotated and thus cannot be matched by chemical identity. Typically, the information available for each feature is retention time (RT), mass-to-charge ratio (m/z), and feature intensity (FI). Pairs of features from the same metabolite in separate datasets can exhibit small but significant differences, making matching very challenging. Current methods to address this issue are too simple or rely on assumptions that cannot be met in all cases. We present a method to find feature correspondence between two similar LC-MS metabolomics experiments or batches using only the features' RT, m/z, and FI. We demonstrate the method on both real and synthetic datasets, using six orthogonal validation strategies to gauge the matching quality. In our main example, 4953 features were uniquely matched, of which 585 (96.8%) of 604 manually annotated features were correct. In a second example, 2324 features could be uniquely matched, with 79 (90.8%) out of 87 annotated features correctly matched. Most of the missed annotated matches are between features that behave very differently from modeled inter-dataset shifts of RT, MZ, and FI. In a third example with simulated data with 4755 features per dataset, 99.6% of the matches were correct. Finally, the results of matching three other dataset pairs using our method are compared with a published alternative method, metabCombiner, showing the advantages of our approach. The method can be applied using M2S (Match 2 Sets), a free, open-source MATLAB toolbox, available at https://github.com/rjdossan/M2S

Metabolome-wide association study on ABCA7 indicates a role of ceramide metabolism in Alzheimer’s disease

Genome-wide association studies (GWASs) have identified genetic loci associated with the risk of Alzheimer’s disease (AD), but the molecular mechanisms by which they confer risk are largely unknown. We conducted a metabolome-wide association study (MWAS) of AD-associated loci from GWASs using untargeted metabolic profiling (metabolomics) by ultraperformance liquid chromatography–mass spectrometry (UPLC-MS). We identified an association of lactosylceramides (LacCer) with AD-related single-nucleotide polymorphisms (SNPs) in ABCA7 (P = 5.0 × 10−5 to 1.3 × 10−44). We showed that plasma LacCer concentrations are associated with cognitive performance and genetically modified levels of LacCer are associated with AD risk. We then showed that concentrations of sphingomyelins, ceramides, and hexosylceramides were altered in brain tissue from Abca7 knockout mice, compared with wild type (WT) (P = 0.049–1.4 × 10−5), but not in a mouse model of amyloidosis. Furthermore, activation of microglia increases intracellular concentrations of hexosylceramides in part through induction in the expression of sphingosine kinase, an enzyme with a high control coefficient for sphingolipid and ceramide synthesis. Our work suggests that the risk for AD arising from functional variations in ABCA7 is mediated at least in part through ceramides. Modulation of their metabolism or downstream signaling may offer new therapeutic opportunities for AD

Metabolome-wide association study on ABCA7 indicates a role of ceramide metabolism in Alzheimer's disease.

Genome-wide association studies (GWASs) have identified genetic loci associated with the risk of Alzheimer's disease (AD), but the molecular mechanisms by which they confer risk are largely unknown. We conducted a metabolome-wide association study (MWAS) of AD-associated loci from GWASs using untargeted metabolic profiling (metabolomics) by ultraperformance liquid chromatography-mass spectrometry (UPLC-MS). We identified an association of lactosylceramides (LacCer) with AD-related single-nucleotide polymorphisms (SNPs) in ABCA7 (P = 5.0 × 10-5 to 1.3 × 10-44). We showed that plasma LacCer concentrations are associated with cognitive performance and genetically modified levels of LacCer are associated with AD risk. We then showed that concentrations of sphingomyelins, ceramides, and hexosylceramides were altered in brain tissue from Abca7 knockout mice, compared with wild type (WT) (P = 0.049-1.4 × 10-5), but not in a mouse model of amyloidosis. Furthermore, activation of microglia increases intracellular concentrations of hexosylceramides in part through induction in the expression of sphingosine kinase, an enzyme with a high control coefficient for sphingolipid and ceramide synthesis. Our work suggests that the risk for AD arising from functional variations in ABCA7 is mediated at least in part through ceramides. Modulation of their metabolism or downstream signaling may offer new therapeutic opportunities for AD

Mapping the internal recognition surface of an octanuclear coordination cage using guest libraries

Size and shape criteria for guest binding inside the cavity of an octanuclear cubic coordination cage in water have been established using a new fluorescence displacement assay to quantify guest binding. For aliphatic cyclic ketones of increasing size (from C5 to C11), there is a linear relationship between ΔG for guest binding and the guest’s surface area: the change in ΔG for binding is 0.3 kJ mol–1 Å–2, corresponding to 5 kJ mol–1 for each additional CH2 group in the guest, in good agreement with expectations based on hydrophobic desolvation. The highest association constant is K = 1.2 × 106 M–1 for cycloundecanone, whose volume is approximately 50% of the cavity volume; for larger C12 and C13 cyclic ketones, the association constant progressively decreases as the guests become too large. For a series of C10 aliphatic ketones differing in shape but not size, ΔG for guest binding showed no correlation with surface area. These guests are close to the volume limit of the cavity (cf. Rebek’s 55% rule), so the association constant is sensitive to shape complementarity, with small changes in guest structure resulting in large changes in binding affinity. The most flexible members of this series (linear aliphatic ketones) did not bind, whereas the more preorganized cyclic ketones all have association constants of 104–105 M–1. A crystal structure of the cage·cycloundecanone complex shows that the guest carbonyl oxygen is directed into a binding pocket defined by a convergent set of CH groups, which act as weak hydrogen-bond donors, and also shows close contacts between the exterior surface of the disc-shaped guest and the interior surface of the pseudospherical cage cavity despite the slight mismatch in shape

SMolESY: An efficient and quantitative alternative to on-instrument macromolecular ¹H-NMR signal suppression

One-dimensional (1D) proton-nuclear magnetic resonance (1H-NMR) spectroscopy is an established technique for measuring small molecules in a wide variety of complex biological sample types. It is demonstrably reproducible, easily automatable and consequently ideal for routine and large-scale application. However, samples containing proteins, lipids, polysaccharides and other macromolecules produce broad signals which overlap and convolute those from small molecules. NMR experiment types designed to suppress macromolecular signals during acquisition may be additionally performed, however these approaches add to the overall sample analysis time and cost, especially for large cohort studies, and fail to produce reliably quantitative data. Here, we propose an alternative way of computationally eliminating macromolecular signals, employing the mathematical differentiation of standard 1H-NMR spectra, producing small molecule-enhanced spectra with preserved quantitative capability and increased resolution. Our approach, presented in its simplest form, was implemented in a cheminformatic toolbox and successfully applied to more than 3000 samples of various biological matrices rich or potentially rich with macromolecules, offering an efficient alternative to on-instrument experimentation, facilitating NMR use in routine and large-scale applications