Samantha E.
Bodman
,
Colum
Breen
,
Felix
Plasser
and
Stephen J.
Butler
*
Department of Chemistry, Loughborough University, Epinal Way, Loughborough, LE11 3TU, UK. E-mail: s.j.butler@lboro.ac.uk
First published on 30th August 2022
The selective recognition of anions in water by artificial receptors remains a significant research challenge. The creation of a receptor selective for adenosine monophosphate (AMP) is particularly difficult due to its similarity in structure with the more negatively charged anions, ADP and ATP. We recently developed a macrocyclic Eu(III) complex that selectively binds AMP in water, by incorporating a sterically demanding quinoline arm that inhibits coordination of commonly interfering anions such as ATP. A phenylboronic acid motif was installed within the ligand to engage the ribose sugar of AMP through reversible covalent bonds. Herein we report two new Eu(III) complexes, [Eu·oBOH2]+ and [Eu·pBOH2]+, to investigate the impact of varying the position of the phenylboronic acid group on the anion binding properties of the Eu(III) receptors. We found that [Eu·pBOH2]+ showed preferential binding to AMP over ATP, but exhibits a lower level of discrimination between AMP and ADP compared with the isomeric complex [Eu·mBOH2]+. Surprisingly, [Eu·oBOH2]+ showed no response to anions but displayed a unique response to pH, ascribed to the direct coordination of the ortho-boronate ester to the Eu(III) centre. Finally, we present first principles computations that offer a promising approach to access the emission spectra of lanthanide complexes, aiding the design of responsive lanthanide probes with specific photophysical properties.
Adenosine monophosphate (AMP) is an important target for detection because of its crucial roles in biological processes, including maintaining energy homeostasis and mediating signal transduction events. However, receptors that selectively bind AMP in water are rare,18–20 due to the tendency of receptors to form stronger electrostatic interactions with substrates with higher negative charges, in particular ADP and ATP. Indeed, the majority of receptors for nucleoside phosphate anions display the selectivity trend ATP > ADP > AMP21–24 and reversing this trend is very difficult. AMP sensing has been reported through water soluble coordination complexes of zinc(II),25 copper(I)26 and europium(III),18,20 as well as gold(II) nanoparticles.27 In a few cases, AMP recognition has been achieved by exploiting the structural aspects of the organic receptors to provide highly complementary sites for AMP binding, including strong electrostatic contacts supported by π–π interactions with the nucleobase.28,29
We recently established a method to selectively bind and sense AMP using a stable Eu(III) complex based on a sterically demanding macrocyclic ligand.20 Two complexes were synthesized, [Eu·Bn]+ and [Eu·mBOH2]+ (Fig. 1a),20 each featuring a bulky 8-(benzyloxy)quinoline pendant arm that coordinates to the Eu(III) ion in a bidentate manner, creating a single coordination site for the monodentate binding of the phosphate component of AMP. Crucially, oxyanions that prefer a chelation binding mode including ATP, bicarbonate and lactate, do not bind to the complexes due to the steric hindrance imposed by the ligand. A boronic acid motif was introduced at the meta-position of the phenyl ring of [Eu·mBOH2]+ to interact reversibly with the 1,2-diol moiety of the ribose sugar. Phenylboronic acids have been utilized in a range of receptors for sensing of saccharides with great success, due to their ability to form reversible covalent bonds with 1,2-diols.30–32
Fig. 1 Family of cationic Eu(III) complexes designed to bind AMP in water. (a) Previously reported complexes [Eu·Bn]+ and [Eu·mBOH2]+.20 (b) New complexes [Eu·oBOH2]+ and [Eu·pBOH2]+ and (c) structures of anions investigated in this work. |
Computations on lanthanide complexes are a powerful means to obtain deeper insight into the structural arrangement of the ligands, the binding modes of the guest molecules, and the properties of the f-electrons. However, these computations can be highly challenging due to the presence of unpaired electrons as well as the importance of solvation effects and conformational flexibility. Well-established protocols for obtaining accurate structures using density functional theory (DFT) are in use, relying on an effective core potential for the f-electrons and implicit solvation models.33–36 Conversely, computing anion binding energies in solution37 is a highly challenging endeavour and is therefore not routinely performed for lanthanide complexes. As opposed to the described DFT simulations, in which the f-electrons are only implicitly considered within the core potential, it would be necessary to include them explicitly for an ab initio computation of emission spectra or magnetic properties. Indeed, magnetic properties of lanthanide complexes have been modelled successfully38–40 using multireference ab initio theory.41 However, to the best of our knowledge, no first principles simulation of lanthanide emission spectra has been reported so far. Emission spectra are usually computed using phenomenological models based on Judd–Ofelt theory.42,43
Herein, we present two new cationic Eu(III) complexes, [Eu·oBOH2]+ and [Eu·pBOH2]+ (Fig. 1b) bearing a boronic acid motif at the ortho- and para-position of the phenyl ring, respectively, allowing us to evaluate the optimal geometry for AMP recognition in water through comparison with the previously reported meta-substituted complex [Eu·mBOH2]+. The ligand design retains sufficient host rigidity and steric hindrance at the metal centre to ensure monodentate binding of AMP, whilst preventing the chelation of common oxyanion interferents. The complexes are studied experimentally using a combination of photophysical techniques, mass spectrometry and NMR spectroscopy. Computations are performed to further investigate binding geometries and energies. We explore the possibility of computing lanthanide emission spectra using a first principles approach.
Complex | λ max/nm | ε/M−1 cm−1 | τ H2O /ms | τ D2O /ms | q | Φ em /% |
---|---|---|---|---|---|---|
Values of hydration state q (calculation error is ±20%) were derived using literature methods.47 Quantum yields were measured using quinine sulfate in 0.05 M H2SO4 as standard (Φem = 59%).48 Calculation error in quantum yields is ±15%.a Mean ± standard deviation for two independent measurements.b Emission lifetimes and hydration state in parentheses were determined at pH 4.0.c Mean ± standard deviation for three independent measurements.d [Eu·mBOH2]+ data from ref. 20 shown for comparison. | ||||||
[Eu·oBOH2]+ | 322 | 3800 | 0.422 ± 0.002 (0.210)b | 0.532 ± 0.004 (0.281)b | 0.3 (1.1)b | 9.0 ± 0.7 |
[Eu·pBOH2]+ | 322 | 2900 | 0.192 ± 0.037 | 0.260 ± 0.069 | 1.2 | 2.0 ± 0.3 |
[Eu·mBOH2]+d | 322 | 2700 | 0.199 ± 0.004 | 0.283 ± 0.002 | 1.5 | 1.5 ± 0.2 |
Each complex displays characteristic Eu(III) emission in the red region of the visible spectrum when excited at 322 nm. The emission spectrum of [Eu·pBOH2]+ (Fig. 2b) displays two distinguishable components in the ΔJ = 1 (585–600 nm) emission band and three components in the ΔJ = 2 (605–630 nm) band. Analysis of the intensity ratio of the ΔJ = 2/ΔJ = 1 emission bands gave the value of approximately 1.4, similar to that seen for complex [Eu·mBOH2]+ (Fig. S4†). A rather different emission spectrum was obtained for [Eu·oBOH2]+ (Fig. 2a), characterised by a dominant ΔJ = 2 emission band and a higher intensity ratio ΔJ = 2/ΔJ = 1 of 3.0, signifying a different coordination environment at the Eu(III) centre.45,46 Excitation spectra collected for each complex were similar to the absorption spectra, confirming that the quinoline antenna is responsible for sensitising the Eu(III) emission (Fig. S5†).
Fig. 2 Emission spectra of (a) [Eu·oBOH2]+ and (b) [Eu·pBOH2] + in 10 mM HEPES at pH 7.0 and 295 K, λexc = 322 nm. |
The quantum yields of the metal-centred luminescence of [Eu·pBOH2]+ and [Eu·oBOH2]+ were found to be 2.0 and 9.0, respectively. The ortho-complex was notably more emissive, and the emission lifetimes in H2O and D2O were also considerably longer, at 0.422 ms and 0.532 ms respectively (Table 1). These differences correspond to a lack of coordinated water (q = 0) in the first coordination sphere of [Eu·oBOH2]+, whereas the number of coordinated water molecules calculated for [Eu·pBOH2]+ is q = 1. The higher emission intensity, distinct spectral shape and longer emission lifetime of [Eu·oBOH2]+ all point to a distinctive conformation that potentially involves an interaction between the ortho-substituted boronic acid and the Eu(III) ion. We corroborate this hypothesis via pH analysis as well as computed binding free energies and emission spectra.
Fig. 3 (a) pH dependent hybridisation of phenylboronic acids, switching from sp2 hybridised boron at low pH to sp3 hybridised at high pH. The pKa value corresponds to the hydroxylation of the boron atom accompanied by the release of a proton;51 (b) structural representation of [Eu·oBOH2]+ showing coordination of the boronic acid with increasing pH, shifting the boron hybridisation from sp2 (left) to tetrahedral sp3 (right); pH profiles of (c) [Eu·pBOH2]+ and (d) [Eu·oBOH2]+ measured in water at 295 K, showing the plot of emission intensity ratio ΔJ = 2/ΔJ = 1 as a function of pH and the pKa determined by nonlinear least squares regression analysis; (e) DFT-optimised structure of [Eu·oBOH3]. |
A much lower pKa of 4.47 ± 0.02 was found for [Eu·oBOH2]+ (Fig. 3d). We postulate that the ortho-substituted boronic acid coordinates to the Eu(III) ion, increasing the Lewis acidity and electrophilicity of the boron centre, shifting the hybridisation from sp2 to tetrahedral sp3 at neutral pH (Fig. 3b). DFT computations support this hypothesis by showing that this molecule has indeed a strong binding affinity for the OH− group, forming a closed [Eu·oBOH3] species with OH− bridging between the Eu(III) and boron atoms (Fig. 3e), occupying a similar coordination site to an isolated water molecule (Fig. 7a). A free energy change of −17.3 kJ mol−1 was obtained for deprotonation of the water bound complex at pH 7, corresponding to a pKa value of 3.96, in good agreement with experimental results.
Analysis of the emission lifetime of [Eu·oBOH2]+ in H2O and D2O at pH 4 corresponded to one water molecule bound to the Eu(III) centre (Table 1), suggesting that coordination of the boronate ester is reversible with variation of pH. The pH response of the [Eu·oBOH2]+ in the range 4.0–5.8 could potentially provide a means to probe the pH of acidic cellular organelles. Previous work in this area includes probes that internalize in the low pH environment of lysosomes in viable cancer cells. These probes contain a BODIPY moiety, with pKa values ranging from 4.4–5.8.52
Complex | Anion | τ H2O /ms | τ D2O /ms | q |
---|---|---|---|---|
a Mean ± standard deviation for two independent measurements. b Values of hydration state q (calculation error is ±20%) were derived using literature methods.47 | ||||
None | 0.192 ± 0.037 | 0.260 ± 0.069 | 1.2 | |
HPO42− | 0.452 ± 0.002 | 0.524 ± 0.016 | 0.1 | |
[Eu·pBOH2]+ | AMP | 0.494 ± 0.003 | 0.581 ± 0.004 | 0.1 |
ADP | 0.489 ± 0.004 | 0.621 ± 0.004 | 0.2 | |
ATP | 0.341 ± 0.004 | 0.508 ± 0.011 | 0.9 | |
Citrate | 0.315 ± 0.044 | 0.412 ± 0.064 | 0.6 |
Next, we assessed the ability of [Eu·pBOH2]+ to discriminate between the nucleotide phosphate anions, AMP, ADP, ATP between the nucleotide phosphate anions, AMP, ADP, ATP and cAMP. Addition of 1 mM AMP to [Eu·pBOH2]+ induced a 4-fold enhancement in overall Eu(III) emission intensity and a striking 6-fold increase in the ΔJ = 2 band (Fig. 4a). Adding 1 mM ADP caused a comparatively smaller (approximately 2-fold increase) in total emission, whereas ATP induced only a minor change in emission and cAMP and adenosine (Fig. S10†) induced no change in emission intensity or shape at all. Similar to [Eu·mBOH2]+, the para-substituted complex exhibits the selectivity profile AMP > ADP > ATP, representing a reversal of the selectivity order observed for most reported nucleoside phosphate receptors, where the binding is governed by electrostatic interactions and thus follows the trend ATP > ADP > AMP.21–24 The remarkable sensitivity of [Eu·pBOH2]+ for AMP was demonstrated by the linear increase in ΔJ = 2/ΔJ = 1 emission intensity ratio over the 0–100 μM range (Fig. S11†).
Based on the close similarities in the emission spectral shape of [Eu·pBOH2]+ in the presence of AMP and hydrogen phosphate (Fig. 4 and S7†), a similar monodentate binding mode to the Eu(III) centre can be assumed. Luminescence lifetimes revealed the displacement of water by AMP and ADP (Table 2), however in the presence of ATP a q value of 1.0 was found, consistent with lack of ATP binding. Our previous NMR spectroscopic studies of a family of conformationally more flexible Eu(III) complexes revealed the preference of ATP to chelate the europium(III) centre via the α and γ-phosphate groups, whereas ADP adopts more than one binding conformation (e.g. both bidentate and monodentate binding).33 In our current systems, we hypothesise that the increased rigidity and steric demand imposed by the 8-(benzyloxy)quinoline pendant arm prevents ATP or ADP from binding in a bidentate manner. Since ADP induces a similar change in emission spectral shape in [Eu·pBOH2]+ to that of AMP (Fig. 4a), it is assumed that ADP has a higher tendency for monodentate binding to the Eu(III) complex via its terminal phosphate group compared with ATP.
In contrast, complex [Eu·oBOH2]+ showed no change in emission intensity or spectral form in the presence of phosphate, bicarbonate, lactate, sulfate, acetate, or the nucleotide phosphate anions, AMP, ADP, ATP and cAMP. This is consistent with the lack of an available binding site due to the coordination of the ortho-substituted boronate ester (Fig. 4b and S8†). The only anion tested that induced a change in emission for both Eu(III) complexes was citrate, which caused a reduction in the overall emission intensity, most notably within the ΔJ = 1 and ΔJ = 4 bands (Fig. S7 and S8†), indicating a common binding mode for this anion. The binding of citrate may cause de-coordination of the quinoline antenna leading to less efficient sensitization and partial hydration of the Eu(III) ion, consistent with q = 0.6, in the presence of citrate (Table 2).20
Apparent binding constants were determined for [Eu·pBOH2]+ with phosphate, AMP and ADP in aqueous buffer at pH 7.0, by following the change in the intensity ratio of the ΔJ = 2/ΔJ = 1 bands as a function of anion concentration, followed by fitting of the data to a 1:1 binding model (Table 3, Fig. S7 and S12†). Stronger binding was observed for AMP (logKa = 4.22) compared with phosphate (logKa = 3.29), whereas a similar binding constant was obtained for ADP (logKa = 4.35). It was not possible to estimate the binding constant for ATP, due to the minor emission changes that took place. The difference in binding affinities between AMP and ADP for [Eu·pBOH2]+ was less favourable compared with the previously reported complex [Eu·mBOH2]+. The presence of the para-substituted boronic acid in [Eu·pBOH2]+ appears to increase affinity for ADP, indicating the meta-substituted boronic acid in [Eu·mBOH2]+ provides a slightly better geometry for AMP binding over ADP. As such [Eu·mBOH2]+ exhibits a better level of discrimination for AMP over ADP (Fig. S9†); a 2.3-fold difference in emission is seen at 615 nm for [Eu·mBOH2]+ compared with a 1.9-fold difference observed for [Eu·pBOH2]+. However, the most striking feature of both these complexes is their large emission enhancement for AMP with almost no interference from ATP.
Anion | [Eu·pBOH2]+a | [Eu·mBOH2]+a,b |
---|---|---|
a Mean ± standard deviation for two independent measurements. b Values previously reported in ref. 20. c Value not determined due to insufficient increases in emission intensity preventing accurate data fitting. | ||
HPO42− | 3.29 ± 0.04 | 3.31 ± 0.04 |
AMP | 4.22 ± 0.01 | 4.07 ± 0.01 |
ADP | 4.35 ± 0.01 | 3.84 ± 0.01 |
ATP | n.d.c | n.d.c |
High resolution mass spectrometric data supported 1:1 binding of AMP, ADP and phosphate to [Eu·pBOH2]+. We observed two overlapping signals at m/z 850.1501 and 854.1217 for the Eu(III) complex with phosphate (Fig. S14a†). The former signal corresponds to the host–guest complex [Eu·pBOH2 + H2PO4 + Na]+, and the latter can be tentatively assigned to loss of one of the boronic acid OH groups [Eu·pBOH2–OH + HPO4 + 2Na]+. In the presence of AMP, the [Eu·pBOH2]+ complex shows three signals (Fig. S15†). The major signal at m/z 1103.2074 can be attributed to the host–guest complex involving a boronate ester-1,2-diol interaction with the ribose sugar of AMP, as depicted in Fig. 5 and Fig. S15.† The smaller signals at m/z 1099.2368 and 1121.2178 are in agreement with the singly charged species [Eu·pBOH2 + AMP + H + Na]+ and [Eu·pBOH2 + AMP + 2Na]+ involving no 1,2-diol interaction. There was only one host–guest species observed for the Eu(III) complex with ADP, [Eu·pBOH2 + ADP + 3Na]+; this signal was very weak compared with that of the free complex at m/z 730.1907 (Fig. S14b†).
The 31P NMR spectrum of [Eu·pBOH2]+ in the presence of 1 equivalent of AMP revealed a signal at 3.4 ppm corresponding to unbound AMP and two broads signals at −207 and −222 ppm, indicating the strong likelihood of exchange between more than one binding mode for AMP (Fig. 6b). The existence of two AMP-receptor complexes may tentatively be explained by the reversible interaction between the phenylboronic acid moiety and the ribose sugar of AMP, consistent with the mass spectral data (Fig. 5 and S15†). In comparison, the addition of 10 equivalents of sodium hydrogen phosphate to [Eu·pBOH2]+ produced 2 resonances in the 31P NMR spectrum (Fig. 6a), corresponding to unbound phosphate (3.5 ppm) and phosphate bound in a single host–guest complex (−187 ppm).
The possibility of the multisite recognition of AMP by complex [Eu·pBOH2]+ was explored, involving Eu(III)–phosphate coordination combined with a boronate ester-1,2-diol interaction with the ribose sugar of AMP. Formation of the cyclic boronic ester was possible (Fig. 7d) and this multisite binding leaves the geometry around europium largely intact with a similar geometry to the monodentate binding mode (Fig. 7c). However, the bond lengths of Eu(III) to the quinoline nitrogen and the ether oxygen are elongated indicating some ring strain involved.
The computational effort required to obtain binding energies with accuracy high enough to reproduce the subtle trends seen in Table 3 is well beyond the scope of this work.37 However, we shall continue with a qualitative discussion of binding energies (Table S1†). First, the computations highlight the strong affinity of these complexes for binding hydrogen phosphate, and in line with experiment, binding to AMP is enhanced. Furthermore, computations indicate that the cyclic structure of AMP with [Eu·pBOH2]+ is bound but that its energy is disfavoured by about 30 kJ mol−1 compared to the open analogue. However, the existence of such a complex with a reasonable geometry and energy makes it plausible that both proposed AMP-receptor complexes exist in aqueous solution, consistent with the 31P NMR spectroscopic data which shows two distinct resonances for the AMP-bound species.
A comparison of experimental and computational emission spectra is shown in Fig. 8. Here, we consider the water bound complex [Eu·pBOH2]+ (as shown in Fig. 7a) and the phosphate bound complex (Fig. 7b) as a model for AMP, considering their similar binding geometries and assuming the electrostatic effect of the phosphate group has the strongest influence on the emission spectrum.
Pleasingly, our computational treatment naturally produces the different emission bands (ΔJ = 0, 1, 2, 3 and 4) along with their relative intensities. Remarkably, it does so without any adjustable parameters aside from three trivial scaling factors (discussed in the ESI†).
Starting with the water-bound complex, the simulations do not reproduce the formally forbidden ΔJ = 0 band. Next, the ΔJ = 1 band is well-reproduced including its crystal field splitting: just as in seen by experiment, two states are close together at lower wavelength with one isolated state at higher wavelength. The presented simulations underestimate the intensity of the hypersensitive ΔJ = 2 band, possibly due to missing vibronic effects. Finally, the low intensity of the ΔJ = 3 band and the high intensity along with its appropriate shape of the ΔJ = 4 band are well reproduced.
Moving from the water-bound complex to the phosphate-bound complex the experimental emission spectrum shows three notable changes: (i) an overall increase in intensity, (ii) an increase in the ΔJ = 2/ΔJ = 1 emission intensity ratio, and (iii) differences in the crystal field splitting. Effect (i) derives from quenching by nearby water molecules and is not covered by the present treatment. However, the other two are well reproduced. The increase in the ΔJ = 2 (605–630 nm) emission intensity is very pronounced showing a 16.1-fold increase when binding phosphate. This rise is more pronounced compared with the experimental value of a 4.2-fold increase but certainly highlights that the correct physics is captured and indicates that the approach can be used to screen for complexes with large ΔJ = 2/ΔJ = 1 changes upon guest binding. Finally, the simulations present appropriate crystal fielding splitting in the phosphate bound complex, specifically the structure of the ΔJ = 1 band changes with respect to the water-bound complex with one state at lower wavelength and two states at higher wavelength (which merge to form one apparent peak in the experimental spectrum, Fig. 8a). The magnitude of the splitting within the ΔJ = 1 band is also reduced. Taken together, the experimental and simulated spectra indicate a change in magnitude and sign of the crystal field coefficient.55 The shape of the ΔJ = 2 band is also well-reproduced with a large peak around 611 nm (deriving from four of the 7F2 states) and a smaller peak at 620 nm (deriving from the fifth state).
These data suggest that first principles computation offers a promising, albeit so far not exploited, approach to access the emission spectra of lanthanide complexes. A qualitatively correct description of the major features of the spectra along with changes upon coordination was readily achieved with the presented approach. We believe that this approach can be utilized as a powerful tool in the rational design of responsive Ln(III) complexes, including host–guest systems with predictable changes in local crystal field and thus specific spectral properties, while also providing a useful way to crosscheck computed structures with experiment.
NMR and mass spectral data supported by DFT calculations indicated the existence of two AMP-receptor complexes in aqueous solution, each stabilised by a phosphate–Eu(III) interaction and with one species involving interaction between the boronate ester and the ribose sugar of AMP. Complex [Eu·mBOH2]+ showed the highest level of discrimination between AMP and ADP, indicating that the meta-substituted boronic acid provides a slightly better geometry for AMP binding over ADP. Both complexes [Eu·pBOH2]+ and [Eu·mBOH2]+ exhibit the sensing selectivity order AMP > ADP > ATP, representing a reversal of that observed for the majority of receptors for nucleoside phosphates. This is attributed primarily to the incorporation of a sterically demanding ligand that caps the ‘axial’ coordination site of the central Eu(III) inhibiting access of the larger chelating polyphosphate anions.
Finally, first principles simulations of lanthanide emission spectra, both in the absence and presence of anionic guests, reproduced all qualitative features well highlighting the power of such an approach in the design of new responsive lanthanide host–guest complexes with specific photophysical properties.
We gratefully acknowledge financial support by the EPSRC [EP/S032339/1].
Footnote |
† Electronic supplementary information (ESI) available: Synthesis, characterisation and photophysical analysis of Eu(III) complexes, spectroscopic and computational studies of host–guest complexes. See DOI: https://doi.org/10.1039/d2qo01067d |
This journal is © the Partner Organisations 2022 |