Landon J.
Kilgallon
a,
Timothy P.
McFadden
b,
Matthew S.
Sigman
b and
Jeremiah A.
Johnson
*acd
aDepartment of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
bDepartment of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
cKoch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA. E-mail: jaj2109@mit.edu
dBroad Institute of MIT and Harvard, Cambridge, MA 02142, USA
First published on 22nd May 2024
Grubbs 3rd-generation (G3) pre-catalyst-initiated ring-opening metathesis polymerization (ROMP) remains an indispensable tool in the polymer chemist's toolbox. Tricyclononenes (TCN) and tricyclononadienes (TCND) represent under-explored classes of monomers for ROMP that have the potential to both advance fundamental knowledge (e.g., structure-polymerization kinetics relationships) and serve as practical tools for the polymer chemist (e.g., post-polymerization functionalization). In this work, a library of TCN and TCND imides, monoesters, and diesters, along with their exo-norbornene counterparts, were synthesized to compare their behaviors in G3-initiated ROMP. Real-time 1H NMR was used to study their polymerization kinetics; propagation rates (kp) were extracted for each monomer. To understand the relationships between monomer structure and ROMP propagation rates, density functional theory methods were used to calculate a variety of electronic and steric parameters for each monomer. While electronic parameters (e.g., HOMO energy levels) correlated positively with the measured kp values, steric parameters generally gave improved correlations, which indicates that monomer size and shape are better predictors for kp than electronic parameters for this data set. Furthermore, the TCND diester—which contains an electron-deficient cyclobutene that is resistant to ROMP—and its polymer p(TCND) are shown to be highly reactive toward DBU-catalyzed conjugate addition reactions with thiols, providing a protecting- and activating-group free strategy for post-polymerization modification.
TCNs and TCNDs are tricyclic, strained olefins with fused cyclobutane and cyclobutene rings, respectively, connected to the two exo positions of the NB skeleton. Compared to other exo-substituted NB monomers, which often require high temperatures9–11 and/or precious metal catalysts8,12 for their synthesis, exo-TCNs and TCNDs can be conveniently assembled by thermal cycloadditions of activated olefins and quadricyclane with perfect exo selectivity.13,14 While TCN and TCND-based polymers have been synthesized using Mo (Schrock-type initiators)15 and Ru (Grubbs 1st and 2nd generation initiators) initiated ROMP16,17 and studied for applications in gas separations,18–25 no studies of their performance under “living” ROMP conditions initiated by G3 or post-polymerization functionalization have been reported. We hypothesized that they may be competitive with commonly used exo-NBs in terms of ROMP propagation rates, but that they could offer advantages of facile synthesis and, in some cases, post-polymerization modification. Moreover, by accessing a training set of monomers, we sought to investigate if the HOMO level was primarily responsible for propagation rate as proposed by Matson and coworkers, or whether other molecular features (e.g., steric parameters) of the monomers may also impact performance.
Here, we synthesize a library of TCNs and TCNDs that incorporates a broad functional group scope, including imides, diesters, and monoesters. We assess their structures and properties, investigate their ROMP behavior using G3, and compare the results to analogous exo-NB derivatives. We identify trends in the relative propagation rates across these monomers, showing that while their HOMO energies are good predictors of reactivity, readily computed steric parameters (e.g., volume, solvent accessible surface area (SASA), and Sterimol B5) offer improved structure–reactivity correlations across all exo-NB, TCN, and TCND monomer classes. Moreover, we demonstrate that the cyclobutene pendants of p(TCND)s can be excellent substrates for conjugate addition reactions with thiol nucleophiles (i.e., “thiol-Michael-type additions”),26–28 offering a straightforward methodology for efficient post-polymerization modification without the need for protecting groups or monomer activation (e.g., activated esters). These studies demonstrate that TCNs and TCNDs are attractive next-generation monomers for ROMP as (a) they are easy to synthesize in high stereochemical purity in as few as two simple steps from commercial starting materials; (b) they display propagation kinetics similar to or faster than many typical exo-NBs; and (c) they offer facile routes to post-polymerization modification that are not possible for common exo-NBs. Finally, this study presents the first application of computed steric parameters to correlate monomer structure with propagation rates in ROMP, providing insights into the predictive design of future ROMP monomers.
We also sought TCN-NBnsyn for comparison to its anti-diastereomer (Fig. 2b). While TCN-NBnsyn and TCN-NBnanti may potentially be separable by chromatography,30 we identified a kinetic resolution method that leveraged the large difference in imidation rates between the syn (slow) and anti (fast) amic acids formed upon initial ring-opening opening of 1, ultimately giving TCN-NBnsyn in 4% yield over 2 steps with >98:2 d.r. (ESI† Section 2). The relative configurations of TCN-NBnsyn and TCN-NBnanti were unambiguously determined by single-crystal X-ray diffraction (Fig. 2b). TCND derivatives TCND-(CO2Me)2 and TCND-CO2Me were synthesized by direct cycloaddition of quadricyclane with dimethylacetylenedicarboxylate and methyl propiolate, respectively (Fig. 2a).13,31
We note that upon storage under ambient atmosphere for extended periods, TCND-(CO2Me)2 yellows and a new product is observed via1H NMR spectroscopy. A small sample of this impurity was isolated chromatographically; its 1H NMR spectrum contained no olefinic resonances and its high-resolution mass spectrum suggested a mass 16 Da (i.e., 1 oxygen atom) larger than TCND-(CO2Me)2. Single-crystal X-ray analysis of the saponified derivative of this impurity revealed its identity to be the exo-epoxidized variant of TCND-(CO2Me)2 (ESI† Section 11). Thus, we recommend storage of uninhibited TCND-(CO2Me)2 under inert atmosphere, as is often done for olefinic monomers, to avoid its oxidative decomposition.
Finally, for comparison of ROMP propagation rates (vide infra), exo-NB analogs of each TCN and TCND were synthesized following previously reported procedures.32–34 All monomers and intermediates were isolated non-chromatographically, and the key intermediates towards TCNs (1) and TCNDs (quadricyclane) were synthesized on the 16 g and 17 g scales, respectively. Notably, the melting temperatures for all but one (TCND-CO2Me) of the TCNs and TCNDs synthesized in this study were greater than their exo-NB analogs (from 10 °C to >150 °C higher), presumably due to their multicyclic, comparably rigid structures. In particular, all of the TCNs studied here have melting points >40 °C above room temperature, making them generally easy to handle and purify (crystallization).
Polymer | M n (theo.) | M n (SEC) | Đ | Yield | E/Z ratiod |
---|---|---|---|---|---|
a Monomer conversion was >98% in each case. b Theoretical Mn values were calculated assuming 100% monomer conversion. c SEC conditions: DMF (0.025 M LiBr) mobile phase, 60 °C, with differential refractive index and multi-angle light scattering detectors. A 100% mass recovery method was used to calculate Mn. d E/Z ratio was determined by the relative integrations of the trans (downfield) and cis (upfield) backbone olefinic protons in the polymers' 1H NMR spectra.35,36 | |||||
p(TCN-NBnanti) | 56.0 kDa | 52.2 kDa | 1.02 | 97% | 37:63 |
p(TCN-NCy) | 54.4 kDa | 54.2 kDa | 1.02 | 95% | 34:66 |
p(TCN-NHex) | 54.8 kDa | 52.2 kDa | 1.02 | 92% | 37:63 |
p(TCN-(CO2Me)2) | 47.4 kDa | 47.8 kDa | 1.02 | 97% | 33:67 |
p(TCND-(CO2Me)2) | 47.0 kDa | 49.5 kDa | 1.04 | 92% | 30:70 |
p(TCND-CO2Me) | 35.3 kDa | 59.3 kDa | 1.97 | 75% | 37:63 |
SEC analysis of p(TCND-(CO2Me)2) prepared under the same conditions (i.e., 15 min reaction time) showed a minor high molecular weight shoulder peak at 2 × Mn indicative of cross-linking, as was also reported for Mo-initiated ROMP of this monomer (Fig. 3e).15 Quenching the ROMP reactions after a shorter time (2 min) using the same batch of TCND-(CO2Me)2 or a freshly distilled sample of TCND-(CO2Me)2 (discarding the first and last ∼10–15% of the distillate) substantially decreased the size of this shoulder peak. Notably, conversion was unaffected by decreasing the reaction time to 2 min, as evidenced by the similar SEC traces, suggesting that this monomer is highly reactive. Moreover, these results indicate that the cyclobutene component of TCND-(CO2Me)2 remains largely intact under these ROMP conditions, providing an opportunity for further modification post-polymerization (vide infra).
While the synthesis of p(TCND-CO2Me) appeared to be no different than for the other samples prepared here, the polymer product did not re-dissolve after the first precipitation. We hypothesized that polymer-associated residual Ru led to cross-linking reactions upon precipitation; the addition of excess pyridine (∼2000 equiv. relative to G3) to the quenched ROMP reaction suppressed this cross-linking (ESI† Section 3), enabling NMR and SEC analysis of p(TCND-CO2Me); however, SEC analysis showed a significant high molecular weight shoulder (Fig. 3f). Notably, while cyclobutene 1-carboxylic esters, like the cyclobutene fragment of TCND-CO2Me, cannot homopolymerize via G3-initiated ROMP,37 they are known to slowly form eneoic carbenes that are kinetically stabilized by intramolecular chelation.38 Thus, we propose that ROMP of TCND-CO2Me occurs first through rapid polymerization of its NB-like alkene followed by slow cross-linking of its cyclobutene pendants, which are more sterically accessible than for TCND-(CO2Me)2 above.
1H NMR spectroscopy was used to characterize the structures of these p(TCN)s and p(TCND)s in more detail (Fig. 4). Resonances associated with each of the monomer protons could be identified in the polymers. Compared to p(exo-NB) analogs (ESI† Section 3), the backbone olefins of p(TCN)s and p(TCND)s were slightly skewed toward the Z configuration: E:Z ranged from 48:52 to 37:63 for p(exo-NB)s while p(TCN)s and p(TCND)s ranged from 37:63 to 30:70 (Table 1). The spectrum of p(TCND-CO2Me) contained an extra downfield resonance at 6.8 ppm that we attribute to the olefinic proton of the cyclobutene. Otherwise, the spectra generally resembled those of their p(exo-NB) counterparts, with the addition of two proton resonances in the carbocyclic region (∼1.4 ppm–3.5 ppm) for the p(TCN)s.
We anticipated that TCN and TCND monomers would display controlled polymerization behavior due to their high ring strain and low likelihood of backbiting/chain transfer. In support of this hypothesis, Mn,SEC values for p(TCN-NBnanti) agreed well with the theoretical values and increased linearly as a function of conversion up to a degree of polymerization (DP) of 500 (Fig. 5a). Moreover, a plot of ln([TCN-NBnanti]0/[TCN-NBnanti]) was linear as a function of time (Fig. 5b), and Đ remained low throughout the polymerization (Fig. 5a), suggesting minimal irreversible termination or chain transfer. Finally, to demonstrate the active nature of the p(TCN-NBnanti) Ru-alkylidene chain ends, a block copolymer p(TCN-NBnanti)100-b-p(TCN-NHex)100 was synthesized via sequential monomer addition (Fig. 5c). SEC analysis showed a clear peak shift to shorter retention time upon addition of TCN-NHex, suggesting complete chain extension; Đ for the block copolymer remained low (1.02).
Monomer | k p (min−1) | t 1/2 (min) |
---|---|---|
a Conditions: 95 mM monomer (200 equiv.), 2.4 mM pyridine (5 equiv.), 0.48 mM G3 (1 equiv.) and internal standard (mesitylene) in CDCl3 at 25 °C. | ||
NB-NBn | 0.157 ± 0.001 | 4.4 |
TCN-NBnanti | 0.117 ± 0.001 | 5.9 |
TCN-NBnsyn | 0.122 ± 0.001 | 5.7 |
NB-NCy | 0.183 ± 0.001 | 3.8 |
TCN-NCy | 0.101 ± 0.001 | 6.9 |
NB-NHex | 0.140 ± 0.001 | 5.0 |
TCN-NHex | 0.094 ± 0.001 | 7.4 |
NB-(CO2Me)2 | 0.225 ± 0.001 | 3.1 |
TCN-(CO2Me)2 | 0.342 ± 0.001 | 2.0 |
TCND-(CO2Me)2 | 0.622 ± 0.007 | 1.1 |
NB-CO2Me | 1.259 ± 0.014 | 0.6 |
TCND-CO2Me | 0.564 ± 0.010 | 1.2 |
The kp values for imides TCN-NBnanti, TCN-NBnsyn, TCN-NCy, and TCN-NHex were within ∼20% of each other and were slightly lower than their exo-NB imide counterparts. A similar trend was observed for the monoesters, but to a larger extent: kp for NB-CO2Me was 2.2-fold higher than for TCND-CO2Me. Notably, however, the opposite trend was observed for the diester monomers; kp for TCND-(CO2Me)2 was greater than TCN-(CO2Me)2, which was greater than NB-(CO2Me)2.
Comparing across the imides, diesters, and monoesters, the imides had the lowest kp values without exception; the fastest imide (NB-NCy) had a ∼19% lower kp than the slowest diester/monoester (NB-(CO2Me)2). This result agrees well with the trends reported for imide versus ester anchor groups for both small molecule monomers and macromonomers.5–7 Comparison of the monoesters and diesters shows that while the monoester TCND-CO2Me and diester TCND-(CO2Me)2 propagate at nearly the same rate, NB-CO2Me propagates 5.6 times faster than its diester counterpart NB-(CO2Me)2. Altogether, these results demonstrate that TCN and TCND monomers possess similar living behavior and G3-initiated ROMP kp values compared to traditional exo-NBs.
Fig. 7 (a) Modeling workflow for steric parameters. AQME39 (rd-kit40) was used to convert monomer SMILES strings into a representative conformer ensemble.44 The ensemble was subjected to DFT computation and conformer properties were obtained. For each conformer the maximum, Boltzmann average, lowest energy (E) conformer, minimum, property value was obtained. Property values for the Vbur-min conformer at C1 were obtained as well to further probe steric effects. (b) Plot of measured monomer ln(kp) versus their HOMO energy levels calculated from their M06-2X/def2-TZVP optimized geometries. (c) Plot of measured monomer ln(kp) versus their calculated Sterimol B5 length from the minimum buried volume conformer of the adjacent active site olefinic carbon C2 (B3LYP/def2SVP/M06-2X/def2TZVP). |
Grubbs3 and Matson also found that monosubstituted exo-NBs polymerize faster than disubstituted variants; Matson attributed this difference either to the possibility that monosubstituted monomers can form regioisomers where their substituents point away from the ruthenium center (a steric effect),7 or to electronic effects (HOMO level and HOMO/LUMO gap).5–7 While the chelating ability of NB monomer substituents has also been shown to significantly affect the activation parameters of polymerization and may slow the polymerization rate by ground-state stabilization in certain cases,3 the chelating ability of exo-NB monomers seems to have little correlation with their propagation rates.5 We note that NB-CO2Me has a ∼4-fold larger kp than TCN-(CO2Me)2 despite having nearly the same calculated HOMO energy (0.6 kcal difference), and that TCND-(CO2Me)2 and TCND-CO2Me both have higher-lying HOMOs (by ∼4 and ∼6 kcal mol−1, respectively) than NB-CO2Me but significantly lower kp values. Overall, while HOMO energy is a useful tool for predicting ROMP reactivity for exo-NB-like scaffolds, other effects, such as steric effects near the propagation site, may impact the rate as well.
These findings led us to consider alternative readily-computed structural descriptors that may provide better correlations with ROMP monomer performance over a wider structure range, and may then also provide greater predictive power for the design of novel monomers. We began by collecting DFT-derived parameters from an AQME derived conformer ensemble39 using RDKit40 for each monomer. From the conformer ensemble, five representative conformers were selected to extract a range of steric and electronic parameters (Fig. 7a).41,42 The ensemble conformer properties were condensed to include the lowest energy conformer (low E), minimum (min), maximum (max), and Boltzmann derived parameter values. Additionally, we also extracted condensed properties from the conformer presenting the minimum percent buried volume (Vbur-min) parameter value at the active site olefinic carbon (C1), as we hypothesized that steric effects may be a driver of the observed rate differences. Percent buried volume probes the proximal steric environment by placing a 3 Å sphere at the site of interest and measuring the resulting overlap.43 The library of parameters was regressed against ln(kp) to reveal a range of good correlations with various parameters describing steric effects (molar volume (MV), solvent accessible surface area (SASA) and Sterimol B5). These correlations all resulted from the Vbur-min conformer values located at the active site olefinic carbon C1. A representative correlation using Sterimol B5 is depicted in Fig. 7c. This parameter describes the display of steric information for how the adjacent olefinic carbon C2 may present itself to the propagating Ru chain end (Fig. 7a). The other two strong correlations, MV and SASA, represent the shape of the monomer while accounting for the local environment defined by the Vbur-min conformer. Notably, the Sterimol B5 and the MV Vbur-min conformers were found to be collinear, and the same parameters derived from other conformers (i.e., low E, min, and max values) do not perform nearly as well (ESI† Section 13). Overall, these data suggest that the overall size and shape of the monomer has a significant impact on polymerization kinetics.
Notably, dimethyl 2,3-dimethylmaleate (DDM), a less-strained electrophile with a similar substitution pattern, did not react under the CDCl3 conditions described above (Table 3, entry 3), and TCND-CO2Me reached only 3% conversion (Table 3, entry 4). These results correlate with the DFT calculated LUMO energies (M06-2X/def-TZVP) for each molecule (Table 3), suggesting that the rigid structure of TCND-(CO2Me)2 paired with its two electron-withdrawing ester substituents make it a better electrophile for conjugate addition. Meanwhile, as discussed above, these structural features are also key for (1) ensuring a high ROMP kp (by avoiding steric hindrance associated with diesters appended directly to the exo-NB) and (2) avoiding undesired cross-linking during ROMP (as observed for TCND-CO2Me), making TCND-(CO2Me)2 a particularly promising monomer for further exploration.
We hypothesized that the electrophilicity of TCND-(CO2Me)2 would not change significantly after ROMP, rendering p(TCND-(CO2Me)2) a similarly useful substrate for post-polymerization conjugate addition reactions. This hypothesis was supported by DFT calculations (M06-2X/def-TZVP); the LUMO energy of a model substrate designed to mimic the p(TCND-(CO2Me)2) backbone—the ethenolysis product of TCND-(CO2Me)2—was very similar (−32 kcal mol−1) to that of TCND-(CO2Me)2 (−30 kcal mol−1). Moreover, nearly quantitative conjugate addition was observed when p(TCND-(CO2Me)2) was exposed to dodecanethiol (5 equiv.) and DBU (20 mol%) in DMF solvent for 1 h under ambient atmosphere (Table 4, entry 1). The reaction solution became cloudy almost immediately after addition of the DBU, suggesting rapid conjugate addition and limited solubility of the functionalized polymer product in DMF. Decreasing the number of equivalents of thiol to 2 (Table 4, entry 2) and lowering the catalyst loading to 5 mol% (Table 4, entry 3) did not affect the conversion. Reduction of the equivalents of thiol to 1.5 (Table 4, entry 4) did not affect conversion either; however, further decreasing the amount of thiol to 1.1 equiv. gave slightly reduced conversion (95%, Table 4, entry 5), which was reduced further with a lower catalyst loading of 2.5 mol% (93%, Table 4, entry 6). Over 95% conversion was achieved with 1.2 equiv. of thiol and 5 mol% DBU (Table 4, entry 7); these conditions were chosen as optimal for further studies.
Entry | H25C12SH equiv. | DBU equiv. | Conversionb (1H NMR) |
---|---|---|---|
a Conditions: p(TCND-(CO2Me)2) in DMF (50 mg mL−1) at room temperature. Dodecanethiol was added, followed by the addition of DBU as a stock solution in DMF (1:9 v/v DBU:DMF) and was allowed to react for 1 hour. b Conversion was assessed by comparing the integration of the backbone olefinic protons (4.8 ppm – 5.7 ppm) versus the integration of the terminal –CH3 of the dodecanethiol side-chain (0.9 ppm). | |||
1 | 5 | 20 mol% | >95% |
2 | 2 | 20 mol% | >95% |
3 | 2 | 5 mol% | >95% |
4 | 1.5 | 5 mol% | >95% |
5 | 1.1 | 5 mol% | 95% |
6 | 1.1 | 2.5 mol% | 93% |
7 | 1.2 | 5 mol% | >95% |
A variety of thiol nucleophiles was explored using these optimized conditions (Table 5); 1H NMR and SEC were used to measure conversion of each reaction, with the two methods showing agreement within ∼15% in most cases. Incomplete conversion was observed for 4-tert-butylthiophenol and 2-naphthalenethiol under ambient (Table 5, entries 1 and 3) or inert (N2) atmosphere (Table 5, entries 2 and 4), which we attribute to the steric bulk of these thiols. Nearly quantitative functionalization was observed for methyl thioglycolate (Table 5, entries 5 and 6) regardless of atmosphere; however, the SEC peak shape for the product showed distinct broadening and tailing compared to the initial p(TCND-(CO2Me)2), which may be due to possible deleterious side-reactions promoted by the relatively acidic α-protons of the α-mercapto ester (ESI† Section 9).49,50 By contrast, 1-methyl-3-mercaptoproprionate, a β-mercapto ester, underwent high conversion to give a polymer product with no significant change in peak shape (Table 5, entry 7). Finally, hexanethiol also gave high conversion under these conditions (Table 5, entry 8).
Overall, this study revealed that TCNs and TCNDs are a synthetically accessible class of ROMP monomers that do not compromise on activity or controlled polymerization behavior in comparison to other commonly used monomers (e.g. exo-NBs). In particular, TCND-(CO2Me)2 is a uniquely promising monomer due to its ease of synthesis (2 steps from commercial materials), fast ROMP kinetics, and its innate activity towards post-polymerization functionalization with thiols; we envision its future widespread use for functional polymer design. Finally, our experimental and computational study revealed that monomer steric parameters correlate most strongly with their measured propagation rates; these results offer the polymer chemist both the intuitive framework and computational tools for developing future novel ROMP monomers.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2326835, 2326836 and 2351614. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc01986e |
This journal is © The Royal Society of Chemistry 2024 |