Articles | Volume 4, issue 1
Research article
16 Feb 2023
Research article |  | 16 Feb 2023

Paramagnetic relaxivity of delocalized long-lived states of protons in chains of CH2 groups

Aiky Razanahoera, Anna Sonnefeld, Geoffrey Bodenhausen, and Kirill Sheberstov

Long-lived states (LLSs) have lifetimes TLLS that can be much longer than longitudinal relaxation times T1. In molecules containing several geminal pairs of protons in neighboring CH2 groups, it has been shown that delocalized LLSs can be excited by converting magnetization into imbalances between the populations of singlet and triplet states of each pair. Since the empirical yield of the conversion and reconversion of observable magnetization into LLSs and back is on the order of 10 % if one uses spin-lock induced crossing (SLIC), it would be desirable to boost the sensitivity by dissolution dynamic nuclear polarization (d-DNP). To enhance the magnetization of nuclear spins by d-DNP, the analytes must be mixed with radicals such as 4-hydroxy-2,2,6,6-tetramethylpiperidin-1-oxyl (TEMPOL). After dissolution, these radicals lead to an undesirable paramagnetic relaxation enhancement (PRE) which shortens not only the longitudinal relaxation times T1 but also the lifetimes TLLS of LLSs. It is shown in this work that PRE by TEMPOL is less deleterious for LLSs than for longitudinal magnetization for four different molecules: 2,2-dimethyl-2-silapentane-5-sulfonate (DSS), homotaurine, taurine, and acetylcholine. The relaxivities rLLS (i.e., the slopes of the relaxation rate constants RLLS as a function of the radical concentration) are 3 to 5 times smaller than the relaxivities r1 of longitudinal magnetization. Partial delocalization of the LLSs across neighboring CH2 groups may decrease this advantage, but in practice, this effect was observed to be small, for example, when comparing taurine containing two CH2 groups and homotaurine with three CH2 groups. Regardless of whether the LLSs are delocalized or not, it is shown that PRE should not be a major problem for experiments combining d-DNP and LLSs, provided the concentration of paramagnetic species after dissolution does not exceed 1 mM, a condition that is readily fulfilled in typical d-DNP experiments. In bullet d-DNP experiments however, it may be necessary to decrease the concentration of TEMPOL or to add ascorbate for chemical reduction.

1 Introduction

The lifetime of spin state populations in nuclear magnetic resonance (NMR) is normally limited by longitudinal relaxation. In certain cases, it is possible to access spin states that have extended lifetimes. In coupled pairs of spins with I=1/2, such long-lived states (LLSs) correspond to population imbalances between singlet and triplet states (Carravetta and Levitt, 2004; Carravetta et al., 2004) that are immune to intra-pair dipole–dipole interactions, which for pairs of protons are normally the dominant cause of longitudinal relaxation. In larger systems, LLSs may involve four, six, or more spins; all these states are weakly affected by dipolar relaxation (Hogben et al., 2011). The relaxation time constants TLLS can be much longer than typical longitudinal relaxation time constants T1. This feature is particularly useful for protein–ligand studies (Salvi et al., 2012; Buratto et al., 2014b, 2016). Applications of LLSs can be combined with different hyperpolarization methods, such as parahydrogen-based methods (Franzoni et al., 2012) or dissolution dynamic nuclear polarization (d-DNP) (Bornet et al., 2014; Kiryutin et al., 2019). The latter, d-DNP, is the most universal method to achieve high spin polarization, and it has found applications in drug screening (Lee et al., 2012; Buratto et al., 2014a; Kim et al., 2016) and in studies of metabolism by in vivo magnetic resonance imaging (MRI) (Nelson et al., 2013). Before dissolution, the saturation of the electron spin transitions by microwave irradiation of a solid sample near 1 K leads to an enhancement of the nuclear spin polarization by up to 4 orders of magnitude, compared to the thermal polarization at room temperature in the same magnetic field. The sample is then quickly dissolved and transferred to a solution-state NMR spectrometer, where the high-resolution spectrum is observed (Ardenkjær-Larsen et al., 2003). In an alternative approach known as “bullet DNP”, the cold solid sample is ejected from the polarizer and rapidly transferred to the NMR spectrometer where it is dissolved (Kouřil et al., 2019). After dissolution, the unpaired electrons of the dilute paramagnetic agent give rise to undesirable paramagnetic relaxation enhancement (PRE). For most molecules of interest, such as metabolites or potential drugs, proton relaxation is so fast that the level of hyperpolarization suffers during dissolution and transfer, which is one of the reasons why d-DNP is more often used for 13C or 15N rather than for protons. Although molecules that are in enriched 13C and 15N offer many possibilities for the excitation of LLSs (Feng et al., 2013; Elliott et al., 2019; Sheberstov et al., 2019), there are several drawbacks of using heteronuclei. Labeled compounds are expensive, and 13C or 15N observation is much less sensitive compared to 1H. After converting proton LLSs back into magnetization, only proton signals of interest are observed, while the background is suppressed. LLSs involving pairs of protons often provide good contrast as protons are often directly exposed to the drug–target interface. On the other hand, the relaxation rate constants of LLSs can be enhanced by mechanisms such as dipolar couplings to solvent nuclei, even with low gyromagnetic ratios, and to paramagnetic species (Kharkov et al., 2022).

Recently, it was discovered that LLSs involving geminal pairs of protons can be readily excited in many molecules containing at least two neighboring CH2 groups (Sonnefeld et al., 2022a, b). Aliphatic chains, which are the focus of this study, are commonly found in potential drugs, so LLSs of CH2 groups could provide a new tool for drug screening using NMR. Hyphenation of the LLS methodology with d-DNP offers promising perspectives, since at very low spin temperatures on the order of 10 mK that are routinely achieved in d-DNP singlet–triplet imbalances can result from a violation of the high-temperature approximation, so LLSs can be excited without any radio-frequency (RF) irradiation (Tayler et al., 2012; Bornet et al., 2014; Kress et al., 2019). LLSs that involve chemically equivalent proton pairs in CH2 groups need not be sustained by RF fields nor protected by shuttling to low fields. Therefore, one can transfer samples with hyperpolarized LLSs to an NMR spectrometer for detection without significant loss of polarization. For small molecules, the ratio TLLS/T1 ranges typically from 2 to 6 for LLSs in CH2 groups in non-degassed samples (Sonnefeld et al., 2022a), although it is possible to achieve a ratio TLLS/T1>30 in some degassed samples containing isolated pairs of protons (Sarkar et al., 2007) or carbon-13 nuclei (Pileio et al., 2012; Stevanato et al., 2015). In this work, we carried out a systematic analysis of relaxivities, i.e., of the dependence of the relaxation rate constants of LLSs and longitudinal magnetization on the concentration of the paramagnetic species 4-hydroxy-2,2,6,6-tetramethylpiperidin-1-oxyl (TEMPOL).

Figure 1(a) Chemical structures of four molecules supporting LLSs of CH2 groups studied in this work: 2,2-dimethyl-2-silapentane-5-sulfonate sodium salt (DSS, I), homotaurine (II), taurine (III), and acetylcholine (IV). CH2 groups supporting LLSs are numbered in each structure and highlighted by green circles. (b) Assignment of the 1H NMR spectrum of a mixture containing all four compounds.


Paramagnetic transition metal ions (Cu2+, Mn2+), lanthanides (Gd3+), and triplet oxygen (O2) have been shown to induce PRE of LLSs, although PRE is not very efficient as the fluctuating external fields at the sites of two closely spaced protons attached to the same carbon atom are strongly correlated (Tayler and Levitt, 2011). The effects of triplet oxygen on LLSs have been investigated in detail (Erriah and Elliott, 2019). The question arises if fluctuating external fields due to the bulky TEMPOL radical are more strongly correlated than for paramagnetic ions or oxygen, in particular when they act on delocalized LLSs involving several neighboring CH2 groups in the molecules shown in Fig. 1. In DSS (I) and homotaurine (II), the LLSs can be delocalized over all six protons of the three CH2 groups, whereas in taurine (III) and acetylcholine (IV) the LLSs always involve all four protons of both CH2 groups. Titration experiments with TEMPOL allowed us to determine to what extent the radical affects the LLS lifetimes and to determine whether it is necessary to quench the radicals after dissolution (Miéville et al., 2010). In low fields, in particular after dissolution during the transfer between the polarizer and the NMR magnet, PRE may be exacerbated by translational diffusion (Borah and Bryant, 1981) of the paramagnetic molecules relative to the analytes (Miéville et al., 2011).

2 Experimental methods

The delocalized LLSs were excited by using spin-lock induced crossing (SLIC) (DeVience et al., 2013) and its polychromatic extension (Sonnefeld et al., 2022b). A generic SLIC pulse sequence is illustrated in Fig. 2a. After a non-selective 90 pulse that rotates the magnetization into the transverse plane, one, two, or three continuous SLIC pulses with a common duration τSLIC are applied to the nuclei of interest, with a common RF amplitude (nutation frequency) ν1 that matches a multiple of the geminal intra-pair J-coupling, i.e., ν1=nJHHintra with n=1 for double-quantum (DQ) SLIC and n=2 for single-quantum (SQ) SLIC. Level anti-crossings (LACs) lead to a transfer of magnetization into LLSs, i.e., into a population imbalance between states with different permutation symmetry. Since pairs of protons in CH2 groups are chemically equivalent in achiral molecules (i.e., have the same chemical shifts), and, in the absence of couplings to heteronuclei, are often nearly magnetically equivalent, there is no need to suppress singlet-to-triplet leakage by transporting the sample into a region of low magnetic field nor by applying an RF field to sustain the imbalance. After allowing the LLSs to relax during a delay τrel, a T00 filter removes short-lived terms (Tayler and Levitt, 2013; Tayler, 2020), and a second SLIC pulse reconverts the remaining LLSs back into observable magnetization for detection. In this work, SLIC experiments with single, double, and triple irradiation (henceforth called single, double, and triple SLIC for simplicity) were carried out to determine TLLS, as shown by wavy arrows in Fig. 2b and c.

Figure 2(a) Generic pulse sequence for single and poly-SLIC where selective RF fields can be applied simultaneously to two or more CH2 groups. (b) Six possible poly-SLIC experiments applied to molecules containing three CH2 groups such as I and II of Fig. 1a. The upper row shows three experiments with irradiation at a single frequency for the creation of LLSs and a single readout pulse applied to the offset of the first, second, or third CH2 group; the lower row shows three experiments using triple irradiation of all three CH2 groups for LLS excitation, combined with a single readout SLIC applied to only one of the three CH2 groups. (c) Two schemes with double SLIC excitation and single SLIC readout for compounds containing only two CH2 groups such as III and IV of Fig. 1a.


Titrations were performed by preparing a set of samples where all compounds except TEMPOL had fixed concentrations. The volume of each sample was 600 µL. A stock solution with 40 mM of each compound was diluted by a factor of 4 to obtain a final concentration of 10 mM for each compound in D2O at pH 7.0 without removing paramagnetic oxygen by degassing. A stock solution of phosphate buffer (70 mM KH2PO4 and 130 mM K2HPO4) was prepared in D2O and diluted by a factor of 4. A 20 mM TEMPOL stock solution was diluted in steps and added to yield final concentrations of 0.5, 1.0, 2.0, 3.0, 4.0, and 6.0 mM. The 1H NMR spectra were obtained by adding 16 transients (for experiments with single SLIC irradiation) or 8 transients (for experiments with multiple SLIC irradiation) using a 500 MHz AVANCE Neo Bruker spectrometer with a 5 mm iProbe at 298 K. Each sample contained a mixture of all four molecules, thus ensuring accurate comparisons of relaxation rate constants of different molecules. The 1H NMR spectrum of the mixture with its assignments is presented in Fig. 1b. Typical signal decays due to LLS relaxation as a function of the TEMPOL concentration are shown in Fig. 3. The intensities of the LLS-derived signals are typically about 5 % for single SLIC experiments and up to 10 % for poly-SLIC experiments. The theoretical maximum efficiency of LLS excitation and reconversion in a four-spin –CH2–CH2– moiety was calculated to be 14 % for single SLIC and 28 % for double SLIC experiments (Sonnefeld et al., 2022a). Simulations of the contributions of different LLS terms to the observed signals were performed using SpinDynamica (Bengs and Levitt, 2018).

Figure 3Decays of LLS-derived signals of DSS (compound I) for different TEMPOL concentrations. The LLSs were excited and reconverted by irradiation with single SLIC pulses applied to CH2(1) with an RF amplitude of 27 Hz to match the condition for single-quantum level anti-crossing (SQ LAC). The solid lines correspond to mono-exponential fits, scaled to begin at 100 %.


3 Results and discussion

3.1 Comparison of relaxivities of long-lived states and of longitudinal magnetization: partly correlated random fields

As apparent in Fig. 4, both the longitudinal relaxation rate constant R1=1/T1 and the long-lived relaxation rate constant RLLS=1/TLLS depend linearly on the concentration of TEMPOL (in units of M or mol L−1):

(1) R 1 = R 1 0 + r 1 [ TEMPOL ] , R LLS = R LLS 0 + r LLS TEMPOL .

The slopes rLLS and r1 are known as relaxivities (in units of M−1 s−1); the intercepts R10 and RLLS0 are the rate constants determined in the absence of TEMPOL. Figure 4 shows that variations of R1 between neighboring CH2 groups within each molecule are much smaller than variations from one molecule to another. Whereas the T1 values of small molecules correlate with the molecular mass – the larger the molecule, the shorter T1 – this is not true for TLLS. In the absence of TEMPOL, the longest TLLS of ca. 15 s was observed for compound III, whereas the shortest TLLS of ca. 5 s was found for compound IV, although their T1 relaxation times and molecular masses are roughly the same, so their correlation times should be similar. The difference in TLLS may be explained by the presence of 12 methyl protons in compound IV, which cause faster relaxation of LLSs.

Wokaun and Ernst famously demonstrated that PRE is less efficient for relaxation of zero-quantum coherences than for single- and double-quantum coherences (Wokaun and Ernst, 1978). Tayler and Levitt demonstrated that a similar logic also applies to LLSs: whereas longitudinal relaxation is enhanced by fluctuations of external local fields induced by unpaired electrons of radicals, a LLS involving two spins I1 and I2 is only relaxed by fluctuating external fields if these are not correlated. In general, the extent of correlation of the two fluctuating fields at the locations of the two spins I1 and I2 can be characterized by the correlation coefficient C=B1B2/B1B2, where Bi=B1B2 is the mean (time-averaged) amplitude. Only the uncorrelated part of the two fluctuating fields given by B1-B22=B12+B22-2B1B2 contributes effectively to LLS relaxation (Tayler and Levitt, 2011). The smaller the radical, the closer it can approach one of the two geminal protons and hence the smaller the correlation coefficient C. It has been shown (Tayler and Levitt, 2011) that the ratio of relaxivities,

(2) κ = r LLS / r 1 ,

is a characteristic measure of the correlation coefficient C; the smaller κ, the larger C. The experimental ratio κ for the (chemically inequivalent) protons of the CH2 group in the (chiral) dipeptide alanine–glycine varied in the range 0.5<κ<0.3, depending on the size of the paramagnetic agent (Tayler and Levitt, 2011). A similar ratio κ=0.36 was observed for the CH2 group in the terminal glycine residue of the tripeptide Ala–Gly–Gly for PRE caused by triplet oxygen (Erriah and Elliott, 2019).

Figure 4Relaxation rate constants R1=1/T1 and RLLS=1/TLLS in CH2 groups of the four molecules I–IV as a function of the TEMPOL concentration. In (a) and (b), the LLSs were excited by triple SLIC, in (c) and (d) by double SLIC, both with an RF amplitude of 13.5 Hz to match the condition for double-quantum level anti-crossing (DQ LAC). In all cases, the LLSs were reconverted into magnetization by single SLIC applied to the CH2(2) group, except for compound IV, where two sets of experiments were performed with reconversion into magnetization of either CH2(1) or CH2(2) groups. The relaxivities r1 and rLLS correspond to the slopes of the linear regressions.


In CH2 chains with chemically equivalent pairs of protons in achiral molecules excited by exploiting magnetic inequivalence, the LLSs can be delocalized over several CH2 groups. Relaxation of a LLS localized within an individual CH2 group will contribute to the decay of a delocalized LLS, so one may expect the relaxivity of delocalized LLS to be more strongly affected by PRE than the relaxivity of a (hypothetical) localized LLS. We must however remain cautious, since one cannot assume that all CH2 groups are equally accessible to radicals. If some of the CH2 groups are more accessible, one may expect delocalized LLSs to have averaged relaxivities. As we shall discuss below, the variations in the observed relaxivities rLLS are not very large for different combinations of excitation and reconversion methods, and intramolecular variations are much smaller than differences between distinct compounds, so one can estimate an average ratio of relaxivities κ=r1/rLLS for all CH2 groups in a given molecule. Compounds I–IV feature average ratios κI〉≈0.22, κII〉≈0.23, κIII〉≈0.18, and κIV〉≈0.32 (see Table 1). To assess the effects of LLS delocalization on the relaxivity, it is useful to compare molecules II and III, as they differ by only one CH2 group. The mean relaxivity rLLS of compound II averaged over three CH2 groups is slightly higher than the relaxivity rLLS measured for compound III (see Table 1). This suggests that a higher degree of delocalization leads to a higher relaxivity. The LLS can be delocalized to a variable extent between all three CH2 groups in I and II, but they are always equally distributed between the two CH2 groups in compounds III and IV.

Table 1Experimentally determined relaxation rate constants (s−1) and relaxivities (M−1 s−1). Standard errors determined from linear regressions are shown in parentheses. For double SLIC, the RF amplitude was chosen to match the condition for double-quantum level anti-crossing (LAC), leading to different imbalances characterized by different rate constants RLLS(0)(SQ) with single SLIC excitation and single SLIC reconversion, and to rate constants RLLS(0)(DQ) with triple SLIC excitation and single SLIC reconversion.

Download Print Version | Download XLSX

3.2 Implications for dissolution DNP

Even though delocalized LLSs are less affected by TEMPOL than longitudinal magnetization, the observed decrease in TLLS is undesirable in the context of d-DNP. Since the use of TEMPOL or other polarizing agents is mandatory for d-DNP experiments, the question arises if it is worth scavenging TEMPOL after dissolution by addition of a reducing agent such as sodium ascorbate (vitamin C) to extend TLLS after dissolution (Miéville et al., 2010, 2011). Note that the preparation of samples comprising two types of beads is rather cumbersome, in particular for bullet DNP. According to Miéville et al. (2010, 2011), the rate of the reduction of TEMPOL by sodium ascorbate may be slow on the timescale of the transfer of the dissolved sample from the polarizer to the NMR magnet. Hence, the reaction may not be entirely completed by the time the sample arrives in the spectrometer. Scavenging by sodium ascorbate may be accelerated ca. 100 times if one uses Frémy's salt instead of TEMPOL (Negroni et al., 2022). Several alternative approaches have been developed to remove radicals once DNP has been achieved. One approach is to use radicals obtained by UV irradiation of frozen pyruvic acid. These radicals are quenched as soon as the temperature increases (Eichhorn et al., 2013). One may also use radicals grafted onto mesostructured silica materials (Gajan et al., 2014) or microporous polymers (Ji et al., 2017; El Daraï et al., 2021). However, the small relaxivities presented in Table 1 suggest that scavenging may not be necessary when using LLSs to preserve the hyperpolarization.

3.3 Experiments and simulations for molecules with three CH2 groups

It was shown (Sonnefeld et al., 2022b) that for the excitation of LLSs in systems with n=3 neighboring CH2 groups, i.e., with 2n=6 spins, there are seven orthogonal LLS product operators that can be created, with seven coefficients λi that depend on the excitation scheme:

(3) σ ^ LLS = - λ A A I ^ A I ^ A - λ M M I ^ M I ^ M - λ X X I ^ X I ^ X - λ A A M M I ^ A I ^ A I ^ M I ^ M - λ A A X X I ^ A I ^ A I ^ X I ^ X - λ M M X X I ^ M I ^ M I ^ X I ^ X - λ A A M M X X I ^ A I ^ A I ^ M I ^ M I ^ X I ^ X .

Here A and A denote the two protons of the CH2(1) group, M and M denote those of the middle CH2(2) group, and X and X correspond to the terminal CH2(3) group. This equation gives a general form of the density operator obtained after poly-SLIC, containing all long-lived terms found by numerical solution of the Liouville–von-Neumann equation. In addition to three bilinear terms, one encounters four higher terms that contain products of four and six spin operators. In principle, each term in Eq. (3) can decay with a different rate constant, so one could distinguish up to seven distinct rate constants RLLS(μ) with μ=AA, MM, XX, AAMM, AAXX, MMXX, and AAMMXX. Each term can be excited with a different amplitude and can contribute with a different weight to the observed signal.

In systems such as compounds III and IV with only two CH2 groups, only one LLS can be excited:

(4) σ ^ LLS = - λ A A I ^ A I ^ A - λ X X I ^ X I ^ X - λ A A X X I ^ A I ^ A I ^ X I ^ X .

The coefficients of the first two bilinear terms are always equal, i.e., λAA=λXX, while the four-spin term is always proportional to the leading bilinear terms, with a weight λAAXX=8/3λAA (Sonnefeld et al., 2022a). This state corresponds to the imbalance between the singlet–singlet state and the triplet–triplet manifold and is therefore expected to decay monoexponentially. In two sets of complementary experiments performed for compound IV, the experimental relaxation rate constants were indeed found to be indistinguishable, as can be seen by comparing the orange triangles and the green inverted triangles in Fig. 4d.

In compounds I and II however, which contain three adjacent CH2 groups, different SLIC excitation schemes lead to populate different LLSs, with different coefficients λLLS(μ) in Eq. (3). There are 9 different ways of exciting miscellaneous LLSs and 9 different ways of reconverting them, giving 81 possible experimental combinations. In order to investigate the relaxivities of these different LLSs which may have different decay rate constants RLLS(μ) and different relaxivities rLLS(μ), we performed six different poly-SLIC experiments with different SLIC pulses for excitation and reconversion, and we indeed found different LLS lifetimes (Fig. 5). Depending on the excitation and reconversion scheme used, there are pronounced differences between the relaxivities rLLS within one and the same molecule.

Figure 5Decay rate constants RLLS=1/TLLS of long-lived states in CH2 groups in (a) DSS (I) and (b) homotaurine (II), each containing three CH2 groups, as a function of the TEMPOL concentration. Six different poly-SLIC experiments with distinct excitation and reconversion methods were performed for each molecule, as indicated by wavy arrows. The relaxivities rLLS correspond to the slopes of the linear regressions.


Figure 6Calculated contributions of the seven different LLS terms P^μ (the 3 two-spin terms are shown in the same color) in the density operator of Eq. (3) to the observed signals for all six different single and poly-SLIC experiments used in this work to determine the relaxivities rLLS(a) in the six-spin systems of DSS (I) and homotaurine (II). The histograms show the products, λμMLLSλ̃μLLSM, of the coefficients of LLS excitation and reconversion methods. The normalization ensures that the sum of all products of coefficients is equal to 1. Experiments with triple SLIC excitation and single SLIC reconversion applied to the middle CH2 group (e) provide LLSs that are almost evenly distributed among all three CH2 groups, whereas the other experiments provide access to LLSs that are in part localized on the group where the reconversion SLIC pulse is applied. The excitation and reconversion of the (yellow) six-spin term is negligible except for case (b).


We calculated the contributions of each of the seven terms to the observable LLS-derived signals, after two consecutive transformations: I^zinσ^LLSI^xobs (see Fig. 6). For each excitation scheme used in this work, we considered all seven coefficients λμMLLS corresponding to the seven terms in Eq. (3), as well as all seven reconversion coefficients λ̃μLLSM. The coefficients were calculated as follows:

(5) λ μ M LLS ( I ^ z in σ ^ LLS ) = Tr P ^ μ σ ^ LLS Tr P ^ μ P ^ μ , λ ̃ μ LLS M ( P ^ µ I ^ x , µ obs ) = Tr I ^ x I ^ x , μ obs Tr I ^ x I ^ x ,

where Tr stands for trace and where the index μ corresponds to one of the seven LLS terms in Eq. (3), the operator P^μ represents the μth LLS term, I^zin is the initial magnetization of the excited spins, I^x is the transverse magnetization of the observed spins after reconversion, and I^x,μobs is the transverse magnetization obtained after reconversion of only the μth term P^μ instead of the full σ^LLS. The observed signal Sμ stemming from the μth term is determined by the product of two coefficients, λμMLLS and λ̃μLLSM, for a given combination of excitation and reconversion SLIC pulses. These contributions are shown in Fig. 6. The sum of all seven amplitudes for each panel in Fig. 6 was normalized to one. These graphs show how the LLSs are delocalized across spin systems comprising n=3 neighboring CH2 groups. We only consider coherent spin dynamics during excitation and reconversion, neglecting possible redistributions of LLSs due to Overhauser-type cross-relaxation effects and neglecting zero-quantum coherences.

Note that a single SLIC pulse applied at the chemical shift of any of the three CH2 groups results in the excitation of a delocalized state, which is predominantly (but not exclusively) associated with the irradiated pair. By using triple SLIC excitation and single SLIC reconversion applied to the middle CH2 group, one can excite a fairly even distribution of LLSs involving all 2n=6 coupled spins. For compound II, the most strongly delocalized state features the largest relaxivity rLLS. For compound I, however, the largest relaxivities were obtained for experiments where the largest contribution to the observed signal came from the terminal group CH2(3) that is closest to the trimethylsilane group. This group has also the largest longitudinal relaxivity r1, as can be seen in Fig. 4a. Detailed calculations of the relaxation superoperator might help to rationalize the experimental results obtained here.

4 Conclusions

The relaxation rate constants of various long-lived states and of the longitudinal magnetization of DSS, homotaurine, taurine, and acetylcholine were measured as a function of the concentration of the radical TEMPOL. In all cases, the relaxivities rLLS are lower by about a factor of 3 compared to the relaxivities r1. This implies that the effects of paramagnetic relaxation enhancement on LLSs due to TEMPOL during sample transfer in dissolution DNP should not be too severe. Furthermore, the LLS relaxivity was studied for different SLIC excitation and reconversion schemes. The results support simulations that show that different LLSs are excited depending on the SLIC sequence and the number of adjacent methylene units. SLIC methods have also been shown to be efficient for other achiral molecules containing neighboring CH2 groups, such as dopamine, γ-aminobutyric acid (GABA), ethanolamine, and β-alanine (Sonnefeld et al., 2022a). All of these molecules contain aliphatic chains, so the effects of paramagnetic polarizing agents like TEMPOL should be similar to what is reported in this work.

Data availability

All original NMR data obtained for this paper are available through the Zenodo repository under (Razanahoera et al., 2022).

Author contributions

Conceptualization: KS and GB. Data collection: AR and KS. Data analysis: AR and KS. Spin dynamic computation: AS and KS. Visualization: AR, AS, and KS. Draft manuscript preparation: AR, AS, GB, and KS. Supervision: GB. Funding acquisition: GB. All authors reviewed the results and approved the final version of the article.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Magnetic Resonance. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


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


We are indebted to the CNRS (Centre National de la Recherche Scientifique) and the ENS (École Normale Supérieure) for support.

Financial support

This research has been supported by the European Research Council (ERC) for the Synergy grant “Highly Informative Drug Screening by Overcoming NMR Restrictions” (HISCORE, grant agreement no. 951459).

Review statement

This paper was edited by Gottfried Otting and reviewed by Alexej Jerschow, Malcolm Levitt, and one anonymous referee.


Ardenkjær-Larsen, J. H., Fridlund, B., Gram, A., Hansson, G., Hansson, L., Lerche, M. H., Servin, R., Thaning, M., and Golman, K.: Increase in signal-to-noise ratio of > 10,000 times in liquid-state NMR, P. Natl. Acad. Sci. USA, 100, 10158–10163,, 2003. 

Bengs, C. and Levitt, M. H.: SpinDynamica: Symbolic and numerical magnetic resonance in a Mathematica environment, Magn. Reson. Chem., 56, 374–414,, 2018. 

Borah, B. and Bryant, R. G.: NMR relaxation dispersion in an aqueous nitroxide system, J. Chem. Phys., 75, 3297–3300,, 1981. 

Bornet, A., Ji, X., Mammoli, D., Vuichoud, B., Milani, J., Bodenhausen, G., and Jannin, S.: Long-Lived States of Magnetically Equivalent Spins Populated by Dissolution-DNP and Revealed by Enzymatic Reactions, Eur. J. Chem., 20, 17113–17118,, 2014. 

Buratto, R., Bornet, A., Milani, J., Mammoli, D., Vuichoud, B., Salvi, N., Singh, M., Laguerre, A., Passemard, S., Gerber-Lemaire, S., Jannin, S., and Bodenhausen, G.: Drug Screening Boosted by Hyperpolarized Long-Lived States in NMR, ChemMedChem, 9, 2509–2515,, 2014a. 

Buratto, R., Mammoli, D., Chiarparin, E., Williams, G., and Bodenhausen, G.: Exploring Weak Ligand–Protein Interactions by Long-Lived NMR States: Improved Contrast in Fragment-Based Drug Screening, Angew. Chem. Int. Edit., 53, 11376–11380,, 2014b. 

Buratto, R., Mammoli, D., Canet, E., and Bodenhausen, G.: Ligand–Protein Affinity Studies Using Long-Lived States of Fluorine-19 Nuclei, J. Med. Chem., 59, 1960–1966,, 2016. 

Carravetta, M. and Levitt, M. H.: Long-Lived Nuclear Spin States in High-Field Solution NMR, J. Am. Chem. Soc., 126, 6228–6229,, 2004. 

Carravetta, M., Johannessen, O. G., and Levitt, M. H.: Beyond the T1 Limit: Singlet Nuclear Spin States in Low Magnetic Fields, Phys. Rev. Lett., 92, 153003,, 2004. 

DeVience, S. J., Walsworth, R. L., and Rosen, M. S.: Preparation of Nuclear Spin Singlet States Using Spin-Lock Induced Crossing, Phys. Rev. Lett., 111, 173002,, 2013. 

Eichhorn, T. R., Takado, Y., Salameh, N., Capozzi, A., Cheng, T., Hyacinthe, J.-N., Mishkovsky, M., Roussel, C., and Comment, A.: Hyperpolarization without persistent radicals for in vivo real-time metabolic imaging, P. Natl. Acad. Sci. USA, 110, 18064–18069,, 2013. 

El Daraï, T., Cousin, S. F., Stern, Q., Ceillier, M., Kempf, J., Eshchenko, D., Melzi, R., Schnell, M., Gremillard, L., Bornet, A., Milani, J., Vuichoud, B., Cala, O., Montarnal, D., and Jannin, S.: Porous functionalized polymers enable generating and transporting hyperpolarized mixtures of metabolites, Nat. Commun., 12, 4695,, 2021. 

Elliott, S. J., Kadeøávek, P., Brown, L. J., Sabba, M., Glöggler, S., O'Leary, D. J., Brown, R. C. D., Ferrage, F., and Levitt, M. H.: Field-cycling long-lived-state NMR of 15N2 spin pairs, Mol. Phys., 117, 861–867,, 2019. 

Erriah, B. and Elliott, S. J.: Experimental evidence for the role of paramagnetic oxygen concentration on the decay of long-lived nuclear spin order, RSC Adv., 9, 23418–23424,, 2019. 

Feng, Y., Theis, T., Liang, X., Wang, Q., Zhou, P., and Warren, W. S.: Storage of Hydrogen Spin Polarization in Long-Lived 13C2 Singlet Order and Implications for Hyperpolarized Magnetic Resonance Imaging, J. Am. Chem. Soc., 135, 9632–9635,, 2013. 

Franzoni, M. B., Buljubasich, L., Spiess, H. W., and Münnemann, K.: Long-Lived 1H Singlet Spin States Originating from Para-Hydrogen in Cs-Symmetric Molecules Stored for Minutes in High Magnetic Fields, J. Am. Chem. Soc., 134, 10393–10396,, 2012. 

Gajan, D., Bornet, A., Vuichoud, B., Milani, J., Melzi, R., van Kalkeren, H. A., Veyre, L., Thieuleux, C., Conley, M. P., Grüning, W. R., Schwarzwälder, M., Lesage, A., Copéret, C., Bodenhausen, G., Emsley, L., and Jannin, S.: Hybrid polarizing solids for pure hyperpolarized liquids through dissolution dynamic nuclear polarization, P. Natl. Acad. Sci. USA, 111, 14693–14697,, 2014. 

Hogben, H. J., Hore, P. J., and Kuprov, I.: Multiple decoherence-free states in multi-spin systems, J. Magn. Reson., 211, 217–220,, 2011. 

Ji, X., Bornet, A., Vuichoud, B., Milani, J., Gajan, D., Rossini, A. J., Emsley, L., Bodenhausen, G., and Jannin, S.: Transportable hyperpolarized metabolites, Nat. Commun., 8, 13975,, 2017. 

Kharkov, B., Duan, X., Rantaharju, J., Sabba, M., Levitt, M. H., Canary, J. W., and Jerschow, A.: Weak nuclear spin singlet relaxation mechanisms revealed by experiment and computation, Phys. Chem. Chem. Phys., 24, 7531–7538,, 2022. 

Kim, Y., Liu, M., and Hilty, C.: Parallelized Ligand Screening Using Dissolution Dynamic Nuclear Polarization, Anal. Chem., 88, 11178–11183,, 2016. 

Kiryutin, A. S., Rodin, B. A., Yurkovskaya, A. V., Ivanov, K. L., Kurzbach, D., Jannin, S., Guarin, D., Abergel, D., and Bodenhausen, G.: Transport of hyperpolarized samples in dissolution-DNP experiments, Phys. Chem. Chem. Phys., 21, 13696–13705,, 2019. 

Kouøil, K., Kouøilová, H., Bartram, S., Levitt, M. H., and Meier, B.: Scalable dissolution-dynamic nuclear polarization with rapid transfer of a polarized solid, Nat. Commun., 10, 1733,, 2019. 

Kress, T., Walrant, A., Bodenhausen, G., and Kurzbach, D.: Long-Lived States in Hyperpolarized Deuterated Methyl Groups Reveal Weak Binding of Small Molecules to Proteins, J. Phys. Chem. Lett., 10, 1523–1529,, 2019. 

Lee, Y., Zeng, H., Ruedisser, S., Gossert, A. D., and Hilty, C.: Nuclear Magnetic Resonance of Hyperpolarized Fluorine for Characterization of Protein–Ligand Interactions, J. Am. Chem. Soc., 134, 17448–17451,, 2012. 

Miéville, P., Ahuja, P., Sarkar, R., Jannin, S., Vasos, P. R., Gerber-Lemaire, S., Mishkovsky, M., Comment, A., Gruetter, R., Ouari, O., Tordo, P., and Bodenhausen, G.: Scavenging Free Radicals To Preserve Enhancement and Extend Relaxation Times in NMR using Dynamic Nuclear Polarization, Angew. Chem. Int. Edit., 49, 6182–6185,, 2010. 

Miéville, P., Jannin, S., and Bodenhausen, G.: Relaxometry of insensitive nuclei: Optimizing dissolution dynamic nuclear polarization, J. Magn. Reson., 210, 137–140,, 2011. 

Negroni, M., Turhan, E., Kress, T., Ceillier, M., Jannin, S., and Kurzbach, D.: Frémy's Salt as a Low-Persistence Hyperpolarization Agent: Efficient Dynamic Nuclear Polarization Plus Rapid Radical Scavenging, J. Am. Chem. Soc., 144, 20680–20686,, 2022. 

Nelson, S. J., Kurhanewicz, J., Vigneron, D. B., Larson, P. E. Z., Harzstark, A. L., Ferrone, M., van Criekinge, M., Chang, J. W., Bok, R., Park, I., Reed, G., Carvajal, L., Small, E. J., Munster, P., Weinberg, V. K., Ardenkjaer-Larsen, J. H., Chen, A. P., Hurd, R. E., Odegardstuen, L.-I., Robb, F. J., Tropp, J., and Murray, J. A.: Metabolic Imaging of Patients with Prostate Cancer Using Hyperpolarized [1-13C]Pyruvate, Sci. Transl. Med., 5, 198ra108,, 2013. 

Pileio, G., Hill-Cousins, J. T., Mitchell, S., Kuprov, I., Brown, L. J., Brown, R. C. D., and Levitt, M. H.: Long-Lived Nuclear Singlet Order in Near-Equivalent 13C Spin Pairs, J. Am. Chem. Soc., 134, 17494–17497,, 2012. 

Razanahoera, A., Sonnefeld, A., Bodenhausen, G., and Sheberstov, K.: Paramagnetic relaxivity of delocalized long-lived states of protons in chains of CH2 groups, Zenodo [data set],, 2022. 

Salvi, N., Buratto, R., Bornet, A., Ulzega, S., Rentero Rebollo, I., Angelini, A., Heinis, C., and Bodenhausen, G.: Boosting the Sensitivity of Ligand–Protein Screening by NMR of Long-Lived States, J. Am. Chem. Soc., 134, 11076–11079,, 2012. 

Sarkar, R., Vasos, P. R., and Bodenhausen, G.: Singlet-State Exchange NMR Spectroscopy for the Study of Very Slow Dynamic Processes, J. Am. Chem. Soc., 129, 328–334,, 2007. 

Sheberstov, K. F., Vieth, H.-M., Zimmermann, H., Rodin, B. A., Ivanov, K. L., Kiryutin, A. S., and Yurkovskaya, A. V.: Generating and sustaining long-lived spin states in 15N, 15N'-azobenzene, Sci. Rep., 9, 20161,, 2019. 

Sonnefeld, A., Razanahoera, A., Pelupessy, P., Bodenhausen, G., and Sheberstov, K.: Long-lived states of methylene protons in achiral molecules, Sci. Adv., 8, eade2113,, 2022a. 

Sonnefeld, A., Bodenhausen, G., and Sheberstov, K.: Polychromatic Excitation of Delocalized Long-Lived Proton Spin States in Aliphatic Chains, Phys. Rev. Lett., 129, 183203,, 2022b. 

Stevanato, G., Hill-Cousins, J. T., Håkansson, P., Roy, S. S., Brown, L. J., Brown, R. C. D., Pileio, G., and Levitt, M. H.: A Nuclear Singlet Lifetime of More than One Hour in Room-Temperature Solution, Angew. Chem. Int. Edit., 54, 3740–3743,, 2015. 

Tayler, M. C. D.: Chapter 10: Filters for Long-lived Spin Order, in: Long-lived Nuclear Spin Order, edited by: Pileio, G., RSC publishing, 188–208, (last access: 14 February 2023), 2020.  

Tayler, M. C. D. and Levitt, M. H.: Paramagnetic relaxation of nuclear singlet states, Phys. Chem. Chem. Phys., 13, 9128–9130,, 2011. 

Tayler, M. C. D. and Levitt, M. H.: Accessing Long-Lived Nuclear Spin Order by Isotope-Induced Symmetry Breaking, J. Am. Chem. Soc., 135, 2120–2123,, 2013. 

Tayler, M. C. D., Marco-Rius, I., Kettunen, M. I., Brindle, K. M., Levitt, M. H., and Pileio, G.: Direct Enhancement of Nuclear Singlet Order by Dynamic Nuclear Polarization, J. Am. Chem. Soc., 134, 7668–7671,, 2012. 

Wokaun, A. and Ernst, R. R.: The use of multiple quantum transitions for relaxation studies in coupled spin systems, Mol. Phys., 36, 317–341,, 1978. 

Short summary
In this work, we study the relaxivity of long-lived states (LLSs) excited across two to three CH2 groups in four different compounds and find that a commonly used polarising agent in dissolution dynamic nuclear polarization (d-DNP) does not enhance relaxation significantly. This result is important for the future implementation of d-DNP hyperpolarization of LLS.