Original Research ARTICLE
In-depth Study of a Novel Class of Ditopic Gadolinium(III)-based MRI Probes Sensitive to Zwitterionic Neurotransmitters
- 1MR Neuroimaging Agents, Max Planck Institute for Biological Cybernetics, Tuebingen, Germany
- 2Centro de Investigacións Científicas Avanzadas, Departamento de Química, Facultade de Ciencias, Universidade da Coruña, A Coruña, Spain
The efficacy of Gd-based low-molecular weight ditopic MRI probes on binding zwitterionic neurotransmitters (ZNTs) relies on their structural compatibility. ZNTs are challenging biomarkers for monitoring chemical neurotransmission due to their intrinsic complexity as target molecules. In this work, we focus on tuning the cyclen- and azacrown ether-based binding sites properties to increase the affinity toward ZNTs. Our approach consisted in performing structural modifications on the binding sites in terms of charge and size, followed by the affinity evaluation through T1-weighted relaxometric titrations. We prepared and investigated six Gd3+ complexes with different structures and thus properties, which were found to be acetylcholine insensitive; moreover, two of them displayed considerably stronger affinity toward glutamate and glycine over hydrogencarbonate and other ZNTs. Complexes with small and non-charged or no substituents on the azacrown moiety displayed the highest affinities toward ZNTs, followed by strong decrease in longitudinal relaxivity r1 of around 70%. In contrast, hosts with negatively charged substituents exhibited lower decrease in r1 of nearly 30%. The thorough investigations involving relaxometric titrations, luminescence, and NMR diffusion experiments, as well as theoretical density functional theory calculations, revealed that the affinity of reported hosts toward ZNTs is greatly affected by the remote pendant on the azacrown derivative.
A step forward into solving still puzzling physiological and pathological brain conditions would not be possible without a thorough understanding of neural activities. These signaling and metabolic processes consist of a chain of complex interactions and biochemical processes regulated by neurotransmitters, whose imbalance results in numerous neuronal pathologies that permanently damage nervous tissue (Lukasiuk and Pitkänen, 2000; Lewerenz and Maher, 2015). Resolving spatial and temporal concentration change patterns of neurotransmitters released into the synaptic cleft will provide information on the type, localization and degree of occurrence of neuropathological events. Such findings would aid establishing an early diagnosis and allow monitoring of a treatment progress. Therefore, to better understand neuronal events, it is necessary to directly profile neurochemical abnormalities by attempting to understand and assess the diversity in their chemical behavior.
Chemical neural signaling that involves zwitterionic neurotransmitters (ZNTs) in living systems occurs at a sub-second time scale yielding concentrations up to the low milimolar range (Kalivas, 2011). For that reason, the crucial requirement for in vivo monitoring besides non-invasiveness and good depth of penetration is high spatio-temporal resolution. These necessities are narrowing down the choice of potential monitoring techniques and due to an effective soft-tissue differentiation, magnetic resonance imaging (MRI) imposes as a modality of choice. An excellent example is the MRI method GluCEST (glutamate-based chemical exchange saturation transfer), which relays on proton exchange between the α-amino group of glutamic acid (Glu) and bulk water (Cai et al., 2012). In practice, this method allows detection of Glu levels across the human brain, which makes it suitable for diagnostic purposes of many pathological states (Roalf et al., 2017). However, it carries intrinsic disadvantages such as signal overlapping with resonances of other metabolites and inability to differentiate between intra and extracellular Glu (Bagga et al., 2018). Thus, one of the ways to overcome these obstacles is imaging living systems more specifically, e.g. via utilization of contrast-enhanced MRI (Lohrke et al., 2016). This approach implies the administration of paramagnetically labeled either extracellular or intracellular sensors able to alter the bulk proton signals upon interaction with targeted metabolites (Cao et al., 2013; Zhang et al., 2017).
ZNTs are rather structurally complex entities as guest molecules. An additional obstacle for their detection is the structural similarity amongst these amino acids, making their differentiation quite difficult. To this end, designing an artificial host capable to selectively recognize and effectively bind ZNTs would require a sophisticated binding platform in terms of structural complexity (Sim and Parker, 2015; Angelovski and Toth, 2017). A notable approach in this direction has been introduced by the Jasanoff lab, which utilized a genetically engineered protein platform that involved a large number of stabilizing interactions (Shapiro et al., 2010). Another approach to sense Glu grounded on the neurotransmitter mimicking strategy has been put forward by the Parker group, aiming to bypass the issues of direct interaction of a probe with Glu (Mishra et al., 2014).
The strategy that allows monitoring the change in longitudinal relaxation time (T1) which occurs due to host-guest association between the small-size molecule sensors and the ZNTs has been established within the common effort of our collaborators and us (Oukhatar et al., 2015a,b). This approach pioneered the ditopic MR hosts that contain a cavity-like binding pocket suited to simultaneously bind the negatively charged carboxylate and the positively charged α-ammonium cation of amino acid NTs (Steed and Atwood, 2009). The association with guests relies on the anticipated interaction between a coordinatively unsaturated Gd3+ with the carboxylate anion and weak hydrogen bonding of a crown ether moiety with the ammonium cation. In these initial studies, it was evidenced that the designed hosts do not distinguish between different ZNTs, making them suitable for functional imaging of neuronal activity that involves excretion of multiple amino acid neurotransmitters.
As a continuation of this study, our investigation focused on the structural optimization of a host binding pocket for sensing the total amount of ZNTs, in specific the most abundant Glu, Gly (glycine), and GABA (γ-aminobutyric acid). The molecular frame of the synthesized complexes consisted of a cyclen-based Gd-chelate as the MR reporting unit covalently linked by a propyl amide linker to the auxiliary 18-azacrown-6 binding site. Moreover, additional efforts were put to optimize the binding sites that would strengthen host-guest interactions. We structurally modified the binding sites, in an attempt to tune the framework to an optimal conformation of the binding pocket and adduct formation. Thus, we varied the size and charge of the substituents on both recognition sites using a set of chemical modifications. We investigated the host-guest binding strength through 1H longitudinal relaxivity measurements. Diffusion NMR and fluorescence spectroscopy techniques, as well as density functional theory (DFT) calculations were employed to obtain information at the molecular level. These were performed to highlight the mechanistic and structural aspects responsible for the recognition process, including the changes of the Ln3+ coordination environment and the effect of diverse remote groups on the solution structure of the free host and their association with structurally different guests.
Materials and Methods
Commercially available reagents and solvents were used without further purification. Procedures for preparation of compounds are provided in the Supplementary Material. Synthesized compounds were purified using silica gel 60 (0.03–0.2 mm) purchased from Carl Roth, Germany. Purification of compounds 11-16 was carried out by preparative TLC using Analtech 500 μm thick silica gel GF plates (20 × 20 cm, Uniplate) and visualized under UV light. For the purpose of 1H DOSY experiments, ligands L3 and L5 were purified using Preparative HPLC on a Varian PrepStar system equipped with the UV–Vis detector model 335 and a binary pump model SD-1 manual injector, controlled by Star chromatography workstation version 6.3 software. 1H, 13C, 31P, DOSY NMR and T1 measurements were performed at 298 K on Bruker AVANCE III 300 MHz spectrometer. The NMR spectra were recorded using either CDCl3 or D2O and referenced to TMS/TSP. Processing was performed using TopSpin 2.1 and the analysis using ACD/SpecManager 9.0 (Advanced Chemistry Development, Inc.). The concentration of the analyzed solutions containing paramagnetic Ln3+ ions was determined using the bulk magnetic susceptibility shift (BMS) method (Corsi et al., 2001). The luminescence emission and lifetime measurements were performed on QuantaMasterTM 3 PH fluorescence spectrometer from Photon Technology International, Inc. (USA). Low resolution mass spectra (ESI-MS) were recorded on ion trap SL 1100 Agilent system with an electrospray ionization source. High resolution mass spectra (ESI-HRMS) were performed on a Bruker BioApex II ESI-FT-ICR, equipped with an Agilent ESI source.
Neurotransmitter Binding Studies
1H longitudinal relaxivity titrations of GdL1−6 were performed at 298 K, pH 7.4 (50 mM HEPES) using standard inversion-recovery pulse sequence. The experiments were performed by stepwise addition of analyte to the solution of the analyzed complex (initial concentration of GdL1−5 and GdL6 were 3.0 and 1.0 mM, respectively) and the 1H longitudinal relaxation time T1, obs was measured after each addition. The longitudinal relaxivity r1 of bulk water protons was calculated according to Equation 1, where T1, obs is the measured value, T1, d represents longitudinal relaxation time (diamagnetic contribution) of the medium and [GdL] is the milimolar concentration of Gd3+.
The calculated r1 values were plotted as a function of added amount of analyte (0–100 equiv.) and fitted to the modified Michaelis-Menten model (Equation 2), where the slope of fitted curve is dissociation constant Kd and [Gd] is the initial concentration of gadolinium complex, while rf and rb are relaxivities of free and bound complexes, respectively. For the better comparison of host-guest binding strength, inverted values for Kd, the affinity constants Ka, were used.
The luminescence measurements with EuL3, 5 (5 mM) were performed in H2O and D2O (298 K, pH 7.4, 50 mM HEPES) before and after addition 50 equiv. of Gly. The emission spectra were obtained after excitation at 395 nm, while setting both the excitation and emission slit widths at 2 nm bandpass. The decays were measured at the emission wavelength 617 nm with a 10 μs resolution. Excitation and emission slits were set to 5 and 15 nm bandpass, respectively. Data sets are an average of 25 scans and each reported value is the mean of three independent measurements. The obtained curves were fitted to a first-order exponential decay. The q values were calculated using modified Horrock's model (Equation 3) (Dickins et al., 1996; Beeby et al., 1999), where n is the number of NH oscillators of coordinated amide groups.
DOSY experiments were performed on a Bruker AVANCE III spectrometer using a BBO probe with z gradients. The samples were dissolved in D2O at pD 7.8 and the data were acquired at 298 K. The viscosity correction for the media was performed using the stebpg1s pulse sequence and the obtained data were processed using the T1/T2 relaxation module. The diffusion measurement on paramagnetic complexes EuL3, 5 were conducted using the standard Bruker pulse sequence stebpg1s19 employing stimulated echo, bipolar gradients and Watergate solvent suppression. The gradient pulse time was set to 2.5 ms with a diffusion time of 200 ms. The gradient strength was increased linearly over 16 experiments from 2 to 95%. Processing was performed using Bruker Topspin software 2.1. The individual slices of pseudo-2D 1H diffusion spectra were phased, base corrected and the diffusion constants were read from the obtained spectra.
The hydrodynamic radii RH were calculated using the Einstein-Stokes relation, which models molecules as hydrodynamic spheres (Equation 4), where k is the Boltzmann constant, T is the absolute temperature, η is solvent viscosity and D is the measured diffusion coefficient.
All calculations were performed using the Gaussian 09 package (Revision E.01) (Frisch et al., 2009) within the hybrid meta-generalized gradient approximation (hybrid meta-GGA) with the TPSSh exchange-correlation functional (Tao et al., 2003). All calculations used the large-core quasi-relativistic effective core potential (LCRECP) approximation developed by Dolg and the corresponding [5s4p3d]-GTO valence basis set for Gd (Dolg et al., 1989), whereas the ligand atoms were described by using te standard 6–31G(d,p) basis set. The stationary points found on the potential-energy surfaces as a result of geometry optimizations were confirmed to correspond to energy minima rather than saddle points by using frequency analysis. Solvent effects were included by using the integral-equation formalism variant of the polarizable continuum model (IEFPCM) (Tomasi et al., 2005). Relative energies include zero point energy (ZPE) corrections obtained with frequency calculations.
Results and Discussion
Design of Ditopic Hosts GdL1−6
We implemented a diversity-oriented synthetic strategy to prepare a set of probes with structurally different pendant arms at both binding sites with the aim to intensify the interaction between the NH group of amino acid NTs and the free electronic pairs of the azacrown ether heteroatoms. To this end, we prepared six different Gd3+ complexes which differ in the structure of the MR reporting unit or the azacrown ether (Figure 1). Thus, L1 was designed by combining the common Gd-chelator based on DO3A (1,4,7,10-tetraazacyclododecane-1,4,7-tricarboxylic acid) with the diazacrown ether that bears a charge-neutral benzyl (Bn) group, whereas the remaining L2−6 contain a DO2A moiety (1,4,7,10-tetraazacyclododecane-1,7-dicarboxylic acid) as Gd-chelator and azacrown ethers with substituents of different charge and size L2−4, 6 or without substituent L5. The rationale behind design of L1, 2 lies on tuning the charge of the fragment containing the Gd3+ complex, which should impact the binding of ZNTs through their carboxylate group. The non-charged and small Bn substituent was expected to allow a facile approach of an ammonium cation. Ligands with negatively charged substituents, an acetate in L3 and a phosphonate in L4, were designed to exchange multiple bifurcated hydrogen bonds of an NH unit and increase the contribution of electrostatic interactions.
In the second approach, we designed L5, 6 to explore the effect of the substituent size on the conformational organization of the binding pocket, which is crucial for the recognition of the guest. Our previous studies evidenced that the carbonyl group of the propyl-amide linker is weakly coordinated to Gd3+, causing cyclen-crown sites to adopt a half opened shell-like conformation, which is well-suited for binding an amino acid functional groups (Oukhatar et al., 2015a,b). Arising from this, we designed L5 with no substituent in an attempt to favor strong amide linker coordination. On the contrary, a sterically bulky substituent methyl-7-methoxycoumarin (Mmc) was attached in the host L6, which could theoretically disturb this interaction due to increased steric constrains. Moreover, the strong fluorescence emission properties of the Mmc group enable additional interaction studies by means of the fluorescence spectroscopy. Last but not least, by combining different types of MR chelators (DO3A or DO2A) and substituents on the azacrown ethers, we set out to investigate the influence of the net charge of the complexes on their affinity toward the ZNTs.
Synthesis of Complexes LnL1−6
The preparation of ligands L1−6 relies on the convergent synthetic strategy, linking cyclen-derived bromides with 18-crown-6 ether fragments via direct N-alkylation. The synthesis of bromides 2 and 5 is summarized in Scheme 1. The bromo fragment 2 was prepared in 4 steps, starting from cyclen. In the first step, cyclen was converted in DO3A tert-butyl ester in a one-step procedure (Machitani et al., 2008), followed by the subsequent N-alkylation with N-(3-bromopropyl)phtalimide and the mild deprotection of Phth with ethylenediamine (EDA) to afford the propyl amine derivative 1. The amine 1 was coupled with bromoacetic acid in the presence of N, N'-dicyclohexylcarbodiimide (DCC) to furnish bromide 2. In a similarly fashion, the DO2A tert-butyl ester was synthesized in 3 steps according to a previously published procedure (Kovacs and Sherry, 1997). This chelator precursor was first functionalized with the Phth-protected propyl amine linker to yield 3. Subsequently, the second aromatic linker was introduced in the cyclen ring in the trans position relative to the propyl linker, and the Phth group was cleaved to yield the amine 4, which was subjected to the same treatment as the DO3A analog to give the desired DO2A-derived bromide 5.
Scheme 1. Synthesis of (A) DO3A-derived and (B) DO2A-derived bromide fragments 2 and 5, respectively. Reagents and conditions: (A) (i) N-(bromopropyl)phthalimide, MeCN, K2CO3, 45°C; (ii) i-PrOH, EDA, 60°C, 75% (over two steps); (iii) CH2Cl2, DCC, bromoacetic acid, 75%; (B) (i) N-(bromopropyl)phthalimide, MeCN, NaHCO3, 3 days, 65%; (ii) 2-bromo-N-(4-nitrophenyl)acetamide, MeCN, K2CO3; (iii) i-PrOH, EDA, 65% (over two steps); (iv) CH2Cl2, DCC, bromoacetic acid, 52%.
The synthesis of the azacrown fragments 6, 7, 9, and 10 commenced from the common precursor 1,10-diaza-18-crown-6-ether, attaching substituents via direct N-alkylation with the corresponding alkyl halides, tert-butyl bromacetate, benzyl bromide, and Br-Mmc, respectively Scheme 2. The adduct 8 was generated under the Kabachnik-Fields reaction conditions, followed by the catalytic debenzylation in a Parr hydrogenator to afford the final mono N-phosphonated derivative 9.
Scheme 2. Synthesis of 18-azacrown-6-ether fragments 6, 9, and 10. Reagents and conditions: (i) tert-butyl bromoacetate (0.95 equiv), CHCl3, Et3N, RT, 48 h, 62%; (ii) BnCl (0.95 equiv), CHCl3, Et3N, RT, 60 h, 58%; (iii) (C2H5O)2P(H)O, CH2O, CH2Cl2, BF3·Et2O, anh. Na2SO4, 65°C, 60%; (iv) H2, Pd(OH)2/C, EtOH, 3 bar, 94%; (v) Br-Mmc, MeCN, NaHCO3, RT, 40%.
Once the precursors of two binding sites were prepared, they were conjugated into the bis-macrocyclic constructs via direct N-alkylation Scheme 3. In general, bridging bromides 2 and 5 with the azacrown ether fragments required over 3 days stirring in dry DMF and N2 atmosphere to acquire 11–16 in acceptable to good yields (25–75%), due to slow kinetics and tendency of bromides to undergo hydrolysis. Subsequently, the obtained tert-butyl esters 11–13, 15, and 16 were hydrolyzed with formic acid at 60°C over 24 h into the final ligands L1−3, 5, 6. Deprotection of 14 required first transesterification with BrSiMe3, followed by hydrolysis of trimethylsilyl esters to give L4. The purity of ligands L1−6 was determined with respect to nitrogen content obtained in CHN analysis and the amount of corresponding LnCl3 salt for complexation was calculated accordingly in stoichiometric ratio of 1:1. Finally, the complexation of L1−6 was performed in water, maintaining pH ~ 7.2 at 60°C for 16 h to ensure formation of the thermodynamically favored cyclen in-cage coordination and yield the final paramagnetic complexes.
Scheme 3. Synthesis of (A) DO3A derived ligand L1, (B) DO2A derived ligands L2−6. Reagents and conditions: (A) (i) dry DMF, N2, NaHCO3, 3 days, 70%; (ii) HCO2H, overnight, RT, 76%. (B) (i) dry DMF, N2, NaHCO3, 3 days, 25–78%; (ii) HCO2H, 32 h, 45°C; (iii) TMBS (16 equiv), dry DCM, N2; (iv) H2O, 50%.
The relaxometric efficacy of the synthesized complexes GdL1−6 was evaluated in T1-titration experiments with the three most abundant ZNTs: Gly, Glu and GABA, the non-zwitterionic neurotransmitter acetylcholine (ACh) and the competitive bicarbonate anion. As well-established, the main contributing factor to the 1H longitudinal relaxivity r1 is the hydration state q of the paramagnetic label (Merbach et al., 2013). Thus, upon addition of the analyte (A), the hydrated chelate GdL(H2O)q associates into the ternary adduct GdL(A), in which the coordinated water molecule(s) are displaced to a different extent. This results in longer longitudinal relaxation times (T1) of bulk water protons, that is, reduced longitudinal relaxivity r1. Furthermore, plotting the r1 values vs. the amount of the added analyte and fitting the obtained curves to a modified Michaelis-Menten mathematical model provides the affinity constants for the formed ternary complexes.
The r1 profiles of GdL1−6 upon the titration with Gly revealed that patterns in relaxometric response are associated with the nature of substituents at the recognition moieties (Figure 2). More specifically, the observed downward trends are controlled by the size and charge of the pendant arm at the auxiliary receptor, according to which complexes were classified in three groups, whilst the dominant factor that determined affinity is found to be the coordination environment on the reporting moiety (Table 1). Specifically, the first group of complexes consists of those which bear the non-charged and sterically small benzyl substituent GdL1, 2 or no substituent GdL5 at the auxiliary moiety. The features that characterize this group are very high initial r1 values that exceeded well over 9.0 mM−1 s−1, and a decrease upon saturation with Gly to 30, 31, and 36% of the starting r1 value, respectively. The second group of complexes are those containing negatively charged pendants GdL3, 4, which displayed relatively high initial relaxivities, 6.88 and 7.14 mM−1 s−1, respectively. However, the relaxivity drop until the saturation point is halved when compared with the first group: addition of Gly lowered the initial r1 by 30 and 37%, respectively. The last complex GdL6, characterized by the sterically bulky and non-charged coumarin group Mmc, exhibited the lowest initial relaxivity of 6.08 mM−1 s−1, but a comparable decrease in r1 with GdL1, 2, 5 of 63% when saturated with Gly.
Figure 2. 1H longitudinal relaxometric titration curves for GdL1−6 with Gly ([GdL1−5] = 3.0 mM, pH 7.4 or [GdL6] = 1.0 mM, pD 7.8, 50 mM HEPES, 298 K, B0 = 7.05 T). The lines represent the fits as explained in the text.
Table 1. Binding affinity constants and relaxivity decrease for complexes GdL1−6 upon saturation with Gly [a][b].
Regarding the affinities toward Gly, the investigated complexes showed slightly different trends with respect to their amplitudes of r1 change (Table 1, Figure 3, and Figure S1). The highest affinities falling in the range 69–75 M−1 were observed for GdL2, 5, 6, the complexes that bear Gd-DO2A as the MR reporting unit and no or non-charged Bn and Mmc substituents on the azacrown ether moiety. The DO3A analog of GdL2, namely complex GdL1, exhibited a lower Ka value, indicating that shifting the net charge of the MR reporting moiety from positive to neutral logically reduced the affinity toward the negatively charged carboxylate group of the ZNTs. On the other hand, the negatively charged substituents on the azacrown ether side in the two DO2A-derived complexes GdL3, 4 did not compensate for their size and potentially additional, undesired interaction with the paramagnetic center (vide infra). Indeed, this group of complexes exhibited the lowest affinities, whereby the more negatively charged phosphonate group influenced GdL4 to show decreased affinity toward Gly to the greater extent and the lowest Ka value of all investigated systems.
Figure 3. 1H longitudinal relaxometric titration curves for GdL2, 3, 6 with different analytes A ([GdL2, 3] = 3.0 mM or [GdL6] = 1.0 mM, 298 K, pH 7.4, 50 mM HEPES, B0 = 7.05 T). The lines represent the fits as explained in the text.
In terms of specificity, it is found that the hosts GdL1, 2 displayed greater affinity toward Glu over other ZNTs and hydrogencarbonate giving values 72 and 311 M−1, respectively, while all the synthesized complexes were ACh-insensitive (Table S2). The improved affinity of GdL1 and GdL2 toward Glu over the other ZNTs is probably related to the longer chain length of the neurotransmitter combined to the favorable effect of the Bn pendant, which introduces low steric hindrance and forces the 18-crown-6 moiety into a conformation that is more suited to accommodate the ammonium cation. It was also noticed that GdL6 exhibited a higher affinity toward sterically less demanding Gly over GABA and Glu. On contrary, the observed diversity in affinity patterns for DO2A derivatives with different substituents on the azacrown ether moiety suggested the synergistic effect of several factors on the recognition process. More specifically, the milder downward trends and lesser change in r1 for GdL3, 4, as well as the lower initial relaxivity for GdL6, indicated that either increase in remote substituent size or the net negative charge have a major impact on binding ZNTs.
We performed luminescence studies on EuL3 and EuL5 to examine their solution structures and the way that auxiliary pendants affect binding with ZNTs (Binnemans, 2015). In the first step, the luminescence emission spectra of both complexes before and upon addition of 50 equiv. Gly were measured. The EuL3(H2O)q complex showed more intense emission than EuL5(H2O)q, but this situation is reversed upon saturation with Gly, with the most intense emission being observed for the ternary adduct EuL5(Gly) (Figure 4). The investigations carried out on similar host-systems evidenced that association of ZNTs with Eu3+ labeled complexes induce hyperchromicity of the hypersensitive 5D0→7F2 transition at ca. 615 nm (Oukhatar et al., 2015b). The intensity increase is related with the coordination of negatively charged carboxylate that displaces the water molecule coordinated directly to Eu3+, thus reducing the quenching efficiency of the Eu3+ 5D0 state by water O-H oscillators (Beeby et al., 1999). Therefore, we next measured the radiative rate constants characterizing the depopulation of the 5D0 excited state in H2O and D2O solutions and the hydration numbers were calculated for both hydrated complexes EuL3, 5(H2O)q and their adducts EuL3, 5(Gly) (Dickins et al., 1996). The calculated q value for EuL5(H2O)q yielded 1.76, which is twice higher than that of EuL3(H2O)q; furthermore, the addition of Gly expectedly decreased the hydration numbers of both complexes (Table 2). The greater drop (Δq = 0.87) occurred for EuL5, pointing out that H2O is displaced more readily than in EuL3, in which change in hydration is almost negligible. The hydration number of EuL5 suggests the presence of two water molecules coordinated to the metal ion, which would imply that the oxygen atom of the propylamide group is not involved in coordination to the metal ion. This is in contrast with previous findings when investigating structurally related systems (Oukhatar et al., 2015b). These results indicate that slight modifications of the ligand structure may have an important impact in the structure of the complexes.
Figure 4. Metal-centered luminescence emission spectra for EuL3, 5 before and after addition 50 equiv. of Gly (λex = 395 nm, [EuL3, 5] = 5 mM, pH 7.4, 50 mM HEPES, 298 K).
Table 2. Luminescence lifetimes of Eu3+ measured at 616 nm from the decay profile of the 5D0 excited state (H2O and D2O) and calculated hydration numbers q [a].
The metal-centered emission spectra revealed diverse coordination environments for the hydrated forms of the complexes (Binnemans, 2015). Namely, EuL3(H2O)q exhibited larger intensity ratio of the magnetic-dipole transition 5D0→7F1, which is insensitive of the local coordination environment, to the hypersensitive 5D0→7F2 transition, which is heavily influenced by the local symmetry of Eu3+ and the nature of ligand atoms (Table S3). These data suggested a solution form for EuL3(H2O)q in which a more polarizable ligand atom than water oxygen is coordinated to Eu3+, partly inhibiting the binding of Gly. Relying on these observations, we hypothesized that the acetate pendant of the lariat ether fragment in EuL3(H2O)q interacts weakly in intra- or intermolecular fashion with the paramagnetic label, favoring a form in which the lanthanide ion is shielded, and thus impedes water coordination. If this is the case, the more polarizable Gly-carboxylate ligand replaces the azacrown ether-acetate pedant by disrupting this interaction, which should induce major conformational changes in EuL3 when Gly is added. However, the exchange of the ligand acetate group by Gly in the coordination sphere of Eu3+ excludes water involvement, resulting in the small changes observed by time-resolved emission decay measurements.
Furthermore, EuL6 appended with the coumarine derived fluorophore, was characterized prior and upon the interaction with Gly by means of fluorescence spectroscopy. However, the complex turned out to be non-emissive when the metal center was excited directly or via the chromophore. Moreover, the fluorophore on EuL6 exhibited strong fluorescence emission intensity at 420 nm when excited at 330 nm, typical for coumarine derivatives (Figure S2). The steady-state fluorescence emission experiments were repeated also in the presence of Gly or Glu, resulting in spectra of similar intensity as with EuL6 alone. Hence, this methodology could not provide further evidences related to interaction of the 18-azacrown-6 ether moiety with NH groups of ZNTs.
Diffusion-ordered 2D NMR spectroscopy (DOSY) allows separation of NMR signals based on the molecular diffusion. Owing to the correlation between self-diffusion and the size/shape of molecular species, represented by Stokes-Einstein model, 1H DOSY-NMR is a readily used technique to disperse 1H resonances along the diffusion axis, relating signals with differently organized forms of a solute (Morris et al., 1994). To verify the hypothesis brought upon luminescence measurements, we performed 1H NMR diffusion measurements on paramagnetic complexes EuL3, 5 and their ternary adducts EuL3, 5(Gly), and compared the calculated hydrodynamic radii (RH) assuming spherical shape of species (Table 3) (Natrajan et al., 2007; Maggini et al., 2012; Denis-Quanquin et al., 2016). Firstly, the 1H DOSY-NMR experiments on aqueous solution of complexes EuL3 and EuL5 were recorded (Figure 5). The obtained pseudo 2D spectra for EuL3(H2O)q identified two distinct sets of 1H resonances with corresponding different diffusion constants (D), while the EuL5(H2O) complex resonated along one diffusion value. Following this step, the diffusion measurements were repeated after 50 equiv. of Gly were added, revealing disappearance of the resonances corresponding to the faster diffusing form in EuL3(Gly).
Table 3. Diffusion coefficients obtained in 1H DOSY NMR experiments and calculated hydrodynamic radii ([EuL3, 5] = 10 mM, D2O as solvent, pD = 7.8, 25°C) [a][b].
Figure 5. 1H DOSY NMR of EuL3 (top) and EuL5 (bottom) in D2O before (left) and after (right) addition of 50 equiv. Gly.
On the contrary, no changes in diffusion of the slower diffusing form in EuL3(Gly) or diffusion of EuL5(Gly) were observed (Figure 5). These results indicate that two solution forms of EuL3(H2O)q are in dynamic equilibria, in which the form with smaller hydrodynamic radius underwent structural changes in the presence of Gly. Transformation of this form into one with the greater hydrodynamic radius, likely eliminates the probability of EuL3(H2O)q intermolecular self-organization into dimeric and oligomeric forms. Instead, the formation of EuL3 species with smaller diameter that disappear with Gly addition rather suggests a weak intramolecular self-assembly, i.e., coordination of the azacrown acetate pendant to Eu3+. This interaction brings closer the two binding moieties toward each other, thus reducing RH. Once Gly approaches the coordination sphere of the paramagnetic metal ion, the acetate pendant is displaced from Eu3+ with the carboxylate of Gly, allowing formation of ternary complex EuL3(Gly) (Figure 6).
Figure 6. Schematic presentation of dynamic equilibrium for the hydrated complex EuL3 and the proposed binding mechanism with ZNTs.
These observations are in good agreement with the results obtained in the relaxometric titrations and luminescence emission experiments (vide supra): existence of the dynamic equilibrium of non-folded and folded forms in EuL3 and GdL3 results in lower initial hydration number q and consequently lower initial r1 value prior any interaction with Gly (or ZNT in general). In such case, the observed r1 changes after addition of Gly occur only due to the displacement of coordinated water on the species without the intramolecular self-assembly. In contrast to the GdL3 system, GdL5 exists only in the non-folded conformation; here, the Gly/ZNTs addition causes greater change in hydration number and hence greater drop in r1.
Theoretical calculations were carried out to rationalize the different behavior of the GdL3 and GdL5 complexes and provide support to the interpretation of relaxometric, luminescence and DOSY experiments. We first analyzed whether the involvement of the carboxylate group attached to the crown moiety in GdL3 in intramolecular coordination to the metal ion is feasible. We thus modeled the GdL3(H2O)2·4H2O system, which included two water molecules coordinated to the metal ion and four additional second-sphere water molecules. The explicit inclusion of a few second-sphere water molecules was done because it was previously found to be essential for providing a good description of the Gd-Owater distances and the spin densities at the 17O nucleus (Esteban-Gomez et al., 2012; Regueiro-Figueroa and Platas-Iglesias, 2015). The calculation provided an energy minimum that contained two water molecules coordinated to the Gd3+ ion, which presents the typical capped square antiprismatic (SAP) coordination environment of this type of complexes. Next, we explored the potential energy surface by approaching the remote carboxylate group to the metal ion, which eventually leads to a second energy minimum in which the carboxylate group coordinates to Gd3+ by replacing a coordinated water molecule (Figure 7). Indeed, these calculations indicated that the coordination of the remote carboxylate group is feasible. Moreover, DFT predicts a very small energy difference between the two energy minima (folded and non-folded species), favoring the open form by only 1.4 kJ mol−1. Calculations performed on the GdL(H2O)·5H2O system (Figure S12) also suggest that coordination of the phosphonate group attached to the crown moiety is likely responsible for the lower binding affinity of this complex toward ZNTs.
Figure 7. Structures of the GdL3(H2O)2·4H2O system obtained with DFT calculations showing the coordination of the carboxylate group of the crown moiety to the metal ion (left) and the open form (right).
We also performed calculations on the GdL5(H2O)2·4H2O and GdL5(Gly)(H2O)·5H2O systems to gain insight into the molecular recognition process. The calculations were done using three model systems with the Gly guest interacting with either the crown moiety or the lanthanide ion, and a third form in which simultaneous recognition by the two binding units is taking place (Figure 8 and Figures S10, S11). The corresponding energy differences provided an estimate of the contributions to the overall recognition process of the coordination to the metal ion and the hydrogen bonding interaction with the crown moiety, which turned out to be 17.3 and 9.7 kJ mol−1. Thus, these calculations indicate that binding to the coordinatively unsaturated GdDO2A unit represents the major driving force of the molecular association process, though the hydrogen-bonding interaction with the crown moiety also provides a sizeable contribution to the overall binding energy.
Figure 8. Structure of the GdL5(Gly)(H2O)·5H2O system obtained with DFT calculations showing the simultaneous binding of Gly to the crown moiety and the Gd3+ ion.
In the present work, we have expanded the scope of Gd3+ based DO3A/DO2A ditopic MR probes sensitive to zwitterionic neurotransmitters (ZNTs) and revealed chemical structure-affinity relationships through the relaxometric response. The results indicated that the lanthanide coordination environment is highly sensitive to influences of the pedant arms on the azacrown ether moieties, resulting in significant differences on both initial relaxivity and the total relaxivity change. Furthermore, the charge neutral DO3A derivative GdL1 displayed reduced affinity over the DO2A analog GdL2, while adding the negative charge on the azacrown ethers reduced the affinity toward guests. Indeed, it is found that the relaxometric response of all derivatives with the same MR reporting unit GdL2−6 is governed by the nature of the auxiliary pendant. More specifically, the substituent at the azacrown ether fragment favors the disposition of binding sites in a charge- and size-dependent manner, thus determining both relaxivity change and binding affinity of the paramagnetic host. The observed r1 patterns revealed that increasing in size of the pendant elevates steric hindrances, resulting in lower initial relaxivities while only slightly diminishing the strength of interactions with guest molecules. Besides steric hindrance, electronic factors in the complexes with negatively charged pendants GdL3, 4 unfavored host-guest association to the greater extent, likely due to the coordination of the remote carboxylate/phosphonate groups to the metal ion. In addition, an achievement of great importance evidenced here is that GdL1−2 displayed higher affinity toward Glu over competitive hydrogencarbonate. The findings reported in this work present a breakthrough in understanding the nature of interactions between the ditopic MR platforms consisted of DO3A-/DO2A-derived units and azacrown ethers with ZNTs. Owning the ability to produce an excellent change in relaxivity, these MR sensors present a great potential for further development of bioresponsive probes that can visualize changes in concentration of the most abundant ZNTs by means of T1-weigthed MRI. Additionally, the results of this study can be used as guideline for development of novel concepts within functional MRI for sensing chemical neurotransmission.
The summary of data is available in the standard reporting form, while the raw data supporting the conclusions of this manuscript will be made available by the authors. Requests to access the datasets should be directed to email@example.com.
GA designed the targeted molecules and supervised research. ÐT synthesized and evaluated the compounds, analyzed the data, wrote the manuscript. CP-I performed DFT calculations and wrote the section concerning computational calculations. All authors critically reviewed the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors gratefully acknowledge the financial support of the German Research Foundation (DFG, grant AN 716/7-1). We thank Dr. Vincent Truffault for assistance with performing diffusion NMR measurements. The authors thank Centro de Supercomputación de Galicia (CESGA) for providing the computer facilities.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00490/full#supplementary-material
Bagga, P., Pickup, S., Crescenzi, R., Martinez, D., Borthakur, A., D'Aquilla, K., et al. (2018). In vivo GluCEST MRI: reproducibility, background contribution and source of glutamate changes in the MPTP model of Parkinson's disease. Sci. Rep. 8:2883. doi: 10.1038/s41598-018-21035-3
Beeby, A., Clarkson, I. M., Dickins, R. S., Faulkner, S., Parker, D., Royle, L., et al. (1999). Non-radiative deactivation of the excited states of europium, terbium and ytterbium complexes by proximate energy-matched OH, NH and CH oscillators: an improved luminescence method for establishing solution hydration states. J. Chem. Soc., Perkin Trans. 2, 493–503. doi: 10.1039/a808692c
Cao, C.-Y., Shen, Y.-Y., Wang, J.-D., Li, L., and Liang, G.-L. (2013). Controlled intracellular self-assembly of gadolinium nanoparticles as smart molecular MR contrast agents. Sci. Rep. 3:1024. doi: 10.1038/srep01024
Corsi, D. M., Platas-Iglesias, C., van Bekkum, H., and Peters, J.A. (2001). Determination of paramagnetic lanthanide(III) concentrations from bulk magnetic susceptibility shifts in NMR spectra. Magn. Reson. Chem. 39, 723–726. doi: 10.1002/mrc.922
Denis-Quanquin, S., Riobe, F., Delsuc, M. A., Maury, O., and Giraud, N. (2016). Paramagnetic DOSY: an accurate tool for the analysis of the supramolecular interactions between lanthanide complexes and proteins. Chem. Eur. J. 22, 18123–18131. doi: 10.1002/chem.201603460
Dickins, R. S., Parker, D., deSousa, A. S., and Williams, J.A.G. (1996). Closely diffusing O-H, amide N-H and methylene C-H oscillators quench the excited state of europium complexes in solution. Chem. Commun. 6, 697–698. doi: 10.1039/cc9960000697
Esteban-Gomez, D., de Blas, A., Rodriguez-Blas, T., Helm, L., and Platas-Iglesias, C. (2012). Hyperfine coupling constants on inner-sphere water molecules of Gd-III-based MRI contrast agents. Chem. Phys. Chem. 13, 3640–3650. doi: 10.1002/cphc.201200417
Lohrke, J., Frenzel, T., Endrikat, J., Alves, F. C., Grist, T. M., Law, M., et al. (2016). 25 Years of contrast-enhanced MRI: developments, current challenges and future perspectives. Adv. Ther. 33, 1–28. doi: 10.1007/s12325-015-0275-4
Machitani, K., Sakamoto, H., Nakahara, Y., and Kimura, K. (2008). Molecular design of tetraazamacrocyclic derivatives bearing a spirobenzopyran and three carboxymethyl moieties and their metal-ion complexing behavior. Anal. Sci. 24, 463–469. doi: 10.2116/analsci.24.463
Maggini, L., Toma, F. M., Feruglio, L., Malicka, J. M., Da Ros, T., Armaroli, N., et al. (2012). Luminescent blooming of dendronic carbon nanotubes through ion-pairing interactions with an EuIII complex. Chem. Eur. J. 18, 5889–5897. doi: 10.1002/chem.201200237
Mishra, A., Mishra, R., Gottschalk, S., Pal, R., Sim, N., Engelmann, J., et al. (2014). Microscopic visualization of metabotropic glutamate receptors on the surface of living cells using bifunctional magnetic resonance imaging probes. ACS Chem. Neurosci. 5, 128–137. doi: 10.1021/cn400175m
Morris, K. F., Stilbs, P., and Johnson, C.S. (1994). Analysis of mixtures based on molecular size and hydrophobicity by means of diffusion-ordered 2D NMR. Anal. Chem. 66, 211–215. doi: 10.1021/ac00074a006
Oukhatar, F., Même, S., Même, W., Szeremeta, F., Logothetis, N. K., Angelovski, G., et al. (2015a). MRI sensing of neurotransmitters with a crown ether appended Gd3+ complex. ACS Chem. Neurosci. 6, 219–225. doi: 10.1021/cn500289y
Oukhatar, F., Meudal, H., Landon, C., Logothetis, N. K., Platas-Iglesias, C., Angelovski, G., et al. (2015b). Macrocyclic Gd3+ Complexes with pendant crown ethers designed for binding zwitterionic neurotransmitters. Chem. Eur. J. 21, 11226–11237. doi: 10.1002/chem.201500542
Regueiro-Figueroa, M., and Platas-Iglesias, C. (2015). Toward the prediction of water exchange rates in magnetic resonance imaging contrast agents: a density functional theory study. J. Phys. Chem. A 119, 6436–6445. doi: 10.1021/acs.jpca.5b01728
Roalf, D. R., Nanga, R. P. R., Rupert, P. E., Hariharan, H., Quarmley, M., Calkins, M. E., et al. (2017). Glutamate imaging (GluCEST) reveals lower brain GluCEST contrast in patients on the psychosis spectrum. Mol. Psychiatry 22:1298. doi: 10.1038/mp.2016.258
Shapiro, M. G., Westmeyer, G. G., Romero, P. A., Szablowski, J. O., Kuster, B., Shah, A., et al. (2010). Directed evolution of a magnetic resonance imaging contrast agent for noninvasive imaging of dopamine. Nat. Biotechnol. 28, 264–270. doi: 10.1038/nbt.1609
Tao, J. M., Perdew, J. P., Staroverov, V. N., and Scuseria, G. E. (2003). Climbing the density functional ladder: Nonempirical meta-generalized gradient approximation designed for molecules and solids. Phys. Rev. Lett. 91:146401. doi: 10.1103/Physrevlett.91.146401
Zhang, L., Martins, A. F., Mai, Y., Zhao, P., Funk, A. M., Clavijo Jordan, M. V., et al. (2017). Imaging extracellular lactate in vitro and in vivo using CEST MRI and a Paramagnetic Shift Reagent. Chem. Eur. J. 23, 1752–1756. doi: 10.1002/chem.201604558
Keywords: crown ethers, gadolinium, magnetic resonance imaging, neurotransmitters, zwitterions
Citation: Toljić Ð, Platas-Iglesias C and Angelovski G (2019) In-depth Study of a Novel Class of Ditopic Gadolinium(III)-based MRI Probes Sensitive to Zwitterionic Neurotransmitters. Front. Chem. 7:490. doi: 10.3389/fchem.2019.00490
Received: 27 March 2019; Accepted: 25 June 2019;
Published: 24 July 2019.
Edited by:Wen-Sheng Chung, National Chiao Tung University, Taiwan
Reviewed by:Apurba Lal Koner, Indian Institute of Science Education and Research, India
Ken Cham-Fai Leung, Hong Kong Baptist University, Hong Kong
Copyright © 2019 Toljić, Platas-Iglesias and Angelovski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Goran Angelovski, firstname.lastname@example.org