Wei
Huang
a,
Kangqiao
Wen
a,
Scott T.
Laughlin
*a and
Jorge
Escorihuela
*b
aDepartment of Chemistry, Stony Brook University, 100 Nicolls Road, Stony Brook, NY 11790, USA. E-mail: scott.laughlin@stonybrook.edu
bDepartamento de Química Orgánica, Universitat de València, Avda. Vicente Andrés Estellés s/n, Burjassot, 46100 Valencia, Spain. E-mail: jorge.escorihuela@uv.es
First published on 18th September 2024
Over the past two decades, click chemistry transformations have revolutionized chemical and biological sciences. Among the different strain-promoted cycloadditions, the inverse electron demand Diels–Alder reaction (IEDDA) has been established as a benchmark reaction. We have theoretically investigated the IEDDA reaction of endo-bicyclo[6.1.0]nonyne (endo-BCN) with 2H-pyran-2-one, 2H-thiopyran-2-one, 2H-pyran-2-thione and 2H-thiopyran-2-thione. These 2H-(thio)pyran-2-(thi)ones have displayed different reactivity towards endo-BCN. Density functional theory (DFT) calculations show, in agreement with experiments, that endo-BCN reacts significantly faster with 2H-thiopyran-2-one compared to other 2H-(thio)pyran-2-(thi)one derivatives because of the lower distortion energy. Experimentally determined second-order rate constants for the reaction of a 2H-pyran-2-thione with different strained derivatives, including a 1-methylcyclopropene derivative and several cycloalkynes (exo-BCN, (1R,8S)-bicyclo[6.1.0]non-4-yne-9,9-diyl)dimethanol, dibenzocycylooctyne and a light activatable silacycloheptyne, were used to validate the computational investigations and shed light on this reaction.
Among the different dienes used in IEDDA reactions, 1,2,4,5-tetrazines have been widely investigated as counterparts for reacting with alkynes/alkenes.13 Conversely, reactivity studies of 2H-pyran-2-ones and analogues, the focus of this report, are scarce. Kranjc and Kočevar reported in 2005 the Diels–Alder reaction of 2H-pyran-2-ones with alkynes under thermal conditions and at high pressures to afford aniline derivatives.14 In 2022, Xian and coworkers reported a click-and-release strategy based on the reaction of 2H-pyran-2-thiones and other analogues with a strained alkyne, (1R,8S,9S)-bicyclo[6.1.0]non-4-yn-9-ylmethanol (endo-BCN, 1).15 In that report, computational studies based on the nucleus-independent chemical shift (NICS) were used to predict the cycloaddition efficiency of endo-BCN with the 2H-(thio)pyran-2-(thi)one analogues. In a recent work, we synthesized substituted pyrones and pyranthiones and studied their biorthogonality, demonstrating their utility in protein labelling and subcellular delivery of hydrogen sulfide in living cells.16 A deeper understanding of these reactions can help in the design of 2H-pyran-2-thiones and heterocyclic analogues with enhanced reactivity for potential applications in chemical biology.
Previous studies have shown that dienes 2–5 can react with strained alkynes through an IEDDA reaction to form an unstable non-aromatic intermediate 6–9viaTS1 (Scheme 1).15 This intermediate then undergoes a retro Diels–Alder reaction through TS2 to provide the aromatic derivative 10 and release a molecule of CO2, COS or CS2 as a byproduct (Scheme 1). Understanding the trends in their reactivity regarding the 2H-(thio)pyran-2-(thi)one derivative will be key for the further development of this click reaction.
![]() | ||
Scheme 1 Schematic representation of the IEDDA reaction of exo-BCN 1 with 2H-pyran-2-one (2), 2H-thiopyran-2-one (3), 2H-pyran-2-thione (4) and 2H-thiopyran-2-thione (5). |
Herein, a systematic theoretical study of cycloadditions of endo-BCN with 2H-pyran-2-one (2), 2H-thiopyran-2-one (3), 2H-pyran-2-thione (4) and 2H-thiopyran-2-thione (5) was performed using the activation strain model. A combined experimental-computational study was conducted to evaluate reactivity of a 2H-pyran-2-thione analogue towards different strained dienophiles including 1-methylcyclopropene (MeCp) and cycloalkynes such as endo-BCN, ((1R,8S,9R)-bicyclo[6.1.0]non-4-yn-9-yl)methanol (exo-BCN), ((1R,8S)-bicyclo[6.1.0]non-4-yne-9,9-diyl)dimethanol (DMBCN), dibenzocycylooctyne (DBCO) and a light activated silacycloheptyne (Si7HT).
The formation of the aromatic derivative 10 can be rationalized in terms of the following general reaction mechanism, which can be divided into two main steps: [4 + 2] cycloaddition and retro Diels–Alder reaction with the elimination COS, CO2 or CS2 depending on the oxygen or sulphur-containing heterocycle. The Gibbs energy profile for the reaction of endo-BCN with the different dienes 2, 3, 4 and 5 is shown in Fig. 1. The first step of the mechanism involves the [4 + 2] cycloaddition and is exergonic. In line with similar Diels–Alder reactions, a pre-reaction complex with higher energy than both separated reactants was located.17 The interaction between the diene and endo-BCN (1) led to the formation of the molecular pre-reaction complex (RC). A decrease in the enthalpy of the reaction system of around 6 kcal mol−1 was computed of this transformation accompanied with substantial reduction in the entropy (Fig. S2, ESI†). Therefore, Gibbs energy (ΔG) of the formation of pre-reaction complex (RC) has a positive value, which excluded the possibility of pre-reaction complexes as stable intermediates. As inferred from Fig. 1, 2H-pyran-2-thione (4) has the lower activation Gibbs energy (21.2 kcal mol−1) of the four heterocyclic derivatives (2–5). For 2H-pyran-2-one (2), the computed activation Gibbs energy was higher (22.7 kcal mol−1), whereas for 2H-thiopyran-2-thione (5) the barrier was up to 25.3 kcal mol−1.
![]() | ||
Fig. 1 Gibbs energy profile (in kcal mol−1) for the cycloaddition of endo-BCN (1) with dienes 2, 3, 4 and 5. |
Finally, 2H-thiopyran-2-one (3) exhibited the highest barrier (26.5 kcal mol−1) for the [4 + 2] cycloaddition among the four derivatives. Fig. 2 shows the transition state structures and activation Gibbs free energies (in kcal mol−1) for the cycloaddition of BCN (1) with heterocyclic structures 2, 3, 4 and 5. Given the lack of symmetry in the heterocyclic dienes, all cycloadditions occurred via non-synchronous transition state structures with C⋯C bond forming distances within the range of 2.13–2.47 Å. The pattern of reactivity was 3 < 5 < 2 < 4. Carbonyl derivatives, i.e., 2H-pyran-2-one (2) and 2H-thiopyran-2-one (3), had activation Gibbs free energies of 22.7 and 26.5 kcal mol−1, respectively, whereas for thiocarbonyl compounds, 2H-pyran-2-thione (4) and 2H-thiopyran-2-thione (5), activation Gibbs free energies of 21.2 and 25.3 kcal mol−1, respectively, were computed. The analysis of transition state structures reveals that 2H-pyran-2-thione (4) has the longest forming bond distance (2.47 Å), and therefore the earliest early transition state leading to the lowest activation energy. The forming bond distances are almost identical for the oxygen and sulphur-containing heterocyclic structures 2, 3, 4 and 5. According to the experimental data, 2H-pyran-2-thione (4) reacts 3.7 times faster than 2H-pyran-2-one (2) in a PBS buffer (pH 7.4 with 10% DMSO).15 Computational data predict a higher reactivity for 2H-pyran-2-thione (4) compared to the other oxygen and sulphur-containing heterocycles.
![]() | ||
Fig. 2 Optimized transition state structures for the [4 + 2] cycloaddition. Activation Gibbs energies (ΔG‡) are given in kcal mol−1, and distances are shown in Angstroms. |
From non-aromatic intermediates 6–9, a retro Diels–Alder reaction affords the aromatic product 10 and releases a molecule of COS for 7 and 8, CO2 for 6 or CS2 for 9. Fig. 3 shows the transition state structures and activation Gibbs free energies (in kcal mol−1) for the retro Diels–Alder reaction to yield the final product 10. The lowest activation Gibbs energy for the retro Diels–Alder reaction was found for intermediate 8 (TS2-4) (coming from 2H-thiopyran-2-thione (4)), with a barrier of 13.4 kcal mol−1, which delivered the final product 10 with an exergonicity of −49.3 kcal mol−1. Thiopyran-derivatives 7 and 9 (coming from 2H-thiopyran-2-one (3) and 2H-thiopyran-2-thione (5)) had barriers of 18.9 and 18.8 kcal mol−1, respectively; whereas intermediate 6 (TS2-2) showed a Gibbs activation barrier of 17.5 kcal mol−1. After the second transition state (TS2), the formation of the aromatic product 10 is highly exergonic, acting as the thermodynamic driving force of this process. The high exergonicity observed from the retro Diels–Alder reaction is due to the recovery of the aromaticity from non-aromatic intermediates 6–9 to deliver the final aromatic product 10 and the release of a molecule of COS, CO2, or CS2. In general, the O⋯C and S⋯C distances for the breaking bonds were around 2.1 and 2.6 Å, respectively, whereas the C⋯C distance for the breaking bond in the case of pyran-derivatives 6 and 8 (TS2-2 and TS2-4) was found to be shorter than that of thiopyran-derivatives 7 and 9 (TS2-3 and TS2-5). Overall, the type of atom at the heterocycle becomes a key factor with a significant effect the thermodynamic reaction feasibility (ΔG°). In this regard, those reactions involving the breaking of a C⋯O bond (TS2-2 and TS2-4) are more exergonic than those in which a C⋯S is being broken (TS2-3 and TS2-5).
Comparing the two transient steps (TS1 and TS2), the [4 + 2] cycloaddition is the rate determining step (RDS) for the IEDDA reaction between endo-BCN and 2H-(thio)pyran-2-(thio)ones, as also observed for IEDDA reactions of 1,2,4,5-tetrazines.18 According to the computed activation free energies for both transition states, and the high exergonicity of the overall process, the cycloadduct (10) would be formed after the reaction with the release of CO2, COS or CS2. Finally, a comparison between M06-2X/6-311+G(d,p), M06-2X/Def2TZVPP and ωB97xD/Def2TZVPP methods was performed to validate the results and each method yielded very similar results in terms of transition state barriers (ΔG‡) (see Table S3, ESI†).
Next, the activation strain model was used to rationalize the trend in reactivity.19Fig. 4a shows the distortion/interaction–activation strain (D/I-AS) analysis of the transition state structures. According to this model, the optimized transition state structure is separated into two fragments (distorted BCN and distorted heterocyclic structure), and the energy difference between the distorted fragment and the corresponding optimized structure is defined as the distortion energy (ΔEdist). The difference between the activation energy and the total distortion energy is the interaction energy (ΔEint). This model has been successfully applied to rationalize the reactivity of cycloadditions, substitution and elimination reactions, and transition metal catalysed reactions.20 As inferred from the analysis, BCN distortion varies slightly (green arrows, 2.6 to 2.9 kcal mol−1), but the variation in distortion energy for each diene is more significant (blue arrows). The variation of geometrical parameters of distorted endo-BCN at the TS was very similar for the four structures (Table S4, ESI†). Small angle deviations (lower than 3°) from the endo-BCN reactant structure agree with the computed low BCN distortion energies. The high distortion energy for the dienes can be associated to the distorted and bended structure at the TS in comparison with the almost planar structure of the reactant. In this regard, 2H-pyran-2-thione (4) displayed the lowest distortion in line with the lowest activation energy. The substitution of oxygen by sulphur at the pyranthione ring (from 4 to 5) was reflected in a significant increase in the distortion energy of the heterocycle, leading to a higher activation energy. A similar effect was observed in the atomic substitution from 2 to 3.
Finally, a non-covalent interaction (NCI) analysis was performed on the transition states (Fig. 4b). Upon inspection of the NCI surfaces, it is clear for all compounds a similarly sized contour was computed, with an attractive interaction in the region where C–C bonds are being formed, but also with a strong repulsive interaction.
Fig. 5 shows the activation strain diagram (ASD) of the cycloaddition of endo-BCN with 2H-pyran-2-one (black line), 2H-thiopyran-2-one (red line), 2H-pyran-2-thione (blue line) and 2H-thiopyran-2-thione (green line). The distortion energies for reactions of 2H-pyran-2-one (black line) and 2H-pyran-2-thione (blue line) are quite similar, and significantly lower than that for 2H-thiopyran-2-one (red line), and 2H-thiopyran-2-thione (green line) along the complete reaction coordinate. This trend is also observed for the activation energy and the interaction energy, and the study suggests the reactivity is controlled by the distortion energy. An earlier transition state indicates that the transition state geometry is closer to that of the reactants, which is reflected in a lower distortion energy for 2H-pyran-2-thione (4). It is evident that the reactivity differences between 2H-pyran-2-one (2) and 2H-pyran-2-thione (4) with endo-BCN comes from the significant differences in the distortion energy between reactants along the reaction coordinate as both reactions show similar interaction energies. Furthermore, these type of reactions are inverse electron-demand Diels–Alder reactions, in which the HOMO of the dienophile interacts with the LUMO of the diene. A more effective orbital overlap interaction was computed between the HOMO of endo-BCN and the LUMO of 2H-pyran-2-thione (4) compared to the other set of compounds.
Given the faster kinetics observed for 2H-pyran-2-thione among the different oxygen and sulphur-containing heterocycles,15,16 we experimentally studied the reactivity of a 2H-pyran-2-thione with several popular dienophiles. For this study, we choose a diester substituted 2H-pyran-2-thione, PyrS2, the current best in class compound based on our previous studies with these reagents in protein labelling experiments.16 The dienophile reaction partners we explored were 1-methylcyclopropene (MeCp) and the cycloalkynes exo-BCN, DMBCN,21 DBCO22 and light activatable Si7HT.23 We conducted kinetics studies for reaction of each of these dienophiles with PyrS2 (see ESI† for further details). Reactivity tests towards exo-BCN, DMBCN, DBCO and MeCp and were carried out at 25 °C in 50% CH3CN/H2O; whereas for Si7HT the kinetic data was determined at 25 °C in CH3CN due to the limited stability of silacycloheptyne after UV illumination (Scheme 2). As expected, PyrS2 exhibited high reactivity with exo-BCN, resulting in a second-order rate constant of 0.36 M−1 s−1. Contrary to other related cycloadditions, endo-BCN and exo-BCN displayed a similar reactivity. Slightly slower kinetics was observed for DMBCN, having a second-order rate constant of 0.15 M−1 s−1. In the case of DBCO, no significant reaction was observed at room temperature after 24 h, whereas the methylcyclopropene derivative (MeCp) had a second-order rate constant of 8 × 10−4 M−1 s−1. Finally, a second-order rate constant of 1.03 M−1 s−1 was determined for Si7HT after light activation. In all cases, the cycloadducts were isolated after the release of COS.
![]() | ||
Scheme 2 Chemical structures of 2H-pyran-2-thione PyrS2, the PyrS2 reaction partners evaluated here, and their second-order rate constants for the reaction with 2H-pyran-2-thione PyrS2. |
To further support our experimental findings, computational investigations using DFT calculations were performed at the M06-2X/6-311+G(d,p) level of theory and including the implicit solvent model SMD. To this purpose, we determined the Gibbs energy barriers for the cycloaddition of 2H-pyran-2-thione PyrS2 with different strained counterparts. The transition state structures and the calculated activation energy (ΔE‡), activation Gibbs energy (ΔG‡), and Gibbs energy of the reaction (ΔGrxn) are displayed in Fig. 6a. It is worth mentioning that activation energy barriers (ΔE‡) are lower than Gibbs energy barriers (ΔG‡). Gibbs energy barriers take into account both enthalpy and entropy changes that occur during a chemical reaction. These differences reveal that entropic contributions on these reactions are significant. The TS structures are all asynchronous with C⋯C bond forming distances ranging from 2.24 to 2.34 Å. The lowest activation Gibbs energy was computed for Si7HT, in agreement with the faster experimentally determined kinetics, having an earlier transition state and a longer forming C⋯C bond distance than other strained counterparts. The computed barriers for the BCN derivatives exo-BCN and DMBCN were found to be 19.0 and 18.8 kcal mol−1, respectively. In contrast, to other related cycloadditions, endo-BCN and exo-BCN displayed a similar reactivity. In the case of DBCO, the high activation Gibbs energy barrier of 25.1 kcal mol−1 is in agreement with the lack of reactivity of this dibenzocycloocytne derivative, as observed for other cycloaddition reactions.12a In agreement with experimental results, the computed rate constants that predict: (1) DMBCN will react faster than endo-BCN and exo-BCN with 2H-pyran-2-thione PyrS2, (2) MeCp has a low reactivity, and that (3) Si7HT will react with PyrS2 around 3 times faster than endo-BCN and exo-BCN. To our surprise, the kinetics observed for DMBCN and PyrS2 were slightly slower than expected from our computational findings, which might imply other factors like solvent effects played some role in reaction kinetics.
The D/I-AS analysis of the transition state structures of 2H-pyran-2-thione PyrS2 with the different studied strained counterparts is shown in Fig. 6b. In this case, the D/I-AS analysis was performed at the transition state as a defined point for each reaction under study. It should be noted that performing the D/I-AS analysis at the state geometry can lead to incorrect interpretations if transition states occur at highly different bond forming distances.24 However, here the forming bond lengths are very similar in all cases, warranting the comparison at the transition state. As shown in Fig. 6b, the distortion energy for PyrS2 (blue arrow) decreased with the smaller ring sizes with relatively smaller energy for MeCp and Si-7HT (14.2 and 13.3 kcal mol−1, respectively). In contrast, biphenyl flanked cyclooctyne, DBCO (19.6 kcal mol−1), showed significantly increased distortion energy for PyrS2 compared with non-fused ring cyclooctynes. The relatively larger distortion energy of PyrS2 was possibly resulting from the large steric effect from DBCO. This was further confirmed by the distortion energy for dienophiles. Another biphenyl flanked dienophile, Si7HT showed much smaller distortion energy compared with DBCO owing to the lower steric interactions after ring size contraction. The non-phenyl flanked cyclooctyne dienophiles had even lower distortion energy, providing more support for the vital role of steric effects in dienophiles. Interestingly, MeCp showed the highest distortion energy (10.5 kcal mol−1) with PyrS2, implying that the interaction between alkynes and PyrS2 were much more favorable than that with alkenes. Additionally, the interaction energy was relevant to steric effects with more favorable interactions with cycloocytnes (exo-BCN, endo-BCN, and DMBCN) with non-fused rings compared to fused ring counterparts (DBCO and Si7HT). The interactions between alkynes were also favorable compared that with alkenes.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ob01263a |
This journal is © The Royal Society of Chemistry 2024 |