Quantitative determination of metabolites in human urine by <Superscript>1</Superscript>H NMR spectroscopy

Conditions for registration of urinary ¹H NMR spectra have been optimized in order to achieve maximal accuracy of quantitative analysis. Urinary samples from patients with acute pancreatitis have been investigated and spectral data of identified urinary metabolites and results of their quantitative determination are given. Employment of ¹H NMR spectra is perspective for the development of new laboratory diagnostic methods.     

Use of quantitative <Superscript>1</Superscript>H NMR chemical shift changes for ligand docking into barnase

¹H NMR complexation-induced changes in chemical shift (CIS) of HN protons have been used to characterize the complexes of barnase with the deoxyoligonucleotides d(GC) and d(CGAC). Quantitative shift changes are used not only to locate the most probable binding site (using ring-current shifts), but also to determine the orientation of the ligand within the binding site, based on a more complete shift calculation including bond magnetic anisotropies and electric field effects. For both ligands, the guanine is in the same binding site cleft, in the same position as identified in the crystal structure of the d(CGAC) complex. By contrast, a previous X-ray crystal structure of the d(GC) complex showed the ligand in the mouth of the active site, rather than at the guanyl-specific site, implying that the location may be an artifact of the crystallisation process.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N backbone NMR chemical shift assignments of the C-terminal P4 domain of Ahnak

Ahnak is a ~?700 kDa polypeptide that was originally identified as a tumour-related nuclear phosphoprotein, but later recognized to play a variety of diverse physiological roles related to cell architecture and migration. A critical function of Ahnak is modulation of Ca²⁺ signaling in cardiomyocytes by interacting with the β subunit of the L-type Ca²⁺ channel (Ca_V1.2). Previous studies have identified the C-terminal region of Ahnak, designated as P3 and P4 domains, as a key mediator of its functional activity. We report here the nearly complete ¹H, ¹³C and ¹⁵N backbone NMR chemical shift assignments of the 11 kDa C-terminal P4 domain of Ahnak. This study lays the foundations for future investigations of functional dynamics, structure determination and interaction site mapping of the Ca_V1.2-Ahnak complex.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N,and <Superscript>13</Superscript>C chemical shift assignments of the regulatory domain of human calcineurin

Calcineurin (CaN) plays an important role in T-cell activation, cardiac system development and nervous system function. Previous studies have demonstrated that the regulatory domain (RD) of CaN binds calmodulin (CaM) towards the N-terminal end. Calcium-loaded CaM activates the serine/threonine phosphatase activity of CaN by binding to the RD, although the mechanistic details of this interaction remain unclear. It is thought that CaM binding at the RD displaces the auto-inhibitory domain (AID) from the active site of CaN, activating phosphatase activity. In the absence of calcium-loaded CaM, the RD is disordered, and binding of CaM induces folding in the RD. In order to provide mechanistic detail about the CaM–CaN interaction, we have undertaken an NMR study of the RD of CaN. Complete ¹³C, ¹⁵N and ¹H assignments of the RD of CaN were obtained using solution NMR spectroscopy. The backbone of RD has been assigned using a combination of ¹³C-detected CON-IPAP experiments as well as traditional HNCO, HNCA, HNCOCA and HNCACB-based 3D NMR spectroscopy. A ¹⁵N-resolved TOCSY experiment has been used to assign Hα and Hβ chemical shifts.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignments of Saccharomyces cerevisiae type 1 thioredoxin in the oxidized state by solution NMR spectroscopy

Thioredoxins (Trx) are ubiquitous proteins that regulate several biochemical processes inside the cell. Trx is an important player, displaying oxidoreductase activity and helping to keep and regulate the oxidative state of the cellular environment. Trx also participates in the regulation of many cellular functions, such as DNA synthesis, protection against oxidative stress, cell cycle and signal transduction. The oxidized Trx is the target for another set of proteins, such as thioredoxin reductase (TrR), which used the reductive potential of NADPH. The oxidized state of Trx also plays important role in regulation of redox state in the cells. In this regard, the oxidized form of Trx is a putative conformer that contributes to the cellular redox environment. Here we report the chemical shift assignments (¹H, ¹³C and ¹⁵N) in solution at 15 °C. We also showed the secondary structure analysis of the oxidized form of yeast thioredoxin (yTrx1) as basis for future NMR studies of protein–target interactions and dynamics. The assignment was done at low concentration (200 µM) because it is important to keep intact the water cavity.     

Backbone <Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignment of full-length human uracil DNA glycosylase UNG2

Human uracil N-glycosylase isoform 2—UNG2 consists of an N-terminal intrinsically disordered regulatory domain (UNG2 residues 1–92, 9.3 kDa) and a C-terminal structured catalytic domain (UNG2 residues 93–313, 25.1 kDa). Here, we report the backbone ¹H, ¹³C, and ¹⁵N chemical shift assignment as well as secondary structure analysis of the N-and C-terminal domains of UNG2 representing the full-length UNG2 protein.     

A comparison of human serum and plasma metabolites using untargeted <Superscript>1</Superscript>H NMR spectroscopy and UPLC-MS

Introduction

Differences in the metabolite profiles between serum and plasma are incompletely understood.

Objectives

To evaluate metabolic profile differences between serum and plasma and among plasma sample subtypes.

Methods

We analyzed serum, platelet rich plasma (PRP), platelet poor plasma (PPP), and platelet free plasma (PFP), collected from 8 non-fasting apparently healthy women, using untargeted standard 1D and CPMG ¹H NMR and reverse phase and hydrophilic (HILIC) UPLC-MS. Differences between metabolic profiles were evaluated using validated principal component and orthogonal partial least squares discriminant analysis.

Results

Explorative analysis showed the main source of variation among samples was due to inter-individual differences with no grouping by sample type. After correcting for inter-individual differences, lipoproteins, lipids in VLDL/LDL, lactate, glutamine, and glucose were found to discriminate serum from plasma in NMR analyses. In UPLC-MS analyses, lysophosphatidylethanolamine (lysoPE)(18:0) and lysophosphatidic acid(20:0) were higher in serum, and phosphatidylcholines (PC)(16:1/18:2, 20:3/18:0, O-20:0/22:4), lysoPC(16:0), PE(O-18:2/20:4), sphingomyelin(18:0/22:0), and linoleic acid were lower. In plasma subtype analyses, isoleucine, leucine, valine, phenylalanine, glutamate, and pyruvate were higher among PRP samples compared with PPP and PFP by NMR while lipids in VLDL/LDL, citrate, and glutamine were lower. By UPLC-MS, PE(18:0/18:2) and PC(P-16:0/20:4) were higher in PRP compared with PFP samples.

Conclusions

Correction for inter-individual variation was required to detect metabolite differences between serum and plasma. Our results suggest the potential importance of inter-individual effects and sample type on the results from serum and plasma metabolic phenotyping studies.

     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N backbone assignment of the human heat-labile enterotoxin B-pentamer and chemical shift mapping of neolactotetraose binding

The major virulence factor of enterotoxigenic Escherichia coli is the heat-labile enterotoxin (LT), an AB₅ toxin closely related to the cholera toxin. LT consists of six subunits, the catalytically active A-subunit and five B-subunits arranged as a pentameric ring (LTB), which enable the toxin to bind to the epithelial cells in the intestinal lumen. LTB has two recognized binding sites; the primary binding site is responsible for anchoring the toxin to its main receptor, the GM₁-ganglioside, while the secondary binding site recognizes blood group antigens. Herein, we report the ¹H, ¹³C, ¹⁵N main chain assignment of LTB from human isolates (hLTB; 103 a.a. per subunit, with a total molecular mass of 58.5 kDa). The secondary structure was predicted based on ¹³C′, ¹³C^{α, 13}C^{β, 1}H^N and ¹⁵N chemical shifts and compared to a published crystal structure of LTB. Neolactotetraose (NEO) was titrated to hLTB and chemical shift perturbations were measured. The chemical shift perturbations were mapped onto the crystal structure, confirming that NEO binds to the primary binding site of hLTB and competes with GM₁-binding. Our new data further lend support to the hypothesis that binding at the primary binding site is transmitted to the secondary binding site of the toxin, where it may influence the binding to blood group antigens.     

<Superscript>31</Superscript>P NMR correlation maps of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O chemical shift isotopic effects for phosphometabolite labeling studies

Intramolecular correlations among the ¹⁸O-labels of metabolic oligophosphates, mapped by J-decoupled ³¹P NMR 2D chemical shift correlation spectroscopy, impart stringent constraints to the ¹⁸O-isotope distributions over the whole oligophosphate moiety. The multiple deduced correlations of isotopic labels enable determination of site-specific fractional isotope enrichments and unravel the isotopologue statistics. This approach ensures accurate determination of ¹⁸O-labeling rates of phosphometabolites, critical in biochemical energy conversion and metabolic flux transmission. The biological usefulness of the J-decoupled ³¹P NMR 2D chemical shift correlation maps was validated on adenosine tri-phosphate fractionally ¹⁸O labeled in perfused mammalian hearts.     

Automated amino acid side-chain NMR assignment of proteins using <Superscript>13</Superscript>C- and <Superscript>15</Superscript>N-resolved 3D [<Superscript>1</Superscript>H,<Superscript>1</Superscript>H]-NOESY

ASCAN is a new algorithm for automatic sequence-specific NMR assignment of amino acid side-chains in proteins, which uses as input the primary structure of the protein, chemical shift lists of (1)H(N), (15)N, (13)C(alpha), (13)C(beta) and possibly (1)H(alpha) from the previous polypeptide backbone assignment, and one or several 3D (13)C- or (15)N-resolved [(1)H,(1)H]-NOESY spectra. ASCAN has also been laid out for the use of TOCSY-type data sets as supplementary input. The program assigns new resonances based on comparison of the NMR signals expected from the chemical structure with the experimentally observed NOESY peak patterns. The core parts of the algorithm are a procedure for generating expected peak positions, which is based on variable combinations of assigned and unassigned resonances that arise for the different amino acid types during the assignment procedure, and a corresponding set of acceptance criteria for assignments based on the NMR experiments used. Expected patterns of NOESY cross peaks involving unassigned resonances are generated using the list of previously assigned resonances, and tentative chemical shift values for the unassigned signals taken from the BMRB statistics for globular proteins. Use of this approach with the 101-amino acid residue protein FimD(25-125) resulted in 84% of the hydrogen atoms and their covalently bound heavy atoms being assigned with a correctness rate of 90%. Use of these side-chain assignments as input for automated NOE assignment and structure calculation with the ATNOS/CANDID/DYANA program suite yielded structure bundles of comparable quality, in terms of precision and accuracy of the atomic coordinates, as those of a reference structure determined with interactive assignment procedures. A rationale for the high quality of the ASCAN-based structure determination results from an analysis of the distribution of the assigned side chains, which revealed near-complete assignments in the core of the protein, with most of the incompletely assigned residues located at or near the protein surface.     

<Superscript>1</Superscript>H NMR metabolic fingerprinting of saffron extracts

The aim of this study was to explore feasibility of ¹H NMR metabolic fingerprinting for discrimination of authenticity of saffron using principal component analysis (PCA) modeling. Authentic reference Iranian saffron (n = 31) and commercial samples (n = 32) were used. Cross-validated PCA models based on ¹H NMR spectra of solutions prepared by direct extraction of grinded saffron with methanol-d ₄ distinguished reference Iranian saffron samples from commercial samples that formed several distinct clusters, some of which represent falsified samples as confirmed by microscopic analysis. The production sites and drying conditions of the authentic reference Iranian samples were not reflected in the current dataset. Picrocrocin and glycosyl esters of crocetin emerged as the most important ¹H NMR markers of authentic saffron by using statistical correlation spectroscopy. In conclusion, ¹H NMR spectra of saffron extracts combined with pattern recognition by PCA provide immediate means of unsupervised classification of saffron samples.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignment of lissencephaly-1 homology (LisH) domain homodimer of human two-hybrid-associated protein 1 with RanBPM (Twa1)

The CTLH complex is a large, highly conserved eukaryotic complex composed of eight proteins that has been associated to several cellular functions, more often described as an E3 ubiquitin ligase complex involved in protein degradation through ubiquitination but also via vacuole-dependent degradation. A common feature observed in several components of this complex is the presence of the domains lissencephaly-1 homology (LisH) and C-terminal to LisH (CTLH). The LisH domain is found in several proteins involved in chromosome segregation, microtubule dynamics, and cell migration. Also, this domain participates in protein dimerization, besides affecting protein half-life, and influencing in specific cellular localization. Among the proteins found in the CTLH complex, Twa1 (Two-hybrid-associated protein 1 with RanBPM), also known as Gid8 (glucose-induced degradation protein 8 homolog) is the smallest, being a good model for structural studies by NMR. In this work we report the chemical shift assignments of the homodimeric LisH domain of Twa1, as a first step to determine its solution structure.     

Critical test of some computational methods for prediction of NMR <Superscript>1</Superscript>H and <Superscript>13</Superscript>C chemical shifts

Performance of 18 DFT functionals (B1B95, B3LYP, B3PW91, B97D, BHandHLYP, BMK, CAM-B3LYP, HSEh1PBE, M06-L, mPW1PW91, O3LYP, OLYP, OPBE, PBE1PBE, tHCTHhyb, TPSSh, wB97xD, VSXC) in combinations with six basis sets (cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, aug-cc-pVTZ, IGLO-II, and IGLO-III) and three methods for calculating magnetic shieldings (GIAO, CSGT, IGAIM) was tested for predicting ¹H and ¹³C chemical shifts for 25 organic compounds, for altogether 86 H and 88 C atoms. Proton shifts varied between 1.03 ppm to 12.00 ppm and carbon shifts between 7.87 ppm to 209.28 ppm. It was found that the best method for calculating ¹³C shifts is PBE1PBE/aug-cc-pVDZ with CSGT or IGAIM approaches (mae?=?1.66 ppm), for ¹H the best results were obtained with HSEh1PBE, mPW1PW91, PBE1PBE, CAM-B3LYP, and B3PW91 functionals with cc-pVTZ basis set and with CSGT or IGAIM approaches (mae?=?0.28 ppm). We found that often larger basis sets do not give better results for chemical shifts. The best basis sets for calculating ¹H and ¹³C chemical shifts were cc-pVTZ and aug-cc-pVDZ, respectively. CSGT and IGAIM NMR approaches can perform really well and are in most cases better than popular GIAO approach.

Open image in new window

Graphical Abstract Mean absolute errors for ¹H and ¹³C chemical shifts and computational times of neutral toluene molecule with aug-cc-pVDZ basis set and CSGT approach

     

A multidimensional <Superscript>1</Superscript>H NMR lipidomics workflow to address chemical food safety issues

Introduction

Although it is still at a very early stage compared to its mass spectrometry (MS) counterpart, proton nuclear magnetic resonance (NMR) lipidomics is worth being investigated as an original and complementary solution for lipidomics. Dedicated sample preparation protocols and adapted data acquisition methods have to be developed to set up an NMR lipidomics workflow; in particular, the considerable overlap observed for lipid signals on 1D spectra may hamper its applicability.

Objectives

The study describes the development of a complete proton NMR lipidomics workflow for application to serum fingerprinting. It includes the assessment of fast 2D NMR strategies, which, besides reducing signal overlap by spreading the signals along a second dimension, offer compatibility with the high-throughput requirements of food quality characterization.

Method

The robustness of the developed sample preparation protocol is assessed in terms of repeatability and ability to provide informative fingerprints; further, different NMR acquisition schemes—including classical 1D, fast 2D based on non-uniform sampling or ultrafast schemes—are evaluated and compared. Finally, as a proof of concept, the developed workflow is applied to characterize lipid profiles disruption in serum from β-agonists diet fed pigs.

Results

Our results show the ability of the workflow to discriminate efficiently sample groups based on their lipidic profile, while using fast 2D NMR methods in an automated acquisition framework.

Conclusion

This work demonstrates the potential of fast multidimensional ¹H NMR—suited with an appropriate sample preparation—for lipidomics fingerprinting as well as its applicability to address chemical food safety issues.

     

Automatic <Superscript>13</Superscript>C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ±?0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ±?0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ±?0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.     

Liquid seaweed extracts identified using <Superscript>1</Superscript>H NMR profiles

Aqueous extracts of Ascophyllum nodosum and several other brown seaweeds are manufactured commercially and widely distributed for use on agricultural crops. The increasingly regulated international trade in such products requires that they be standardized and defined to a degree not previously required. We examined commercially available extracts using quantitative ¹H NMR and principal components analysis (PCA) techniques. Extracts manufactured over a 4-year period using the same process exhibited characteristic profiles that, on PCA, clustered as a discrete group distinct from the other commercial products examined. In addition to recognizing extracts made from different seaweeds, analysis of the ¹H spectra in the 0.35–4.70 ppm region allowed us to distinguish amongst extracts produced from the same algal species by different manufacturers. This result established that the process used to make an extract is an important variable in defining its composition. A comparison of the ¹H NMR integrals for the regions 1.0–3.0 ppm and 3.0–4.38 ppm revealed small but significant changes in the A. nodosum spectra that we attribute to seasonal variation in gross composition of the harvested seaweed. Such changes are reflected in the PCA scores plots and contribute to the scatter observed within the data point cluster observed for Acadian soluble extracts when all data are pooled. Quantitative analysis using ¹H NMR (qNMR) with a certified external standard (caffeine) showed a linear relationship with extract concentration over at least an order of magnitude (2.5–33 mg/mL; R ² > 0.97) for both spectral regions integrated. We conclude that qNMR can be used to profile (or “fingerprint”) commercial seaweed extracts and to quantify the amount of extract present relative to a suitably chosen standard. Issued as NRCC no. 42,652.     

Use of EDTA to minimize ionic strength dependent frequency shifts in the <Superscript>1</Superscript>H NMR spectra of urine

The ¹H NMR spectrum of urine exhibits a large number of detectable metabolites and is, therefore, highly suitable for the study of perturbations caused by disease, toxicity, nutrition or environmental factors in humans and animals. However, variations in the chemical shifts and intensities due to altered pH and ionic strength present a challenge in NMR-based studies. With a view towards understanding and minimizing the effects of these variations, we have extensively studied the effects of ionic strength and pH on the chemical shifts of common urine metabolites and their possible reduction using EDTA (ethylenediaminetetraacetic acid). ¹H NMR chemical shifts for alanine, citrate, creatinine, dimethylamine, glycine, histidine, hippurate, formate and the internal reference, TSP (trimethylsilylpropionic acid-d₄, sodium salt) obtained under different conditions were used to assess each effect individually. EDTA minimizes the frequency shifts of the metabolites that have a propensity for metal binding. Chelation of such metal ions is evident from the appearance of signals from EDTA complexed to divalent metal ions such as calcium and magnesium. Not surprisingly, increasing the buffer concentration or buffer volume also minimizes pH dependent frequency shifts. The combination of EDTA and an appropriate buffer effectively minimizes both pH dependent frequency shifts and ionic strength dependent intensity variations in urine NMR spectra. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

4D prediction of protein <Superscript>1</Superscript>H chemical shifts

A 4D approach for protein ¹H chemical shift prediction was explored. The 4th dimension is the molecular flexibility, mapped using molecular dynamics simulations. The chemical shifts were predicted with a principal component model based on atom coordinates from a database of 40 protein structures. When compared to the corresponding non-dynamic (3D) model, the 4th dimension improved prediction by 6–7%. The prediction method achieved RMS errors of 0.29 and 0.50 ppm for Hα and HN shifts, respectively. However, for individual proteins the RMS errors were 0.17–0.34 and 0.34–0.65 ppm for the Hα and HN shifts, respectively. X-ray structures gave better predictions than the corresponding NMR structures, indicating that chemical shifts contain invaluable information about local structures. The ¹H chemical shift prediction tool 4DSPOT is available from .