Open Access Article
Malwina J. Niedźwiedź
,
Wojciech Ignaczak,
Peter Sobolewski
,
Agata Goszczyńska,
Gokhan Demirci* and
Miroslawa El Fray
*
Department of Polymer and Biomaterials Science, Faculty of Chemical Technology and Engineering, West Pomeranian University of Technology in Szczecin, Al. Piastów 45, 70-311 Szczecin, Poland. E-mail: mirfray@zut.edu.pl; gdemirci@zut.edu.pl
First published on 19th June 2023
Injectable and in situ photocurable biomaterials are receiving a lot of attention due to their ease of application via syringe or dedicated applicator and ability to be used in laparoscopic and robotic minimally invasive procedures. The aim of this work was to synthesize photocurable ester-urethane macromonomers using a heterometallic magnesium-titanium catalyst, magnesium-titanium(IV) butoxide for elastomeric polymer networks. The progress of the two-step synthesis of macromonomers was monitored using infrared spectroscopy. The chemical structure and molecular weight of the obtained macromonomers were characterized using nuclear magnetic resonance spectroscopy and gel permeation chromatography. The dynamic viscosity of the obtained macromonomers was evaluated by a rheometer. Next, the photocuring process was studied under both air and argon atmospheres. Both the thermal and dynamic mechanical thermal properties of the photocured soft and elastomeric networks were investigated. Finally, in vitro cytotoxicity screening of polymer networks based on ISO10993-5 revealed high cell viability (over 77%) regardless of curing atmosphere. Overall, our results indicate that this heterometallic magnesium-titanium butoxide catalyst can be an attractive alternative to commonly used homometallic catalysts for the synthesis of injectable and photocurable materials for medical applications.
In our previous work, we demonstrated that bismuth and zinc-based homometallic catalysts can be attractive, less toxic alternatives to DBTDL for the synthesis of photocurable ester-urethane macromonomers for elastomeric networks preparation. In that work, we discussed there two possible catalysis mechanisms: insertion and Lewis acid.17 Importantly, because isophorone diisocyanate (IPDI) is an asymmetric cycloaliphatic diisocyanate with primary and secondary isocyanate groups, the choice of catalyst plays an important role. Typically, the secondary cycloaliphatic group is more reactive than the primary aliphatic group.18 For the case of Lewis acid catalysts (such as DBTDL), the reactivity of the secondary isocyanate group is further increased. Meanwhile, for the case of Lewis base catalyst (such as DABCO), the primary isocyanate group becomes slightly more reactive, leading to similar reactivity of both NCO groups.
Recently there has been growing interest in heterometallic complexes which can operate via multiple interactions and can result in higher catalytic activity compared to their homometallic analogues.19 For example, heterometallic complexes of Ti/Zn,20 La/Mg,21 and Li/Mg22 have been shown to have superior catalytic activities to their homometallic analogues for the ring-opening polymerization. Most of these complexes feature a M–O–M′ framework, enabling intermetallic electronic communication and/or “ate”-type activation (vide supra), which can lead to enhanced nucleophilicity of the M–R bond (e.g. R = alkoxide). Thus far, experimental observations presented in the above, as well as many other works suggest that metals in multi-metallic catalytic systems have different functions: (1) larger and more electropositive metals (e.g. potassium, lithium, sodium) act as monomer coordination sites; (2) smaller and less electropositive metals (e.g. magnesium, zinc, aluminium, indium, yttrium, germanium, tin) act as the source of metal-alkoxide nucleophiles.19
In this context, a potentially promising heterometallic catalytic system is magnesium-titanium(IV) butoxide. This catalytic system has been already used for different type of reactions, including polycondensation23 and polyesterification24–26 of polyamides,23,25 poly(ester-ether)s,27,28 poly(butylene terephthalate) (PBT)-based,23,29–31 and poly(butylene-succinate) (PBS)-based copolymers.32 Unfortunately, the organo-bimetallic magnesium-titanium catalytic mechanism is not well described. However, it can be assumed that the isocyanate is activated by coordination to the Lewis acid, while the alcohol is activated through a hydrogen bond with a basic tert-butyl alcohol ligand. The alcohol is thus delivered through a six-membered transition state. The methalated carbamate then captures a proton from tertbutyl alcohol or from the added alcohol (Fig. 1).33
Building upon our previous work, here we investigated synthesis of the telechelic macromonomers based on polyester diol using magnesium-titanium(IV) butoxide as catalyst. We hypothesized that this heterometallic catalytic system would enhance selectivity toward the secondary isocyanate function of IPDI and would result in shorter reaction times than our previous work with bismuth and zinc-based homometallic catalysts,17 while at the same time obtaining one dominant macromonomer that would remain non-toxic. Towards this aim, we examined chemical structure and physio-chemical properties of the obtained macromonomers and we performed the cytotoxicity screening of fabricated elastomeric networks to assess their suitability for medical applications, especially for delivery of a precursor in minimally invasive medical procedures.
C
O groups and forming urethane bonds (bending N–H vibrations), as calculated from FT-IR spectra, reached approx. 5. Based on NCO titration, this value corresponds to approx. 50% of reacted groups. For the case of the synthesis with use of the 0.1 mol% of catalyst, the ratio stabilized at approx. 10, therefore the second step was started with this higher ratio. Then, 6 mg of PTZ (inhibitor), a second dose of catalyst (the same amount as in the first step), and 13.2 mL of HEMA were added to the reaction mixture. The reaction was considered finished when all of the isocyanate groups were converted, as determined using FT-IR (absence of the band at 2262 cm−1). The flask was then removed from the oil bath and allowed to cool down to room temperature. The product was purified from unreacted reagents by triple precipitation in four-fold excess of ice-cold methanol and any residual solvents were evaporated under reduced pressure.
Nuclear magnetic resonance spectroscopy (1H and 13C NMR) was performed to confirm the chemical structure of the obtained telechelic macromonomers. Spectra were recorded using Bruker DPX HD-400 MHz at 25 °C. CDCl3 was used as solvent and all shifts were determined with respect to tetramethylsilane (TMS). Additionally, the degree of acrylation34 was calculated based on proton NMR according to eqn (1):
| Degree of acrylation = δu,w × (4/δa) | (1) |
Dynamic viscosity of the macromonomers was assessed using a DV3TRV rotary cone-plate rheometer (Brookfield AMETEK, USA). The parameters for the measurement were as follows: measuring head in the cone-plate system with a diameter of ϕ = 40 mm, distance between cone and plate h = 1 mm, deformation of 30%, constant shear rate
(s−1) 0.200, and temperature 25 °C.
Differential scanning calorimetry analysis (DSC) was performed using Q2500 DSC (TA Instruments, New Castle, Delaware, USA) calorimeter to examine phase change behaviour of macromonomers. Samples were weighed (∼5 mg) into aluminium pans and hermetically sealed before the analysis. Samples were cooled down to −90 °C, held isothermally for 3 min, and then heated up to 120 °C. The heating/cooling rate was 5 °C min−1. The data were analysed using TRIOS software.
Gel permeation chromatography (GPC) was used to determine average molar (apparent) mass (
n and
w) and molar mass dispersity, Ð (
w/
n) of obtained macromonomers. The GPC Wyatt (Germany) instrument was composed of a guard column and four Perfect Separation Solutions (PSS) columns (50, 100, 1000, and 100
000 Å), Dn 2010 WGE Dr Bures differential refractometer (RI), and Wyatt MALLS DAWN HELEOS multi-angle light scattering (LS) detectors. 8 different linear polystyrene standards (molecular weight of 970, 1990, 5030, 10
680, 19
760, 34
800, 70
950, and 126
700 g mol−1 (Polymer Laboratories and Solutions LLC, USA)) were used for calibration of the equipment. The tetrahydrofuran was used as the mobile phase at a flow rate of 1 mL min−1 at 35 °C. PSS WinGPC Unity software was used for data collection and molecular mass calculations.
In order to assess gel fraction of photocured elastomeric networks, samples were refluxed for 6 hours in a Soxhlet extraction apparatus (Behr Labor-Technik, Germany) in EtOAc. The samples were weighed before extraction (Winitial) and after drying (under reduced pressure until a constant mass was achieved) (Wfinal). The gel fraction was calculated according to eqn (2):
| Gel fraction (%) = Wfinal/Winitial × 100 | (2) |
Swelling tests of photocured materials were performed according to EN ISO:175. Discs (1.6 cm in diameter) were placed in chlorobenzene for 24 h at 40 °C. After that time, samples in closed laboratory glass vessels were removed from the incubator and left for 20 min at room temperature to cool down (15–30 min). The discs were then washed with THF, dried with filter paper, and weighed. The weighing after removal from the solvent took no longer than 30 s. The swelling ratio was calculated according to eqn (3), where ms indicates for swollen mass and md is dry mass of the sample.
| Swelling (%) = (ms − md)/md × 100 | (3) |
Dynamic mechanical thermal analysis (DMTA) measurements of photocrosslinked films were performed using a Q800 DMTA instrument (TA Instruments). Samples were measured at a frequency of 1 Hz and amplitude of 1 μm. Samples were first cooled to −90 °C and then heated to 120 °C. The rate of heating was set to 3 °C per minute.
Finally, elastomeric polymer networks were screened for cytotoxicity according to ISO10993-5 (ref. 35) using murine fibroblast cell line L929, as described previously.17 Extracts were prepared in duplicate, by incubating samples of each photocured material (3 cm2 area, thickness 0.5 mm, cut into small pieces) in 1 mL of complete growth medium (Dulbecco's Modified Eagle Medium (DMEM), containing 10% Fetal Bovine Serum (FBS), 2 mM L-glutamine, 100 U mL−1 penicillin, and 100 μg mL−1 streptomycin) for 24 hours at 37 °C, in a 24-well plate. At the same time, 10 × 103 L929 cells (passages 10–25) were plated per well in a 96-well plate in complete growth medium and incubated for 24 hours in a cell culture incubator. After 24 hours, the medium was aspirated and replaced with 100 μL of extract, followed by an additional 24 hours of incubation (5 technical replicates). Cell viability was then assessed using light microscopy (Delta Optical IB-100, Mińsk Mazowiecki, Poland) and the resazurin viability assay.36 A multi-functional plate reader (Biotek Synergy HTX, Winooski, VT, USA) was used to measure the fluorescence signal (excitation 540 nm, emission 590 nm). The experiments were repeated with similar results. The cell viability was expressed as percent of sham-extract treated cells.
Therefore, the observed differences in kinetics relative to our previous work may be explained by differences in selectivity and activity of a catalyst. Further, the higher ratio (∼11) NCO/NH observed for the case of synthesis using 0.1 mol% of catalyst may be also related to the pot life of Mg–Ti butoxide catalyst: at this low concentration the reaction proceeds slowly and the catalytic activity of the catalyst may begin to decrease before the lower ratio can be reached. All of the FT-IR spectra for all reactions are available in the ESI, Fig. S2–S5.† The kinetics of all reactions, shown as ratio of absorbance bands (A2262) to (A1526) as a function of time, are presented in Fig. 3. Meanwhile, Table 1 presents the reaction times and yields.
![]() | ||
| Fig. 3 Kinetics of macromonomers synthesis performed with use of different concentrations of Mg–Ti butoxide catalyst. | ||
| Material | Catalyst concentration (per step) [mol%] | Reaction time (I step) [hours] | A2262/A1525 ratio after I step | Reaction time (II step) [hours] | Total reaction time [hours] | Reaction yield [%] |
|---|---|---|---|---|---|---|
| PrMg-Ti_0.1 | 0.1 | 11 | 10.43 | 47.5 | 58.5 | 53 |
| PrMg-Ti_0.25 | 0.25 | 4.5 | 6.53 | 23.5 | 28 | 67 |
| PrMg-Ti_0.5 | 0.5 | 3.5 | 5.21 | 5.5 | 9 | 57 |
| PrMg-Ti_1 | 1 | 0.7 | 4.91 | 3.3 | 4 | 63 |
As it can be seen from Fig. 3 and Table 1, as can be expected, higher concentrations of catalyst resulted in shorter total reaction time (58.5 h for 0.1 mol% and 4 h for 1 mol%) and lower ratio of (A2262)/(A1525) after the first step (10.43 for 0.1 mol% and 4.91 for 1 mol%). The reaction yields are similar, ranging from 53% to 67%, which is consistent with our previous work where we observed a similar range (62–70%).17
Importantly, compared to our previous work,17 here the total reaction time for the highest concentration of Mg–Ti butoxide catalyst was markedly shorter, as compared to homometallic bismuth and zinc catalysts (7 hours and 25 hours, respectively)—despite a lower catalyst concentration (1 mol% vs. 2 mol%). This confirms that a different catalytic mechanism was involved for In the previous work,17 the homometallic catalysts were acting via either insertion mechanism or as Lewis acids.17,38 Meanwhile, in the heterometallic system used here, we hypothesize that the catalyst acts similarly to titanium alkoxide catalyst via insertion of the isocyanate into the Ti–O bond of the alkoxy ligand, but that due to the higher stability of the heterometallic catalyst, the pot life of the catalyst is elongated.33,37,38
![]() | ||
| Fig. 4 FT-IR spectra of macromonomers obtained with use of various concentrations of Mg–Ti catalyst. | ||
The presence of absorbance bands at 3368 and 1525 cm−1 confirmed IPDI incorporation (corresponding to stretching vibrations of N–H in urethane groups). The band at 1643 cm−1 corresponding to stretching vibration of C
C, confirmed the presence of HEMA. Additionally, the bands at 2921 and 2853 cm−1 correspond to stretching vibrations of C–H in –CH2 and –CH3. The band 1732 cm−1 is characteristic for stretching of C
O in ester bonds of the polyester diol, urethane and methacrylic groups. The band at 1460 cm−1 is characteristic for R–CH2–CH3 and –C–C– in rings, while the band at 1369 cm−1 reflects rocking vibrations of C–H in fatty acid chains of Priplast 1838. At 1302 and 1049 cm−1, the presence of RCOOR′ vibration in esters was detected, while the bands at 1236 cm−1 and 1168 cm−1 reflect the –C–O– stretching vibration and –C–O–C(
O) stretching vibration, respectively.
1H NMR and 13C NMR spectroscopy were used as more in-depth methods of verifying the chemical structure of synthesized macromonomers. Overall, the analysis of NMR spectra confirmed the FT-IR findings. All of the obtained macromonomers had similar structure, consistent with our expectations based on prior work.17
As a representative example, Fig. 5 presents 1H and 13C NMR spectra of macromonomer obtained using 0.1 mol% Mg–Ti butoxide catalyst per step with peak assignments. The 1H NMR spectra for the other catalyst concentrations are available in the ESI (see Fig. S6†).
![]() | ||
| Fig. 5 1H (A) and 13C NMR (B) spectra of macromonomers synthesized with of 0.1 mol% Mg–Ti butoxide catalyst. | ||
The degree of acrylation (Fig. 6 and Table S1†) was calculated based on NMR spectra which influences the crosslinking density of the obtained elastomeric networks. It can be observed that increasing the concentration of catalyst led to a decrease in the degree of acrylation, with the exception of sample PrMg-Ti_0.5 (0.5 mol% of catalyst). Higher catalyst concentration increases the reaction rate to achieve the equilibrium faster to the product which results in lower number of accessible NCO groups for methacrylate (HEMA) attachment. Lower catalyst concentration may result in reduced IPDI addition to polyester diol due to which some free terminal OH groups are still in the molecule which also leads to lower HEMA attachment. These data are also correlated with gel fraction of films after photocuring, which will be discussed in Section 3.7.
Unfortunately, obtained spectra did not allow us to explicitly support the presence one predominant macromonomer. We expected to note a shift of carbons and protons corresponding to IPDI, that would confirm that the Mg–Ti butoxide catalyst first increased the reactivity of the secondary isocyanate group.
n and
w values were lower with increasing catalyst concentration. The higher catalyst concentrations shorten the total reaction time substantially, because the reaction proceeds simultaneously at a greater number of active centers of the catalyst. At the same time, this effect favours the formation of shorter macrochains. Overall, the dispersity indexes of the macromonomers are low (1.5–1.7), which is consistent with our hypothesis that the Mg–Ti butoxide catalyst would favour formation of one predominant macromonomer.
| Material | Catalyst concentration (mol%) | Mn [g mol−1] | Mw [g mol−1] | Ð [Mw/Mn] | Dynamic viscosity [Pa s] |
|---|---|---|---|---|---|
| Priplast 1838 | — | 3700 | 7800 | 2.12 | 18 ± 2 |
| PrMg-Ti_0.1 | 0.1 | 10 000 |
17 500 |
1.74 | 172 ± 12 |
| PrMg-Ti_0.25 | 0.25 | 8200 | 14 100 |
1.71 | 162 ± 28 |
| PrMg-Ti_0.5 | 0.5 | 8200 | 13 000 |
1.58 | 136 ± 15 |
| PrMg-Ti_1 | 1 | 7500 | 11 400 |
1.53 | 119 ± 38 |
As expected, the dynamic viscosity of the macromonomers reflects the same trend (Fig. 7 and Table 2), with a decrease from 172 Pa s for the material with the lowest concentration of the catalyst to 119 Pa s for material synthesized with the highest catalyst concentration.
![]() | ||
| Fig. 7 Dynamic viscosity for macromonomers synthesized with different concentrations of Mg–Ti butoxide catalyst. | ||
![]() | ||
| Fig. 8 DSC thermograms of macromonomers synthesized with different concentrations of Mg–Ti butoxide catalyst. | ||
In all of the spectra presented in Fig. 9, similar bands related to the relevant groups were present. The decrease in the band intensity at 1643 cm−1, corresponding to the stretching vibration of C
C in methacrylate group, is evidence for the reaction of the unsaturated functional groups during the photopolymerization. Additionally, due to the reaction of methacrylate functional groups, a new band appears at 1634 cm−1, corresponding to the resultant C–C bond. For the case of materials that were cured under air atmosphere, we note that the signal at 1643 cm−1 is more pronounced, which is likely due to oxygen inhibition.
To further assess the obtained photocured polymer networks, we performed swelling tests according to EN ISO:175 using chlorobenzene. The obtained results confirmed our hypothesis regarding the differences in the density of polymer networks. Further, in this case, the results (Fig. 11 and ESI Table S4†) clearly showed differences between the films cured under air and under argon, indicating an influence of oxygen inhibition on the curing process. The differences between the results arise from differences between the methods. In assessing gel fraction, the loss of material in solvent was measured, whereas in the swelling test the solvent uptake was determined.
δ) are presented in Fig. 12 (the numerical values are presented in Table S5†). Interestingly, the values of the storage modulus, E′ (Fig. 12a), representing elastic behavior in viscoelastic polymers, indicated that within the glassy state region (below −20 °C), there are some significant differences in polymers stiffness (values of E′ differ from ∼600 MPa to ∼1100 MPa), independent of the curing atmosphere and the catalyst concentration. After the phase transition, all of the polymer networks photocured under the same atmosphere exhibit comparable elasticity (more details in Table S5†). However, overall samples cured under argon were less elastic, as compared to those cured under air atmosphere (i.e. 0.2 versus 0.078 MPa for 0.1 mol% catalyst concentration). This effect is most likely due to the lower conversion of macromonomers for the case of samples photocured in air atmosphere, due the negative effect of oxygen inhibition. The resulting lower crosslinking density would result in more space for macrochains to arrange. Additionally, possible residual unreacted macroradicals, captured in the polymer network, can act as lubricants. Combined these two factors are likely responsible for the observed increase in elasticity. For the case of the sample prepared from macromonomers synthesized with 1 mol% of catalyst and photocured under argon atmosphere, fracture occurred during the measurement. This is likely due to its relatively low crosslinking density, related to a low degree of macromonomer acrylation. At high catalyst concentration (1 mol%), the isocyanate groups reacted with the hydroxyl termini very rapidly (stage I lasted less than 1 hour), facilitating the formation of shorter macrochains, which was confirmed by significant changes in the molecular weight and dynamic viscosity of this macromonomer. Additionally, at the same molar ratio of P1838 and IPDI maintained within the 1st step of the reaction, the obtained shorter chains will statistically have less NCO end-groups. This was likely to significantly influence the 2nd step of the reaction, resulting in lowest values of macromonomer acrylation degree.
Overall, the storage modulus values indicate (0.02–0.45 MPa) that the obtained elastomeric polymer networks can be suitable for soft tissue applications (the numerical values are presented in Table S5†). For reference, the values of storage modulus of human abdominal wall is in the range of 0.1–10 MPa, while the stiffness of the heart muscle at the end of diastole is between 0.2–0.3 MPa.39
In terms of the loss modulus (E′′) (a measure of viscous behavior in polymers), the results (Fig. 12b), showed only one relaxation maximum at ∼35 °C, corresponding to the amorphous phase of the polymer networks. Differences in relaxation temperatures were subtle (up to 5 °C), but a shift towards lower temperatures with increasing concentration of the catalyst used during the synthesis was noticed for both photocrosslinking atmospheres.
The changes in tan
δ, a dimensionless parameter that does not represent any physical quantity directly, but is a measure the ratio of energy dissipated to energy stored with periodically variable deformation, are presented in Fig. 12c. The analysis of damping properties (relation of E′′ to E′) revealed that samples photocured under argon atmosphere more effectively suppress mechanical vibrations, most likely due to higher cross-linking density and their overall viscoelastic properties. Analyzing the influence of catalyst concentration on the damping properties of polymer networks, that the highest values of tan
δ (at the reference temperature point of 37 °C) were noted for samples where 0.5 mol% catalyst was used during the macromonomer synthesis. This result is correlated with these samples having the highest values of gel fraction, lowest swelling degrees, and the macromonomers having the highest acrylation degrees.
![]() | ||
| Fig. 13 Cell viability L929 cell viability after 24 hours of exposure to extracts of materials obtained with use of various concentrations of Mg–Ti butoxide catalysts. | ||
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra02157b |
| This journal is © The Royal Society of Chemistry 2023 |