Articles | Volume 2, issue 2
Magn. Reson., 2, 619–627, 2021

Special issue: Geoffrey Bodenhausen Festschrift

Magn. Reson., 2, 619–627, 2021

Research article 10 Aug 2021

Research article | 10 Aug 2021

Virtual decoupling to break the simplification versus resolution trade-off in nuclear magnetic resonance of complex metabolic mixtures

Virtual decoupling to break the simplification versus resolution trade-off in nuclear magnetic resonance of complex metabolic mixtures
Cyril Charlier1, Neil Cox1, Sophie Martine Prud'homme2,a, Alain Geffard2, Jean-Marc Nuzillard3, Burkhard Luy4, and Guy Lippens1 Cyril Charlier et al.
  • 1Toulouse Biotechnology Institute (TBI), Université de Toulouse, CNRS, INRAE, INSA, Toulouse, France
  • 2Université de Reims Champagne-Ardenne (URCA), UMR-I 02 SEBIO (Stress Environnementaux et Biosurveillance des milieux aquatiques), Moulin de la Housse, Reims, France
  • 3Université de Reims Champagne Ardenne, CNRS, ICMR UMR 7312, 51097 Reims, France
  • 4Institute for Biological Interfaces 4 – Magnetic Resonance, Karlsruhe Institute of Technology (KIT), Herrmann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
  • apresent address: Université de Lorraine, CNRS, LIEC, 57000, Metz, France

Correspondence: Cyril Charlier ( and Guy Lippens (


The heteronuclear single quantum correlation (HSQC) experiment developed by Bodenhausen and Ruben (1980) in the early days of modern nuclear magnetic resonance (NMR) is without a doubt one of the most widely used experiments, with applications in almost every aspect of NMR including metabolomics. Acquiring this experiment, however, always implies a trade-off: simplification versus resolution. Here, we present a method that artificially lifts this barrier and demonstrate its application towards metabolite identification in a complex mixture. Based on the measurement of clean in-phase and clean anti-phase (CLIP/CLAP) HSQC spectra (Enthart et al., 2008), we construct a virtually decoupled HSQC (vd-HSQC) spectrum that maintains the highest possible resolution in the proton dimension. Combining this vd-HSQC spectrum with a J-resolved spectrum (Pell and Keeler, 2007) provides useful information for the one-dimensional proton spectrum assignment and for the identification of metabolites in Dreissena polymorpha (Prud'homme et al., 2020).

1 Introduction

The recognition that a given nucleus was characterized by a specific chemical shift value depending on its exact environment in a molecule (Proctor and Yu, 1950; Dickinson, 1950) ushered nuclear magnetic resonance (NMR) from its initial discovery in a nuclear physics department (Purcell et al., 1946; Bloch et al., 1946) into the chemistry sphere. Magnetic shielding of nuclei was identified at the origin of the phenomenon (Ramsey, 1950) and representative chemical shift values of the different protons in organic molecules were rapidly established (Arnold et al., 1951; Bernstein and Schneider, 1956). The further recognition of homonuclear scalar coupling (the “J-coupling”) patterns (Gutowsky and McCall, 1951) increased the information content, and notions as singlet, doublet, triplet, etc. were rapidly used to characterize the individual proton lines in the NMR spectra and to assist the identification of the molecule under study. Realizing that these same J couplings could be used to transfer magnetization from one spin to another, in 1971, Jeener proposed an indirect acquisition scheme to reconstruct a 2-D map (Jeener and Alewaeters, 2016), whereby the off-diagonal peaks connect different protons through their J coupling. The homonuclear correlation spectroscopy (COSY) experiment was born (Aue et al., 1976), and many two- and higher-dimensional homonuclear pulse sequences would adopt a similar principle.

One problem with proton (1H) NMR is the limited range of chemical shift values (typically 10 ppm) with line widths on the order of magnitude of 1 Hz, leading to a crowded spectrum especially when complex mixtures are studied. The carbon (13C) chemical shift range is significantly larger ( 200 ppm), making it an attractive alternative to characterize a (mixture of) molecule(s). However, its low natural abundance ( 1.1 %), inherently poorer sensitivity with its 4-fold lower gyromagnetic ratio and long relaxation times increase the required acquisition times by a factor of 106 or higher when a similar S/N as for the proton spectrum is desired. The introduction of a two-dimensional (2-D) heteronuclear single quantum correlation (HSQC) spectrum by Bodenhausen and Ruben in 1980 has literally been a game changer in the NMR field (Bodenhausen and Ruben, 1980). In this pioneering publication, cited 3000 times as of the time of writing, “the detection of NMR spectra of less sensitive nuclei coupled to protons was shown to be significantly improved by a two-dimensional Fourier transform technique involving a double transfer of polarization”. Building on the ideas of 2-D NMR (Jeener and Alewaeters, 2016; Aue et al., 1976) and magnetization transfer between a heteronucleus and proton by their scalar coupling (Maudsley and Ernst, 1977; Morris and Freeman, 1979), polarization transfer from the proton to the insensitive heteronucleus in the excitation step and back to proton for the detection led to a tremendous gain in sensitivity. The HSQC experiment was born, and although initially described for a proton directly linked to its amide nitrogen, it was readily applied to other nuclei such as 199Hg (Roberts et al., 1980), 113Cd (Live et al., 1985), 31P (Bolton and Bodenhausen, 1982) and 13C (Bendall et al., 1983). The HSQC experiment applied to the latter 1H–13C pair has not only solved sensitivity problems associated with the direct 13C 1-D spectrum but also greatly increased the information content of both individual 1H and 13C spectra as it connects a carbon signal directly to its proton binding partner. The success story of the HSQC spectrum was born and formed the basis for a myriad of 2- to n-D spectra used from protein NMR to the analysis of metabolomic samples.

Indeed, beyond synthetic chemistry, NMR has found its place in the realm of metabolomics (Emwas et al., 2019; Wishart, 2019), and together with mass spectrometry, it is the method of choice for characterizing complex mixtures such as biofluids. If throughput is required, the 1-D 1H spectrum remains the general workhorse, whose assignment can be based on database searches, on the spiking of relevant standards within the sample or on analysis of the corresponding 2-D spectra acquired on one or more samples. The 1H–13C HSQC spectrum proposed by Bodenhausen and Ruben is again a crucial element in the latter approach, as it spreads out the proton signals according to the 13C chemical shift of its attached carbons and allows connection of both in the same time. The use of 2-D experiments in the field of metabolomics, despite inherent longer experimental times, has grown extensively, in part due to the introduction of ultra-fast acquisition (Tal and Frydman, 2010), ASAP (acceleration by sharing adjacent polarization) techniques (Schätzlein et al., 2018) and non-uniformly sampling (NUS) acquisition scheme (Guennec et al., 2014; Schlippenbach et al., 2018; Zhang et al., 2020). Based on these developments, the use of 2-D experiments for larger number of samples has become within reach and will probably place the 1H–13C HSQC spectrum as a cornerstone of metabolomic studies by NMR. However, for quantitative studies, the 1-D proton spectrum with a long relaxation delay will remain the reference, as both limits on the relaxation time and variable J-coupling values make the HSQC spectrum inherently non-quantitative. Assignment of the 1-D spectrum on the basis of the HSQC spectrum and extracting quantitative information from the former is however one way to proceed.

In the initial HSQC paper, it was recognized that obtaining a single peak for the 15N–1H pair requires decoupling both during the indirect detection (obtained by a single proton π pulse) and the direct proton detection (through a broadband decoupling scheme). However, for technical reasons related to the large spread of the 13C resonances, decoupling during the direct acquisition period can lead to sample heating and/or probe arcing. One typically limits acquisition times to 100 ms, and especially with current high S/N cryogenically cooled probes, constructors limit the power delivery during this decoupling period (Bahadoor et al., 2021). As a result, the resolution in the 1H dimension is limited and the homonuclear coupling pattern of the protons cannot be observed, despite the fact that they carry valuable information.

Here, we explore the capacity to virtually decouple a HSQC spectrum without physical decoupling during the acquisition time, thereby removing all physical limits on the acquisition time. We reconstruct the decoupled spectrum by recording two coupled HSQC spectra: one with the direct 1H–13C doublet in phase and one with the doublet in anti-phase. Comparable sequences aimed at measuring the 1J coupling constant have been published before, especially in the framework of determining the residual dipolar couplings for samples partially oriented in the liquid phase (Andersson et al., 1998). However, here we use both spectra only to distinguish the up- and downfield components of the doublet, thereby providing the necessary information for shifting one component with its complex homonuclear coupling pattern to the central position. As a result, we obtain a HSQC spectrum with a high resolution in the proton dimension. Traces from it can be immediately superimposed on the 1-D 1H spectrum or on a high-resolution J-resolved (J-res) spectrum, thereby helping the assignment of the latter. We first demonstrate the approach on the spectrum of an isolated oligosaccharide and show its use on a zebra mussel (Dreissena polymorpha) hydrophilic extract that we recently studied in the framework of an ecotoxicological study (Prud'homme et al., 2020).

2 Methods

All experiments were performed on an Avance III 800 MHz spectrometer (Bruker BioSpin GmbH, Karlsruhe, Germany), equipped with a 1.7 mm triple-resonance HCP or a 5 mm quadruple resonance QCI-P (H/P-C/N/D) cryoprobe and were recorded at 298 K. Data were processed with TopSpin 4.0.8 (Bruker BioSpin GmbH, Karlsruhe, Germany). All proton spectra were referenced to the methyl proton of 3-trimethylsilylpropionic acid-d4 (TSP-d4). 13C chemical shifts were determined by indirect referencing (Markley et al., 1998).

Pulse sequences were initially tested on a sample of 1 mg of dextran oligosaccharide dissolved in 40 µL of D2O plus 2 µL of TSP-d4 on the 1.7 mm cryoprobe. Clean in-phase and clean anti-phase (CLIP/CLAP) experiments were acquired with 16 384 (1H) × 256 (13C) time points, a spectral width of 12.0172 ppm (1H) × 60 ppm (13C) and NS = 8, DS = 16, RG = 512, d1 = 1.5 s, experimental time = 1 h 23 min 2 s (CLIP) and 1 h 22 min 39 s (CLAP). Data were transformed to a matrix of 16 384 (1H) × 1024 (13C) frequency points. 1H–13C HSQC with decoupling was recorded with 4096 (1H) × 256 (13C) time points, a spectral width of 12.0172 ppm (1H) × 60 ppm (13C) and NS = 8, DS = 32, RG = 512, d1 = 1 s, experimental time = 43 min 2 s. Data were transformed to a matrix of 4096 (1H) × 1024 (13C) frequency points.

Hydrophilic extracts of zebra mussel were prepared a described in (Prud'homme et al., 2020). Briefly, 50 mg of biomass were dried prior to resuspension in 50 µL of 100 mM potassium phosphate buffer dissolved in D2O supplemented with 1 mM sodium azide and 0.5 mM TSP-d4. From this solution, 40 µL were then transferred into the NMR sample. CLIP/CLAP experiments were acquired with 16 384 (1H) × 256 (13C) time points, a spectral width of 13.9486 ppm (1H) × 100 ppm (13C) and NS = 64, DS = 64, RG =  512, d1 = 1.5 s, experimental time = 10 h 28 min 51 s (CLIP) and 10 h 25 min 48 s (CLAP). Data were transformed to a matrix of 16 384 (1H) × 1024 (13C) frequency points. The 1-D 1H, 1H–13C HSQC with decoupling and the J-resolved experiments were recorded on the 5 mm cryoprobe on a sample of 200 µL. A high-resolution 1-D 1H spectrum was acquired with 32 768 time points, a spectral width of 12.0172 ppm and NS = 64, DS = 4, RG = 8, d1 = 5.0 s. The 1H–13C HSQC was recorded with 2048 (1H) × 256 (13C) time points, a spectral width of 13.9486 ppm (1H) × 100 ppm (13C) and NS = 64, DS = 64, RG =  912, d1 = 1 s, experimental time = 8 h 2 min 26 s. Data were transformed to a matrix of 2048 (1H) × 512 (13C) frequency points. The J-resolved experiment was recorded with 16 384 time points, a spectral width of 12.0172 ppm and NS = 128, DS = 64, RG = 912, d1 = 1.32 s for an experimental time of 18 h 1 min 29 s.

3 Results and discussion

As pointed out in the original paper of Bodenhausen and Ruben, without physical decoupling of the direct 1H–13C interaction, the theoretical pattern of a peak in the direct dimension (F2) of the HSQC spectrum is a doublet, due to the heteronuclear coupling constant (1JCH for a 1H /13C spectrum). Each component of the doublet can moreover display a more or less complex coupling pattern due to homo- (JHH) or heteronuclear (for example, proton–phosphorus coupling, JHP) coupling interactions involving the observed proton. The active 1JCH coupling can be removed by heteronuclear spin decoupling during signal acquisition and renders more easily interpretable 2-D maps with a single cross peak per 1H /13C pair. Decoupling, however, has its drawbacks, especially when the bandwidth to be decoupled increases. Although numerous approaches have been developed to decrease the delivered power (Kobzar et al., 2004; Kupče, 2020) and hence reduce potential sample heating, limiting the acquisition time remains the first option, be it at the detriment of resolution and hence loss of intrinsic coupling pattern. Ideally, one would want both – a decoupled 1H spectrum with a single peak per 1H /13C pair with a high resolution leaving the intricate coupling pattern of the proton intact.

Omitting the decoupling during the acquisition is the obvious solution to the heating problem but leads to a doubling of the number of peaks in the spectrum and thereby potentially increases spectral overlap. Moreover, the peaks resonate at ±J/2 Hz from their true chemical shift, whereby J is the active coupling at the origin of the cross peak and hence cannot directly be superimposed on the 1-D proton spectrum dominated by the contribution of the 12C-linked protons. For isolated peaks, however, this doublet structure can be easily recognized, and shifting both components back to their central position solves the problem.

This reconstruction of the decoupled spectrum meets problems, however, when peak overlap becomes important. One issue is whether a given peak corresponds to the down- or upfield component of the doublet, necessitating to look to the left or right for its corresponding component. For this, however, solutions have long existed, notably in the endeavour to measure (small) coupling constants. Combining absorption/dispersion non-decoupled spectra has been used to measure coupling constants down to few hertz for well-resolved signals (Kessler et al., 1985; Oschkinat and Freeman, 1984). A 2-D X-filtered total correlation spectroscopy (TOCSY) (Wollborn and Leibfritz, 1992) or its corresponding 1-D version (Nuzillard and Bernassau, 1994) was proposed to improve the accuracy of the coupling constant measurement. Other experiments such as α/β-HSQC/heteronuclear multiple quantum coherence (HMQC) (Andersson et al., 1998) or α/β-HSQC-TOCSY (Koźmiński, 1999) were developed based on the concept of spin-state selection (Meissner et al., 1997) to selectively observe in each subspectrum a single component of the coupling multiplet. However, as these methods are not applicable without directly linked protons to the carbon, workarounds have been proposed using heteronuclear single quantum multiple bond correlation (HSQMBC) (Williamson et al., 2000) or HSQC spectra (Titman et al., 1989) supplemented with Carr–Purcell–Meiboom–Gill (CPMG) train pulses (Boros and Kövér, 2011; Kövér et al., 2006) to improve the intrinsic twisted linewidth due to the evolution of homonuclear proton–proton couplings. New pulse sequences were later developed to facilitate measurement of one-bond coupling constants by eliminating sources of line-shape distortion, such as the clean in-phase (CLIP)-HSQC experiment (Enthart et al., 2008). Acquiring the anti-phase magnetization in a separate clean anti-phase (CLAP)-HSQC experiment followed by addition/subtraction of both spectra yields high-resolution α and β-state subspectra. BEBOP (broadband excitation by optimized pulses)/BIBOP (broadband inversion by optimized pulses) (Kobzar et al., 2004; Luy et al., 2005) pulses on 13C can be used to obtain uniform performance of the experiments across the large carbon spectral width. Long-range proton–carbon coupling constants have been measured with the CLIP-HSQMBC (Saurí et al., 2013), PIP-HSQMBC (Castañar et al., 2014) and CSSF-CLIP-HSQMBC (Moreno et al., 2019).

Proton multiplet patterns are often seen as a source of increased peak overlap that needs to be addressed by a myriad of pure-shift methods (Zangger and Sterk, 1997; Foroozandeh et al., 2014; Castañar and Parella, 2015). These have been combined with HSQC sequences to alleviate extraction of coupling constants (Timári et al., 2016). However, we believe that the information contained in these patterns can be crucial towards the assignment and identification of metabolites in complex mixtures. Herein, we combine the original CLIP/CLAP pulse sequences (Enthart et al., 2008) with virtual decoupling of a 1H–13C spectrum based on automated recognition of the α and β states of the multiplets and subsequent back-shifting of the high-resolution lines to their true chemical shift. We thereby obtain both a single peak per 1H /13C moiety and the required high resolution for identification of its J-coupling pattern.

We first tested the potential of virtual decoupling on an isolated dextran oligosaccharide. Both CLIP/CLAP experiments acquired under identical conditions led to two high-quality spectra, as shown in Figs. 1a, b and S1 in the Supplement. The CLIP experiment in which the anti-phase proton magnetization is refocused with a 180 13C pulse during the final INEPT followed by a 90 13C pulse prior to detection to remove any dispersive component shows two purely in-phase signals separated by the heteronuclear coupling constant 1JCH. In the CLAP experiment, these two pulses are omitted, leading to the observation of the anti-phase doublet. The absence of decoupling during the acquisition in the CLIP/CLAP experiment allows the observation of each component of the doublets with high resolution (Fig. S1). We developed a semi-automatic Python procedure for the virtual decoupling that is based on the nmrglue package (Helmus and Jaroniec, 2013) and is usable within TopSpin 4. It can be downloaded from (last access: 4 August 2021). While launching the script from TopSpin, the user will have several options that are detailed in the tutorial available online. First, taking advantage of the sign difference between the α and β components of the doublet in the CLAP experiment (Fig. S2), we reduce the resolution and hence the fine structure of both components by processing the CLAP spectrum with only 1024 points in the 1H dimension, thereby intentionally destroying the high-resolution J-coupling pattern and obtaining just a single peak per component. Automatic peak detection is then performed twice: once on the positive components and once on the negative components of the spectrum and can be performed on the full spectral width or only for a user-defined region with a threshold that can be adapted if necessary. A clustering step based on 13C chemical shifts obtains a pairwise selection of the signals. The central position of this pair is used to define the J/2 value for this particular 1H /13C pair, thereby individualizing the back shift and avoiding the impossible use of a common value for all peaks. If exactly two signals of opposite sign are found with the same 13C frequency (in digital points), the script will automatically bring the upfield component extracted from CLIP-HSQC spectrum with full resolution (processed with 16 384 points in the time domain) back to the centre of the doublet and thereby create a virtually decoupled HSQC (vd-HSQC) spectrum (Figs. 1c and S3). In the case of more than two signals identified at the same 13C frequency value, a user interface will pop up and allow manual selection and clustering of the signals.

Figure 1Application of virtual decoupling isolated dextran oligosaccharide. CLIP (a) and CLAP (b) are shown with positive contours in blue and negative contours in green. (c) Virtually decoupled HSQC. (d) Decoupled HSQC. The traces shown above panels (c) and (d) are extracted at the centre of the peak of the 2-D spectra.


Figure 2Application of virtual decoupling to a zebra mussel sample. (a) Overlay of the virtually decoupled HSQC (red) and decoupled HSQC (grey) and (b) 1-D 1H trace extracted from the vd-HSQC at 55.623 ppm. (c) J-resolved spectrum with (d) the 1H projections from the coloured areas on the 2-D spectrum and named A/B (top) and 1/2/3 (bottom). Coloured areas were extracted using the zerter command in TopSpin (Cox et al., 2021). (e) Virtually decoupled HSQC (red) and decoupled HSQC (grey) and (f) 1-D 1H trace extracted from the vd-HSQC at 38.426 ppm (top) and 43.298 ppm (bottom). 1/2/3 correspond to the three doublets of doublets identified in areas 1/2/3 of the J-res experiment.


After developing the procedure on this simple oligosaccharide sample, we turned to a zebra mussel (Dreissena polymorpha) extract as representative of a more complex sample. We recently annotated the spectrum in the framework of a NMR-based metabolomic ecotoxicological study and identified and assigned a large number of metabolites that can be used as reporters of environmental health (Prud'homme et al., 2020). While the initial metabolomic study of the zebra mussel (Watanabe et al., 2015) was based on 1-D 1H and 2-D 1H–13C HSQC spectrum, our more recent work (Prud'homme et al., 2020) also included 2-D 1H–1H, 2-D 1H–31P and J-resolved spectra and led to the assignment of  50 % of the signals observed on the 1-D 1H spectrum. To evaluate the value of virtual decoupling on such a complex sample, we measured the CLIP/CLAP spectra of the D. polymorpha whole-body polar extract in order to reconstruct the vd-HSQC and combined it with a high-resolution J-res spectrum. Here, we focus on the region of the trimethylamines from 3.1 to 3.3 ppm (Fig. 2). The J-resolved experiment shows the presence of two singlets around 3.12 ppm (peaks A and B in Fig. 2c). Because they almost overlap and moreover differ by a factor of 10 in intensity, the decoupled HSQC (Fig. 2a grey spectrum) only shows a single peak in which the two signals cannot be resolved. However, in the corresponding region of the virtually decoupled HSQC (Fig. 2a), the full proton pattern containing two signals is observed, confirming that both protons are linked to a carbon with identical chemical shift. The trace extracted at 55.623 ppm through the vd-HSQC (Fig. 2b) perfectly matches the trace extracted at 0 Hz in the J-res spectrum (Fig. 2d), in agreement with a negligible isotope shift of the 13C nucleus on the proton chemical shift (Tiainen et al., 2010). This first example illustrates the ability of the vd-HSQC to separate two signals that cannot be resolved with a conventional decoupled HSQC (d-HSQC).

A second example of the vd-HSQC's capacity to assist and enhance the assignment is given by the three doublets of doublets (dd) highlighted by the numbers on the J-res spectrum (Fig. 2c). Using the zerter programme recently developed in our group (Cox et al., 2021) first to isolate and zero the singlet at 3.27 ppm and then extract the individual pseudo-triplets, each extracted trace shows the expected 1 : 2 : 1 pattern (lines 1, 2 and 3 in Fig. 2d). These 1H 1-D traces can be used to search through the vd-HSQC spectrum to identify the corresponding 13C resonance frequencies. Because of the high resolution of the vd-HSQC, the scan is not only based on the chemical shift value but equally on the proton–proton coupling pattern, thereby significantly enhancing the information content. Two patterns were identified at 13C resonance frequencies of 38.426 and 43.298 ppm (Fig. 2e). In both cases, the vd-HSQC shows the correct 1H frequency and the correct coupling pattern. Indeed, the trace extracted from the signal at the higher 13C frequency (Fig. 2f) shows the presence of a 1 : 2 : 1 coupling pattern which fits that of the projected spectrum of the J-res. The peak identified at the lower 13C frequency in the vd-HSQC shows two 1 : 2 : 1 patterns for which the 1H trace can be matched with the projection of the J-res spectrum. This combined information enhances our confidence that the three dd signals, although closely together in proton chemical shift and coupling pattern, actually represent organic moieties that differ by 5 ppm in carbon chemical shift. If we compare this with a scan through the d-HSQC, assigning the broad peak at 38.426 ppm to the two dd signals seems more hazardous, especially as the singlet at 3.271 ppm could also be erroneously assigned to this proton frequency. The coupling pattern hence increases the information content and, when combined with a high-resolution J-resolved spectrum, can lead to accurate information on individual proton resonances.

Figure 3(a) 1H traces extracted from the virtually decoupling at different 13C chemical shifts and (b) 1-D 1H spectrum with high resolution (acquired with 16 384​​​​​​​ points). Signals are labelled as in Fig. 2.


As stated before, presently most studies aiming at the analysis of a large number of samples are based on the high-resolution 1-D proton spectrum. To evaluate how well the information of our vd-HSQC spectrum can be transferred to the 1-D spectrum, we show in Fig. 3 the traces extracted from the vd-HSQC spectrum on top of the high-resolution 1-D spectrum of the D. polymorpha extract (Fig. 3). The two singlets at the 13C frequency of 55.623 ppm can immediately be identified in the 1-D spectrum, and their corresponding intensities represent the concentrations of both entities if the 1-D spectrum is acquired with a sufficiently long relaxation delay. But even in the crowded region around 3.25 ppm, we can use the traces of the vd-HSQC to assign the individual pseudo-triplets and associate them with their corresponding 13C values. Interestingly, the 1-D spectrum also contains the singlet at 3.271 ppm (indicated by an asterisk (*) in Fig. 3), but the combined information from the high-resolution vd-HSQC and J-res spectra immediately singles it out as such.

4 Conclusion

Beyond the time constraints associated with the 2-D nature of the spectrum, the use of the 1H–13C HSQC in metabolomics is always associated with an NMR dilemma: broadband decoupling during acquisition greatly simplifies the spectra, but due to limitations in the duty cycle it leads to the loss of all information about proton multiplicities; no decoupling during acquisition leads to the observation of the complete proton multiplet pattern but doubles the number of resonances. Virtual decoupling of the spectrum is of great interest for this issue because it provides a single resonance that maintains the highest achievable resolution. Here, we explore the benefits of this approach and demonstrate that the virtually decoupled spectrum can be directly compared with the high-resolution 1-D proton spectrum, thereby helping in the process of metabolite identification over many samples. The methodology presented here, while being robust, is rather simplistic and will benefit from the rise of deep learning and artificial intelligence applications to NMR spectra (Karunanithy et al., 2021) to provide a new entry point into metabolite databases.

Code and data availability

The code written in Python for virtual decoupling and compatible within TopSpin 4 can be found on GitHub: (, Charlier et al., 2021). A readme file guides potential users through the data used in the paper.


The supplement related to this article is available online at:

Author contributions

CC was responsible for software, conceptualization and writing (review and editing); NC, AG, JMN and SMP were responsible for resources, review and editing; BL was responsible for software, review and editing; GL was responsible for conceptualization, project administration, supervision and writing (review and editing).​​​​​​​

Competing interests

The authors declare that they have no conflict of interest.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Geoffrey Bodenhausen Festschrift”. It is not associated with a conference.


The authors thanks Pierre Millard​​​​​​​ (INRAE, TBI, Toulouse) for insightful discussions. We thank Edern Cahoreau and Lindsay Peyriga for expert support in the MetaToul (Toulouse metabolomics & fluxomics facilities,, last access: 5 August 2021) NMR facility. MetaToul is part of the French National Infrastructure for Metabolomics and Fluxomics MetaboHUB-AR-11-INBS-0010 (, last access: 5 August 2021) and is supported by the Région Midi-Pyrénées, the ERDF, the SICOVAL and the French Minister of Education & Research, who are all gratefully acknowledged.

Financial support

This research has been supported by the CNRS.

Review statement

This paper was edited by Fabien Ferrage and reviewed by Oscar Millet and two anonymous referees.


Andersson, P., Nordstrand, K., Sunnerhagen, M., Liepinsh, E., Turovskis, I., and Otting, G.: Heteronuclear correlation experiments for the determination of one-bond coupling constants, J. Biomol. NMR, 11, 445–450,, 1998. 

Arnold, J. T., Dharmatti, S. S., and Packard, M. E.: Chemical Effects on Nuclear Induction Signals from Organic Compounds, J. Chem. Phys., 19, 507–507,, 1951. 

Aue, W. P., Bartholdi, E., and Ernst, R. R.: Two-dimensional spectroscopy. Application to nuclear magnetic resonance, J. Chem. Phys., 64, 2229–2246,, 1976. 

Bahadoor, A., Brinkmann, A., and Melanson, J. E.: 13C-Satellite Decoupling Strategies for Improving Accuracy in Quantitative Nuclear Magnetic Resonance, Anal. Chem., 93, 851–858,, 2021. 

Bendall, M. R., Pegg, D. T., and Doddrell, D. M.: Pulse sequences utilizing the correlated motion of coupled heteronuclei in the transverse plane of the doubly rotating frame, J. Magn. Reson., 52, 81–117,, 1983. 

Bernstein, H. J. and Schneider, W. G.: Nuclear Magnetic Resonance Spectra of Pyridine and Some Deuterated and Methylated Pyridines, J. Chem. Phys., 24, 469–470,, 1956. 

Bloch, F., Hansen, W. W., and Packard, M.: Nuclear Induction, Phys. Rev., 69, 127–127,, 1946. 

Bodenhausen, G. and Ruben, D. J.: Natural abundance nitrogen-15 NMR by enhanced heteronuclear spectroscopy, Chem. Phys. Lett., 69, 185–189,, 1980. 

Bolton, P. H. and Bodenhausen, G.: Resolution enhancement in heteronuclear two-dimensional spectroscopy by realignment of coherence transfer echoes, J. Magn. Reson., 46, 306–318,, 1982. 

Boros, S. and Kövér, K. E.: Low-power composite CPMG HSQMBC experiment for accurate measurement of long-range heteronuclear coupling constants: Low-power composite CPMG HSQMBC experiment, Magn. Reson. Chem., 49, 106–110,, 2011. 

Castañar, L. and Parella, T.: Broadband 1H homodecoupled NMR experiments: recent developments, methods and applications, Magn. Reson. Chem., 53, 399–426,, 2015. 

Castañar, L., Saurí, J., Williamson, R. T., Virgili, A., and Parella, T.: Pure In-Phase Heteronuclear Correlation NMR Experiments, Angew. Chem. Int. Ed., 53, 8379–8382,, 2014. 

Charlier, C., NMRTeamTBI, and Millard, P.: NMRTeamTBI/VirtualDecoupling: Virtual Decoupling v1.0 (main), Zenodo [code],, 2021. 

Cox, N., Millard, P., Charlier, C., and Lippens, G.: Improved NMR Detection of Phospho-Metabolites in a Complex Mixture, Anal. Chem., 93, 4818–4824,, 2021. 

Dickinson, W. C.: Dependence of the F19 Nuclear Resonance Position on Chemical Compound, Phys. Rev., 77, 736–737,, 1950. 

Emwas, A.-H., Roy, R., McKay, R. T., Tenori, L., Saccenti, E., Gowda, G. A. N., Raftery, D., Alahmari, F., Jaremko, L., Jaremko, M., and Wishart, D. S.: NMR Spectroscopy for Metabolomics Research, Metabolites, 9, 123,, 2019. 

Enthart, A., Freudenberger, J. C., Furrer, J., Kessler, H., and Luy, B.: The CLIP/CLAP-HSQC: Pure absorptive spectra for the measurement of one-bond couplings, J. Magn. Reson., 192, 314–322,, 2008. 

Foroozandeh, M., Adams, R. W., Meharry, N. J., Jeannerat, D., Nilsson, M., and Morris, G. A.: Ultrahigh-resolution NMR spectroscopy, Angew. Chem. Int. Ed. Engl., 53, 6990–6992,, 2014. 

Guennec, A. L., Giraudeau, P., and Caldarelli, S.: Evaluation of Fast 2D NMR for Metabolomics, Anal. Chem., 86, 5946–5954,, 2014. 

Gutowsky, H. S. and McCall, D. W.: Nuclear Magnetic Resonance Fine Structure in Liquids, Phys. Rev., 82, 748–749,, 1951. 

Helmus, J. J. and Jaroniec, C. P.: Nmrglue: an open source Python package for the analysis of multidimensional NMR data, J. Biomol. NMR, 55, 355–367,, 2013. 

Jeener, J. and Alewaeters, G.: “Pulse pair technique in high resolution NMR” a reprint of the historical 1971 lecture notes on two-dimensional spectroscopy, Prog. Nucl. Mag. Res. Sp., 94–95, 75–80,, 2016. 

Karunanithy, G., Mackenzie, H., and Hansen, F.: Virtual Homonuclear Decoupling in Direct Detection NMR Experiments using Deep Neural Networks, ChemRxiv, Cambridge Open Engage, Cambridge, UK,, in review, 2021. 

Kessler, H., Müller, A., and Oschkinat, H.: Differences andsums of traces within,COSY spectra (DISCO) for the extraction of coupling constants: “Decoupling” after the measurement, Magn. Reson. Chem., 23, 844–852,, 1985. 

Kobzar, K., Skinner, T. E., Khaneja, N., Glaser, S. J., and Luy, B.: Exploring the limits of broadband excitation and inversion pulses, J. Magn. Reson., 170, 236–243,, 2004. 

Kövér, K. E., Batta, G., and Fehér, K.: Accurate measurement of long-range heteronuclear coupling constants from undistorted multiplets of an enhanced CPMG-HSQMBC experiment, J. Magn. Reson., 181, 89–97,, 2006. 

Koźmiński, W.: Simplified Multiplet Pattern HSQC-TOCSY Experiment for Accurate Determination of Long-Range Heteronuclear Coupling Constants, J. Magn. Reson., 137, 408–412,, 1999. 

Kupče, Ē.: Perspectives of adiabatic decoupling in liquids, J. Magn. Reson., 318, 106799,, 2020. 

Live, D., Armitage, I. M., Dalgarno, D. C., and Cowburn, D.: Two-Dimensional 1H-113Cd Chemical-Shift Correlation Maps by 1H-Detected Multiple-Quantum NMR in Metal Complexes and Metalloproteins, J. Am. Chem. Soc., 107, 1775–1777,, 1985. 

Luy, B., Kobzar, K., Skinner, T. E., Khaneja, N., and Glaser, S. J.: Construction of universal rotations from point-to-point transformations, J. Magn. Reson., 176, 179–186,, 2005. 

Markley, J. L., Bax, A., Arata, Y., Hilbers, C. W., Kaptein, R., Sykes, B. D., Wright, P. E., and Wüthrich, K.: Recommendations for the presentation of NMR structures of proteins and nucleic acids–IUPAC-IUBMB-IUPAB Inter-Union Task Group on the standardization of data bases of protein and nucleic acid structures determined by NMR spectroscopy, Eur. J. Biochem., 256, 1–15, 1998. 

Maudsley, A. A. and Ernst, R. R.: Indirect detection of magnetic resonance by heteronuclear two-dimensional spectroscopy, Chem. Phys. Lett., 50, 368–372,, 1977. 

Meissner, A., Duus, J. ø., and Sørensen, O. W.: Spin-State-Selective Excitation. Application for E.COSY-Type Measurement ofJHHCoupling Constants, J. Magn. Reson., 128, 92–97,, 1997. 

Moreno, A., Hansen, K. Ø., and Isaksson, J.: CSSF-CLIP-HSQMBC: measurement of heteronuclear coupling constants in severely crowded spectral regions, RSC Adv., 9, 36082–36087,, 2019. 

Morris, G. A. and Freeman, R.: Enhancement of nuclear magnetic resonance signals by polarization transfer, J. Am. Chem. Soc., 101, 760–762, 1979. 

Nuzillard, J. M. and Bernassau, J. M.: A DISCO Approach to the Measurement of Heteronuclear Long-Range Coupling Constants, J. Magn. Reson. B, 103, 284–287,, 1994. 

Oschkinat, H. and Freeman, R.: Fine structure in two-dimensional NMR correlation spectroscopy, J. Magn. Reson., 60, 164–169,, 1984. 

Pell, A. J. and Keeler, J.: Two-dimensional J-spectra with absorption-mode lineshapes, J. Magn. Reson., 189, 293–299,, 2007. 

Proctor, W. G. and Yu, F. C.: The Dependence of a Nuclear Magnetic Resonance Frequency upon Chemical Compound, Phys. Rev., 77, 717–717,, 1950. 

Prud'homme, S. M., Hani, Y. M. I., Cox, N., Lippens, G., Nuzillard, J.-M., and Geffard, A.: The Zebra Mussel (Dreissena polymorpha) as a Model Organism for Ecotoxicological Studies: A Prior 1H NMR Spectrum Interpretation of a Whole Body Extract for Metabolism Monitoring, Metabolites, 10, 256,, 2020. 

Purcell, E. M., Torrey, H. C., and Pound, R. V.: Resonance Absorption by Nuclear Magnetic Moments in a Solid, Phys. Rev., 69, 37–38,, 1946. 

Ramsey, N. F.: Magnetic Shielding of Nuclei in Molecules, Phys. Rev., 78, 699–703,, 1950. 

Roberts, M. F., Vidusek, D. A., and Bodenhausen, G.: Adducts of ethylmercury phosphate with amino acids studied by indirect detection of 199Hg NMR, FEBS Lett., 117, 311–314,, 1980. 

Saurí, J., Parella, T., and Espinosa, J. F.: CLIP-HSQMBC: easy measurement of small proton–carbon coupling constants in organic molecules, Org. Biomol. Chem., 11, 4473,, 2013. 

Schätzlein, M. P., Becker, J., Schulze-Sünninghausen, D., Pineda-Lucena, A., Herance, J. R., and 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., 410, 2793–2804,, 2018. 

Schlippenbach, T. von, Oefner, P. J., and Gronwald, W.: Systematic Evaluation of Non-Uniform Sampling Parameters in the Targeted Analysis of Urine Metabolites by 1H,1H 2D NMR Spectroscopy, Sci. Rep., 8, 4249,, 2018. 

Tal, A. and Frydman, L.: Single-scan multidimensional magnetic resonance, Prog. Nucl. Mag. Res. Sp., 57, 241–292,, 2010. 

Tiainen, M., Maaheimo, H., Soininen, P., and Laatikainen, R.: 13 C isotope effects on 1 H chemical shifts: NMR spectral analysis of 13 C-labelled D-glucose and some 13 C-labelled amino acids: 13 C isotope effects on 1 H chemical shifts, Magn. Reson. Chem., 48, 117–122,, 2010. 

Timári, I., Kaltschnee, L., Raics, M. H., Roth, F., Bell, N. G. A., Adams, R. W., Nilsson, M., Uhrín, D., Morris, G. A., Thiele, C. M., and Kövér, K. E.: Real-time broadband proton-homodecoupled CLIP/CLAP-HSQC for automated measurement of heteronuclear one-bond coupling constants, RSC Adv., 6, 87848–87855,, 2016.  

Titman, J. J., Neuhaus, D., and Keeler, J.: Measurement of long-range heteronuclear coupling constants, J. Magn. Reson., 85, 111–131,, 1989. 

Watanabe, M., Meyer, K. A., Jackson, T. M., Schock, T. B., Johnson, W. E., and Bearden, D. W.: Application of NMR-based metabolomics for environmental assessment in the Great Lakes using zebra mussel (Dreissena polymorpha), Metabolomics, 11, 1302–1315,, 2015. 

Williamson, R. T., Márquez, B. L., Gerwick, W. H., and Kövér, K. E.: One- and two-dimensional gradient-selected HSQMBC NMR experiments for the efficient analysis of long-range heteronuclear coupling constants, 38, 265–273​​​​​​​, 2000. 

Wishart, D. S.: NMR metabolomics: A look ahead, J. Magn. Reson., 306, 155–161,, 2019. 

Wollborn, U. and Leibfritz, D.: Measurements of heteronuclear long-range coupling constants from inverse homonuclear 2D NMR spectra, J. Magn. Reson., 98, 142–146,, 1992. 

Zangger, K. and Sterk, H.: Homonuclear Broadband-Decoupled NMR Spectra, J. Magn. Reson., 124, 486–489,, 1997. 

Zhang, B., Powers, R., and O'Day, E. M.: Evaluation of Non-Uniform Sampling 2D 1H–13C HSQC Spectra for Semi-Quantitative Metabolomics, Metabolites, 10, 203,, 2020. 

Short summary
The HSQC experiment developed by Bodenhausen and Ruben is a cornerstone for modern NMR. When used in the field of metabolomics, the common practice of decoupling in the proton dimension limits the acquisition time and hence the resolution. Here, we present a virtual decoupling method to maintain both spectral simplicity and resolution, and demonstrate how it increases information content with the zebra mussel metabolome as an example.