Đorđe
Toljić
and
Goran
Angelovski
*
MR Neuroimaging Agents, Max Planck Institute for Biological Cybernetics, Tübingen, Germany. E-mail: goran.angelovski@tuebingen.mpg.de
First published on 17th September 2019
We report a novel ditopic Gd(III)-based probe selective to zwitterionic amino acid neurotransmitters (ZNTs) crafted for ratiometric MRI imaging. The probe displayed increased binding affinity to ZNTs and non-synchronized concentration-dependent changes of the r1- and r2-relaxivity. Through the application of a T2/T1 weighted MRI strategy, we demonstrated signal enhancement for cooperatively bound glutamate and γ-aminobutyric acid ZNTs over competitive hydrogencarbonate, which remained MR silent.
Molecular recognition between magnetic host sensors and guest molecules occurs via the formation of a ternary adduct, and is evaluated via its binding affinity. In the case of a ‘turn-off’ response, the interaction of biomarkers with the coordination cage of an MRI probe restricts inner-sphere water accessibility to the paramagnetic center resulting in a signal decrease.5 So far, the direct sensing of NTs has been approached with two distinct designs, via genetically engineered metalloproteins and by small-molecule ditopic paramagnetic probes.6–8 Although large molecule CAs tailored for monoamine neurotransmitters were shown to be capable of detecting target metabolites in the μM range, further improvement of small-sized probes is imperative due to issues associated with translation into biological systems. In terms of host–guest interaction, the low-molecular-weight sensors are bismacrocyclic Gd(III)-based complexes with a 1,4,7,10-tetraazacyclododecane-1,7-dicarboxylic acid (DO2A) chelator, whose geometry allows for easy access to small anions.9 In particular, the selectivity of the current sensors is greatly hindered by hydrogencarbonate (HCO3−), which binds in a similar fashion as the ZNTs. Owing to high extracellular concentration (∼25 mM),10 small molecular size and consequently low steric hindrances, HCO3− predominantly binds to the paramagnetic label, thus preventing selective interaction with ZNTs.7,8 On the other hand, the advantageous feature of ZNTs is their existence in the form of an ion-pair entity at physiological pH, which opens the possibility of exploiting a ditopic platform to strengthen binding through cooperativity.11 In accordance with this, different binding modes may unequally affect parameters such as inner-sphere water exchange rate or molecular tumbling (τm and τR, respectively), which determine both longitudinal (r1) and transverse (r2) relaxivities. In turn, the occurrence of non-synchronised changes of r1 and r2 would allow for the utilization of the ratiometric r2/r1 MRI approach.12–14
In order to develop such a sensor, it is necessary to evoke non-synchronized stimuli-induced changes of parameters that determine r1 and r2. Therefore, we set out to exploit the cooperative binding potential of a ditopic low-molecular weight Gd-chelate to improve the sensing of ZNTs. The probe design consisted of a bismacrocyclic framework, providing receptor sites to favour cooperative ion-pair binding of ZNTs. The host–guest interactions involve simultaneous coordination of the carboxylate to Gd(III) and multiple hydrogen bonds established between the 18-crown-6 derived moiety (18C6) and the NH3+ terminal of ZNTs (Fig. 1).7,8 The rationale for choosing 1,4,7,10-tetraazacyclododecane-1,4,7-tricarboxylic acid (DO3A) as the MR responsive unit over the coordinatively less saturated DO2A macrocycle is due to its higher selectivity towards carboxylates from ZNTs over HCO3−, owing to stronger electronic repulsions from the acetate pendant arms. To reinforce the affinity towards the ammonium cation, we incorporated a formanilide pendant that bears two amide groups, one proximate and one distant to the 18C6 ring-frame, expecting to engage as the H-bond donor and additionally stabilize the binding of the ammonium cation.15 The role of the flexible amide linker is to preserve the high degree of conformational mobility of the binding sites, providing an adaptable geometry for cooperative binding of the zwitterions (Fig. 2).
![]() | ||
Fig. 2 Binding modes of structurally different guests with GdL. (a) Monotopically bound HCO3− and (b) ditopically bound ZNT. |
The proposed sensor GdL was synthesized by bridging two macrocyclic components, followed by complexation with Gd(III) (Scheme 1). The 18C6 fragment 1 was obtained by direct mono-N alkylation of the 1,10-diaza-18-crown-6 ether with 2-bromo-N-(4-nitrophenethyl)acetamide. Next, the macrocyclic precursors 1 and 2 were covalently linked to form bismacrocyclic derivative 3. As the preparation of bismacrocycles of this type is often limited by a bridging step, the yield of 62% can be considered as excellent. In the following step, the aromatic nitro group of 3 was subjected to Pd(OH)2/C-catalysed hydrogenation to afford the amine 4. Simultaneous hydrolysis of the tert-butyl esters and conversion of the amine into formamide was conducted in formic acid to obtain the final ligand L with a yield of 64%, and a ratio of isomers of 72:
28.16 The complexation was performed in water media at neutral pH with GdCl3 to give the corresponding complex GdL. To avoid entrapment of a Gd(III) ion within the 18C6 component, the complexation was performed using a ligand to Gd(III) ratio of 1.0
:
0.8 and the excess of ligand was effectively removed by the HPLC. The absence of free Gd(III) ions was additionally confirmed by the xylenol orange test.
The bulk 1H T1- and T2-relaxation time dependencies of GdL on the concentration of the ZNTs: Glu, glycine (Gly), GABA, the non-zwitterionic transmitter acetylcholine (ACh) and competitive metabolite HCO3− were assessed in relaxometric titrations (Fig. 3). The initial values calculated for r1 and r2 were 10.14 and 13.46 mM−1 s−1, respectively. The displayed downward trends revealed that GdL is ACh-insensitive while HCO3− exhibited a comparable drop in r1 (69%) to the ZNTs that induced changes, with −Δr1 ranging from 61 to 69%. However, notably diverse r2 changes between HCO3− and ZNTs were observed. Upon saturation with the guests (60 equiv.), the r2 of ternary adduct GdL-HCO3− was 3.69 mM−1 s−1, while the values for the -Gly, -Glu and -GABA adducts were 5.64, 5.78 and 5.89 mM−1 s−1, respectively.12 The gathered results are in agreement with the simulations performed for low-molecular-weight Gd(III)-based CAs at high magnetic fields (>1.5 T).12 Namely, T1 is relatively constant over a broad range of values for the τm and τR parameters. On the contrary, for the same parameters, T2 undergoes a rather drastic change.12 Hence, the occurrence of non-synchronized changes of T1 and T2 within the same MR reporter, opens the possibility to implement a r2/r1 ratiometric MRI strategy. When presented as the r2/r1 ratio, there is an increase from the initial value of 1.32 for GdL to 1.7, 1.9 and 2.1 for Gly, Glu and GABA, respectively. In the case of HCO3−, the decrease in r1 is followed by a proportional decrease in r2, resulting in insignificant fluctuations for the r2/r1 values over the entire titration range.
![]() | ||
Fig. 3
1H relaxometric study of GdL. (a and b) the experimentally measured (symbols) and fitted data (curves) for (a) r1 and (b) r2 relaxivities for GdL (2.5 mM) with different analytes (0–100 equiv.) in 50 mM HEPES (pH 7.4, 7.05 T, 298 K). Lines in (a) correspond to the fit described in the ESI,† while the dashed lines in (b) show the mismatch of the experimental to those fitted in the same way, while assuming Ka values obtained from fitting the r1 data. (c) Decreases in r1 and r2 for the host–guest molar ratio of 1![]() ![]() |
The binding affinities (Ka) were obtained by fitting the r1 data to eqn (S2) (ESI†) (Table 1). Comparison of Ka values revealed higher affinities of 116 and 106 M−1 for GABA and Glu, respectively, over HCO3− and Gly, which were 86 and 64 M−1, respectively. The observed trends suggested greater structural complementarity of the receptor-binding pocket to the long-chain ZNTs, Glu and GABA. Comparing HCO3− and Gly, the former binds stronger due to its potential to interact in the bidentate binding mode. The decreased affinity of Gly compared to GABA and Glu is likely due to the absence of cooperative host–guest binding as the carboxylate– amine terminals in GABA and Glu are more suited to the spatial disposition of binding sites within the host. The ambient concentrations of Glu, GABA, Gly and Gln in physiological brain states of ∼6 μM, 0.23 μM, 16.7 μM and 1.738 mM, respectively, are still below the detection level of the prepared probe.17 However, in pathological states the concentrations of ZNTs drastically increase, e.g., in the case of ischemia, Gly reaches a value of 513 μM. Therefore, the affinity constants for GdL, which are on the level of ∼10 mM, still require slight improvements to allow for the detection of the relevant ZNTs, whose concentrations are closer to low mM values.
The conformational mobility of the binding sites in the host molecule is affected by the interactions with guest molecules. Therefore, to gain insight into the host–guest binding modes, various NMR studies were performed on the diamagnetic yttrium(III) analogue, YL.18 High-resolution 1H NMR spectra were recorded with YL in the presence of 20 equiv. of Glu or HCO3− (Fig. 4). The resonances of both ternary complexes shifted upfield, with a more pronounced effect for Glu, which suggests its stronger host–guest interaction. Furthermore, splitting of the broad signal at 1.69 ppm in ternary adducts, corresponding to the methylene CH22CH2 of the amide linker that connects the two binding sites, indicated the existence of both bound and non-bound forms. Overall, this revealed host–guest binding with a slow exchange rate and a greater magnitude for YL-Glu interactions.
Anion binding to the paramagnetic label coordinated by DO3A-based chelates is well documented, and thus, we focused on investigating interactions that involve the 18C6 binding site. 2D EXSY NMR spectra displayed conformational diversity for the 18C6 moiety within the ternary adducts (Fig. 4 and Fig. S6–S9 in ESI†). The higher density of cross-peaks in YL-Glu pointed towards a binding mode that involves multiple interactions between the 18C6 and the ammonium cation. Conversely, the more defined cross-resonances for 18C6 in YL-HCO3− were assigned to Na+ inclusion (bicarbonate was added as the monosodium salt). Moreover, no change was observed in the ratio of the isomers in the 1H NMR spectra (see above), which led us to exclude the involvement of the formanilide moiety in binding.
Another parameter sensitive to molecular motions of bound protons is the T1 relaxation time.19 Thus, to quantify the observed conformational flipping, we measured T1 of protons in the NCH22O groups within the 18C6 fragment and the meta-H in the remote pendant (Table S2 and Fig. S10 in ESI†). The acquired T1 values for the 18C6 moiety in YL-Glu and YL-HCO3− were 692 ms and 651 ms, respectively, whereas the one of binary YL was measured to be 628 ms. The elevated T1 value in the presence of Glu indicated a higher rate of conformational exchange of the 18C6 binding site. This was rationalized as the impact of reversible interactions of multiple hydrogen bonds between NH3+ of the zwitterion and free electronic pairs on the heteroatoms in 18C6. Slightly increased mobility of 18C6 in the monotopically bound YL-HCO3− is likely related to stronger interactions with entrapped Na+. The identical trend of mobility was observed for the 18C6 pendant, which is in alignment with the higher degree of freedom for linear compounds.
Finally, 1H DOSY NMR measurements provided information on the host–guest association in solution.20 The obtained self-diffusion constants displayed faster translational motion for the ternary adducts, YL-Glu and YL-HCO3−, compared to the binary YL (Fig. S11 and Table S3 in ESI†). The corresponding hydrodynamic radii shrunk from 1.54 nm to 1.29 nm upon moving from the binary complex to the ternary species, suggesting a different spatial arrangement of YL prior to its interaction with guests. These values excluded the formation of aggregates or intermolecular self-assembly. The rotational correlation times (τR), calculated according to the Stokes–Einstein-Debye model (see ESI†), were identical for the ternary complexes (234 ps), which ruled out molecular tumbling as the explanation for the observed diversity in T2 trends.21
To capitalize on this effect, we tested the selectivity of the GdL probe by comparing the signal enhancements of T1w and T2w with T2/T1w images of MRI tube phantoms in a 7.05 T MRI scanner (Fig. 5 and Fig. S12, S13 and Tables S4, S5 in ESI†). The contrast enhancements were calculated as differences in obtained signal-to-noise ratio (SNR). The T1w images gave negligible ΔSNRs for Gly, GABA and Glu with respect to the signal enhancement produced by HCO3−, which are in the range 1–4%. On the other hand, HCO3− caused greater MRI signal alteration in the T2w images, thus masking the effects of other ZNTs. Specifically, Gly, GABA and Glu only reached 70, 47 and 57% of the signal enhancement observed for HCO3−, respectively. In contrast to the results obtained using T1w and T2w MRI protocols, the T2/T1w phantoms gave more intense signals for Gly, GABA, and Glu compared to HCO3− (4, 11 and 10%, respectively). Moreover, HCO3− remained completely MR silent with the ΔSNR equal to that of GdL. In addition, comparing the greater potential of the T2/T1w MRI protocol to the commonly used T2w analogue to produce images with greater SNR,14 it should be noted that the effective discrimination of HCO3− was successfully achieved under physiologically-relevant conditions.
![]() | ||
Fig. 5
T
2/T1w phantom image of GdL (1.5 mM) with ZNTs and competitive metabolites in a ratio of 1![]() ![]() |
In summary, we have developed a novel ditopic bismacrocyclic MRI probe with an increased binding affinity towards the major inhibitory (GABA) and excitatory (Glu) neurotransmitters and employed an r2/r1 imaging strategy to discriminate from their physiological competitor, HCO3−. We demonstrated that the bound guests induced uneven changes of r1 and r2 in the supramolecular adducts. Unlike the long-chain ZNTs, the monotopically bound HCO3− induced a proportional decrease in both r1 and r2, whose ratio remained equal to that of unbound GdL. Hence, greater SNR changes were observed for GABA and Glu in MRI phantoms, and were less emphasized for Gly. Overall, these findings establish an efficient strategy to sense polydentate metabolites by employing low-molecular-weight polytopic MR probes. We anticipate that the introduced concept of utilizing a polytopic host with a ratiometric imaging methodology will be beneficial in the future development of MR probes for sensing polycharged guest molecules.
The authors are thankful to Dr Vincent Truffault for measuring NMR spectra and Dr Tanja Savić for recording the MRI phantom images. The financial support of the German Research Foundation (DFG, grant AN 716/7-1) is gratefully acknowledged. Open Access funding provided by the Max Planck Society.
Footnote |
† Electronic supplementary information (ESI) available: Methods, synthetic procedures, and selected spectroscopic and MRI data. See DOI: 10.1039/c9cc06463j |
This journal is © The Royal Society of Chemistry 2019 |