Open Access Article
N.
Broftová
a,
T.
Landovský
a,
H.
Dvořáková
b,
V.
Eigner
c,
M.
Krupička
a and
P.
Lhoták
*a
aDepartment of Organic Chemistry, University of Chemistry and Technology, Prague (UCTP), Technická 5, 166 28 Prague 6, Czech Republic. E-mail: lhotakp@vscht.cz
bLaboratory of NMR Spectroscopy, UCTP, Technická 5, 166 28 Prague 6, Czech Republic
cDepartment of Solid State Chemistry, UCTP, Technická 5, 166 28 Prague 6, Czech Republic
First published on 15th May 2023
A phenoxathiin-based macrocycle represents an inherently chiral building block, well accessible in two steps from the starting thiacalix[4]arene. The oxidized derivatives bearing one sulfoxide group and three sulfonyl groups were found to exhibit unexpected stereochemical preferences of the sulfoxide group during transformations. The sulfoxide moiety is always pointing out of the cavity (S
O out), while the opposite (S
O in) configuration was never obtained by direct oxidation. In order to achieve full oxidation to sulfone, the configuration of the sulfoxide group must first be changed by a photochemical inversion before the final oxidation occurs. The phenomenon of stereomutation of the sulfoxide group in the thiacalixarene series was studied using a combination of experimental (NMR and single crystal X-ray analysis) and theoretical (DFT) approaches.
Although the formation of spirodienone derivatives4 has long been known in the classical calixarene series and is used for the synthesis of some unusual derivatisation patterns,5 the application of this approach to thiacalixarenes led to the discovery of a completely new type of chemistry. As shown in Scheme 1, the starting thiacalix[4]arene 1 can be smoothly transformed into spirodienone 2
,6 which is under acidic conditions unexpectedly rearranged into phenoxathiin-based macrocycle 3 in a very high yield.7 Due to the presence of the condensed heterocyclic moiety, compound 3 represents a highly rigidified inherently chiral system with potential application in the design of chiral receptors. During our attempts at the shaping of this skeleton (formation of various conformations) and the oxidation of sulfur bridges, we found8 that sulfur bridges can be oxidised to yield sulfone 4, which was alkylated to form intermediate 5.
Recently, rigid calix[4]arenes9 with high internal strain have been recognized as surprisingly reactive systems enabling unusual transformations.10 Similarly, it was found that compound 5 can be used for further unexpected reactions. In the presence of alkoxides, it is subject to a nucleophilic attack11 leading to a hitherto unexplored type of macrocycle 6, or it undergoes an unprecedented cascade of nucleophilic cleavage/rearrangements (with SH−) to form linear oligomers 7 (Scheme 1).12
In this paper, we focused on some other interesting aspects of oxidised compound 5 and related systems concerning their unusual stereochemical preferences during the oxidation of sulfur bridges. We hope that a better understanding of this phenomenon can be used in further derivatization of these systems and contribute to the application of these compounds in supramolecular chemistry.
Based on the above results, we can therefore postulate the following unexpected assumption: for the oxidation of the sulfoxide group in compound 8, a change of the configuration to the opposite stereoisomer is first required, and only then can the oxidation occur by itself (Scheme 2).
The stereochemistry of the sulfoxide group in compound 8 was already assigned by single crystal X-ray measurement8 which revealed its “out of cavity” configuration (Scheme 2). The fact that the reaction stopped spontaneously at the 3 × SO2, 1 × SO stage (compound 8) indicates that the access of the oxidizing agent from inside the cavity is greatly disadvantaged. Even to such an extent that a change in configuration is required first for oxidation (to compound 4) to occur at all.
This is reminiscent of our earlier observation13 that the oxidation of sulfur to sulfoxide in thiacalix[4]arenes in the cone conformation always occurs on the side opposite to the lower rim of calixarene.14 In other words, the oxidizing agent always approaches from the side of the upper rim, which is probably sterically more accessible (Fig. 1). In this case, however, the selectivity is not so pronounced and after the oxidation of all bridges to sulfoxide, further oxidation to sulfone normally follows.
Since the non-alkylated derivatives 4 and 8 are difficult to handle (they form long tails in thin-layer chromatography), we decided to study the last oxidation step on alkylated derivatives, which are easily chromatographically separable. Compound 8 was alkylated using an excess of Cs2CO3 as a base and alkyl iodide (MeI, EtI, and PrI) as an alkylating agent in acetone. The corresponding peralkylated compounds 9a–c were isolated in good yields (Scheme 3).
![]() | ||
| Scheme 3 Stereomutation of the sulfoxide group in phenoxathiin-based thiacalixarenes and subsequent oxidation. | ||
In order to study the possible photochemical stereomutation15 of the sulfoxide moiety, we first chose an appropriate wavelength of light. The UV/vis spectrum of the starting sulfoxide 9a in CHCl3 (see the ESI Fig. 43†) showed two maxima at 295 nm (ε295 = 9.02 × 103 l mol−1 cm−1) and 328 nm (ε328 = 7.77 × 103 l mol−1 cm−1), and the maxima for 9b and 9c were almost identical (ESI Fig. 44 and 45†). The sample of 9a was irradiated with light (commercially available LED with a maximum at a wavelength of 311 nm, which was close to the measured maximum of the compound). The reaction was carried out at room temperature and quenched after one day of irradiation. However, we found that under these conditions, the molecule was degraded to some unidentified products and the expected thiacalixarene derivatives were not found in the reaction mixture.16
To avoid harsh conditions, we used light with a maximum at 365 nm (commercially available LED). After 24 h, derivative 10a was isolated from the reaction mixture in 40% yield, together with 31% of the unreacted starting material 9a (Scheme 3). A similar reaction of 9b gave product 10b with the opposite configuration in 58% yield (11% of 9b was isolated), while a propoxy derivative provided 10c in 66% yield (14% of 9c was isolated).
The structure of 10a was confirmed by HRMS ESI+ analysis, which showed a peak at m/z = 895.2283, which was in very good agreement with the value of the [M + Na]+ cation (895.2285) predicted for 10a. The number of peaks and the splitting patterns of the 1H NMR spectrum (CDCl3, 400 MHz, 298 K) of 10a are the same as those of the starting compound. However, some differences are clearly visible here, such as the shift of the singlet of the MeO group from 3.55 (9a) to 3.17 (10a) or the change in the chemical shift of the tert-Bu group from 0.73 (9a) to 0.97 (10a). In both cases, these are groups that come into close proximity to the S
O function facing the cavity of compound 10a.
The unequivocal evidence of the inversion of sulfoxide stereochemistry was obtained by single crystal X-ray analysis. Slow crystallization from a CHCl3–EtOH mixture provided 10a as a solvate with ethanol in the partial cone conformation (monoclinic system, space group P21/c). As shown in Fig. 2a, the sulfur atom of the phenoxathiin system is tilted into the cavity, and the sulfoxide group points toward the neighbouring But and MeO groups (see above for changes in the NMR spectrum). The inversion of configuration is best highlighted by a comparison with the starting compound 9a (Fig. 2b). At the same time, the sulfoxide group is responsible for trapping ethanol, which is bound through a hydrogen bond (S
O⋯H–O = 1.977 Å).
![]() | ||
| Fig. 2 X-ray structure of: (a) 10a – side view (EtOH not shown for better clarity); (b) the starting 9a – side view (ref. 8a); (c) 10b – side view; (d) 10c – side view, (sulfoxide group shown as balls). | ||
Compound 10b crystallized in a triclinic system, space group P
, adopting the partial cone conformation. The solvent molecule (CH3OH) is held by HB interactions with the sulfoxide group (S
O⋯H–O = 2.051 Å), resembling thus the solvate of 10a (Fig. 2c). The arrangement is further stabilised by intramolecular close contacts between the sulfonyl oxygen atoms and the ether oxygen atoms of the calixarene skeleton (3.034 Å, 2.861 Å and 3.024 Å). These distances are shorter than the corresponding sum of the van der Waals radii (ca.3.04 Å),17 indicating attractive interactions between the oxygen atoms which can be classified as chalcogen bonding.18 The same features can be also found within the above 10a·EtOH complex. Compared to Me or Et derivatives, the introduction of propyl groups had only a minimal effect on the basic geometrical parameters of the system; see the structure of 10c (Fig. 2d).
Having in hand the corresponding stereoisomers 9a–c and 10a–c possessing the out and in configuration of the sulfoxide group, we were interested in how their subsequent oxidation would proceed. The oxidation of 9a to tetrasulfone 11a carried out in the dark under the standard conditions (30% aq. H2O2/TFA/CHCl3, room temp) showed only 17% conversion of the starting compound after 28 h. Similar results were also obtained for ethyl and propyl derivatives 9b and 9c, where we observed 13% and 17% conversion (1H NMR), respectively (Scheme 3).
In sharp contrast, oxidations of derivatives 10a–c with the reversed configuration proceeded smoothly in all cases with 100% conversion, and the respective sulfones 11a–c were isolated by preparative TLC in 91%, 87%, and 89% yields, respectively. The above result confirms our assumptions that the oxidation of the “in” configuration is much easier than that of the “out” configuration.
In order to confirm these conclusions also for non-alkylated derivatives, compound 8 was irradiated in the same way (light 365 nm, CHCl3, rt) as was done for compounds 9a–c to achieve the stereomutation of the sulfoxide group.
In addition to the unreacted starting compound 8, the crude reaction mixture showed the formation of a small amount of sulfone 4, identified by a comparison of the NMR signals with an authentic sample, while the main product 12 remained unknown (Scheme 3). Since the chromatographic separation of non-alkylated derivatives is very difficult, both products were converted by a reaction with methyl iodide (MeI, K2CO3/acetone/reflux) to the corresponding methoxy derivatives and then separated by preparative TLC on silica gel. The formation of sulfone 4 was confirmed by its alkylated product 11a. Interestingly, the alkylation of 12 led to a mixture of two products 12a and 12b possessing the same molecular mass but different 1H NMR spectra, which indicated the formation of two different conformers. Surprisingly, the HRMS ESI+ analysis of the methoxy derivative 12a showed a molecular peak at m/z = 879.2335, which was in perfect agreement with the [M + Na]+ cation (879.2336) predicted for a compound lacking one oxygen atom. Thus, the main product is not sulfoxide with altered stereochemistry, but compounds 12a and 12b possessing a sulfide bridge.
While methylation of tetrasulfone 4 or trisulfone monosulfoxide 8 leads to a single conformation (partial cone), the formation of two conformers 12a and 12b suggests that the conformational preferences of the sulfur-bridged derivative 12 are significantly altered. To reveal the possible mobility of alkylated derivatives 12a and 12b, the VT 1H NMR studies of both products were carried out. While heating of 12a (500 MHz, CDCl2–CDCl2) did not lead to any change of the spectrum (see the ESI†), compound 12b assigned as the 1,2-alternate was gradually transformed into 12a at 110 °C (see Fig. 4). This suggests that (i) the sulfur bridge in 12 presents less steric hindrance than sulfone 4 or sulfoxide 8, where a similar change is not possible; (ii) the partial cone conformer 12a represents a thermodynamically more stable conformation than 12b. This was further confirmed by heating a solution of compound 12b in CDCl2–CDCl2 overnight at 120 °C, which led to complete conversion to compound 12a.
The final proof of the structure was obtained by the single crystal X-ray analysis of 12a. This revealed that compound 12a crystallized as a solvate with cyclohexane (12a·C6H12) in a monoclinic system, space group P21/c. As shown in Fig. 3a, the molecule adopts the partial cone conformation and the missing oxygen atom was removed in the phenoxathiin moiety, thus forming a sulfide bridge. The conformation was stabilised by an array of intramolecular close contacts between the oxygen atoms from SO2 and OR moieties (2.890 Å, 2.883 Å, 2.869 Å and 3.039 Å), which can be classified as chalcogen bonds (Fig. 3b).
The crystal packing of 12a (Fig. 3c) was characterised by the infinite chains of calixarenes held together by the intermolecular chalcogen bonds between the sulfide S atoms and one of the S
O bonds of the opposite sulfone bridge. An S⋯O
S distance of 3.186 Å was reasonably shorter than the corresponding sum of vdW radii (3.32 Å).17 The individual chains are interconnected (see Fig. 3d) by non-classical hydrogen bonds between S
O groups and hydrogens in tert-butyl groups (S
O⋯H = 2.658 Å) and on the aromatic subunits (S
O⋯Harom = 2.666 Å).
The formation of both reduced and oxidised compounds during the attempted stereomutation of 8 was rather unexpected. To avoid possible confusion about the oxygen source for oxidation, the reaction was performed under argon in carefully degassed chloroform, but the result was the same. It indicates that the sulfoxide S
O bond in 8 was cleaved during the reaction, yielding 12. Such reactions (deoxygenations) are known from the chemistry of aromatic sulfoxides,19 and the resulting highly reactive oxygen species20 usually react with the environment (solvent) to form oxidized products. In our case, the only oxidizable molecule in the system (besides chloroform) is the starting sulfoxide 8, so the formation of a small amount of by-product 4 seems to be a more or less logical consequence.
To further support this idea, we repeated the reaction under identical conditions but with the addition of diphenyl sulfide as the presumed oxygen acceptor. After 11 hours of irradiation (364 nm) and subsequent methylation of the crude reaction mixture with MeI, the 1H NMR analysis showed three calixarene compounds 12a, 12b and 9a in 14%, 7% and 79% yields. The isomers 12a and 12b were obtained after column chromatography in 12% and 5% yields.
In any case, compound 12 represents the only example of its kind, where there is an unoxidized sulfur bridge together with the sulfone bridges. To the best of our knowledge, a similar type of compound has not yet been described for thiacalixarenes, because the oxidation of sulfur atoms always proceeds in such a way that sulfones begin to form only after all the sulfur bridges have been oxidized to sulfoxides.
Important parts of the reaction pathway were modelled using DFT. The detailed sequence of the multiple oxidation steps is beyond the scope of this research, hence we assumed that the sulphur bridges are fully oxidized before phenoxathiine oxidation takes place. Conformational screening of all structures was performed using CREST, and 30 lowest-lying conformers were optimized at the B3LYP/def2-SVP level with D3BJ corrections. The energy of the most favourable conformer is reported. The reactivity is expected to be guided by steric repulsion, which is usually well described at the DFT level, but due to a wide range of functional groups present in the molecule, final DLPNO-CCSD(T)/def2-QZVPP single point energy evaluation was performed. The relative energies at the DFT level differ quite substantially from CCSD(T) and both values are reported, with the DFT value (in italics).
From the thermodynamic point of view, the structures of 9a and 10a, possessing the out and in configuration of the sulfoxide group are close in energy, with the isomer 9a (out) being preferred by 3.3 (0.5) kcal mol−1. The reaction energy of the second oxidation step from sulfoxide 9a to sulfone 11a is −68 (−51) kcal mol−1, significantly more favoured than the first step, sulfide 12a to sulfoxide (9a or 10a), which is −41 (−28) kcal mol−1.
The kinetics of the oxidation steps was examined with respect to the approach of the oxidant, modelled using trifluoroperoxoacetic acid (TFPA). Both transition states leading to 9a and 10a were found, and the attack from outside was preferred by about 2.1 (−1.7) kcal mol−1. The second oxidation step, from sulfoxide to sulfone, proceeds from above the molecule, rather across the calixarene ring. The activation energy of the oxidation of 9a is 1.1 (−0.3) kcal mol−1 higher than that of the 10a isomer. The whole reaction sequence is shown in Fig. 5.
The difference in activation energy for the isomers agrees with the experimental results of the first oxidation step. In the second step, the calculated energy difference favors the oxidation of 10a, with S
O pointing inwards, which is impossible without the isomerization of sulfoxide. The intuitive explanation using the argument of the steric hindrance of the calixarene ring is probably justified in this case. In addition, the extent of HOMO delocalization to aromatic rings is smaller in 10a, making the lone pair on sulphur more reactive in the S
O in isomer.
The insufficiency of a computational model may come into play. Nevertheless, steric interactions are usually well described at the applied DFT-D3 level of theory, providing good geometries of intermediates and transition states. The role of solvent is not expected to be very important, namely, when comparing two very similar stereoisomers. The steric effects are presumably not the only contribution, for example in Fig. 5, phenoxathiin-oxide is deformed inwards during the oxidation of 10a, contributing to the increased activation energy.
O out), while the opposite (S
O in) configuration has never been obtained by direct oxidation. In order to achieve full oxidation, the configuration of the sulfoxide group must first be changed by a photochemical reaction before the final oxidation step occurs. The phenomenon of stereomutation of the sulfoxide group in this thiacalix[4]arene series was studied using a combination of experimental (NMR and single crystal X-ray analysis) and theoretical (DFT) approaches.
Data for 4 were in accordance with the previously published data.8
Data for 8 were in accordance with the previously published data.8 Compounds 9a–c were prepared according to the published procedure.8
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 111.2 mg, 40%).
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 117.3 mg, 58%).
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 73.7 mg, 66%).
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 23.2 mg, 91%).
Data for 11a were in accordance with the previously published data.8
Compound 11a was obtained in 17% yield under identical reaction conditions using derivative 9a as the starting compound.
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 22.6 mg, 89%).
Data for 11b were in accordance with the previously published data.8
Compound 11b was obtained in 13% yield under identical reaction conditions using derivative 9b as the starting compound.
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 22.1 mg, 87%).
Data for 11c were in accordance with the previously published data.8
Compound 11c was obtained in 17% yield under identical reaction conditions using derivative 9c as the starting compound.
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 4.3 mg, 17%).
Data for 11a were in accordance with the previously published data.8
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 3.3 mg, 13%).
Data for 11b were in accordance with the previously published data.8
:
cyclohexane 1
:
10 v/v) as a racemic mixture (white solid, 4.3 mg, 17%).
Data for 11c were in accordance with the previously published data.8
The residue was suspended in dry acetonitrile (10 mL) and caesium carbonate (800 mg, 2.45 mmol) was added to the mixture. The suspension was stirred for 30 minutes. After that, methyl iodide (0.18 mL, 2.9 mmol) was added and the reaction mixture was heated under reflux overnight. Then, the solvent was removed under reduced pressure. The residue was dissolved in chloroform (10 mL) and water (10 mL). The layers were separated and the aqueous phase was extracted with chloroform (3 × 20 mL). Combined organic layers were dried over MgSO4. The solvent was evaporated and the residue was purified by silica gel chromatography (eluent ethyl acetate
:
cyclohexane 1
:
6 v/v) to give compound 12a (12 mg, 12%) and compound 12b (5 mg, 5%) as white solids together with 55 mg of compound 11a (from unreacted 8).
, a = 13.7722 (5) Å, b = 14.5337 (5) Å, c = 25.6391 (9) Å, α = 79.646 (2)°, β = 88.001 (2)°, γ = 78.359 (2)°, Z = 4, V = 4944.5 (3) Å3, Dc = 1.280 g.cm−3, μ(Cu-Kα) = 2.25 mm−1, crystal dimensions of 0.42 × 0.27 × 0.07 mm. Data were collected at 180 (2) K on a Bruker D8 Venture Photon CMOS diffractometer with Incoatec microfocus sealed tube Cu-Kα radiation. Data reduction, scaling and absorption correction were performed in Apex4.25 The structure was solved by charge flipping methods22 and anisotropically refined by full matrix least squares on F squared using CRYSTALS23 to final value R = 0.075 and wR = 0.191 using 18
737 independent reflections (θmax = 70.4°), 1232 parameters and 66 restraints. The hydrogen atoms bonded to carbon atoms were placed in calculated positions refined with riding constraints, while hydrogen atoms bonded to oxygen were refined using soft restraints, after which they were refined using riding constraints. The disordered functional group positions were found in different electron density maps and refined with restrained geometry. MCE24 was used for the visualization of electron density maps. The occupancy of disordered functional groups was constrained to full. The structure has been deposited into the Cambridge Structural Database under number CCDC 2227441.†
, a = 12.4947 (5) Å, b = 12.7653 (6) Å, c = 19.6418 (8) Å, α = 92.4959 (17)°, β = 90.0704 (16)°, γ = 112.6404 (15)°, Z = 2, V = 2888.1 (2) Å3, Dc = 1.311 g cm−3, μ(Mo-Kα) = 0.38 mm−1, crystal dimensions of 0.53 × 0.34 × 0.23 mm. Data were collected at 180 (2) K on a Bruker D8 Venture Photon CMOS diffractometer with Incoatec microfocus sealed tube Mo-Kα radiation. Data reduction, scaling and absorption correction were performed in Apex4.25 The structure was solved by charge flipping methods22 and anisotropically refined by full matrix least squares on F squared using CRYSTALS23 to the final value R = 0.057 and wR = 0.169 using 12
169 independent reflections (θmax = 26.7°), 736 parameters and 89 restrains. The hydrogen atoms bonded to carbon atoms were placed in calculated positions refined with riding constraints, while hydrogen atoms bonded to oxygen were refined using soft restraints after which they were refined using riding constraints. The disordered functional group positions were found in different electron density maps and refined with restrained geometry. MCE24 was used for the visualization of electron density maps. The occupancy of disordered functional groups was constrained to full. The structure has been deposited into the Cambridge Structural Database under number CCDC 2227442.†
, a = 14.1073 (2) Å, b = 17.6934 (5) Å, c = 19.5984 (4) Å, α = 70.921 (2)°, β = 89.7620 (14)°, γ = 83.5928 (17)°, Z = 4, V = 4591.58 (3) Å3, Dc = 1.286 g cm−3, μ(Cu-Kα) = 2.39 mm−1, crystal dimensions of 0.42 × 0.29 × 0.07 mm. Data were collected at 95 (2) K on a Rigaku OD Supernova using an Atlas S2 CCD detector with Cu-Kα radiation from a microfocus sealed X-ray tube. Data reduction, scaling and absorption correction were performed in CrysAlis PRO.26 The structure was solved by charge flipping methods22 and anisotropically refined by full matrix least squares on F squared using CRYSTALS23 to the final value R = 0.051 and wR = 0.136 using 18
212 independent reflections (θmax = 74.6°), 1167 parameters and 116 restraints. The hydrogen atoms were placed in calculated positions refined with riding constraints. The disordered functional group positions were found in different electron density maps and refined with restrained geometry. MCE24 was used for the visualization of electron density maps. The occupancy of disordered functional groups was constrained to full. Despite significant efforts, a suitable solvent disorder model could not be found; therefore, the solvent contribution to the structure factors was accounted for using PLATON squeeze,27 resulting in a total solvent void volume of 374.6 Å3. The structure has been deposited into the Cambridge Structural Database under number CCDC 2227443.†
265 independent reflections (θmax = 26.4°), 710 parameters and 255 restraints. The hydrogen atoms were placed in calculated positions refined with riding constraints. The disordered functional group positions were found in different electron density maps and refined with restrained geometry. MCE24 was used for the visualization of electron density maps. The occupancy of disordered functional groups was constrained to full. The structure has been deposited into Cambridge Structural Database under number CCDC 2227444.†
Final single point energies were obtained at DLPNO-CCSD(T) with a def2-QZVPP31 basis set, def2/J auxiliary basis set33 and def-QZVPP/C auxiliary basis set for CC,35 using the RIJCOSX approximation.32
Geometries of the most favourable isomers of all stationary points at the DFT level are reported in the ESI.†
Footnote |
| † Electronic supplementary information (ESI) available: Spectral characterization of new compounds, temperature-dependent NMR spectra of 12a and 12b, and DFT calculations. CCDC 2227440–2227444. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3ob00530e |
| This journal is © The Royal Society of Chemistry 2023 |