Massimo Christian
D'Alterio
a,
Serena
Moccia
a,
Yolanda
Rusconi
b,
Claudio
De Rosa
a and
Giovanni
Talarico
*ab
aDepartment of Chemical Sciences, Università Degli Studi Di Napoli Federico II, Via Cintia, 80126 Napoli, Italy. E-mail: talarico@unina.it
bScuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy
First published on 9th August 2024
Stereoselective ring opening polymerization (ROP) of racemic lactide (rac-LA) is a challenging goal because a rationale connecting the catalyst structure and polymer microstructure (as has been established for α-olefin polymerization) is still missing. In this work, we reveal the origin of the stereoselective preference for D and L-lactide with two enantiopure salen–Al complexes, which have so far been claimed as the most efficient in enantiomorphic site control, using Density Functional Theory calculations. We introduce active site reorganization and monomer/chain switching throughout the reaction pathway, unconventional aspects necessitating careful consideration when confronting the intricacies associated with chiral catalyst recognition. We show how the catalytic pocket easily rearranges in the reaction path establishing a novel concept of the ligand coordination controlled by monomer binding. The resulting final picture of PLA stereoselectivity is much more complex than that of α-olefin polymerization catalysis, and a “complete” prediction by brute-force is (currently) hard, but the principles evolving should – even in their incomplete form – be useful in the design of new selective catalysts.
The lactide (LA) monomer has two chiral centers and, as a consequence, three diastereoisomers, namely D-LA (RR), L-LA (SS) and meso-LA (RS). The homopolymerization of a single enantiomer by ring opening polymerization (ROP) results in an isotactic polymer whereas the stereocontrolled ROP of rac- and meso-lactide utilizing metal complexes as initiators may lead to a wider range of polymer microstructures (Scheme 1).5,6 Analogous to α-olefin polymerization transition metal (TM) catalysis,7–10 the ROP stereochemistry can be dictated by either the enantiomorphic site control (ESC) mechanism, where catalyst chirality dictates homochiral preference, or the chain-end control (CEC) mechanism, where the preference for homochirality or heterochirality depends on the last inserted monomer unit.11 Indeed, the stereochemistry influence on the physical and mechanical properties of PLA materials has been carefully assessed, and manipulating the microstructure becomes pivotal for tailoring PLA properties with the aim to replace polyolefins in target applications.12 However, simple relationships between ligand catalysts/PLA microstructures are not yet achieved13–20 diversely from the α-olefin counterpart where the polymer microstructure can be finely tuned by suitable ligand catalyst modification a priori.21
In a seminal study of Spassky et al.,22 the enantiomerically pure binaphthyl–salen aluminum complex (system (R)-1, Scheme 2) was employed for the synthesis of a block copolymer poly(D-lactide)–poly(L-lactide) [(PDLA–PLLA)n] using rac-LA monomer feed. Indeed, system (R)-1 revealed, for the first time, a clear ESC preference toward D-LA rather than L-LA (kRR/kSS = 20, Scheme 2), with the formation of a stereocopolymer due to the incorporation of L-LA at a high conversion rate, when the reaction mixture is poor in D-LA. Relevant insights have been later reported by Coates et al. by using rac-1 for rac-LA polymerization still obtaining isotactic-enriched PLA with a high melting temperature (Tm = 179–191 °C) even at low conversion.23 The formation of a stereoblock microstructure at low conversion was explained by assuming a polymer chain exchange among the two enantiomeric forms of rac-1 catalysts.23 Interestingly, Feijen et al.24 demonstrated that the enantiopure aluminum complexes incorporating a Jacobsen ligand ((R,R)-2A, Scheme 2) also exhibited significant isoselectivity in rac-LA ROP by ESC but with a clear preference toward L-LA (kSS/kRR = 25, Scheme 2). As mentioned for system (R)-1, the progressive increase in the relative concentration of the D-LA isomer within the monomer pool, ascribed to the initial L-LA transformation into PLLA, results into a stereosequence gradient, transitioning from extended isotactic L-lactyl units to extended isotactic D-lactyl units. Finally, the enantioselectivity dramatically drops when the tBu moieties in the ortho position are substituted with H atoms (Scheme 2), as recently reported by Wang and coworkers (kSS/kRR = 5, Scheme 2 for ((R,R)-2B).25
Inspired by these results, the enantiomeric resolution and stereoselective polymerization of rac-LA have become a challenge in the last two decades with the final target to achieve high levels of isoselectivity combined with an enhanced productivity.26–30
However, factors governing ESC and CEC for ROP mechanisms, as well as the interplay between those two stereocontrol mechanisms remain obscure.31,32 In our recent communication,33 we anticipated that the origin of the stereoselectivity in rac-LA ROP achieved by the enantiopure system (R)-1 is different from the classical ES reported so far. Novel mechanistic steps have been used to trace the experimental stereoselectivity as summarized in Scheme 3. This (remarkable) complexity is due to several factors: 1) multiple ligand wrapping modes at the aluminum center; 2) variations of the wrapping mode configurations during the ROP catalytic cycle; 3) rate limiting step (RLS) energetics dependent on the monomer (D and L) LA as well as on the chirality of the last inserted unit in the growing chain.34 This complicated scenario makes the prediction of the stereoselective ROP of rac-LA a challenging game for chemists31 and justifies the lack of ligand structure/polymer microstructure relationships as consolidated in stereoselective TM catalyzed olefin polymerization.7,10,21 Going into more detail, we found that the ligand wrapping mode for (R)-1 may assume fac–mer 1 (fm1), fac–mer 2 (fm2) and fac–fac (ff) configurations and they may interconvert going from the nucleophilic addition (TS1) to the ring opening (TS2) steps.
We identified several reaction paths under the assumption that for the first mechanism (M1), both TS1 and TS2 occur with the same wrapping mode (fm1, fm2 or ff, M1-A, M1-B, M1-C, respectively). For the second mechanism (M2), TS1 occurs at fm1 and TS2 at fm2 and vice versa (M2-A and M2-B, respectively). The third mechanism, (M3), based on the interconversion among fac-mer and fac-fac, is also possible but it requires an additional TS that involves active site reorganization via the modification of the OL–Al–OL valence angle (defined as α) (see later).
The aim of this work is to furnish, by using computational methods rooted into Density Functional Theory (DFT), a comprehensive analysis of the ESC mechanism(s) for the stereoselective ROP of rac-LA promoted by (R)-1 and (R,R)-2A to understand the key factors leading to asymmetric kinetic resolution polymerization (AKRP) processes35 for D and L-LA. We will elucidate how the two enantiopure salen–Al complexes, despite the ligand/metal similarity, act inversely in the AKRP of rac-LA (Scheme 2). The structural features of the ligands, the mechanistic ROP intricacies and the role of the chirality of the growing polymer chain will be analyzed. We will also highlight the microstructure differences derived from our DFT calculations for the so-called stereoblock PLAs promoted by enantiopure (R)-1 and (R,R)-2A complexes. Finally, a novel concept for the stereoselective ROP of rac-LA will be defined, aiming to establish ligand catalyst/PLA microstructure relationships.
Fig. 1 Preferred Gibbs energy paths for RR-LA (right) and SS-LA (left) propagations promoted by (R)-1 following the M1–M2 mechanisms. The RLSs for the reaction paths (see text) are reported in bold. |
In Scheme 3, we draw the (easy) interconversion between fm1 and fm2 wrapping modes going from TS1 to TS2. In detail, this interconversion is possible through a barrierless rotation at INT1 around the dihedral angle indicated as θ1 and involves a single bond (Scheme 3). This mechanism is particularly noteworthy since it enables the catalysts to find the lowest energy paths among the various accessible on the potential energy surface, particularly when the ring opening step (TS2) shows remarkable differences in the Gibbs energies depending on the ligand wrapping mode (Fig. 1). The DFT geometries of these two TSs are reported in Fig. 2 and the main steric interactions leading to the TS2 energetic preference for fm2 (Fig. 2A) vs. fm1 (Fig. 2B) configurations are indicated with arrows.
The steric maps obtained with a molecular descriptor (% VBur)42–44 reported in Fig. 3 may help to visualize the (R)-1 ligand features leading to the energetic differences reported in Fig. 2. Indeed, the key aspect of the Spassky ligand is the binaphthyl linker (see the red zone in Fig. 3) that is closer to the monomer moiety in TS2 having an fm1 configuration (Fig. 3B) than to the more flexible polymer chain in fm2 (Fig. 3A). Indeed, at fm2, the couple RR-chain and RR-LA monomer finds a better accommodation rather than at fm1, where the monomer points toward the hindered area of the NE quadrant. On the other side, at fm2, the monomer lies in the free SW quadrant and the chain avoids clashing, positioning itself on the z-axis out of the plane. In the previous considerations, the chirality of the chain also plays an important role in the orientation assumed during the TSs. Indeed, in its conformational minimum, the RR-chain at the C2–O2–C3–C4 dihedral angle is G+ (≈+60°), whereas the SS-chain at the same dihedral angle is G− (≈−60°).
The interconversion among fac–mer and fac–fac is also possible but it requires the energetic evaluation of an additional TS that involves active site reorganization via the modification of the OL–Al–OL valence angle (defined as α) from ≈90 to ≈180° (Fig. 4). In Fig. 4A, we report the preferred Gibbs energy path for the mechanism moving from TS1 (in the preferred fm1) to TS2 (ff) as a function of α for (RR)-LA polymerization promoted by (R)-1 and we located the TS for such ligand interconversion at α = 141°, with a barrier of 14.6 kcal mol−1 very close to the RLS found for M2-A (14.4 kcal mol−1, see Fig. 1). These computations suggest that fm–ff transformation has to be considered in the stereoselectivity computation, increasing the number of mechanistic paths to be computed. However, we noted that such a feature is strictly related to the catalyst flexibility. In our computational screening looking for a simple method to evaluate the ligand wrapping mode energetics, we found that the use of a bidentate ligand such as acetylacetonate (acac−) to simulate the monomer and the growing chain binding appears promising. In fact, the computations with acac− (Fig. 4B) to mimic the polylactide growing chain avoided the complications derived from the several degrees of freedom of the chain (Fig. 4A), still showing similar energetic profile and geometrical features of TSα.
These findings achieved for (R)-1 need to be better discussed and inserted in a more general framework for a comparison with system (R,R)-2. Four different diastereoisomers are derived from the various wrapping modes of the salen ligand around the Al metal center (Scheme 4). These geometries arise from the combination of fac and mer with Λ and Δ isomerism in the octahedral geometries adopted by the ligand in the ROP TSs (Scheme 3). The mer–mer isomer was ruled out since the chain and the monomer would be in trans of each other, forming an angle of ≈180° (not feasible with the ROP mechanism). Complexes with ff or fm configurations are inherently chiral and the chirality of the bridge (indicated as X in Scheme 4) “forces” the catalyst to adopt a specific octahedral wrapping mode.
Scheme 4 Illustrations of the combination of fac and mer with Λ and Δ isomerism for octahedral chiral salen-type catalysts. |
It is worth to note that (R)-1, bearing X = (R)-binaphthyl as the bridge, presents two stable diastereoisomers, fm Δ and ff Λ which differ only by 1 kcal mol−1 (Table 1) due to the flexible binaphthyl linker, constituted by four carbon atoms.
Diastereoisomers | Catalyst | |
---|---|---|
(R)-1 | (R,R)-2A | |
fac–mer Δ | 0.0 | 0.0 |
fac–mer Λ | 29.2 | 8.5 |
fac–fac Δ | 40.2 | 35.9 |
fac–fac Λ | 1.0 | 13.7 |
This explains why the SS-LA monomer dictates the ff configuration of the catalyst (Fig. 1) with a Λ octahedral chirality that is preferred to the fm with a Δ octahedral chirality. Indeed, the SS-LA polymerization is more easily accommodated in the ff Λ wrapping mode (Fig. 5A) than in the fm Δ wrapping mode (Fig. 5B) with an energetic preference of 4.0 kcal mol−1 (Table S1†). In Fig. 5A, where the favored situation is depicted, the monomer is located in the southern region of the map and the chain lies parallel to the z-axis. Steric interactions between the ligand and one of the methyl groups of SS-LA are present, explaining the energetic difference in TS1 with RR-LA (16.2 vs. 12.7 kcal mol−1, Fig. 1). However, the situation is additionally disadvantaged with fm1 (Fig. 5B), by the fact that the chain is directed toward the hindered SE zone of the map.
It is worth to remark that the ff wrapping mode is precluded for (R,R)-2A having X = (R,R)-cyclohexyl bridge, which is far less flexible being constituted by two carbon atoms. Indeed, in this case, only the fac–mer Δ isomer is energetically feasible (Table 1). This important aspect that we will take into account in the following sections sheds a new light on the understanding of chiral catalyst recognition toward rac-LA.
Once a reasonable agreement with the experimental data is achieved, the origin of the stereoselectivity can be rationalized. In particular, the high steric hindrance exerted by bulky tBu groups in the ortho position of the (R,R)-2A system is crucial for inferring a stereoselective outcome to this reaction. This feature can be visualized by the % VBur steric maps reported in Fig. 7. Indeed, TS1 for SS propagation accommodates both the monomer and the growing chain in the least hindered area of the ligand (Fig. 7A). For RR propagation, the orientation of the RR-chain points towards the SW occupied by tBu groups (Fig. 7B), justifying the energetic difference between the two RLSs determining the enantioselectivity. As a matter of fact, the conformation of the chain with the steric hindrance exerted by the ligand pushes the disfavored axial positioning of the lactide methyl substituents of the nucleophilic addition TS (Fig. S2†). Lastly, the % VBur steric maps also explain the energetic reason for the switching of the monomer and growing chain positions in the ring opening step (TS2) with fm1 (Fig. 7C) having less unfavorable interactions with respect to fm2 conformations (Fig. 7D).
The importance of tBu groups in the ortho position has been further confirmed by computing the stereoselectivity of the (R,R)-2B system (replacing tBu with H atoms) in the rac-LA polymerization. The calculated enantioselectivity of (R,R)-2B decreases to 0.5 kcal mol−1 (RLS structures are reported in Fig. S3†), in good agreement with recent experimental data by Wang et al. who claimed a poor stereoselectivity in the ROP of rac-LA toward SS-LA promoted by (R,R)-2B.25 The % VBur steric maps of the RLS for SS-LA (Fig. 8A) and RR-LA (Fig. 8B) obtained with (R,R)-2B show that the occupancies of NW and SW quadrants drop respectively to 57.1% and 43.9% from 76.5% and 49.2% calculated for (R,R)-2A (Fig. 7A and B).
Fig. 8 % VBur steric maps of the RLS corresponding to TS1 for SS-LA (A) and RR-LA (B) propagation with catalyst (R,R)-2B. |
The DFT energetic plots in Fig. 1 and 6, combined with steric maps in Fig. 3 and 5, rationalized the origin of ESC control in the rac-LA ROP promoted by (R)-1 (RR preference) and by (R,R)-2A (SS preference) including the intermediate situation of (R,R)-2B. The multiple mechanistic paths we proposed are the only way to explain the experimental data (Table S4†). However, additional insights are disclosed by computing the activation Gibbs energies for all enchainments, including the RR and SS insertions on SS- and RR-chains, respectively. The main results summarized in Table 2 (complete lists reported in Tables S5 and S6†) suggest that the occasional insertion of SS-LA is corrected by the ESC of (R)-1 being the heterotactic enchainment intermediate between the isotactic ones (Table 2-left and Fig. S4†). Occasional RR insertion at (R,R)-2A are, instead, favoured by the last inserted RR-chain (Table 2-right and Fig. S5†). The additional CEC interplay we found might explain the isotactic stereosequences reported by Nomura et al.29 in the rac-LA ROP promoted by an achiral salen catalyst (X = Et; R2 = tBu, Scheme 4) as well as the subtle microstructure differences of PLA stereoblocks obtained by (R)-1 and (R,R)-2A.
System (R)-1 | System (R,R)-2A | ||||
---|---|---|---|---|---|
Assembly mode | ΔG‡RLS (kcal mol−1) | ΔΔG (kcal mol−1) | Assembly mode | ΔG‡RLS (kcal mol−1) | ΔΔG (kcal mol−1) |
RR-LA + RR chain | 14.4 | 0.0 | SS-LA + SS chain | 11.3 | 0.0 |
RR-LA + SS chain | 14.9 | 0.5 | RR-LA + RR chain | 13.5 | 2.2 |
SS-LA + RR chain | 15.2 | 0.8 | RR-LA + SS chain | 13.9 | 2.6 |
SS-LA + SS chain | 16.2 | 1.8 | SS-LA + RR chain | 14.5 | 3.2 |
In contrast, with a less flexible chiral bridge, such as (R,R)-cyclohexyl in (R,R)-2A, access to ff chiral octahedral wrapping modes is precluded. However, multiple monomer-controlled reaction paths remain. The stereoselective outcome, in this case, is ensured by a fixed chiral octahedral environment enforced by the steric hindrance of the tBu groups in the ortho position of the salen ligand, leading to the preferred formation of PLLA from rac-LA. This feature was further confirmed by modeling the ROP of rac-LA promoted by (R,R)-2B, which has H atoms in the ortho position, resulting in a lower AKRP.25
Finally, the DFT calculations of the overall microstructures revealed the interplay between ESC and CEC, opening the way to rationalize the CEC promoted by achiral salen systems.29 We believe that the insights from this work can be useful for both the theoretical and experimental communities in the rational development of ROP stereoselective catalysts,45 thereby bridging the gap with TM-catalyzed olefin polymerization.46
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cy00937a |
This journal is © The Royal Society of Chemistry 2024 |