Open Access Article
Thomas R.
Berki
ab,
Jonathan
Martinelli
c,
Lorenzo
Tei
c,
Helen
Willcock
*b and
Stephen J.
Butler
*a
aDepartment of Chemistry, Loughborough University, Leicestershire LE11 3TU, UK. E-mail: s.j.butler@lboro.ac.uk
bDepartment of Materials, Loughborough University, Leicestershire LE11 3TU, UK. E-mail: h.willcock2@lboro.ac.uk
cDepartment of Science and Technological Innovation, Università del Piemonte Orientale, I15121, Alessandria, Italy
First published on 1st February 2021
A new synthetic strategy for the preparation of macromolecular MRI contrast agents (CAs) is reported. Four gadolinium(III) complexes bearing either one or two polymerizable methacrylamide groups were synthesized, serving as monomers or crosslinkers for the preparation of water-soluble, polymeric CAs using Reversible Addition–Fragmentation Chain Transfer (RAFT) polymerization. Using this approach, macromolecular CAs were synthesized with different architectures, including linear, hyperbranched polymers and gels. The relaxivities of the polymeric CAs were determined by NMR relaxometry, revealing an up to 5-fold increase in relaxivity (60 MHz, 310 K) for the linear polymers compared with the clinically used CA, Gd-DOTA. Moreover, hyperbranched polymers obtained from Gd(III) crosslinkers, displayed even higher relaxivities up to 22.8 mM−1 s−1, approximately 8 times higher than that of Gd-DOTA (60 MHz, 310 K). A detailed NMRD study revealed that the enhanced relaxivities of the hyperbranched polymers were obtained by limiting the local motion of the crosslinked Gd(III) chelate. The versatility of RAFT polymerization of Gd(III) monomers and crosslinkers opens the doors to more advanced polymeric CAs capable of multimodal, bioresponsive or targeting properties.
An exciting prospect in MRI CA design is the development of macromolecular systems that possess significantly higher relaxivities. Numerous strategies have been pursued,5–8 including the conjugation of one or several Gd(III) complexes to polymers,9–12 dendrimers,13–18 micelles19–23 or nanoparticles,24–31 or via the non-covalent association with biomolecules (e.g. serum albumin protein)32,33 or nano-assembled capsules.34,35 Some macromolecular systems have been shown to possess significantly higher relaxivities compared with commercial CAs (up to 40 mM−1 s−1), attributed predominantly to the slower tumbling of the macromolecules, and the incorporation of several Gd(III) complexes within a single system.4,6
Despite these advances, the challenge still remains to develop high molecular weight CAs wherein the global motion of the macromolecule is effectively coupled to the motion of the paramagnetic centre.34,36 One way to approach this challenge involves positioning a single Gd(III) chelate at the barycentre of the macromolecule (e.g. dendritic systems).17,37 The rotational correlation time of these high molecular weight Gd(III) chelates is thus defined by the motion of the macromolecule. Polymeric CAs have been developed which incorporate multiple Ln(III) chelates, linked via a single flexible arm to the polymeric scaffold. However, the use of flexible linkers permits local motion of the paramagnetic centre, limiting the degree of coupling with the more slowly tumbling polymer.9–12,38
A further challenge involves the development of new methods to control the macromolecular size, structure and shape, as this could lead to well-defined second spheres of hydration, whilst allowing fast-inner sphere water exchange. The majority of polymeric CAs are prepared by classical conjugation of appropriate ligands to a polymer bearing reactive pendant groups (e.g. maleimide group, ester activated monomers). Typically, a protected ligand is covalently attached to a synthesized polymer and the resulting conjugated ligand is deprotected, followed by complexation with Gd(III).9,10,39 In related work, Sherry and co-workers have presented a strategy for polymeric PARACEST agents in which non-metal containing ligands are directly incorporated by the free radical polymerisation of acrylamide functionalised ligands, followed by complexation with Eu(III).38 However, this approach has limitations, because it is difficult to quantify both the extent of ligand conjugation and the degree of lanthanide(III) complexation. Furthermore, the removal of residual lanthanide(III) ions can be challenging, but is crucial if such macromolecules are to be considered for biomedical or clinical applications.
Herein we present a new synthetic approach to macromolecular MRI CAs, involving the single-step incorporation of kinetically stable, monomeric Gd(III) complexes within well-defined macromolecular CAs. We have synthesized four DOTA-like Gd(III) complexes (Gd·L1–4, Fig. 1), bearing one or two pendant methacrylamide arms, which serve as monomers and crosslinkers respectively, for the synthesis of linear and hyperbranched polymers using reversible addition fragmentation chain transfer (RAFT) polymerization. Gd·L1 contains a single methacrylamide arm, capable of forming linear polymers, whereas complexes Gd·L2–4 contain two methacrylamide arms, serving as crosslinkers to create hyperbranched polymers, in which the motion of the Gd(III) centre is effectively coupled to the slowly tumbling macromolecule. Complex Gd·L2 possesses two trans-related polymerizable arms, whereas for Gd·L3 and Gd·L4 the polymerizable arms are in cis-geometry, with Gd·L3 possessing shorter arms. The impact of these structural and geometric modifications on the relaxivity and tumbling motion of the resulting polymers was evaluated. Each complex possesses a DOTA-like core, to confer maximal thermodynamic and kinetic stability.3,40 Moreover, Gd·L1–4 possess a single coordinated water molecule (q = 1) and are negatively charged overall, thereby allowing relatively fast water exchange to overcome the limitations of some previously reported macromolecular CAs.7,8
For complexes Gd·L1–3, mono or bis-alkylation of the macrocyclic free amines with α-bromoester 1, prepared from diazotization and bromination of Cbz-protected L-ornithine (Scheme S2†), led to formation of the fully protected macrocyclic ligands (e.g.3 for Gd·L2, Scheme 1). The yields of the alkylations were significantly improved by the addition of potassium iodide to the reaction mixture (K2CO3/acetonitrile), allowing iodide/bromide exchange.44 It is possible that partial racemisation of the alkylating agent 1 or racemisation during the alkylation reaction occurred,45 leading to the formation of the protected ligands as a mixture of stereoisomers. The methacrylamide arms were introduced by Cbz deprotection of the ornithine sidechains to give the bis-amine (e.g. macrocycle 4), followed by coupling with N-hydroxysuccinimide methacrylate ester 5. Next, the tert-butyl esters were deprotected using trifluoroacetic acid, followed by the addition of GdCl3 in water at pH 7, to afford the water soluble Gd(III) complexes Gd·L1–3 after purification by preparative reverse-phase HPLC (Fig. S1–S3†). The Gd(III) complexes of a given isomer of ligand L1–3 will have further elements of chirality arising from the sign and torsion angles of the cyclen NCCN chelate rings, and the NCCO chelates defining the helicity of the pendant arms.3,36 As such, the Gd(III) complexes will exist as a mixture of stereoisomers in solution, which may interconvert by either cyclen ring inversion or arm rotation. The separation of stereoisomers was not attempted in this work.
The synthesis of Gd·L4 involved initial bis-alkylation of cis-DO2A(OtBu)2 (8) with α-bromoester 7, prepared from L-glutamic acid (Scheme S2†) to give protected ligand 9. Again, it is possible that partial racemisation of α-bromoester 7 occurred, or racemisation during the alkylation reaction, resulting in a mixture of stereoisomers of protected ligand 9. Subsequent deprotection of the tert-butyl esters of 9 using trifluoroacetic acid, followed by the addition of slight excess of GdCl3 in water at pH 7, gave the precursor Gd(III) complex 10. Finally, the methacrylamide groups were introduced via coupling the terminal carboxylic acids of 10 to N-(3-aminopropyl)-methacrylamide, using the coupling reagent TNTU (2-(5-norborene-2,3-dicarboximido)-1,1,3,3-tetramethyl-uronium tetrafluoroborate), to give the water soluble complex Gd·L4 after purification by reverse-phase HPLC. Analysis of the purified complexes Gd·L1–4 by analytical reverse-phase HPLC revealed a single peak in each case, and high-resolution mass spectral data confirmed formation of the desired complexes (Fig. 2 and S1–S4†). A major signal corresponding to the negatively charged molecular ion, [M]−, was observed in each case, and the isotopic distribution was in excellent agreement with the theoretical data.
| Entry | N NAM0,Chain | N Gd·L1 0,Chain | NAMc conv., % | Expected Mnd, g mol−1 | M n SEC, g mol−1 | Ð | [Gd]g, μg L−1 | N GdChain |
|---|---|---|---|---|---|---|---|---|
a RAFT polymerization in the presence of Gd-DOTA.
b Theoretical number of NAM and Gd·L1 units per chain depending of the initial polymerization reaction composition.
c Conversion determined by 1H NMR spectroscopy.
d Expected Mn = ([NAM]0MNAMConvNAM + [Gd·L1]0MGd·L1ConvGd·L1)/[CTA]0 + MCTA, with ConvGd·L1 = 1.
e Obtained by SEC analysis (CH3Cl/triethylamine 98 : 2 v/v, RID detectors).
f Obtained by SEC analysis (H2O/MeOH 80/20 v/v with 0.1 M NaNO3, RID detector).
g Gd(III) concentration determined by ICP-MS based on mass spectral signal of 157Gd isotope.
h
N
GdChain = number of Gd(III) ions per polymer chain, estimated from ICP-MS data (ESI, Section 2, Table S2).
i Not determined due to the high content of Gd·L1.
|
||||||||
| 1a | 100 | 0 + 4 Gd·DOTA | 97.6 | 14 181 |
13 900e |
1.08e | 0 | 0 |
| 2 | 99.1 | 0.9 | 98.2 | 15 632 |
11 800e |
1.18e | 7308 | 0.37 |
| 3 | 98.1 | 1.9 | 96.3 | 15 894 |
10 900e |
1.14e | 14 680 |
1.2 |
| 4 | 97.2 | 2.8 | 97.5 | 16 619 |
10 200e |
1.12e | 16 680 |
1.37 |
| 5 | 96.2 | 3.8 | 96.2 | 16 946 |
7,600f | 1.19f | 18 260 |
1.59 |
| 6 | 95.5 | 4.5 | 98.3 | 17 863 |
9,100f | 1.22f | 19 560 |
0.53 |
| 7 | 91.4 | 8.6 | 98.1 | 20 532 |
12.100f | 1.20f | 57 840 |
2.98 |
| 8 | 83.7 | 16.7 | N/Di | 25 335 |
14 200f |
1.25f | 62 680 |
8.02 |
Next, the synthesis of linear P(NAM-r-Gd·L1) copolymers was investigated using different molar proportions of Gd·L1, ranging from 1 to nearly 17 mol%. The polymerizations were conducted in a mixture of DMSO/water (80
:
20) at 80 °C (Scheme 2), and the polymers were purified by dialysis through a semi-permeable membrane against distilled water (15 MΩ cm−1). Successful synthesis of the target linear P(NAM-r-Gd·L1) copolymers was confirmed by SEC analysis (Fig. 3 and S7†): polymers with number average molecular weights (Mn) between 7600 and 14
200 g mol−1 were formed with low dispersity (Ð) values, ranging from 1.08 to 1.25 (Table 1, entries 2–8). Due to the charged nature of the Gd(III) complex and the difference between the PNAM and the standard used (PS), the polymers displayed lower molecular weights than the theoretical values. Standard 1H NMR analysis to confirm Mn for the polymers was not possible due to severe line broadening imposed by the paramagnetic Gd(III) ion; however, polymerization of NAM in the presence of free Gd-DOTA resulted in polymers with the expected Mn by 1H NMR spectroscopy.
![]() | ||
| Scheme 2 Synthesis of linear (DP 100 and m between 1 to 17) and crosslinked (DP 100 and m = 0.9) polymeric MRI CAs using NAM and Gd·L1–4. | ||
Analysis of the bulk magnetic susceptibility (BMS) shift of each sample and ICP-MS analysis of the copolymers after dialysis enabled an approximation of the number of Gd(III) ions per polymer chain (Tables 1, S2 and S3†). As expected, the number of Gd·L1 units incorporated into the polymer increased with increasing molar ratio of Gd·L1/NAM monomers used for the polymerization. In the highest case, eight Gd·L1 units were incorporated per polymer chain (Table 1, entry 8).
Having demonstrated that linear copolymers can be prepared in a controlled manner using the monomers Gd·L1 and NAM, we turned our attention to the crosslinker complexes Gd·L2–4. In order to prevent gelation from occurring and based on some of our previous work and that of Flynn et al.,50,53 the crosslinked polymers were obtained by fixing the concentration ratio of the crosslinker (Gd·L2–4) relative to the chain transfer agent (CTA), such that [Gd·L2–4]/[CTA] = 0.9. Initially, a range of water soluble hyperbranched polymers were synthesized by RAFT polymerization of NAM and the crosslinker Gd·L2 at different initial concentrations (1–2 M) in DMSO, using poly(ethylene glycol) methyl ether 2-(dodecylthiocarbono-thioylthio)-2-methyl-propionate as the CTA (Table 2, entries 1–4). For each reaction, SEC analysis of the molecular weight distribution indicated the formation of hyperbranched polymers, evident from the high molecular weight shoulder (Fig. 3, centre). The broad molecular weight distribution was reflected in the high dispersities, Ð, ranging between 4.2 and 19.9 (Table 2), consistent with crosslinking of the growing polymer chains during RAFT polymerization. Further evidence for crosslinking of Gd·L2 was given by analysis of the BMS shift and ICP-MS analysis (Tables 2 and S3†), which confirmed an increasing number of crosslinked polymer chains at higher initial monomer concentration. For example, when [NAM]0 was 1.00 M, the average number of crosslinker Gd·L2 units was estimated to be 7.0 per hyperbranched polymer chain, whereas this increased to approximately 34 units per chain when [NAM]0 was increased to 2.00 M. A higher initial monomer concentration led to higher molecular weight and dispersity, with values up to Mn = 744
000 g mol−1 and Ð = 19.9 when [NAM]0 = 2.00 M. However, polymerization conducted at higher initial monomer concentration, such that [Gd·L2]/[CTA] > 1, consistently led to the formation of gels (Table 2, entry 5). This is in accordance with the work of Perrier and co-workers, who showed that EGDMA/CTA ratios of over 1 led to gelation in a RAFT system.54
| Entry | Complex | Gd·L2–4 equiv. | [NAM]0 , mol L−1 | Expected Mna, g mol−1 | NAMb conv., % | M n SECc, g mol−1 | M w SECc, g mol−1 | Ð | [Gd]d, μg L−1 | N GdChain |
|---|---|---|---|---|---|---|---|---|---|---|
| a Expected Mn = (DP × n/100 × MNAM) + (DP × m/100 × MGd·LX) + MCTA, with ConvGd·L = 1. b Determined by 1H NMR spectroscopy. c Determined by SEC (DMF with 5 mM NH4BF4, RID/UV/LS detector). d Estimated number of Gd(III) ions per polymer chain, determined by ICP-MS (ESI, Table S3). | ||||||||||
| 1 | Gd·L2 | 0.9 | 1.00 | 15 817 |
99.1 | 33 200 |
138 600 |
4.2 | 7128 | 7.0 |
| 2 | Gd·L2 | 0.9 | 1.25 | 15 817 |
99.7 | 38 900 |
378 600 |
9.7 | 4920 | 21.9 |
| 3 | Gd·L2 | 0.9 | 1.50 | 15 817 |
98.8 | 30 800 |
592 700 |
19.2 | 6276 | 30.5 |
| 4 | Gd·L2 | 0.9 | 2.00 | 15 817 |
99.6 | 37 400 |
744 700 |
19.9 | 6008 | 33.6 |
| 5 | Gd·L2 | >1.0 | 2.00 | 15 884 |
>95 | Gel | Gel | Gel | N/D | N/D |
| 6 | Gd·L3 | 0.9 | 2.00 | 15 817 |
100 | 26 600 |
302 800 |
11.4 | N/D | N/D |
| 7 | Gd·L4 | 0.9 | 2.00 | 16 945 |
99.6 | 24 800 |
119 900 |
4.8 | N/D | N/D |
The optimal conditions found for the synthesis of the hyperbranched P(NAM-r-Gd·L2) polymers were [NAM]0 = 2.00 M, [Gd·L2]/[CTA] = 0.9 (entry 4). These parameters were applied to the synthesis of hyperbranched polymers using crosslinkers Gd·L3 and Gd·L4, bearing cis-related polymerizable arms (entries 6 and 7). Hyperbranched polymers P(NAM-r-Gd·L3) and P(NAM-r-Gd·L4) were successfully formed. Notably, they displayed lower molecular weights and dispersities compared with those obtained using Gd·L2 under similar conditions ([NAM]0 = 2.00 M, [Gd]/[CTA] = 0.9), suggesting that these polymerizations could be achieved at higher initial monomer concentration (>2.00 M).
The average hydrodynamic diameter Dh and dispersity of representative examples of linear and hyperbranched polymeric CAs were estimated by diffraction light scattering (DLS) measurements. From the number weighted particle size distribution, the linear polymer P(NAM-2%-Gd·L1) has a Dh of 4.3 ± 0.9 nm. Crosslinked systems with Gd·L3 or Gd·L4 display Dh of 13.2 ± 4.6 nm and 15.1 ± 4.3 nm, respectively (Fig. 4). This indicates that the polymers exist primarily as unimers in solution. The sizes of the hyperbranched polymers are comparable to previously reported hydrophilic and charged hyperbranched polymers (Dh between 10 to 20 nm for MW between 100 to 500 kDa).50 The Gd(III) monomers are charged and therefore hydrophilic, and as expected, do not direct the assembly of these polymers into higher order structures.
![]() | ||
| Fig. 4 Representative examples of DLS number weighted particle size distribution of linear and hyperbranched polymeric CAs obtained in deionized water at 25 °C after filtration with a syringe filter (220 nm cut-off). Associated correlograms and detailed DLS results are presented in the ESI (Fig. S8 and Table S4†). | ||
| Parameters | Gd·L1 | Gd·L2 | Gd·L3 | Gd·L4 |
|---|---|---|---|---|
| a To fit 1H NMRD data at 298 K, the following parameters were fixed in the fitting procedure: q = 1, rGd-H = 3.0 Å, aGd-H = 4 Å, 298DGd-H = 2.25 × 10−5 cm2 s−1. b From the fitting of the 17O NMR data, with the fixed value EV = 1 kJ mol−1 and ER = 20 kJ mol−1. | ||||
| 298 r 1 20MHz/mM−1 s−1 | 6.5 ± 0.1 | 8.4 ± 0.1 | 7.6 ± 0.1 | 9.1 ± 0.1 |
| 310 r 1 20MHz/mM−1 s−1 | 5.1 ± 0.1 | 6.6 ± 0.1 | 5.9 ± 0.1 | 7.1 ± 0.1 |
| 298 r 1 60MHz/mM−1 s−1 | 5.6 ± 0.1 | 7.6 ± 0.1 | 7.1 ± 0.1 | 8.6 ± 0.1 |
| 310 r 1 60MHz/mM−1 s−1 | 4.5 ± 0.1 | 6.1 ± 0.1 | 5.5 ± 0.1 | 6.6 ± 0.1 |
| Δ 2/1018 s−2 | 9.4 ± 0.1 | 8.1 ± 0.1 | 8.2 ± 0.1 | 7.0 ± 0.1 |
| 298 τ V/ps | 32.5 ± 0.3 | 48.8 ± 0.4 | 29.6 ± 0.4 | 35.1 ± 0.4 |
| 298 τ R/ps | 106 ± 1 | 161 ± 1 | 142 ± 1 | 187 ± 1 |
| 298 τ M/ns | 154 ± 2 | 166 ± 2 | 125 ± 1 | 119 ± 3 |
| ΔH#m/kJ mol−1b | 43.8 ± 0.4 | 43.1 ± 0.4 | 43.5 ± 0.3 | 39 ± 1 |
Upon copolymerization of Gd·L1 with NAM by RAFT, the resulting linear polymers P(NAM-r-Gd·L1) possessed significantly higher relaxivities (Table 4) in the range 12.6–13.5 mM−1 s−1 at 60 MHz, and 14.3–15.4 mM−1 s−1 at 20 MHz (310 K, pH 7.4), up to ca. 5 times higher than Gd-DOTA. These relaxivity values are similar to those obtained previously by Davis, Boyer and coworkers,10,57 for both discrete core crosslinked star polymers and hyperbranched polymers, each containing Gd(III) chelates attached via a single pendant arm.
| Parameters | Linear polymers P(NAM-r-X% Gd·L1) | Hyperbranched polymers P(NAM-r-0.9% Gd·L2–4) | |||||
|---|---|---|---|---|---|---|---|
| 1% | 5% | 9% | 17% | Gd·L2 | Gd·L3 | Gd·L4 | |
| a To fit the 1H NMRD data at 298 K, the following parameters were fixed: q = 1, rGd-H = 3.0 Å, aGd-H = 4 Å, 298DGd-H = 2.25 × 10−5 cm2 s−1. b Estimated from VT NMR relaxivity profiles (Fig. 5) at fixed magnetic field (20 MHz). | |||||||
| 298 r 1 20MHz/mM−1 s−1 | 18.6 ± 0.2 | 17.0 ± 0.2 | 17.1 ± 0.2 | 17.6 ± 0.2 | 30.4 ± 0.4 | 32.7 ± 0.5 | 26.0 ± 0.4 |
| 310 r 1 20MHz/mM−1 s−1 | 15.4 ± 0.2 | 14.3 ± 0.1 | 14.4 ± 0.2 | 15.0 ± 0.1 | 30.7 ± 0.4 | 33.5 ± 0.5 | 23.0 ± 0.4 |
| 298 r 1 60MHz/mM−1 s−1 | 15.8 ± 0.2 | 14.8 ± 0.1 | 14.8 ± 0.1 | 15.2 ± 0.1 | 20.7 ± 0.3 | 22.1 ± 0.2 | 20.9 ± 0.2 |
| 310 r 1 60MHz/mM−1 s−1 | 13.5 ± 0.2 | 12.6 ± 0.1 | 12.9 ± 0.1 | 13.2 ± 0.1 | 20.7 ± 0.3 | 22.8 ± 0.2 | 18.6 ± 0.2 |
| Δ 2/1018 s−2 | 4.37 ± 0.06 | 5.94 ± 0.07 | 5.61 ± 0.04 | 5.31 ± 0.05 | 6.2 ± 0.2 | 6.3 ± 0.4 | 5.7 ± 0.3 |
| 298 τ V/ps | 42.1 ± 0.6 | 38.1 ± 0.5 | 38.2 ± 0.3 | 41.8 ± 0.4 | 16.1 ± 1 | 10.4 ± 0.5 | 12.6 ± 0.3 |
| 298 τ RL/ps | 342 ± 6 | 369 ± 5 | 313 ± 3 | 346 ± 4 | 501 ± 45 | 418 ± 99 | 447 ± 34 |
| 298 τ RG/ps | 2680 ± 90 | 2870 ± 110 | 2420 ± 50 | 2500 ± 60 | 6366 ± 486 | 6668 ± 812 | 3810 ± 256 |
| S 2 | 0.175 | 0.118 | 0.163 | 0.159 | 0.500 | 0.601 | 0.403 |
| 298 τ M/nsb | 330 ± 5 | 330 ± 5 | 330 ± 3 | 330 ± 4 | 326 ± 8 | 308 ± 12 | 315 ± 9 |
The hyperbranched polymers P(NAM-r-Gd·L2–4), obtained using crosslinkers Gd·L2–4,possessed even higher relaxivities in the range 18.6–22.8 mM−1 s−1 at 60 MHz, and 23.0–33.5 mM−1 s−1 at 20 MHz (310 K, pH 7.4) (Table 4). The enhancements in relaxivity of the polymeric CAs relative to the reference agent Gd-DOTA are shown in Fig. 6. Notably, we observe a substantial 9 to 10-fold increase in relaxivity for the crosslinked polymers based on Gd·L2 and Gd·L3, relative to Gd-DOTA at 20 MHz, and a 2-fold increase relative to the linear polymers. Such high gains in relaxivity can be ascribed to the role of Gd·L2–3 crosslinkers, which reduce the rate of tumbling of the Gd(III) chelate in the resulting hyperbranched polymers. This limited rotational flexibility leads to much higher relaxivity.4 It is also clear from Fig. 6 that the gains in relaxivity for the linear polymers are essentially the same at 20 and 60 MHz, whereas for the crosslinked polymers based on Gd·L2 and Gd·L3, the relaxivity gains are greater when measured at 20 MHz. In contrast, for the crosslinked polymer based on Gd·L4 the relaxivity enhancement at 20 MHz is less substantial. This can be explained by the faster local tumbling motion of the paramagnetic centre of Gd·L4, due to the longer and more flexible crosslinking arms.
![]() | ||
| Fig. 6 Relaxivity gains observed for our linear or hyperbranched polymers, relative to Gd·DOTA (3.4 and 2.9 mM−1 s−1 at 20 MHz and 60 MHz, respectively), determined at 310 K at 20 MHz and 60 MHz. | ||
Nuclear Magnetic Resonance Dispersion (NMRD) profiles, i.e. the variation of relaxivity (r1) as a function of the applied magnetic field strength, were measured at 298 and 310 K and at pH 7.4 in the proton Larmor frequency range 0.01–70 MHz (0.000234–1.64 T). Representative examples of the NMRD profiles of the monomeric complexes, and of the linear and crosslinked polymers are presented in Fig. 5. The monomeric complexes Gd·L1–4 displayed profiles typical for fast tumbling, low molecular weight complexes, each characterized by a steady decrease in their relaxivity at low magnetic field (<1 MHz), a drop in relaxivity between 1 MHz to 10 MHz, followed by a second plateau in the high magnetic field region (>10 MHz), governed by the rotational correlation time, τR. In the high field region, the relaxivity is similar for the four monomers Gd·L1–4 as expected, since their similar molecular weight, size and charge results in similar rotational dynamics. Raw data and fitted NMRD profiles are reported in the ESI (Tables S5–S8 and Fig. S9†).
The NMRD profiles of Gd·L1–4 were fitted according to the established theory of paramagnetic relaxation expressed by the Solomon–Bloembergen–Morgan (SBM) and Freed's equations for the inner- (IS) and outer-sphere (OS) proton relaxation mechanisms, respectively (see ESI, eqn (S1)–(S11)†).58–60 Certain parameters were fixed to reasonable values according to previously reported examples:25,56,61 a hydration number of q = 1 was assumed, and the distance between the Gd(III) ion and a bound water molecule proton, rGdH, was set to 3.0 Å, based on crystallographic data for Gd-DOTA. The closest approach of the bulk water molecules, aGdH, was set to 4.0 Å, and the diffusion coefficient of a water proton away from the Gd(III) centre was assumed equal to DGdH = 2.25 × 10−5 cm2 s−1 at 298 K, or 3.1 × 10−5 cm2 s−1 at 310 K.
To provide better estimates of the rotational correlation times of Gd·L1–4, by fitting of the NMRD profiles, the water residence lifetimes (τM = 1/kex) were determined by 17O NMR relaxometry. Thus, the temperature dependence of the transverse relaxation rate (R2) and chemical shifts (Δωr) were determined at 11.75 T at neutral pH using relatively concentrated solutions of Gd·L1–4 (Fig. 5 and S10, Table S9–S12†) and the profiles were fitted according to the Swift–Connick theory for 17O relaxation.62,63 For complexes Gd·L1–4 the τM values were similar (τM = 119–166 ns) and in line with the values determined for similar Gd-DOTAGA (DOTAGA = 2-(4,7,10-tris-carboxymethyl-1,4,7,10-tetraazacyclododecan-1-yl)pentanoic acid) derivatives under the same conditions.25 However, it appears that the complexes Gd·L3 and Gd·L4, bearing cis-methacrylamide arms, exhibit slightly faster water exchange.
The NMRD profiles obtained for the linear and crosslinked polymers revealed a significant increase in relaxivity over the entire proton Larmor frequency range (Fig. 5), compared with monomeric complexes Gd·L1–4 (raw data and fitted NMRD profiles are reported in Tables S13–S19 and Fig. S11 and S12†). The profile shapes were also distinctly different, indicating successful incorporation of the monomeric complexes into higher molecular weight macromolecules. The profiles of the linear copolymers P(NAM-r-Gd·L1) displayed a broad relaxivity peak in the Larmor frequency range of 10–60 MHz, whereas in the case of the hyperbranched (crosslinked) polymers, the relaxivity peak was sharper in the 20–50 MHz region, indicating slower tumbling of the crosslinked polymers compared with the linear polymers.
To obtain more accurate fitting of the NMRD profiles of the polymeric systems, τM values were estimated by fitting of the variable temperature 1H NMR profiles at 20 MHz for the linear polymers containing between 1–17 mol% Gd·L1 theoretically, and for all hyperbranched polymers (Table S22–S24†). In fact, the relatively low concentration of Gd(III) in the polymeric systems prevented the acquisition of 17O NMR data, which typically require Gd(III) concentrations in excess of 5 mM. The water exchange rate of the coordinated water molecule of the polymerized Gd(III) complexes were determined to be two times slower than the corresponding monomeric complexes (τM values between 300–330 ns). Slower water exchange kinetics has been observed previously for other macromolecular Gd(III) systems,25,27 and in the current study this may be tentatively ascribed to weak non-covalent interactions between the Gd(III) chelate and polymer backbone. Alternatively, a reduced rate of water diffusion through the polymer could also contribute to the slower water exchange rate.10,64 For linear polymers with 5 and 9 mol% of Gd·L1, the water residence lifetime was assumed to be 330 ns, since these polymers have similar molecular weights, dimensions and relaxivities to those containing 1 and 17 mol% Gd·L1.
The NMRD profiles for the linear and crosslinked polymers were fitted based on SBM theory and modified with the Lipari–Szabo approach for the description of the rotational dynamics of Gd(III) chelates covalently linked to macromolecules (see equations in Section 5).65,66 In particular, the contributions of the fast local tumbling motion of the Gd-chelate (τRL) were separated from the slower global tumbling of the macromolecule (τRG). τRL and τRG are associated with the order parameter, S2, which describes the level of interconnectivity between the local and the global motions (i.e., if S2 = 0 the motions are independent, if S2 = 1 the motions are fully linked).
The NMRD profiles of the linear polymers were fitted over the entire range of magnetic fields investigated (0.01 to 70 MHz, Fig. 5-B2 and S11†), whereas for the hyperbranched systems, the profiles were better fitted using only the high field data, i.e. above 1 MHz (Fig. 5-B3 and S12†),‡ as commonly performed for large macromolecular MRI CAs.1
Compared with the discrete complexes Gd·L1–4, the linear polymers displayed slower local (τRL greater than 2-fold longer) and slow global (τRG greater than 20-fold longer) reorientation correlation times (Table 4). This can be ascribed to the incorporation of Gd·L1 into polymer chains via the pendant arm of the macrocyclic ligand. The order parameter, S2, obtained for the linear polymers ranged between 0.12 and 0.18, which is reasonable for macromolecules containing Gd-chelates conjugated via a single flexible linker, which allows relatively fast local tumbling. The hyperbranched polymers containing Gd·L2–4 showed a large increase in relaxivity (18.6–22.8 mM−1 s−1 at 60 MHz), primarily attributed to the slower tumbling of the crosslinked Gd(III) complexes. In fact, both τRL and τRG increased relative to the linear polymers (Table 4), and the order parameter S2, was also much higher than for the linear polymers (S2 ≈ 0.60 for P(NAM) containing 0.9% Gd·L3), consistent with a more restricted motion of the Gd(III) chelate within the crosslinked systems.
Further inspection of the relaxivity data revealed that for the crosslinked polymers prepared from Gd·L2 and Gd·L3, there is very little difference in relaxivity at 298 K and 310 K (Fig. 5-A3). This cannot be explained by differences in water exchange rate, since this parameter is very similar for the different polymers synthesised (Table 4). Rather, this can be attributed to more effective coupling of the local and global tumbling motion (S2 up to 0.60) in the crosslinked systems, hence we do not lose any relaxivity gains at higher temperature. In the case of the crosslinked polymer synthesised from Gd·L4, the local and global motion is less effectively coupled (S2 = 0.40), and consequently the relaxivity at 310 K is approximately 11% lower than at 298 K. A similar decrease in relaxivity is observed for the linear polymers at 310 K, where the motion of the Gd(III) centre is less efficiently coupled to the more slowly tumbling polymer (S2 = 0.12–0.18).
| Comparison with linear systems | Comparison with hyperbranched systems | |||||
|---|---|---|---|---|---|---|
| This work | Ref. 25 | Ref. 27 | This work | Ref. 25 | Ref. 67 | |
| a r 1 values are given at 60 MHz. Relaxivities of literature examples were estimated from NMRD profiles. For ref. 25, a range of relaxivities is given, accounting for differences observed depending on the type of assembly formed (e.g. micelles or liposomes). b Water residence times (τM) were estimated from VT NMR relaxivity profiles at fixed magnetic field (20 or 40 MHz). | ||||||
| Description | Linear copolymer formed using monomer Gd·L1 | Self-assembled structures from an amphiphilic Gd complex | Fibril-shaped nanoparticles from block copolymers | Hyperbranched polymer made from crosslinker Gd·L3 | Self-assembled structures made from a rigid amphiphilic Gd complex | Hyperbranched (dendrimer like) amino-functionalized polyglycerol |
| Name | P(NAM-r-Gd·L1) | Gd-DOTAGAC12 | FMN | P(NAM-r-Gd·L3) | Gd-DOTA(GAC12)2 | HB-PG-Gd(DOTA-pBn) |
| 298 r 1/mM−1 s−1 | 15.8 | ≈ 13–14.5 | — | 22.1 | ≈ 27–29 | ≈ 25 |
| 310 r/mM−1 s−1 | 13.5 | ≈ 11–11.5 | ≈ 13 | 22.8 | ≈ 23–25 | — |
| Δ 2/1018 s−2 | 4.37 ± 0.06 | 4.9 | 7 | 6.3 ± 0.4 | 5.2 | — |
| 298 τ V/ps | 42.1 ± 0.6 | 11 | 53 | 10.4 ± 0.5 | 13 | — |
| 298 τ RL/ps | 342 ± 6 | 210 | 150 | 418 ± 99 | 820 | 530 |
| 298 τ RG/ps | 2680 ± 90 | 2900 | 2800 | 6668 ± 812 | 4700 | 4000 |
| S 2 | 0.175 | 0.14 | 0.25 | 0.601 | 0.7 | 0.36 |
| 298 τ M/nsb | 330 ± 5 | 220 | 350 | 308 ± 12 | 297 | 152 |
The hyperbranched polymers containing Gd·L2–4 showed a large increase in relaxivity (18.6–22.8 mM−1 s−1 at 60 MHz). Interestingly, the values of the rotational correlation times and of S2 are comparable to those reported for micellar aggregates obtained by self-assembly of a Gd·DOTAGA2 complex bearing two C12 aliphatic chains in cis-position,25 as seen in Table 5. Also in that example, the restricted local motion of the Gd(III) complex was responsible for a strong relaxivity enhancement with respect to analogous micelles embedding a Gd(III) complex bearing only one aliphatic chain. Our hyperbranched polymers also displayed similar relaxivities and NMRD parameters to those obtained for hyperbranched dendrimers, conjugated with Gd·DOTA-pBn via a single arm (r1 ≈ 25 mM−1 s−1, 298 K).67
Comparing our hyperbranched polymers at 60 MHz and 310 K (Table 4), the system prepared from crosslinker Gd·L3, bearing the shortest pendant arms in a cis orientation, displayed a higher relaxivity (22.8 mM−1 s−1), than systems prepared from the trans-oriented crosslinker Gd·L2 (20.7 mM−1 s−1) or the cis-oriented crosslinker Gd·L4 with longer arms (18.6 mM−1 s−1). This indicates that the combination of shorter crosslinker arms in a cis-geometry is most ideal for limiting the local motion of the Gd(III) complex within the hyperbranched macromolecules.
Of the very few reported macromolecular CAs wherein the Gd(III) complex behaves as a crosslinker,68–70 our system is the only example which takes advantages of the crosslinking to efficiently reduce the local rotational tumbling. This leads to higher relaxivity than those previously reported (e.g. crosslinked acrylamide nanogels68 bearing Gd-DOTA or DTPA like ligands show r1 = 9.7–17.6 mM−1 s−1, 60 MHz, 310 K). To the best of our knowledge, only one other crosslinked system displays a slightly higher relaxivity (24.1 mM−1 s−1 at 60 MHz, 310 K),69 which we propose is due to the slower global tumbling of the nanoparticles (average diameter of 65 nm), and faster water exchange (since all the Gd(III) complexes are on the outside of the nanoparticle), despite exhibiting faster local tumbling.
Copolymerization of Gd·L1, bearing a single methacrylamide arm, with NAM led to the formation of linear polymers with higher relaxivities (12.6–13.5 mM−1 s−1 at 60 MHz) and slower tumbling compared with the discrete monomeric complexes Gd·L1–4. Moreover, hyperbranched polymers prepared via the incorporation of crosslinked Gd(III) chelates Gd·L2–4 displayed significantly higher relaxivities and slower tumbling compared with the linear polymers. Analysis of the NMRD profiles revealed that the higher relaxivities of the hyperbranched polymers is due to the restricted motion of the crosslinked Gd(III) chelates, which approximately doubles the local and global reorientation correlation times relative to the linear polymers containing Gd·L1. Crucially, the global motion of the hyperbranched polymers was more effectively coupled to the motion of the paramagnetic centre. This is apparent from an increase in the order parameter, S2, relative to the linear polymers, thus showing that the rotational flexibility was significantly reduced for polymers containing crosslinkers Gd·L2–4.
Hyperbranched polymers prepared from Gd·L3 displayed the highest relaxivity (22.8 mM−1 s−1 at 60 MHz), suggesting that the combination of shorter polymerizable arms in a cis-orientation is optimal for limiting the local motion of the Gd(III) complex within a hyperbranched polymer. In comparison, hyperbranched polymers prepared from Gd·L4, bearing more flexible arms in a cis-geometry, or from Gd·L2 bearing two trans-related polymerizable arms, showed lower relaxivities (18.6–20.7 mM−1 s−1 at 60 MHz). The results obtained herein will guide the design of second generation Gd(III) monomers possessing three or four polymerizable arms, in order to access macromolecules with even higher relaxivities.
Our synthetic approach to macromolecular CAs combines the simplicity of a single polymerization step (with no post-polymerization modification) and scalability. The monomeric complexes Gd·L1–4 serve as building blocks for the construction of more complex polymeric MRI CAs possessing responsive or theragnostic properties.71 Further, polymeric CAs capable of in vivo targeting may be addressed through the copolymerization of Gd·L1–4 with monomers containing water-soluble sugar moieties or small peptide sequences, which modulate the second sphere of hydration. Work in this regard is ongoing in our laboratories.
:
1 v/v) to give compound 1 as a colourless oil (3.32 g, 8.60 mmol, 46% over 2 steps). IR (νmax/cm−1, neat): 3341, 2940, 1705, 1528, 1250, 1141, 741, 694. Rf (Pet. ether/EtOAc v/v 70
:
30): 0.42. 1H NMR (400 MHz, CDCl3) δ (ppm): 7.49–7.28 (m, 5H, CH aromatic), 5.08 (s, 2H, CH2Ph), 4.12 (t, 3J = 7.2 Hz, 1H, CH), 3.22 (app. q, 3J = 6.7 Hz, 2H, NHCH2), 2.15–1.80 (m, 2H, CH2CHBr), 1.75–1.50 (m, 2H, CH2CH2CH2), 1.46 (s, 9H, CH3). N–H signals are not observed. 13C NMR (101 MHz, CDCl3) δ (ppm): 168.8 (CO2tBu), 156.5 (NHCO2CH2), 136.6 (CIV aromatic), 128.7, 128.3 (CH aromatic), 82.7 (C(CH3)3), 66.9 (CH2Ph), 47.2 (CH), 40.3 (NHCH2), 32.1 (CH2CHBr), 27.9 (CH3 and CH2CH2CH2). HRMS (ESI+, m/z) calculated for M = C17H2479BrNO4, [M + Na]+: 408.0781, found 408.0779.
:
8 v/v, with an increment of 2%) to yield compound 3 (1.76 mg, 1.74 mmol, 87%) as a yellow solid. IR (νmax/cm−1, neat): 3254, 2975, 2933, 2837, 1711, 1517, 1227, 1154, 1113. Rf (CH2Cl2/CH3OH, v/v 90
:
10): 0.32. 1H NMR (400 MHz, CDCl3) δ (ppm): 7.40–7.26 (m, 10H, CH aromatic), 5.06 (s, 4H, CH2Ph), 3.42 (d, 2J = 17.1 Hz, 2H, CO2CH2N), 3.32 (m, 2H, CH), 3.19 (m, 4H, NHCH2CH2), 3.07 (app. t, Japp = 13.3 Hz, 2H, CH2 cyclen), 2.94 (app. t, Japp = 13.3 Hz, 2H, CH2 cyclen), 2.76 (d, 2J = 17.1 Hz, 2H, CO2CH2N), 2.65–2.25 (m, 8H, 4 × CH2 cyclen), 2.20 (app. d, Japp = 13.3 Hz, 2H, CH2 cyclen), 2.08 (app.d, Japp = 13.3 Hz, 2H, CH2 cyclen), 1.86 (m, 2H, CH2CH2CH2), 1.69 (m, 2H, CH2CHBr), 1.53–1.65 (m, 4H, overlap CH2CHBr + CH2CH2CH2), 1.53–1.30 (m, 36H, CH3). N–H signals are not observed. 13C NMR (101 MHz, CDCl3) δ (ppm): 174.9, 172.9 (CO2tBu), 156.6 (PhCH2CO2), 136.8 (CVI aromatic), 128.5, 128.0 (CH aromatic), 82.1, 81.9 (C(CH3)3), 66.4 (CH2Ph), 61.0 (CH), 56.0 (NCH2CO2), 52.7, 48.8, 47.3, 44.6 (CH2 cyclen), 40.8, 40.7 (NHCH2CH2), 29.8 (NHCH2CH2), 28.3, 28.2, 27.9 (CH3), 22.0 (CHCH2). HRMS (ESI+, m/z) calculated for M = C54H86N6O12, [M + Na]+: 1033.6196, found 1033.6189.
:
10 v/v) to give compound 7 (3.32 g, 10.3 mmol, 17% over 2 steps) as a clear oil. IR (νmax/cm−1, neat): 2978, 1726, 1367, 1255, 1136, 843, 755. Rf (Pet.ether/EtOAc, v/v 90
:
10): 0.34. 1H NMR (400 MHz, CDCl3) δ (ppm): 4.24 (dd, 3J = 8.4, 5.8 Hz, 1H, CH), 2.46–2.36 (m, 2H, CH2CO2), 2.27–2.16 (m, 2H, CH2CHBr), 1.49 (s, 9H, CHBrCO2tBu), 1.46 (s, 9H, CH2CO2tBu). 13C NMR (101 MHz, CDCl3) δ (ppm): 171.5 (CHBrCO2), 168.6 (CO2CH2), 82.6 (CH2CO2C(CH3)3), 80.9 (BrCHCO2C(CH3)3), 47.0 (CHBr), 32.9 (CH2CHBr), 30.0 (CH2CO2tBu), 28.2 (BrCHCO2tBu), 27.8 (CH2CO2C(CH3)3). HRMS (ESI+, m/z): calculated for M = C13H2379BrO4, [M + H]+: 345.0672, found 345.0671. Spectral data were in accordance with that reported in the literature.25
:
10 v/v, with an increment of 2%) to give compound 7 (823 mg, 0.93 mmol, 75%) as a yellow amorphous solid. IR (νmax/cm−1, neat): 3322, 2974, 2932, 2832, 1718, 1366, 1227, 1147. 1H NMR (400 MHz, CDCl3) δ (ppm): 3.41 (app. d, 3J = 10.0 Hz, 2H, 2 × CH), 3.34 (d, 2J = 17.3 Hz, 2H, NCH2CO2), 3.25–2.85 (m, 4H, 2 × CH2 cyclen), 2.79 (d, 2J = 17.3 Hz, 2H, NCH2CO2), 2.60–2.00 (m, 12H, 6 × CH2 cyclen), 2.00–1.80 (m, 4H, CO2CH2CH2CH), 1.65–1.50 (m, 4H, CO2CH2CH2CH), 1.50–1.35 (m, 54H, 18 × CH3). 13C NMR (101 MHz, CDCl3) δ (ppm): 174.9, 172.9, 172.8, 172.4, 172.3, 172.0 (CO2tBu), 82.4, 82.3, 81.9, 81.8, 80.6, 80.5 (C(CH3)3), 60.2, 60.1 (CH), 55.8, 55.7, 55.0 (NCH2CO2tBu), 52.6, 48.6(2), 48.5(7), 47.4, 47.3, 44.7, 44.5 (CH2 cyclen), 33.7 (CO2CH2CH2CH), 31.8 (CO2CH2CH2CH), 28.3, 28.2, 28.1, 28.0, 27.9(3), 27.8(6) (CH3). HRMS (ESI+, m/z) calculated for M = C46H84N4O12, [M + H]+: 885.6159, found 885.6158.
:
Gd·L1
:
CTA equal to 95
:
5
:
1.
4-Acryloylmorpholine (67.1 mg, 60.0 μL, 0.5 mmol, 95.0 equiv.), Gd·L1 (17.1 mg, 25.0 μmol, 5 equiv.), trioxane (3.8 mg, 41.7 μmol, 8.3 equiv.), cyano-4-[(dodecylsulfanyl-thiocarbonyl)-sulfanyl]pentanoic acid (2.0 mg, 5.0 μmol, 1.00 equiv.) and AIBN (0.2 mg, 1.0 μmol, 0.20 equiv.) were stirred in DMSO/H2O (80
:
20 v/v, 440 μL), the overall reaction volume being 0.5 mL, in a Schlenk tube equipped with a stirring bar. The reaction mixture was stirred at room temperature for 5 minutes and degassed with N2 through 3 successive freeze–pump–thaw cycles. The reaction mixture was then stirred at 80 °C overnight (15 h). The resulting polymer was purified by extensive dialysis against water (6 × 4 hours) and lyophilization to give pure polymer as an amorphous white solid. IR (νmax/cm−1, neat): 3433, 2963, 2856, 1626, 1439, 1233, 1111. All P(NAM-r-Gd·L1) copolymers displayed similar IR spectra and showed only small differences (±30 cm−1) in their stretches νmax (ESI Section 2). 1H NMR (400 MHz, CDCl3) δ (ppm): 4.25–3.10 (8H, NCH2CH2O), 2.30–2.9 (1H, CH backbone), 2.00–1.00 (2H, CH2 backbone). All P(NAM-r-Gd·L1) copolymers displayed similar NMR spectra with identical peaks numbers and peak shifts. The only noticeable difference is an increase in peak broadness with an increase of the Gd·L1 percentage in copolymer formulation (ESI Section 2). 13C NMR (101 MHz, CDCl3) δ (ppm): 172.9 (CONH), 66.9, 66.4 (CH2O), 46.1, 42.3 (CH2N), 35.5 (CH backbone), 34.8 (CH2 backbone). All linear P(NAM-r-Gd·L1) copolymers displayed similar 13C NMR spectra, with identical numbers of resonances and similar chemical shifts.
:
20 v/v, 0.5 mL overall) in a Schlenk tube equipped with a stirring bar. The reaction mixture was stirred at room temperature for 5 minutes and degassed with N2 through 3 successive freeze–pump–thaw cycles. The reaction mixture was then stirred at 80 °C overnight (15 h). The resulting polymer was purified by extensive dialysis against water (6 × 4 hours) to give, after lyophilization, pure polymer as an amorphous white (to slightly yellow) solid. IR (νmax/cm−1, neat): 3433, 2963, 2856, 1630, 1436, 1231, 1111, 1029. All hyperbranched P(NAM-r-Gd·L2–4) polymers displayed similar IR spectra and showed only small differences (±30 cm−1) in their stretches νmax. 1H NMR (400 MHz, CDCl3) δ (ppm): 4.25–3.10 (8H, NCH2CH2O), 2.30–2.9 (1H, CH backbone), 2.00–1.00 (2H, CH2 backbone). All hyperbranched P(NAM-r-Gd·L2–4) copolymers displayed similar NMR spectra with identical peaks numbers and peak shifts. 13C NMR (101 MHz, CDCl3) δ (ppm): 172.9 (CONH), 66.9, 66.4 (CH2O), 46.1, 42.3 (CH2N), 35.5 (CH backbone), 34.8 (CH2 backbone). All linear P(NAM-r-Gd·L2–4) copolymers displayed similar 13C NMR spectra, with identical numbers of resonances and similar chemical shifts.
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc04750c |
| ‡ NMRD profiles for the hyperbranched systems were also fitted over the entire Larmor frequency range (0.01–100 MHz) and are provided in the ESI (Fig. S13 and Table S21). The resulting NMRD parameters are very similar to those obtained from the high field fitting, presented in Fig. 5. |
| This journal is © The Royal Society of Chemistry 2021 |