Long-lived states of magnetically equivalent spins populated by dissolution-DNP and revealed by enzymatic reactions.

Hyperpolarization by dissolution dynamic nuclear polarization (D-DNP) offers a way of enhancing NMR signals by up to five orders of magnitude in metabolites and other small molecules. Nevertheless, the lifetime of hyperpolarization is inexorably limited, as it decays toward thermal equilibrium with the nuclear spin-lattice relaxation time. This lifetime can be extended by storing the hyperpolarization in the form of long-lived states (LLS) that are immune to most dominant relaxation mechanisms. Levitt and co-workers have shown how LLS can be prepared for a pair of inequivalent spins by D-DNP. Here, we demonstrate that this approach can also be applied to magnetically equivalent pairs of spins such as the two protons of fumarate, which can have very long LLS lifetimes. As in the case of para-hydrogen, these hyperpolarized equivalent LLS (HELLS) are not magnetically active. However, a chemical reaction such as the enzymatic conversion of fumarate into malate can break the magnetic equivalence and reveal intense NMR signals.

Introduction

Spin hyperpolarization by dissolution dynamic nuclear polarization (d-DNP) has become a major area of research in NMR. This emerging method provides a way of boosting the sensitivity of NMR experiments by enhancing the intrinsically low nuclear spin polarization dictated by Boltzmann’s law. Thus, 13C NMR signals of small molecules have been enhanced by up to four orders of magnitude.1 In a typical d-DNP experiment, the frozen sample is initially polarized at low temperatures and moderate fields, and the signals are subsequently measured in solution in a separate detection apparatus operating at room temperature. This implies that the polarized frozen sample needs to be dissolved and transferred rapidly. This transfer, sometimes poetically called voyage, can be performed either manually or by means of a pneumatic system. During the voyage, the hyperpolarized molecules experience low magnetic fields (sometimes as low as the earth’s field, or even lower), which have detrimental effects on the enhanced polarization. This is one of the reasons why D-DNP has been most useful for nuclear spins with long T1, such as the isolated low-gamma quaternary 13C spin in 1-13C pyruvic acid.2 On the other hand, apart from some exotic experiments,3 1H spins have hardly been exploited by d-DNP, as their short T1 relaxation times mean that the hyperpolarization is driven back rapidly toward Boltzmann equilibrium. However, we have shown recently that 1H can be polarized very efficiently and rapidly up to P(1H)=91 % with a buildup time constant as short as τDNP(1H)=150 s at B0=6.7 T and T=1.2 K.4 One possible strategy for taking advantage of this large 1H hyperpolarization consists in storing the magnetization in the form of long-lived states (LLS)5 with extended lifetimes.

In recent years, several successful studies combining d-DNP with LLS6 have shown that hyperpolarized magnetization can be converted into LLS with extended lifetimes TLLS≫T1. In a pair of equivalent spins 1/2, the singlet state S0=(|αβ〉-|βα〉)/√2 is largely disconnected from the triplet states T+1=|αα〉, T0=(|αβ〉+|βα〉)/√2 and T−1=|ββ〉 because relaxation mechanisms that are symmetric with respect to spin exchange (such as the dipole–dipole interaction between the two spins) cannot induce singlet–triplet transitions.7 Therefore, if a triplet–singlet population imbalance (TSI) is prepared by any means, it is likely to be long-lived. We use the expression TSI in analogy to the A/E imbalance (AEI) recently described for methyl groups by Benno Meier et al.8 Both TSI and AEI refer to a difference between the average populations of spin states belonging to different irreducible representations of the spin permutation group, that is, ΓA and ΓE in methyl groups and Γg and Γu (or triplet and singlet states) in pairs of equivalent spins. The excitation and detection of an LLS involving a pair of equivalent spins is challenging because the magnetic equivalence needs to be lifted during both excitation and detection but preserved during storage. In this context, two possible scenarios are: 1) in most experiments described so far, the symmetry is imposed on an otherwise inequivalent two-spin system during the storage period only, or 2) in this work, the symmetry of an inherently equivalent two-spin system is broken during both excitation and detection.

Para-hydrogen9 offers the best example of nuclear singlet order in a molecule with two equivalent spins. The singlet state of H2 can be produced at low temperatures (typically 40 K) in the presence of a paramagnetic catalyst, which allows singlet–triplet interconversion by lifting the symmetry of H2 near the catalytic surface. The singlet spin state of H2 has the lowest energy, primarily determined by the quantization of its rotational state, and therefore, is predominantly populated at low temperatures. This leads to the creation of a large TSI compared to H2 in Boltzmann equilibrium at room temperature. Para-H2 is not magnetically active, and therefore cannot be observed directly by NMR, but it can be converted into observable signals through an asymmetric hydrogenation reaction by which the two protons stemming from para-H2 become inequivalent. On the other hand, a symmetric hydrogenation process can generate a molecule in which the equivalence of the two protons is preserved. This is the case, for example, for the hydrogenation of acetylene to form ethylene,10 or of dimethyl acetylene dicarboxylate to produce dimethyl maleate.11 The preserved para state can subsequently be rendered accessible to NMR by another chemical reaction that lifts the symmetry of the molecule.

If one starts with an inequivalent two-spin system, a precursor state, that is, a state that acquires a long-lived property as soon as the two spins are made equivalent during the storage interval, can be prepared by using suitable rf pulse sequences,12 by adiabatic transport to low fields,13 or by chemical reactions.6b Alternatively, a compromise can be found by using systems containing nearly equivalent spins14 in which the singlet and triplet states are only weakly mixed, but with an admixture that can be augmented by suitable pulse sequences to induce a singlet-to-magnetization (S2M) conversion.

Most experiments in which d-DNP is combined with LLS6a–e rely on rf pulse sequences to prepare the LLS, usually after the transfer of the hyperpolarized sample to the detection magnet. As a result, extensive relaxation occurs during the transfer. However, Tayler and co-workers6f have shown that LLS order can be populated directly before the transfer by d-DNP for the two inequivalent 13C spins in 1,2-13C2-pyruvic acid. They pointed out that the polarization of the singlet state Ps (=PTSI, vide supra) is proportional to the square of the spin Zeeman polarization P (i.e., PTSI=Ps=−1/3P2). Therefore, provided a high spin polarization can be reached by d-DNP, say PZ=50 %, a significant amount of singlet order, in this example P=8.33 %, can be created directly without any rf pulses. In the case where PZ=91 % can be attained, one obtains P=28 %. Such high levels of polarization can indeed be prepared directly by DNP for 1H at B0=6.7 T and T=1.2 K and indirectly for 13C or other nuclei through cross polarization from protons.4

In this work, we demonstrate that a TSI can be efficiently populated by d-DNP for the pair of magnetically equivalent 1H spins in fumarate, and that this is preserved in the liquid state after dissolution for a long time TTSI, which was estimated to be of the order of 50 s. We refer to this type of LLS as hyperpolarized equivalent long-lived states (HELLS). We show how HELLS can be readily “revealed” by allowing fumarate to undergo a biologically relevant enzymatic conversion into malate.

Experiments

For efficient d-DNP, the samples usually consist of frozen glassy solids containing typically 10–50 mM polarizing agents such as TEMPOL in addition to the molecules of interest. In our experiment, the molecule of interest shall possess two spins I and S that are magnetically equivalent in the liquid phase, but inequivalent in the frozen state and in moderate magnetic fields because they are exposed to slightly different environments and therefore experience different chemical shifts because of chemical shift anisotropies (CSAs) and different internuclear and electron–nuclear dipolar couplings. Given that freezing to low temperatures lifts the equivalence, the energy levels are better expressed in the product basis (PB). At T=1.2 K and B0=6.7 T, the proton Boltzmann polarization without DNP is PZ=0.57 %. Therefore, the deviations of diagonal elements from the demagnetized state Δσ=σ−E will be (Δnαα, Δnαβ, Δnβα, Δnββ)=1/4(2P+PZ2, −PZ2, −PZ2, −2PZ+PZ2)=(0.003, 0, 0, −0.003). Assuming, for simplicity, that DNP could confer a Zeeman polarization PZ=100 %, only the lowest energy level |αα〉 would be populated by hyperpolarization, so that (Δnαα, Δnαβ, Δnβα, Δnββ)=(0.75, −0.25, −0.25, −0.25) (See Figure 1 a: TSI Preparation). As soon as the polarized sample is heated and dissolved to the liquid state, CSAs and dipolar couplings are averaged out, so that the spins I and S become magnetically equivalent. The density operator can therefore better be expressed in the singlet–triplet basis (STB). In our case, as nα=nβα, (and hence Δnα=Δnβα), it is easily seen that σ(PB)=σ(STB) [and hence Δσ(PB)=Δσ(STB)], with the following diagonal elements: (Δnαα, Δnαβ, Δnβα, Δnββ)=(ΔnT+1, ΔnT0, ΔnS0, ΔnT−1)=1/4(2P+PZ2, −PZ2, −PZ2, −2PZ+PZ2). Hence, PTSI=ΔnS0−1/3(ΔnT+1+ΔnT0+ΔnT−1)=−Pz2/3. The TSI will thus result from the depletion of nαβ and nβα by hyperpolarization (Figure 1 b). The spins are equivalent, so the TSI can be stored indifferently in a low or high magnetic field (in a magnetic path for example). During the storage period, the populations of the three triplet states will equilibrate, that is, the deviations of the population of the three triplet levels will average out to give (ΔnT+1)′=(ΔnT0)′=(ΔnT−1)′=1/3(ΔnT+1+ΔnT0+ΔnT−1)=1/4PZ2/3. The singlet should not be affected by dipole–dipole relaxation, so the TSI in principle remains equal to PTSI=−PZ2/3. (See Figure 1 b: TSI storage). The sample is then transferred to the NMR or MRI magnet for detection. The system of two equivalent spins can then be transformed (chemically or enzymatically) into a system of two inequivalent spins, so that the “sealed” hyperpolarization can be “revealed” by conversion into observable magnetization. If the reaction is fast and goes to completion, one can convert Δσ from the STB back to the PB by using a suitable base transformation (see Ref. 15). Δσ(PB) resulting from this transformation can be expressed as a superposition of longitudinal two-spin order and zero-quantum coherence since Δσ(PB)=Pz2/6 (2IzSz+2ZQx). (See Figure 1 c: TSI revelation).

Figure 1

Schematic deviations of the populations from the fully saturated state Δσ=σ−E among the energy levels of the two protons of fumarate a) in the polarizer at 6.7 T and 1.2 K without DNP at Boltzmann equilibrium (PZ(1H)=0.57 %) and after DNP polarization to the theoretical limit PZ(1H)=100 %, b) during the transfer, which may go through low magnetic fields or through a magnetic tunnel to sustain a higher field, and c) in the detection magnet, typically at 7 T and 300 K, where the spins are made inequivalent by an enzymatic conversion. In each scheme, the deviations of the diagonal elements from the demagnetized state Δσ=σ−E are given as a function of the polarization PZ. In (c), the full density matrix is given to show the off-diagonal elements.

Enzymatic reactions are not instantaneous, and do not necessarily lead to complete conversion into the product. Figure 2 shows an example of the conversion of fumarate into malate by fumarase under conditions that can be combined with d-DNP. The steady-state concentrations are only reached after 25 min. This has important implications for our experiment. In fact, a highly polarized state Δσ=2IzSz+2ZQx is indeed produced instantaneously in malate whenever fumarate molecules carrying a TSI undergo an enzymatic conversion, but the ZQx term immediately starts evolving under the difference of chemical shifts, and therefore rapidly dephases and averages to zero as the reaction proceeds. Furthermore, the hyperpolarized TSI of fumarate, once it is transferred to malate, will tend to relax to thermal Boltzmann equilibrium.

Figure 2

Enzymatic conversion of fumarate into malate by fumarase at 300 K monitored by integration of the conventional 1H NMR signals of the two species. The nuclear polarizations are in thermal Boltzmann equilibrium, without resorting to DNP. A solution of fumarase (4 μL, 5.8 mg mL−1, i.e., 10 units) was injected into a fumarate solution (500 μL, 50 mM) at pH 8 in a buffer of 25 mM TRIS and 200 mM NaCl.

It is, however, possible to “sustain” the LLS of malate by so-called “high-field” methods,7c, 12, 15b for example, by applying an rf irradiation halfway between the two chemical shifts (either continuous-wave (CW), or, if desired, by applying a WALTZ-16 pulse train),16 thus preserving the full Δσ=2IzSz+2ZQx state. This strategy allows one to slow down relaxation of 2IzSz and prevent dephasing of ZQx. For the two inequivalent protons in malate, we thus determined TLLS=6 s at B0=7 T and T=298 K. Moreover, the use of WALTZ-16 pulse trains has the advantage of wiping out any single-quantum magnetization that would not arise from HELLS. A conventional LLS detection sequence, for example, the second half of the “Sarkar sequence”15b(Figure 3), can then be used to transform Δσ=2IzSz+2ZQx into observable magnetization.

Figure 3

Timing of a HELLS experiment. After dissolution, the hyperpolarized solution containing fumarate that carries the TSI is transferred to a holding chamber just above the NMR tube to determine its lifetime TTSI during a variable preinjection delay τTSI. The fumarate solution is then injected into a solution containing fumarase to start the conversion of fumarate into malate, accompanied by a conversion of the TSI on fumarate into an LLS on malate. A WALTZ-16 pulse train is applied during the delay τLLS with the carrier halfway between the chemical shifts of the two protons of malate to make these two protons effectively equivalent. The remainder of the pulse sequence is identical to the second half of the “Sarkar sequence”.15b The conversion of the LLS into observable magnetization is most efficient when τ1=1/(4 J)and τ2=1/(2Δν)(J=10.4 Hz and Δν=960 Hz at 300 MHz). The detection scheme can be repeated n times, bearing in mind that the LLS on malate is replenished during each sustaining interval τLLS by enzymatic conversion of fumarate that carries a slowly relaxing TSI.

The lifetime of the LLS of malate (TLLSM=6 s at 300 MHz if the rf amplitude of the CW field is ν1=3 kHz) is short compared to the enzymatic transformation, so the time τLLS (see Figure 3) allocated for the LLS to accumulate in malate before it is converted into observable signals needs to be optimized carefully. The concentrations [F] and [M] of fumarate and malate can be described by pseudo first-order kinetics as shown in Equation (1), in which [F](t) and [M](t) are the concentrations of fumarate and malate, kFM and kMF are the apparent kinetic constants of the overall enzymatic conversion of fumarate into malate and vice versa, without considering the details of the Michaelis–Menten mechanism.(1)

The temporal evolution of the expectation value PLLSM in malate arising from the conversion of fumarate can be obtained by solving numerically the rate equations [Eq. (2)], in which PTSIF and PLLSM are the expectation values of the TSI in fumarate and of the LLS in malate, and RTSIF and RLLSM are their relaxation rates.(2)

The “apparent” rate constants kFM and kMF can be obtained by fitting the signal amplitudes in Figure 2 to the rate equations in Equation (1). One can then calculate the temporal evolution of PTSIF in fumarate in the presence of ten units of enzyme, as well as PLLSM of malate obtained by the conversion of the TSI of fumarate into an LLS of malate that relaxes with TLLSM (Figure 4). These curves were obtained by assuming that TTSIF=60 s for fumarate (on the basis of preliminary observations as discussed below), and using the experimentally determined time constant TLLSM=6 s for malate. According to Figure 4, the optimal delay to maximize the conversion of the TSI of fumarate into the LLS of malate is 10 s. Thus, one should wait τLLS=10 s while sustaining the LLS by a suitable rf field before attempting to convert the LLS of malate into observable magnetization. The alternation of rf irradiation and signal observation can be repeated n times. During each interval τLLS, the LLS on malate will be replenished by the enzymatic conversion of the slowly relaxing TSI of fumarate. The decay of the magnetically silent TSI of fumarate will be reflected indirectly in the decay of the malate signal as n increases. Moreover, it can be seen in Figure 4 b that only around 1 % of the HELLS of fumarate is transferred to malate during each loop n=1, 2, …, N.

Figure 4

a) Temporal evolution of the (unobservable) PTSI of fumarate and b) signal of malate obtained by numerical solution of Equation (2) with TTSIF=60 s and TLLSM=6 s, and the apparent forward and backward rate constants kMF=4.5×10−4 s−1 and kFM=3.7×10−4 s−1=0.825 kMF optimized by fitting the curves in Figure 2. The vertical scale was increased 100 times in (b) to show the malate signal, which is barely visible as (-•-) in (a) because of the slow rate of the enzymatic conversion.

Results

A sample comprising ten frozen pellets of 10 μL each of 0.5 M fumarate with 50 mM TEMPOL was hyperpolarized by microwave irradiation at B0=6.7 T and T=1.2 K for about 20 min. The sample was then dissolved, together with ten frozen pellets of 10 μL each of 3 M sodium ascorbate in D2O,17 with 5 mL D2O at 400 K and 1.0 MPa, and transferred in 4.5 s to a holding chamber just above the magnetic center of a 7 T NMR (300 MHz) spectrometer, where the static field is Bhold>6.5 T. After a preinjection delay 1<τTSI<60 s, which allows one to assess the lifetime TTSI of the TSI (PTSI) of hyperpolarized fumarate in the holding chamber, the solution was injected into a 5 mm NMR tube containing fumarase to start the conversion of fumarate into malate, and to transfer concomitantly the TSI of fumarate into an LLS on malate. The latter was sustained by a WALTZ-16 pulse train with an rf amplitude ν1=3 kHz. The sequence of Figure 3 was then used to convert the LLS of malate into observable magnetization.

Figure 5 d shows four spectra of malate acquired at 7 s intervals (N=4 loops, each comprising a sustaining interval τLLS=6 s and an acquisition time of 1 s) after the injection of hyperpolarized fumarate into the NMR tube containing fumarase. In this case, the preinjection delay τTSI=1 s during which the fumarate was kept in the holding chamber was negligible compared with TTSIF. The enzymatic conversion is relatively slow, so the signals in Figure 5 d arise from the conversion of a small fraction of fumarate into malate (≈1 % every 7 s, according to Figure 4 b). The decay of the malate signal with increasing n reflects 1) the decay of the inaccessible TSI of fumarate with a time constant TTSIF owing to its relaxation (believed to be very slow), 2) the consumption of fumarate with a time constant 1/kFM owing to its enzymatic conversion into malate, and 3) the decay of the LLS of malate with a time constant TLLSM=6 s (Figure 4 a). However, it is risky to extract a reliable estimate of TTSIF from numerical fits of a single decay.

Figure 5

a) Conventional NMR spectrum excited by a 90° pulse 25 min after injection into a solution containing fumarase, when the enzymatic reaction has reached a steady state and the hyperpolarization (both PTSIF in fumarate and PLLSM in malate) has decayed to thermal equilibrium. Note the signals of fumarate, malate, ethanol, and buffer. The HDO peak was attenuated by presaturation with a selective pulse with an rf amplitude of 75 Hz and a duration of 5 s. b) Spectrum of malate (without significant stopover in the holding chamber since τTSI=1 s≪TTSIF) recorded with the sequence of Figure 3, shortly after injection (n=1) into a solution containing fumarase in the 7 T NMR system. c) Spectrum of malate recorded after keeping the hyperpolarized fumarate for τTSI=60 s in the holding chamber at B0>6.5 T prior to injection into the fumarase solution. d) The first four spectra of malate acquired with n=1, 2, 3, and 4 at intervals of 7 s using the sequence in Figure 3 (τTSI=1 s, τLLS=6 s, acquisition time 1 s) showing that PLLSM is replenished through the enzymatic reaction.

The lifetime TTSIF can be estimated more accurately by repeating the entire experiment such that the fumarate that carries the hyperpolarized TSI is kept in the holding chamber during a longer preinjection delay τTSI=60 s. Although it is challenging to reproduce the experiment under identical conditions, we observed that the remaining signal of malate after τTSI=60 s was reduced by a factor of approximately 3.5 (Figure 5 b,c), implying that TTSIF≈50 s. This is somewhat shorter than the lifetime TTSI=270 s reported by Zhang et al.11 for deuterated dimethyl maleate produced by the addition of para-H2 to deuterated dimethyl acetylene dicarboxylate. This discrepancy may be caused by the presence of dissolved paramagnetic triplet oxygen in the superheated water used in our dissolution experiments, or to the presence of some residual TEMPOL radicals, as the reduction by ascorbate may not be quantitative. Because the WALTZ-16 pulse train destroys magnetization arising from any sources other than HELLS, single-quantum terms arising either from d-DNP or from a partial return to thermal Boltzmann equilibrium are wiped out (compare Figure 5 b,c with Figure 5 a). The detected signals can therefore unambiguously be traced back to the TSI of fumarate prepared by d-DNP, stored in the holding chamber for a time τTSI, and converted into LLS on malate by the enzyme. The remaining peaks of fumarate and HDO in the spectrum of Figure 5 b probably stem from hyperpolarized single-quantum magnetization that was not fully saturated by the WALTZ-16 pulse train and was brought into the active volume of the rf coil by convection.

Conclusion

We have shown that a pure TSI can be created readily by d-DNP in a system that contains two magnetically equivalent spins in solution. Once dissolved, this imbalance displays a lifetime TTSI that is much longer than the longitudinal relaxation time T1. We believe this to be the first proof of principle of the creation of hyperpolarized long-lived states for equivalent spins (HELLS) by d-DNP. Such a long-lived spin order can be used readily to monitor a slow enzymatic process of biochemical relevance, but may find applications in other areas of magnetic resonance such as imaging (MRI), for which hyperpolarization by d-DNP has become a technique of choice to enable metabolic imaging, and in which short lifetimes of hyperpolarized molecules are usually a major limitation. The HELLS methodology will be applied to more challenging molecules containing magnetically equivalent pairs of spins, such as CH2RR′, CH2Cl2, and possibly H2O. We are currently investigating molecules with interesting lifetimes and interesting chemical or biochemical properties that can be addressed by HELLS. As fumarate plays a crucial role in the Krebs cycle, it may be of interest for in vivo studies as it has been demonstrated to be a probe for cellular necrosis.18

Experimental Section

DNP samples

Solutions of 0.5 M dibasic sodium fumarate (Sigma–Aldrich) in the glass-forming mixture D2O:[D6]ethanol (60:40 v/v) were doped with 50 mM TEMPOL (Sigma–Aldrich). Ethanol was added drop by drop to avoid precipitation. The solution was then sonicated for 10 min. Ten frozen pellets of 10 μL each of this mixture were inserted in the polarizer, along with ten frozen pellets of 10 μL each containing 3 M ascorbate (Sigma–Aldrich) in D2O to scavenge the radicals after dissolution.17

DNP polarization and dissolution

DNP was performed at 1.2 K and 6.7 T in a home-built polarizer by applying frequency-modulated microwave irradiation19 at fμW=188.3 GHz and PμW=100 mW, with a modulation frequency of 10 kHz and modulation amplitude of 50 MHz. The polarized pellets were dissolved in 0.7 s with 5 mL D2O, preheated to T=400 K at P=1.0 MPa, and transferred in 4.5 s to a 7 T (300 MHz) magnet by pushing with helium gas at 0.6 MPa through a PTFE tube (1.5 mm inner diameter) running through a magnetic tunnel (3 m length).

Enzymatic detection

The hyperpolarized solution was kept at B0>6.5 T in a holding chamber just above the NMR sample tube for a variable delay τTSI to monitor the relaxation of the TSI of fumarate. The sample was then injected in 2 s into an NMR tube containing D2O (200 μL) for field-frequency locking, NaCl (200 mM) and TRIS buffer (25 mM), and fumarase (5 μL, 5.8 mg mL−1, 12.5 units) from porcine heart (Sigma–Aldrich). Finally, the LLS detection sequence described in Figure 3 was applied with n sustaining delays of τLLS=6 s each with WALTZ-16 irradiation.