Ben. J.
Tickner
a,
Jennifer S.
Lewis
a,
Richard O.
John
a,
Adrian C.
Whitwood
b and
Simon B.
Duckett
*a
aCenter for Hyperpolarization in Magnetic Resonance (CHyM), University of York, Heslington, York YO10 5NY, UK. E-mail: simon.duckett@york.ac.uk
bDepartment of Chemistry, University of York, Heslington, York YO10 5DD, UK
First published on 24th September 2019
Signal Amplification By Reversible Exchange (SABRE) is a hyperpolarisation technique that commonly uses [Ir(H)2(carbene)(substrate)3]Cl complexes to catalytically transfer magnetisation from para-hydrogen derived hydride ligands to coordinated substrates. Here, we explore the reactivity of a novel class of such catalysts based on sulfoxide containing [IrCl(H)2(carbene)(DMSO)2], which are involved in the hyperpolarisation of pyruvate using SABRE. We probe the reactivity of this species by NMR and DFT and upon reaction with sodium pyruvate establish the formation of two isomers of [Ir(H)2(η2-pyruvate)(DMSO)(IMes)]. Studies with related disodium oxalate yield [Ir2(H)4(IMes)2(DMSO)2(η2-κ2-Oxalate)] that is characterised by NMR and X-ray diffraction.
Many catalytic processes exploit the oxidative addition of H2 to a transition metal centre by enabling the subsequent transfer of the two hydrogen atoms into an unsaturated centre13,14 However, recent examples in main group chemistry involving frustrated Lewis pairs mean that direct addition is also possible.15 Understanding of these reactions can be improved by exploiting a technique known as para-Hydrogen Induced Polarisation (PHIP) to detect reaction intermediates.16,17 This process incorporates the two protons of a single para-hydrogen (p-H2) molecule into a substrate via a hydrogenation reaction to see the PHIP effect. p-H2 is a spin isomer of H2 and exists as a singlet state with a nuclear spin order of zero. H2 gas can be enriched (>98%) in its para state by simply cooling it in the presence of a paramagnetic catalyst such as Fe2O3 or activated charcoal.18–20 If the spin orientation of p-H2 survives the ensuing hydrogenation reaction the NMR properties of the product can be harvested to increase product detectability. This effect was exemplified in 1986 by Bowers and Weitekamp21,22 with Eisenberg and Bargon23,24 producing similar independent observations around the same time. As the net nuclear spin of p-H2 is zero it is invisible to an NMR experiment but once the symmetry of its two protons is broken the resulting spin order can be detected. In fact, the NMR signal intensity of the now NMR visible, and hyperpolarised, product is derived from the resulting large non-Boltzmann populations that lie across its nuclear spin energy levels. Since its introduction PHIP has been used to detect low concentration analytes and true intermediates in the field of catalysis.22,25–27 Excitement has resulted from the hydrogenation of some unsaturated clinical agents which has led to their in vivo MRI detection.28–30
Since 2009, the p-H2 based method Signal Amplification By Reversible Exchange (SABRE) has been used to harness signal gains from p-H2 without the need for the direct hydrogenation of a substrate (sub).31 It achieves the catalytic transfer of magnetisation into a substrate through the formation of a J-coupled network within the associated catalyst.32 The first step in SABRE typically involves the conversion of a stable 16 electron precursor such as [IrCl(COD)(IMes)] (1) (where IMes = 1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene and COD = cis,cis-1,5-cyclooctadiene) into an intermediate [Ir(COD)(IMes)(sub)]Cl which then undergoes hydrogen addition to form highly reactive [Ir(H)2(COD)(IMes)(sub)]Cl.33–35 Subsequent hydrogenation of the COD ligand then leads to a SABRE catalyst such as [Ir(H)2(IMes)(sub)3]Cl which importantly reversibly binds both H2 and sub. In SABRE, H2 addition to square planar [Ir(COD)(IMes)(Sub)]Cl typically occurs over the COD-Ir-Sub axis rather than the COD-Ir-IMes axis.36 This observation has been explained in terms of both reduced steric interactions with the bulky carbene ligand and electronic effects in terms of orbital interactions.36 It has been suggested that H2 addition is favoured over axis containing ligands with π acceptor orbitals which provide additional stabilising interactions as the geometry of the complex changes during the addition step.37
The most common substrate molecules used in SABRE have proven to feature N-donor sites that readily coordinate to iridium. N-Heterocycles such as pyridine,31,38–41 nicotinamides,31,42,43 pyrazines,39,40 and pyrazoles41 reflect common examples although other N-functionalities found in nitriles44 and amines45 have been used. It has recently been reported that in the presence of a stabilising sulfoxide coligand, the reversible coordination and subsequent hyperpolarisation of oxygen ligating pyruvate can occur.46
The direct hyperpolarisation of pyruvate was made possible in this case by the formation of the novel polarization transfer catalyst [Ir(H)2(η2-pyruvate)(DMSO)(IMes)]. In addition to this active catalyst, [IrCl(H)2(DMSO)2(IMes)], is also present in solution and is expected to be critical in improving the efficiency of this important hyperpolarisation transfer process. In this work we investigate the formation, behaviour, and ligand exchange processes exhibited by this sulfoxide containing polarisation transfer catalyst. We do this in solvents that were purchased from Sigma and used without further purification as our aim is to study reactivity without taking any special precautions. This is reflective of the most likely scenario when used by the non-specialist. We extend this method to include an examination of the behaviour of oxalate with the aim of developing the range of materials that can be hyperpolarised with sulfoxide containing SABRE polarisation transfer catalysts.
DFT calculations have been used to further confirm these product assignments, with their predicted relative energies detailed in Table 1 (see ESI‡ for details). These calculations used full models of the complexes, [Ir(COD)(IMes)(X)], and are relative to a zero point which includes all species (indicated complex, DMSO, H2O and CH3OH) in an equimolar amount. They do not reflect transition state barriers and excess reagents are not accounted for. These calculations confirm that 1, where X = Cl, is most stable in agreement with experiment. [Ir(COD)(IMes)(OH2)]Cl, 2 is predicted to be more stable than the corresponding methanol or sulfoxide adducts, again supporting our NMR observations.
X | H/kJ mol−1 | G/kJ mol−1 |
---|---|---|
Cl | 0 | 0 |
OH2 | 52.8 | 57.0 |
S-DMSO | 42.3 | 68.9 |
HOCH3 | 56.1 | 69.4 |
The identity of 3 was further confirmed by repeating this experiment using [IrBr(COD)(IMes)], as the corresponding reaction product [IrBr(H)2(DMSO)2(IMes)] yields hydride signals at δ −15.67 and δ −20.45 in methanol-d4 which are clearly different to those of 3. Furthermore, when [Ir(CH3CN)(COD)(IMes)]PF6 is used instead, [Ir(H)2(CH3CN)(DMSO)2(IMes)]PF6 forms which yields signals at δ −15.76 and δ −19.83 at 255 K; we note there was no evidence for the displacement of CH3CN by DMSO in the associated NMR spectra. Hence, the hydride resonance in these complexes shift according to the identity of the ligand that is trans to it thereby confirming the indicated product identities. The identity of 3 was also studied by DFT as detailed in the results shown in Table 2. [IrCl(H)2(DMSO)2(IMes)] was found to be lower in energy than the corresponding species where chloride is replaced by water, methanol or sulfoxide.
X | H/kJ mol−1 | G/kJ mol−1 |
---|---|---|
Cl | 0 | 0 |
OH2 | 24.4 | 34.3 |
S-DMSO | 39.1 | 70.8 |
ODCH3 | 29.3 | 45.4 |
When the reaction of H2 with an equilibrium mixture of 1 and 2 in methanol-d4 containing DMSO was monitored at 245 K (Fig. 2) no evidence for the formation of any H2 addition products except 3 is observed. This suggests that H2 addition is slow, and subsequent COD hydrogenation in intermediate 1-H2 of Fig. 2 is fast. Interestingly, the proportion of 2 remains roughly constant as this conversion proceeds thereby suggesting any equilibrations involving it are also slow. The route to 3 is therefore most likely to involve direct H2 addition to 1 rather than 2. DFT confirms that H2 addition to 1 proceeds over its COD-Ir–Cl axis rather than the COD-Ir-IMes axis according to the relative energies of the corresponding products (see ESI‡). This is supported by the fact that when this H2 addition reaction is monitored in CDCl3 at 245 K, resonances corresponding to 1-H2, at δ −13.39 and δ −18.42, are detected in addition to those of 3. The resonances for 1-H2 rapidly disappear upon warming this solution to 298 K where 3 then forms.
Measuring the change in these resonance's signal intensities as a function of time after the initial r.f. excitation step allowed the associated ligand loss rates to be calculated. The rate of reversible hydrogen loss, kH2, was found to be 3.31 ± 0.25 s−1 at 263 K and is faster than that seen in many related systems.52 Loss, and hence exchange of the DMSO ligand cis to IMes, was found to occur on a similar timescale with kDMSO 3.35 ± 0.01 s−1 at 263 K. Repeating these measurements at different temperatures allowed the transition state barriers for these two processes to be determined as ΔH±(H2) 78 ± 6 kJ mol−1, ΔS±(H2) 60 ± 25 J K−1 mol−1 and ΔH±(DMSO) 79 ± 6 kJ mol−1, ΔS±(DMSO) 66 ± 8 J K−1 mol−1 respectively. This enthalpic barrier to hydrogen loss is comparable to those for similar complexes, overall the entropy change suggests that a dissociative pathway is adopted.53
The effects of H2 pressure and DMSO concentration on these exchange processes were also studied in methanol-d4 at 243 K for differing reagent concentrations. kH2 proved to increase at higher H2 concentrations before reaching a plateau as detailed in the ESI.‡ The rate of DMSO exchange is unaffected by increasing DMSO or H2 concentration. We therefore propose that DMSO and H2 loss proceed via formation of the common 16-electron intermediate [IrCl(H)2(DMSO)(IMes)] in what is a dissociative first step. H2 exchange then proceeds via the formation of [IrCl(H)2(H2)(DMSO)(IMes)]; such indirect H2 exchange, rather than direct H2 loss, has been reported for several similar systems.51,54 When similar EXSY data were collected in CD2Cl2, the rate of hydrogen and DMSO loss proved to be slower than those in methanol-d4, as summarised in Table 3.
Solvent | Process | k (263 K)/s−1 | ΔH≠/kJ mol−1 | ΔS≠/J K−1 mol−1 |
---|---|---|---|---|
Methanol-d4 | Hydrogen exchange (kH2) | 3.31 ± 0.25 | 78 ± 6 | 60 ± 25 |
DMSO exchange (kDMSO) | 3.35 ± 0.04 | 79 ± 6 | 66 ± 8 | |
Dichloromethane-d2 | Hydrogen exchange (kH2) | 1.16 ± 0.04 | 83 ± 8 | 73 ± 29 |
DMSO exchange (kDMSO) | 1.56 ± 0.01 | 84 ± 2 | 78 ± 8 |
These deductions were confirmed through further DFT calculations that revealed the products of direct H2 loss to form a 16 electron product as being highly energetically unfavourable, as shown in Table 4. The five coordinate product formed by loss of DMSO trans to carbene proved highly unstable, undergoing spontaneous rearrangement to form an intermediate with a vacancy in the equatorial plane; this is consistent with the EXSY data. In contrast, the predicted ligand dissociation pathway involves the formation of five-coordinate [IrCl(H)2(DMSO)(IMes)] via loss of the DMSO ligand that lies cis to carbene.
Loss of | H/kJ mol−1 | G/kJ mol−1 |
---|---|---|
H2 | 89.0 | 58.5 |
Cl | 39.6 | 3.9 |
DMSO cis carbene | 32.7 | −10.1 |
Even at low temperature (245–265 K), H/D exchange leading to [IrCl(D)(H)(DMSO)2(IMes)] (3-d) and [IrCl(D)2(DMSO)2(IMes)] was evident. Binding of solvent methanol to [IrCl(H)2(DMSO)(IMes)] is an obvious route to deuterium exchange and the formation of these species alongside HD and D2 gas. While we do not directly observe methanol bound adducts, they are suggested to form in similar systems and indeed proposed to account for the hydrogen isotope exchange reaction that is often observed to run in parallel to SABRE.50,55 Experimentally, at low temperature it proved possible to reliably and selectively excite the hydride resonance of 3 or 3-d. The selective excitation of the hydride resonances for these two species revealed their selective chemical exchange into H2 and HD, from 3 and 3-d respectively. Therefore, kH2 and kHD can be determined as previously described. These values proved to be the same within error and hence, there is little or no kinetic isotope effect which is consistent with other reports.56,57 We note that exchange between 3 and 3-d is not observed in the associated EXSY data which provides confirmation that the underlying deuterium exchange processes involving methanol-d4 are slow.
Complex | H/kJ mol−1 | G/kJ mol−1 |
---|---|---|
4a | 8.5 | 14.8 |
4b | 0 | 0 |
4c | 22.6 | 28.1 |
When this solution was examined with a 32 scan 1H NMR measurement at 298 K, unlike the 245 K data of Fig. 3b, the main hydride containing complex observed is 4b and resonances for 3 and 4a cannot be discerned. 4a and 4b have previously been implicated in pyruvate polarisation using SABRE and their NMR characterisation data has been reported.464 is also detected immediately after H2 addition to a solution of 1 containing pyruvate (5 equivalents relative to iridium) and DMSO (10 equivalents) at 298 K. When an equilibrium mixture of 3 and 4 is shaken with 3 bar p-H2 at 65 G enhanced hydride resonances are observed in the 1H NMR spectrum as shown in Fig. 3c, which are strongest for 4b.
If sodium-1,2-pyruvate-[13C2] is used as the substrate and the p-H2 shaking process is performed under SABRE-SHEATH conditions in a mu metal shield (see Experimental) an enhanced 13C response for the free material at δ 169 and δ 203 can also be readily detected.46 Two enhanced signals can be also seen for the bound pyruvate ligand in 4b at δ 168 and δ 207 in these NMR spectra. It is also possible to observe extremely weak hyperpolarised signals for the bound pyruvate ligand in 4a at δ 161 and δ 198, as shown in Fig. 4d. Interestingly, when EXSY measurements are used to probe these signals, the selective excitation of the hydride resonances of 4b reveals no exchange into H2 on the NMR timescale. Furthermore, when 13C-EXSY is used to probe the bound signals of 4b, no exchange is seen into free pyruvate on this timescale. Hence 4b appears to be relatively stable in agreement with the DFT study, but this is not consistent with the strong hydride polarisation that is evident for 4b in Fig. 3. Consequently, a role for 3 which is present and undergoes rapid H2 exchange as discussed earlier in its formation is indicated. It is by this route that the delivery of fresh p-H2 into 4 and the resulting enhancement of hydride and pyruvate ligand signals is achieved. As the NMR relaxation times of the 13C resonances in pyruvate are much longer than those of 1H, slow exchange is still commensurate with the build-up of pyruvate polarisation in solution. Related relayed p-H2 enhancement effects have been described recently.58,59
An equilibrium mixture containing 4 can also be formed if sodium-1,2-pyruvate-[13C2] is added to a preformed methanol-d4 solution of [IrBr(H)2(DMSO)2(IMes)] (3-Br) or [Ir(H)2(CH3CN)(DMSO)2(IMes)]PF6 (3-Acn). When these samples are shaken with p-H2 for 10 seconds at 65 G strongly hyperpolarised hydride resonances can again be observed for 4b in the corresponding high field NMR measurement. These resonances appear with significantly lower intensity than those achieved with 3 derived systems (35% and 31% of this signal intensity respectively). Furthermore, when the p-H2 shaking step is repeated in a mu metal shield, hyperpolarised resonances for 13C-pyruvate are again visible. The corresponding 13C-pyruvate signal enhancements are typically 580-fold and 120-fold for the 3-Br and 3-Acn derived systems respectively. These are lower than those achieved for the analogous 3 system (1070-fold) which is consistent with the reduced amount of 4b which is present at equilibrium when these precursors are used (22 and 23% respectively when compared to 92% with 3).
We note that the linewidths of the hydride resonances of 4b at 298 K are similar in each of these samples (43–45 Hz). Similarly, the linewidths of 3, 3-Br and 3-Acn are comparable at 245 K (3–5 Hz). Therefore, we do not expect that differences in hydride polarisation levels seen for 4b in these solutions are due to different hydrogen exchange rates in 3, but rather the different binding strengths of the coligands that must be displaced by pyruvate to form 4 from 3. Hence, whilst there appears to be a link between pyruvate polarisation level and 4b concentration, the identity of L in [Ir(H)2(DMSO)2(IMes)L] must play a large effect on the level of pyruvate signal enhancement. These data further confirm that 3 is important in mediating efficient H2 exchange within the 3/4 hyperpolarisation mixture.
While the complexation dynamics of oxalate are complex,62–64 based on the pyruvate observations described earlier it should be possible to form a species such as [Ir(H)2(η2-oxalate)(DMSO)(IMes)], 5 which would account for this observation. In fact, the 13C2 NMR signal profile shown in Fig. 4 is indicative of Zeeman magnetisation in such a product rather than singlet spin order. This indicates that while 13C polarisation transfer from a set of p-H2 derived hydride ligands initially results this must be associated with a [AA′BB′] spin system, which converts to the [AMSX] type with inequivalent hydride ligands in high field. A product with ligand arrangement analogous to 4a would match with this hypothesis.46 Unfortunately, examination of the hydride region of the resulting 1H NMR spectra reveal signals for 3 at δ−15.45 and δ −21.53 alongside over seven minor hydride containing complexes. None of the hydride resonances of these species correlated to these two 13C signals in an overnight 2D measurement at 243 K. Studying this reaction again at 315 K, proved to be of no benefit as conversion of 3 into the same range of hydride containing species results.
In contrast, when this experiment was repeated in DCM-d2 the hydride resonances of 3 are again observed at δ −15.71 and δ −21.27 but no 13C oxalate derived polarisation is seen. There is very low solubility of sodium oxalate in this medium. As a consequence, samples of 3 in CD2Cl2 or CH3OD were prepared and reacted with 50 μL D2O solutions of sodium-1,2-oxalate-[13C2]. Now a new product forms cleanly, that yields a hydride signal at δ −27.1 alongside diagnostic resonances for the IMes at δ 2.10, δ 2.34, δ 6.99, δ 7.15 and DMSO at δ 2.92 whose relative signal intensities suggest the presence of a [Ir(H)2(IMes)(DMSO)(L)n]x species, 6. Upon shaking with p-H2 the hydride resonances for 6 do not enhance, however, the addition of 3 mL of degassed hexane results in the growth of single crystals. Subsequent X-ray diffraction studies confirmed that the product 6 corresponds to the dimer [Ir2(H)4(IMes)2(DMSO)2(η2-κ2-Oxalate)] with structure shown in Fig. 5 (full X-ray and NMR characterisation data are given in the in the ESI‡). The ability of oxalate to form such dimers restricts its dissociation and thereby limits its use in these hyperpolarisation studies. Nonetheless, the structure of 6 reflects a common binding mode of oxalate with transition metals.63,65,66
Fig. 5 Structure of [Ir2(H)4(IMes)2(DMSO)2(η2-κ2-Oxalate)], 6, determined by X-ray crystallography. All non-hydride hydrogen atoms and solvent of crystallisation have been omitted for clarity. |
Samples were prepared containing 2 mg [IrCl(COD)(IMes)] precatalyst (where IMes = 1,3-bis(2,4,6-trimethyl-phenyl)imidazole-2-ylidene and COD = cis,cis-1,5-cyclooctadiene) in 0.6 mL of deuterated methanol-d4 unless otherwise stated in a 5 mm NMR tube that was fitted with a J. Young's tap. All commercial compounds were purchased from Sigma-Aldrich, Fluorochem, or Alfa-Aesar and used as supplied. [IrCl(COD)(IMes)] was synthesized according to a literature procedure.67 The resulting solutions were degassed by two freeze–pump–thaw cycles before the addition of 3-bar H2.
The shake and drop method was employed for recording hyperpolarised NMR spectra. This involves filling NMR tubes with p-H2 at 3 bar pressure and shaking them vigorously for 10 seconds in a 65 G magnetic field if 1H NMR spectra are to be recorded at 9.4 T. A mu metal shield is used if 13C NMR spectra are to be recorded. Unless otherwise stated multiple shake and drop measurements are undertaken and average signal enhancement values quoted. Signal enhancements and exchanges rates were calculated according to literature procedures.68
Footnotes |
† Dedicated to Prof. Robin N. Perutz on the occasion of his 70th birthday. |
‡ Electronic supplementary information (ESI) available. CCDC 1941589. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c9dt02951f |
This journal is © The Royal Society of Chemistry 2019 |