Open Access Article
David
Bauer
ab,
Markus
Blumberg
ab,
Martin
Köckerling
c and
Constantin
Mamat
*ab
aInstitut für Radiopharmazeutische Krebsforschung, Helmholtz-Zentrum Dresden-Rossendorf, Bautzner Landstraße 400, D-013228 Dresden, Germany. E-mail: c.mamat@hzdr.de
bFakultät Chemie und Lebensmittelchemie, TU Dresden, D-01062 Dresden, Germany
cInstitut für Chemie-Anorganische Festkörperchemie, Universität Rostock, Albert-Einstein-Straße 3a, D-18059 Rostock, Germany
First published on 10th October 2019
Metals, which form divalent cations, including the alkaline earth metals offer radionuclides like lead-203, lead-212, barium-131, and strontium-89, which are promising candidates for radiopharmaceutical applications. Besides, the heavy homologous nuclides radium-223 and radium-224 – with similar properties to barium – are suitable alpha-emitters for targeted alpha-particle therapy. However, there is a lack of suitable chelation agents, especially for heavy group 2 metals. The macrocycle calix[4]arene-1,3-crown-6 seems to interact with these metals strongly. Therefore, this ligand and its coordination to the divalent cations of barium, strontium, and lead have been investigated. The complex formation was analyzed by NMR and UV/Vis titration experiments in acetonitrile, and stability constants were determined to be >4 with both methods. It was found that the stability of these complexes increase in the order strontium, barium, and lead. Additional to these investigations, X-ray crystallography, solvent-dependent 1H NMR, and 207Pb NMR measurements were performed to deliver deeper insight into the coordination chemistry of this ligand.
For this work, the divalent cations Sr2+, Ba2+, and Pb2+ have been selected, since they provide radionuclides, which possess useful properties for application in radiopharmacy as well as nuclear medicine and are therefore of particular interest.
In the field of targeted and pre-targeted radioimmunotherapy, 212Pb has demonstrated significant utility in both in vitro and in vivo systems.4 The radionuclide 212Pb is not only a promising β−-emitter but serves as an in vivo generator for the alpha-emitter 212Bi.5 Alpha-emitters combine a short range with a high linear energy transfer, which results in the relatively high biological effect and cytotoxicity.6 Furthermore, the combination of 212Pb and 203Pb is a matched radionuclide pair for image-guided radionuclide therapy. 203Pb decays by electron capture with a γ-emission of 279 keV (80%) which makes it a promising SPECT (single photon emission computer tomography) imaging agent.7
Alkaline earth metals belong to the most prominent group of divalent ions, which offer a variety of interesting radionuclides. The utility of 89Sr in the treatment of advanced metastatic prostate cancer has been examined in numerous clinical trials.8–1189Sr in the form of its dichloride (Metastron) is FDA-approved for radionuclide therapy since 1993 and is used for the palliation of painful osseous metastases.12 The element barium offers the radionuclide 131Ba, which decays by electron capture while emitting suitable γ-rays for diagnostic use. The radionuclide 131Ba could be a promising bone-scanning agent in scintigraphy.13,14 The group of heavy alkaline earth metals can also be expanded to its heaviest homolog. The element radium offers the α-emitters 223Ra and 224Ra. Due to their suitable half-lives of 11.4 d, respectively 3.6 d, and their excellent decay properties (a cascade of 4 × α and 2 × β decays), both are ideal nuclides for the targeted alpha-particle therapy.15 Radium, in the form of [223Ra]RaCl2 (Xofigo), is indicated for the treatment of metastatic castration-resistant prostate cancer (mCRPC) and received marketing approval by the EMA and FDA in 2013.16 Since there is no stable radium isotope, it is not a trivial challenge to investigate its chemistry. However, due to its related chemical behavior, barium can be seen as its non-radioactive surrogate and findings derived from the barium chemistry will most likely apply for radium as well.17–19 All mentioned radionuclides and their properties are listed in Table 1.
| Nuclide | Radii (Å) | t 1/2 | Decay mode | E (keV) | Production method | Comments | Application |
|---|---|---|---|---|---|---|---|
| 203Pb | 1.29 | 51.9 h | ε (100%) | 279 (γ) | 203Tl(p,n)203Pb | SPECT | |
| 212Pb | 10.6 h | β− (100%) | 570 | 224Ra/212Pb generator | Daughter is 212Bi (α) | β−/α therapy | |
| 89Sr | 1.26 | 50.6 d | β− (100%) | 1.5 × 103 | 88Sr(n,γ)89Sr | β− therapy | |
| 131Ba | 1.42 | 11.5 d | ε (100%) | 124, 216 (γ) | 131Cs(p,n)131Ba | SPECT | |
| 223Ra | 1.48 | 11.4 d | α (100%) | 6.0 × 103 | 227Ac/223Ra generator | Decay chain with 4 × α and 2 × β | α therapy |
| 224Ra | 3.6 d | α (100%) | 5.8 × 103 | 228Th/224Ra generator |
Although these nuclides possess suitable properties, the absence of stable bifunctional chelating agents for these radiometals hampers their use for targeted radiotherapy and radioimaging.1 The essential requirement for a good chelator is the thermodynamic and kinetic stability of its metal–chelate complex in vivo; this necessitates the design of ligands appropriate for the particular radionuclide.2,20 To compare and evaluate chelators with high complexation ability, the knowledge of the stability constant (log
K) is mandatory.
There is a constant search for ligands perfectly adapted to the particular application. Noteworthy, the recently developed DOTA-derivative TCMC (1,4,7,10-tetrakis(carbamoylmethyl)-1,4,7,10-tetraazacyclododecane) seems to be promising for Pb.28,29 Though, there are no remotely suitable chelating agents known for heavy alkaline earth metals.
Calix-crown ethers are a widely investigated class of ligands for cations based on the calixarene backbone. They provide a suitable coordination chemistry towards alkali metal and alkaline earth metal cations.30–34 Depending on the metal, the coordination appears to involve not just the crown ether-oxygen donors, but also the aromatic phenol-units of the calixarene.351H NMR studies can provide further experimental evidence for the influence of cation–π-interactions and are a key part of this work.
The calix[4]arene-1,3-crown-6 bears two hydroxyl groups located on opposite sides. These groups act as ionizable units by the loss of their protons under basic conditions and thus to form a neutral complex witch a divalent cation. Additionally, the hydroxyl functions can straightforwardly be modified by attaching other ionizable units.
The objective of this research was to evaluate calix[4]arene-1,3-crown-6 (1) as a lead compound that could, upon further modifications, yield a viable ligand for the above mentioned divalent cations. Literature about alkaline earth metal ligands, especially about radium, are mostly focused on extraction studies without information about comparable stability constants.32,36,37 Therefore, UV/Vis and NMR titration experiments were established as reliable and constant methods for the determination of the stability constants of calix 1 with Sr2+, Ba2+, and Pb2+. Due to an activity-related concentration limit, Ra2+ could not be investigated.
The second route is a one-pot reaction starting directly with the calix[4]arene and pentaethylene glycol ditosylate. Potassium tert-butoxide is used as a base to deprotonate the hydroxyl groups; additionally the potassium cations serve the purpose of a template effect and arrange the formation of the cone conformation. However, this reaction takes 3 days and obtains a moderate yield.38 For this work, ligand 1 was obtained with a higher yield of 50% using a modified reaction according to the second route.39
When slowly evaporating the solvent of a saturated solution of 1 in dichloromethane, the formation of single crystals was observed. These crystals were analyzed by single-crystal X-ray diffraction, which unfolded an inclusion complex with dichloromethane located in the upper rim of the aromatic calix[4]arene part (Fig. 1). Compound 1 crystallizes in the monoclinic crystal system with the centrosymmetric space group P2/n. The n glide plane of the space group cuts through the calixarene molecule as well as through the carbon atom of the dichloromethane molecule. Therefore, the asymmetric unit contains only half of both molecules. The part of the crown ether ring starting from O3 up to O3′ (
− x, y, ½ − z) is disordered with two different orientations as well as the chlorine atom of the dichloromethane molecule. This disorder is refined using split positions (A and B) for each affected atom (group). The different orientations of the ether part and the dichloromethane molecule are shown in Fig. 1B.
The four aromatic rings of the calixarene form a cavity, which is expanded at the upper rim. The opposite rings are inclined to each other with an angle of 34.5° of the two rings, which are attached to the crown ether moiety, compared to the other two unsubstituted rings, where the angle between the two main planes through the ring carbon atoms was determined with 84.2°. A little above the upper rim, the (disordered) dichloromethane molecule is located. The two hydrogen atoms, whose positions were initially located from the difference Fourier map and later refined using riding models, are oriented such that they have relatively close distances to the centers of the phenyl groups, which are attached to the crown ether residue. The distance of these protons to the center of the rings is 2.485 Å. This is close enough to be discussed as weak C–H–π interaction.40 The oxygen atoms of the crown ether moiety (O2, O3, O4(A, B) and symmetry equivalents) are not primarily directed toward the inner part of the cavity, as they would do whenever a metal cation is coordinated. The diameter of this cavity, as defined by the distances of opposite oxygen atoms, is 5.887 Å (O2–O4A′ (−½ + x, 1 − y, ½ + z) 2 times) and 6.849 Å (O3–O3′ (−½ + x, 1 − y, ½ + z)). The shortest hydrogen bonds are found intermolecularly as O1–H1A⋯O2 contact at a D⋯A distance of 2.767(1) Å and a D–H⋯A angle of 167.7°. All other possible proton contacts have D⋯A distance longer than 3.0 Å.
| 1H signala | Multiplicity [Hz] | Integral | Shift [ppm] | Assignment | ||||
|---|---|---|---|---|---|---|---|---|
| CDCl3 | CD3CN | C6D6 | (CD3)2SO | (CD3)2CO | ||||
| a Is not representing the IUPAC assignment. Compare the assignment with Fig. 2. | ||||||||
| H1 | s | 2 | 7.50 | 7.68 | 7.83 | 8.02 | 8.00 | OH |
| H2 | d# | 4 | 6.84 | 6.94 | 6.69 | 6.99 | 6.67 | Ar–H |
| H3 | t# | 2 | 6.73–6.67 | 6.77 | 6.26 | 6.75 | 6.95 | Ar–H |
| H4 | d§ | 4 | 7.07 | 7.13 | 7.09 | 7.13 | 7.12 | Ar–H |
| H5 | t§ | 2 | 6.67–6.63 | 6.67 | 6.85 | 6.58 | 6.60 | Ar–H |
| H6 | d* | 4 | 3.37 | 3.42 | 3.27 | 3.41 | 3.42 | Hendo |
| H7 | d* | 4 | 4.42 | 4.40 | 4.58 | 4.30 | 4.47 | Hexo |
| H8 | m | 4 | 4.18–4.12 | 4.14–4.11 | 3.82–3.78 | 4.10–4.06 | 4.20–4.15 | CH2O |
| H9 | m | 4 | 4.04–4.00 | 3.93–3.89 | 3.70–3.63 | 3.98–3.94 | 4.09–4.05 | CH2O |
| H10 | m | 4 | 3.96–3.91 | 3.82–3.78 | 3.70–3.63 | 3.81–3.77 | 3.92–3.87 | CH2O |
| H11 | m | 4 | 3.87–3.83 | 3.73–3.68 | 3.70–3.63 | 3.74–3.70 | 3.80–3.76 | CH2O |
| H12 | s | 4 | 3.71 | 3.57 | 3.54 | 3.58 | 3.64 | CH2O |
| #2 J [Hz] = | 7.6 | 7.6 | 7.6 | 7.6 | 7.5 | |||
| §2 J [Hz] = | 7.5 | 7.5 | 7.5 | 7.5 | 7.5 | |||
| *2J [Hz] = | 13.1 | 13.1 | 13.1 | 12.8 | 13.0 | |||
No significant influence of the solvent was observed when comparing the coupling constants from Table 2 (2J), and except of benzene-d6, the signals in the spectra are similar arranged for the different solvents. As expected, there is an influence of the solvent on the hydroxyls (H1), depending on its polarity and amount of acceptor sites. The stronger the hydroxyl oxygen can coordinate to the solvent molecules, the less electron density remains at the H1 proton (downfield shift). However, this does not apply for the solvent benzene, which has the lowest polarity in comparison, but a significant H1-downfield shift. A likely explanation for the solvent effect of benzene on the calix 1 is illustrated in Fig. 3. The solvent interacts strongly with the phenolic rings (π–π-interaction), which will cause a flattening of the cone (Fig. 3-1b) and a contraction of the upper-rim of the calix 1. Hence, the intramolecular hydrogen bonds break and the H1 signal shifts downfield. The protons H2 and H3 are facing into the cavity. Due to the aromatic ring current effect induced by the magnetic field of the NMR spectrometer, these protons perceive a lower magnetic field, and this will be noticed as a highfield shift in the spectrum. Since the aromatic rings are facing the crown ether, the closest protons (H8–H10) will also perceive a lower magnetic field (highfield shift of their signals).
The other smaller solvents, as represented in the crystal structure (Fig. 1), instead form an inclusion complex with calix 1, which leads to an expansion of the upper rim (Fig. 3-1a).
Comparing the different spectra, the signals received from the acetonitrile-measurement, are sharp and well separated, and there does not seem to be significant interaction with the hydroxyl functions. Additionally, we use metal perchlorate salts to study the ligand–metal-interaction; those are highly soluble in acetonitrile. The use of acetonitrile-d3 is superior over DMSO-d6, because DMSO is known to form complexes with divalent cations.39,46 For these reasons and since the compound is not soluble in water, the stability constants are determined in acetonitrile. Two different methods have been used to determine the stability constants: NMR and UV/Vis titration with the respective metal perchlorate salts.
![]() | ||
| Fig. 4 1H NMR spectra of the 1:1 complex of ligand 1 with Pb2+, Sr2+, and Ba2+ in acetonitrile-d3. The signal at 7.58 represent residual CHCl3 and do not affect the method. | ||
Stability constants of these complexes were determined by the relationship between the shifts of a selected proton signal and the equivalents of the cation. A reliable method to determine the stability constants via1H NMR spectroscopy was developed in the past,39 using Pb2+, Sr2+, and Ba2+ as their perchlorate salts and acetonitrile-d3 as the solvent. For the calculation of the stability constants, the chemical shift has to be precisely determined. To ensure this, the changes in the chemical shifts of the 1H signals between ligand 1 and the cation-containing complexes Pb-1, Sr-1, and Ba-1 must be sufficiently pronounced. Considering all 1H NMR signals and their shifts, especially H4 is a suitable signal for these calculations, since it remains sharp and does not overlap with other signals upon complexation. This signal is shifted downfield, as the H4 proton moves out of the aromatic ring current zone caused by the contraction of the lower rim. The dependence of chemical shifts on the metal concentration is illustrated in Fig. 5 by the example of Ba2+ (see ESI† for Sr2+and Pb2+).
![]() | ||
| Fig. 5 1H NMR spectra of ligand 1 at different Ba(ClO4)2 concentrations measured in acetonitrile-d3. Signal H4 is sufficiently pronounced to be used for the stability constant calculations. | ||
When plotting the 1H NMR chemical shifts of the signal H4 against the equivalents of Ba2+ (Fig. 6), a change of the slope is obtained at the ratio of 1:1, indicating the formation of a 1:1-complex. The obtained data for the shift of the signal H4 was evaluated by using the WinEQNMR2 software,47 calculating a stability constant of 4.6 ± 0.4.
![]() | ||
| Fig. 6 Shifts for the 1H NMR H4-signal of compound 1 at different Ba(ClO4)2 concentrations measured in acetonitrile-d3. | ||
Comparable plots and calculations were obtained for titrations with Sr2+ and Pb2+ (see ESI† and Table 3). No changes in the chemical shifts were observed when using the perchlorates of Na+ and Bu4N+, indicating that ligand 1 shows a size selectivity for the cations of interest.
| Method | log K for the selected cations |
||||
|---|---|---|---|---|---|
| Pb2+ | Sr2+ | Ba2+ | Na+ | Bu4N+ | |
| 1H NMR | 5.5 ± 0.2 | 4.3 ± 0.2 | 4.6 ± 0.4 | — | — |
| UV/Vis | 6.36 ± 0.05 | 4.07 ± 0.08 | 4.66 ± 0.02 | — | — |
Additionally, the formation of the 1:1-Pb2+-calix-complex Pb-1 was confirmed by 207Pb NMR. The isotope lead-207 (natural abundance: 22.6%) has a medium sensitivity NMR spin-½ nucleus.48,49 For this experiment, lead(II) perchlorate trihydrate was dissolved in acetonitrile-d3 and titrated with ligand 1. The difference between the 207Pb signal of Pb(ClO4)2 and Pb-1 resulted in a chemical shift of Δδ = 296 ppm (Fig. 7). Since the trihydrate salt was used, it has to be assumed that upon complexation, the hydration shell of Pb2+ was released to interact with the calix 1. This exchange resulted in a downfield shifted signal.
![]() | ||
| Fig. 8 (A) UV/Vis spectra of compound 1 at different Ba(ClO4)2 concentrations measured in acetonitrile. (B) Absorption at the wavelength 279 nm of compound 1 at different Ba(ClO4)2 concentrations. | ||
K values for all divalent metals of >4 (Table 3). It was found that the stability of these complexes increase in the order of strontium, barium, and lead. No complexation was observed for Na+ and Bu4N+. Thus, high selectivity of the heavy alkaline earth metal ions and of Pb2+ was found over the control cations Na+ and Bu4N+. The log
K values of both methods are comparable under the given conditions. However, both of these methods were performed close to their limits of validation,51 therefore, values >5 may show the right tendency but are not of high accuracy. The limits of these methods are related to the log
K value itself and to the concentration of the host and guest molecule that is demanded to adequately detect 1H NMR signals or the UV/Vis absorption, respectively.52,53
The tendency that Ba2+ is showing a higher log
K than Sr2+, the cation with the smaller radius, suggest that radium might provide even higher stability. Though, alternative methods have to be developed to determine higher log
K values and to test for the radium stability as well.
:
metal ratio of 2
:
3, 30 μL portions of a 1.0 M perchlorate salt solution were used, and stepwise additions were continued until a ligand
:
metal ratio of 1
:
6 was reached to exclude the formation of a complex with another stoichiometry. The displacements of selected 1H NMR signals of ligand 1 upon addition of the perchlorate salt were used to calculate the complex stability constants. All calculations were performed using the WinEQNMR2 software.47 The advised range for the data input covers the addition of metal to ligand from 0.1 to 0.9 equivalents. This instruction was followed and 9 points in this range were measured (steps of 0.1 equiv.) and used for the calculation. The formation of a 1:1 complex was proven by plotting the changes of selected signals against the cation to chelate ratio, observing the change of the slope at a ligand
:
metal ratio of 1
:
1.
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 1950592. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c9ra07293d |
| This journal is © The Royal Society of Chemistry 2019 |