Georgina L.
Gregory
a,
Elizabeth M.
Hierons
a,
Gabriele
Kociok-Köhn
a,
Ram I.
Sharma
b and
Antoine
Buchard
*a
aDepartment of Chemistry, University of Bath, Bath BA2 7AY, UK. E-mail: a.buchard@bath.ac.uk
bDepartment of Chemical Engineering, University of Bath, Bath BA2 7AY, UK
First published on 13th February 2017
The development of biodegradable polymers from renewable resources is vital in addressing the dependence of plastics on petroleum-based feedstocks and growing ocean and landfill waste. Herein, both CO2 and natural sugar diols are utilised as abundant, safe and renewable building blocks for the synthesis of degradable and biocompatible aliphatic polycarbonates. Despite a strong potential for advanced polymer properties, inspired by Nature's supramolecular base-pairing, polycarbonates from the sugar components of DNA, 2′-deoxyribonucleosides have been limited by the inability of phosgene derivatives to form the cyclic carbonate monomers that would allow for controlled ring-opening polymerisation. CO2 insertion at 1 atm pressure into methylated thymidine 2′-deoxyribonucleoside, facilitated by organic base 1,8-diazabicyclo-[5.4.0]-undec-7-ene, affected an intramolecular SN2-like displacement of a tosyl leaving group to yield the cyclic carbonate by stereochemical inversion. Organocatalytic ring-opening polymerisation proceeded rapidly in solution resulting in high monomer conversions of 93% and number-average molecular weights, substantially greater and more controlled than via polycondensation routes. The thermodynamic parameters of the polymerisation (ΔHp = −12.3 ± 0.4 kJ mol−1 and ΔSp = −29 ± 1.1 J mol−1 K−1) were determined from the equilibrium monomer conversions over a temperature range of 0 to 80 °C and pseudo-first order kinetics demonstrated. The amorphous thymidine-based polycarbonates exhibited high glass transition temperatures of 156 °C and were found to be highly degradable to the constituent diol under basic aqueous conditions. Static water contact angle measurements and cell studies with MG-63 cell line indicated slightly hydrophilic and biocompatible materials, promising for tissue-engineering applications. The novel, CO2-driven approach to cyclic carbonate synthesis represents a means of expanding the scope of sugar-based monomers for tailored material properties derived from natural products.
Polycarbonates have been widely prepared by traditional polycondensation with phosgene, phosgene derivatives or dialkyl carbonates.18 This however overlooks the direct utilisation of CO2, involves hazardous reagents and often presents challenges in molecular weight control. The copolymerisation of epoxides with CO2 has been extensively studied as a greener alternative but can be limited by the substrate scope particularly with regards to carbohydrate feedstocks.19 Enabling larger and well-defined number average molecular weights (Mn) with narrow dispersities (Đ), the controlled ring-opening polymerisation (ROP) of typically 6-membered cyclic carbonates has emerged as an attractive choice in polycarbonate production and often proceeds under mild reaction conditions.20 Both for polymer synthesis and various other applications (e.g. battery electrolytes),21 a growing research effort has focused on the improved synthesis of cyclic carbonates that avoids the use of phosgene or phosgene derivatives and directly incorporates CO2 (ideally at ambient pressure and temperature).22
Currently, sugar-based cyclic carbonates for ROP are limited to those derived from D-glucose,23D-mannose24 and D-xylose25 (Fig. 1). Cyclocarbonation of trans-configured 1,3-diols in 5-membered furanose sugars such as DNA 2′-deoxyribonucleosides has never been achieved to the best of our knowledge. In 1972, Tittensor and Mellish26 reported the attempted synthesis of the 3′,5′-cyclic carbonate of thymidine and suggested the need for an active ester at the 3′-position rather than 5′-hydroxyl group for cyclisation to occur. Subsequently in 2011, Suzuki et al.10b selectively placed a carbonyldimidazole (CDI) component at the 3′-position of thymidine and adenosine nucleosides. However, only polycondensation resulted, leading to thymidine-based polycarbonates of low Mn (4700 g mol−1, Đ 1.86).
![]() | ||
Fig. 1 Unlike D-mannose, D-glucose and D-xylose derived monomers, the cyclic carbonate of thymidine cannot be prepared by traditional routes. |
Herein, the synthesis and successful ROP of a thymidine-based cyclic carbonate incorporating CO2 at 1 atm pressure is reported. Motivated by the work of Jessop et al.,27 using 1,8-diazabicyclo-[5.4.0]-undec-7-ene (DBU) reagent to capture CO2 into alcohols, we employ the resulting carbonate as a nucleophile to affect cyclisation via an intramolecular SN2-type reaction. We present DFT calculations that rationalise the monomer synthesis and its ability to be polymerised. The thermodynamic and kinetic parameters of the polymerisation are reported, as well as the resulting new polymer microstructure and thermal properties. A preliminary tissue engineering study is finally presented.
However, DFT calculations28 of the ring-closing kinetics (transition state ΔΔG‡ = +13.3 kcal mol−1) and thermodynamic driving force (ΔΔG = −20.8 kcal mol−1) suggested that this pathway should be accessible at ambient reaction conditions (Scheme S2†). Calculations also indicated that the ring-strain of the hypothetical monomer, consisting of a 6-membered cyclic carbonate trans-fused to a furanose sugar ring, was much higher compared to when fused to the 6-membered pyranose ring in reported D-mannose and D-glucose based monomers (Schemes S4 and S5†). For example, ring-opening with methanol is thermodynamically favored with ΔΔG = −6.7 kcal mol−1 compared to −0.9 and −0.2 kcal mol−1 for the mannose and glucose monomers, respectively. The calculated enthalpy of isodesmic ring-opening with dimethyl carbonate (ΔΔH) is also 5.5 kcal mol−1 more favorable for the trans-fused furanose monomer compared to the trans-fused D-mannopyranose cyclic carbonate. The pyranose ring can adopt a chair conformation that can accommodate a trans-fused ring whereas a furanose-cored monomer may be too highly strained to be isolated (especially in the presence of alcohol). The cis-3′,5′-thymidine analogue (cis-1, Fig. 2) however, lies 9.2 kcal mol−1 lower in energy compared to the trans- and shows a similar isodesmic ring-strain and in fact greater thermodynamic driving force for ring-opening compared to the D-xylofuranose derived monomer prepared by Gross and co-workers.25 In that report, cyclocarbonation with ethylchloroformate of the already cis-3,5-diol in D-xylose yielded a cis-fused cyclic carbonate that successfully underwent ROP. Guided by this insight and further DFT calculations, suggesting favorable cyclisation kinetics and thermodynamics (Scheme S3†), we endeavored to invert the sugar stereochemistry of thymidine and prepare the cis-cyclic carbonate by nucleophilic displacement of a tosyl leaving group with a CO2-generated carbonate.
As shown in Scheme 1, selective tosylation of the 3′-hydroxyl position in thymidine was achieved in good yield by a sequential one-pot silylation followed by tosylation reaction. Subsequent deprotection of the 5′-silyl protecting group was carried out under mild reaction conditions with I2 in methanol following the procedure outlined by Vaino and Szarek.29 Unfortunately, protection of the 3-N position of the nucleobase was found to be necessary to avoid undesired side reactions during the cyclisation step and potentially during the ROP. Bases such as DBU are known to lead to formation of 2,3′-anhydrothymidine by intramolecular ring-closing of the nucleobase with 3′-leaving groups.30 As an initial study and proof of concept of the monomer synthesis and ROP potential, the position was methylated in a simple, high yielding step. 3-N-Methyl-3′-tosyl-thymidine 2 was then reacted with DBU under 1 atm CO2 pressure to form a transient amidinium alkyl 5′-carbonate salt, which underwent intramolecular nucleophilic cyclisation with 71% conversion. Subsequent isolation by column chromatography and purification by recrystallisation from toluene gave monomer cis-1 (now referred to as 1 for simplicity), with polymerisation-grade purity, in reasonable yield of 52% (cf. 41% for D-xylose, 25–36% for glucose and 57% for mannose-based cyclic carbonates).
The strategy represents a generally applicable method to cyclic carbonate synthesis from natural furanose sugars bearing a trans-1,3-diol. It extends the range of sugar-based cyclic carbonate monomers beyond those possible with known sugars and classical reagents such as phosgene and its derivatives, which require a nucleophilic addition-elimination mechanism. Although the need to protect the 3-N position of the nucleobase is somewhat limiting, the functionalisation potential of the pyrimidine arm is tremendous, ranging from Re(CO)3 complexes31 to alkyne functionalities for click chemistry.5c,32
As a further demonstration, we also prepared the 3-N-benzoyl derivative of 1 (1-Bz). Benzoylation could be carried out in tandem with the silylation and tosylation, reducing the number of overall synthetic steps (Scheme S1 and Fig. S12–S14†). Literature precedence has furthermore shown post-polymerisation removal of benzoyl protecting groups.10b
The structure of 1 was confirmed by 1H NMR (Fig. 6), 13C{1H} NMR and FTIR spectroscopies as well as electrospray ionisation mass spectrometry (ESI-MS), elemental analysis and single crystal X-ray diffraction. Full NMR assignment was aided by DEPT135, HSQC and COSY experiments (Fig. S1–S5†). In the 13C NMR spectrum, further to the carbonyl environments of the nucleobase, an additional CO resonance was observed at 147 ppm, typical of a 6-membered cyclic carbonate, and a vicinal 3J3′–4′ coupling constant in the 1H NMR spectrum of just 3.9 Hz was consistent with a cis-configuration. X-ray analysis of crystals of 1 grown from toluene (Fig. 3) further corroborated stereochemical inversion at the 3′-carbon to form the 3′,5′-cis-cyclic carbonate in which the furanose ring adopts a strained 4′-endo-3′-exo twist (4T3) conformation. This is likely retained in solution as the ring-methylene H-2′ protons display very different chemical shifts and J-coupling constants consistent with the Karplus equation. In contrast the unconstrained furanose ring in the preceding trans-3′,5′-diol 2 adopts a 2-endo (2E) envelope conformation (see crystal structure in ESI†).
![]() | ||
Fig. 3 ORTEP33 drawing of the crystal structure of 1 with thermal ellipsoids at the 50% probability level (see ESI†). Selected bond lengths (Å) and dihedral angles (°): C(2)–O(2) 1.466 (2), O(2)–C(1) 1.326 (3), O(1)–C(1) 1.202 (3), C(1)–O(3) 1.326 (4), O(3)–C(6) 1.456 (3); C(4)–C(3)–C(2)–C(5) −29.61 (19), C(3)–C(2)–C(5)–O(4) −36.87 (19), C(2)–C(5)–O(4)–C(4) 30.01 (2), C(5)–O(4)–C(4)–C(3) 10.5 (2), O(4)–C(4)–C(3)–C(2) 12.6 (2). |
A linear kinetic plot is consistent with pseudo-first order kinetics, typical for ROP, from which a kapp of 1.39 ± 0.005 h−1 was determined (inset Fig. 4 and S15†). The high melting point of the monomer (204–205 °C) and low thermal stability (178–205 °C, 88% mass loss, Tinf 248 °C) prevented polymerisations from being carried out in the melt. The temperature dependence of the equilibrium polymerization was also investigated at fixed [M]0, catalyst loading and initiator concentration in 1,2-dichloroethane solvent to access a wider temperature range of 0–80 °C. A plot of ln[M]eqversus the reciprocal of the absolute temperature gave an estimate of the polymerisation thermodynamic parameters (Fig. 5): ΔHp = −12.3 ± 0.4 kJ mol−1 and ΔSp = −29 ± 1.1 J mol−1 K−1. For comparison, Endo and co-workers34 reported an enthalpic driving force of −26.4 kJ mol−1 for the anionic polymerisation in THF of unsubstituted trimethylene carbonate and accompanying −44.8 J mol−1 K−1 entropy decrease ([M]0 = 1 mol L−1). A decrease in ΔHp with increased steric bulk of substituted derivatives was attributed to distortion of the polymer backbone as opposed to a conformational effect on the monomer. Interestingly, for polymer recyclability, the polymerisation mixture at 0 °C, which showed a 72% conversion to polymer at equilibrium, resulted in near quantitative recovery of the monomer upon heating to 80 °C.
A linear increase in Mn with percentage monomer conversion whilst maintaining narrow dispersities indicated that the polymerisation proceeded in a controlled fashion up to 60% monomer conversion (Fig. 6). Initially excellent agreement was observed between the calculated Mn and those estimated by size exclusion chromatography (SEC) versus polystyrene standards with CHCl3 eluent. For example, the polymer at 54% monomer conversion had a Mn of 7210 g mol−1 (Đ 1.12) as estimated by SEC (Fig. S17†), only slightly lower than the predicted35 7740 g mol−1. Once the thermodynamic equilibrium was reached and with longer reaction times, a broadening of the polymer distribution was observed due to transesterification reactions.
Matrix-assisted laser desorption/ionisation time-of-flight (MALDI-ToF) mass spectrometry of this polymer showed a major linear polymer series of Mn 7120 g mol−1 and expected 4-methyl benzyl alcohol and O–H end groups (Fig. 7 and S19†). A minor (19%) cyclic species with m/z consistent with no end-groups and an integer number of repeat units (m/z ∼282.25), was also observed due to backbiting of the polymer chain. As the polymer chain grows and monomer concentration decreases, a greater rate of chain backbiting relative to chain propagation may account for deviations observed between the calculated and SEC estimated Mn at higher monomer conversions. Dilution of the reaction conditions to 0.5 mol L−1 in an attempt to favour macrocycle formation resulted in lower Mn (Table 1, entry 1) and in more cyclic species observed by MALDI-ToF analysis (Fig. S21†). The lack of chain end in macrocycles should prevent reversion of the polymer to the monomer at elevated temperatures, however, difficulty in separating the cyclic and linear species prevented us from rigorously investigating the issue. Polymeric macrocycles are also interesting in their own right, including those derived from nucleosides, e.g. as supramolecular complexes with metal coordination characteristics.12
![]() | ||
Fig. 7 MALDI-ToF mass spectrum of polymer (Mn,SEC 7210 g mol−1, Đ 1.12) showing major linear and minor cyclic polymeric species (distribution assigned based on molecular masses and corresponding formulas, see Fig. S19† for details). |
Entry | Temp. (°C) | Cat. | [M]0![]() ![]() ![]() ![]() |
Conv.b (%) | Time (h) |
M
n![]() |
---|---|---|---|---|---|---|
a All polymerisations were carried out at [M]0 = 2.5 mol L−1 in CH2Cl2 unless stated otherwise. b Determined by 1H NMR spectroscopy in CDCl3. c g mol−1, calculated by SEC with CHCl3 eluent versus polystyrene standards (RI detector). d [M]0 = 0.25 mol L−1. e L1 = 2,4-di-tert-butyl-6-{[(2′-dimethylaminoethyl)methylamino]-methyl}phenolate. f BDI-1 = 2-((2,6-diisopropylphenyl)amido)-4-((2,6-diisopropylphenyl)-imino)-2-pentene. | ||||||
1d | 0 | TBD | 100![]() ![]() ![]() ![]() |
68 | 4 | 4260 [1.09] |
2 | 0 | TBD | 150![]() ![]() ![]() ![]() |
90 | 3 | 17![]() |
3e | 25 | L1ZnOEt | 100![]() ![]() |
71 | 20 | 7250 [1.22] |
4f | 25 | (BDI-1) | ||||
ZnEt | 100![]() ![]() ![]() ![]() |
60 | 24 | 5240 [1.15] | ||
5 | 25 | TBD | 100![]() ![]() ![]() ![]() |
90 | 3 | 11![]() |
6 | 0 | TBD | 100![]() ![]() ![]() ![]() |
93 | 3 | 15![]() |
The largest molecular weights, up to 17000 g mol−1 (Đ 1.3) were achieved at 0 °C and [M]0 of 2.5 mol L−1 (Table 1, entry 2). A greater disparity between the theoretical Mn and those determined by SEC was observed at higher [M]0/[I]0 ratios, most likely due to chain back-biting reactions. To explore the influence of the catalytic system on the control of Mn, metal complexes were tested for the ROP of 1. Sn(Oct)2 catalyst with 4-MeBnOH initiator at 120 °C in toluene did not yield homopolymer. Y(OiPr)3 or Al(OTf)3 with BnOH initiator resulted in no control with only low Mn (1800–3910 g mol−1, Đ 1.35–1.39) in DCM, dioxane and toluene solvents at 0, 25, 70 and 120 °C. Although highly active for lactide ROP in CH2Cl2 at rt, the zinc(II) alkoxide catalyst (L1ZnOEt) reported by Williams et al.36 was slower compared to TBD for the ROP of 1, as was the (BDI-1)ZnEt catalyst37 with 4-MeBnOH initiator (Table 1, entries 3 and 4).
![]() | ||
Fig. 8 Comparison of the 1H NMR spectra (400 MHz, CDCl3) of monomer 1 (bottom) and its polymer (top, Table 1, entry 6). |
The thermal properties of the polycarbonates were investigated by thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). The onset of thermal degradation occurred at relatively low temperatures of ∼170 °C, reaching a maximum rate of mass loss (Tinf) around 246 °C for polymers of Mn 11000 and 15
400 g mol−1 (Table 1, entries 5 and 6). High mass losses of 91 and 94%, respectively were achieved by ∼300 °C (Fig. S22†). Analysis of the gases evolved by tandem mass spectrometry confirmed the presence of m/z 44 attributed to the loss of CO2+ (Fig. S24†). The polymers exhibited high glass transition temperatures (Fig. S25–S27†); a Tg of 156 °C was measured for polycarbonate of Mn 15
400 g mol−1 and is attributed to restricted rotation about the furanose ring backbone and to the bulky thymine group. For comparison, furanose-cored and isopropylidene-protected D-xylose polycarbonate exhibited a Tg of 128 °C (for Mn = 13
200 g mol−1).25 Pyranose backboned polycarbonates derived from isopropylidene-protected D-mannose and methyl protected D-glucose showed Tg of 152 °C (13
600 g mol−1)23 and 122 °C (14
700 g mol−1),24 respectively. No transitions associated with a melting domain were observed implying an amorphous polymer, which was further confirmed by powder X-ray diffraction (Fig. S28†).
The hydrolytic stability of the water insoluble polymer was investigated at 0.25 wt% at rt in water, 1 mol L−1 HCl and 1 mol L−1 NaOH, as well as in PBS buffer at 37 °C in presence of Candida antartica lipase B. After 1 week, no dissolution or loss of mass by SEC was observed in the case of water and 1 mol L−1 HCl. NMR analysis of the supernatant in D2O revealed no environments. Under the alkali conditions, the polymer was observed to dissolve after 4 hours. 1H NMR spectroscopy in D2O after 24 h at 4 wt% revealed proton environments consistent with the cis-diol of 3-N-methyl thymidine (Fig. S29 and S30†). Mass spectrometry (HR-MS ESI) after 24 h further confirmed the presence of the diol suggesting complete hydrolytic degradation of the polymer to its constitute components. In the case of lipase B, unknown degradation products were detected in solution after one week. Static water contact angle measurements of polymer thin films drop cast onto glass slides revealed an average contact angle of 48 ± 1°. This was more hydrophobic than the monomer but less hydrophilic compared to the glass control (Fig. S32†).
Cell attachment studies were carried out with a human osteoblast cancer cell line, MG-63, on polymer films drop cast onto glass slides. Cell attachment after a 1 h incubation period (37 °C, 5% CO2) was comparable to the empty tissue culture polystyrene well-plate and glass controls, and indicated the possibility of thymidine-based materials to support cells (Fig. S33–S35†). After a longer 24 h incubation period, cell attachment dropped compared to the empty well-plate control suggesting further modifications to the polymer surface are required for it to serve as an effective tissue engineering scaffold. Nevertheless, the choice of the 3-N protecting group and general applicability of the strategy present wide-scope and future opportunities for novel biomaterials with applications in the tissue engineering and regenerative medicine fields. Other degradable polymers attractive in this area, such as PLA and PLGA, are notoriously difficult to functionalise.38 Therefore, the successful preliminary ROP of benzoyl-protected monomer (1-Bz, see ESI†) offers promise for more elaborate functionality as well as post-polymerisation modifications, including self-assembly with base-paired nucleoside polycarbonates.
Footnote |
† Electronic supplementary information (ESI) available: Experimental and computational details; spectroscopic and crystallographic data for 1 and 2; SEC traces, MALDI-ToF spectra, spectroscopic and thermal (TGA, DSC) data for polymers; DFT calculations data and associated digital repositories; tissue engineering protocol and results. CCDC 1523960 and 1523961. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c7py00118e |
This journal is © The Royal Society of Chemistry 2017 |