Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

A computational triage approach to the synthesis of novel difluorocyclopentenes and fluorinated cycloheptadienes using thermal rearrangements

David Orr a, Jonathan M. Percy *a and Zoë A. Harrison b
aWestCHEM Department of Pure and Applied Chemistry, University of Strathclyde, Thomas Graham Building, 295 Cathedral Street, Glasgow, G1 1XL, UK. E-mail: jonathan.percy@strath.ac.uk; Fax: +44 (0)141 548 4898
bRefractory Respiratory Inflammation DPU, GlaxoSmithKline Medicines Research Centre, Gunnels Wood Road, Stevenage, SG1 2NY, UK

Received 22nd March 2016 , Accepted 15th June 2016

First published on 16th June 2016


Abstract

Electronic structure calculations have been used for the effective triage of substituent effects on difluorinated vinylcyclopropane precursors and their ability to undergo vinyl cyclopropane rearrangements (VCPR). Groups which effectively stabilised radicals, specifically heteroarenes, were found to result in the lowest energy barriers. Ten novel precursors were synthesised to test the accuracy of computational predictions; the most reactive species which contained heteroarenes underwent thermal rearrangements at room temperature to afford novel difluorocyclopentenes and fluorinated benzocycloheptadienes through competing VCPR and [3,3]-rearrangement pathways, respectively. More controlled rearrangement of ethyl 3-(1′(2′2′-difluoro-3′-benzo[d][1,3]dioxol-5-yl)cyclopropyl)propenoate (22) allowed these competing pathways to be monitored at the same time and activation energies for both reactions were determined; Ea(VCPR) = (23.4 ± 0.2) kcal mol−1 and Ea([3,3]) = (24.9 ± 0.3) kcal mol−1. Comparing our calculated activation energies with these parameters showed that no single method stood out as the most accurate for predicting barrier heights; (U)M05-2X/6-31+G* methodology remained the best for VCPR but M06-2X/6-31G* was better for the [3,3]-rearrangement. The consistency observed with (U)B3LYP/6-31G* calculations meant that it came closest to a universal method for dealing with these systems. The developed computational design model correctly predicted the observed selectivity of rearrangement pathways for both our system and literature compounds.


Introduction

The continual rise in the processing power of modern computer systems has allowed the computational chemist to perform ever larger and more accurate calculations over shorter and shorter periods of time. Calculations are routinely performed for synthetically interesting reactions in order to characterise pathways in detail and rationalise experimental observations.1 Typically, these calculations are performed after considerable synthetic optimisation, which can be time consuming and expensive. Computational evaluation of reactions prior to the commitment of experimental resource is now becoming less rare having been shown to streamline investigations into a range of organic transformations.2

We found that electronic structure calculations were extremely helpful in developing our understanding of the rearrangement of 1a to difluorocyclopentene 2 (Scheme 1).3


image file: c6sc01289b-s1.tif
Scheme 1 Correlation between experimental and calculated rates for ordering the rearrangement pathways of 1aG values (blue) calculated using UB3LYP/6-31G*, gas phase, 298 K, Gaussian'09, all energy values are in kcal mol−1).

We observed good agreement between calculated (29.8 kcal mol−1, UM05-2X/6-31G*, with the conformationally simpler Me ester) and experimental (28.6 kcal mol−1) VCPR activation energies, but only the less expensive UB3LYP/6-31G* calculations predicted the order of reactivity between competing cyclopropane stereoisomerisation and [3,3]-sigmatropic rearrangement correctly. The system is most unusual in that reactions involving open shell singlets and triplet species are concurrent with a more classically concerted pericyclic rearrangement to 3. Calculations also revealed the unexpected role that the alkene stereochemistry played in controlling rearrangement pathways; E-alkenoates underwent VCPR whereas Z-alkenoates reacted via initial [3,3]-rearrangement.

The two ring expansion rearrangements observed in our system have both been used successfully by other groups in total synthesis projects. VCPRs were utilised as key steps in securing the 5-membered cyclic cores for a range of complex structures,4 including α-vetispirene 4 from the sesquiterpene family of compounds,5 as well as the prostaglandin E2 methyl ester 5 (Fig. 1a).6 There have also been reports of related [3,3]-sigmatropic rearrangements,7 including divinylcyclopropane to cycloheptadiene rearrangements,8 being used in a similar fashion to generate molecular complexity. However, despite being utilised in the successful synthesis of marine sponge derived natural product frondosin B 69 and the highly potent SIRT-inhibitor 7,10 the [3,3]-rearrangements from our system are much less common (Fig. 1b).


image file: c6sc01289b-f1.tif
Fig. 1 Natural products and active pharmaceutical compounds synthesised using (a) VCPR and (b) aromatic-vinylcyclopropane rearrangements (key ring structures formed during the rearrangements are highlighted in blue).

The key rearrangements in the synthesis of 4 and 5 required high and very high temperatures of 190 °C and 560 °C, respectively. Substrate modifications11 and the utilisation of transition metal catalysts12 have started to allow access to room temperature rearrangements. Fluorine atom substitution on the cyclopropane ring has also been shown to have a beneficial effect.13 The ring strain energy is ∼10 kcal mol−1 higher for gem-difluorinated VCP 8a than that for the parent hydrocarbon 8b (Scheme 2).14 This induces regiospecific ring opening in 8a, due to the weakened distal C3–C5 bond, and lowers the activation energy for VCPR (ΔΔG = 10.2 kcal mol−1 between 8a and 8b).13b,15


image file: c6sc01289b-s2.tif
Scheme 2 Effect of fluorine atom substitution on VCPR activation energies (ΔG experimental values, 298 K, re-calculated from activation parameters, kcal mol−1).16

Our VCPR of precursor 1a exploited this effect, alongside stabilisation of the open shell singlet transition state through benzylic resonance, to lower the activation energy further (ΔΔG = 20.6 kcal mol−1 from 1a to 8b). This system offers a highly atom-efficient route into a novel difluorocyclopentene at practical rearrangement temperatures (100 °C).3 Despite recent synthetic advances for securing difluorinated 5-membered ring structures,17 difluorinated cyclopentenes are still rare motifs in the literature.18 To our knowledge, there are currently no reported routes for the synthesis of structures similar to fluorinated benzocycloheptadiene 3. Since both of these desirable fluorinated compounds were accessible from difluoro-VCP 1a, our system represents an exciting building block class for accessing new fluorinated chemical space.19

We now look to understand more fully what effect changing the functional groups around our difluorinated precursors has on the balance between these two rearrangement pathways; we wish to be able to design precursors which rearrange at relatively low temperatures. By securing accurate transition structures for all three possible rearrangement pathways (Fig. 2), we have ensured a relatively unique opportunity to start our study by assessing the scope and limitations of our system using electronic structure calculations before committing to any synthetic chemistry.


image file: c6sc01289b-f2.tif
Fig. 2 Geometries found representing the open shell singlet-VCPR, closed-shell [3,3]-sigmatropic rearrangement, and triplet cyclopropane stereoisomerisation transition states derived from VCP-1a/1b.

The computational screening would start by separately predicting the effect different substituents, attached either to the difluorocyclopropane or the alkene portions of the precursor, would have on the barrier for VCPR. A selection of interesting compounds based on the predicted ease of rearrangement would then be synthesised to assess the accuracy of the theoretical predictions. As the overall goal is to identify computational models that can be easily implemented by synthetic chemists, off-the-peg and less computationally-intensive methods would be preferred over calculations which required bespoke methods and were more demanding in terms of computational resource.20

Results and discussion

Effect of difluorocyclopropane substitution

We reported that the most cost effective method for calculating VCPR activation energies was an unrestricted B3LYP21 (UB3LYP) method with the 6-31G* basis set.3 This method was used to assess the impact of modification of the left-hand side of difluoro-VCP 9, with calculated activation barriers ranked alongside phenyl-substituted precursor 9j (Table 1 and Fig. 3).
Table 1 Effects of substitution on difluorocyclopropane precursors 9 (UB3LYP/6-31G*, gas phase, 298 K, Spartan'10)

image file: c6sc01289b-u1.tif

VCP R1 R2 ΔGB3LYP ΔΔGB3LYPa
a ΔΔGB3LYP = (ΔGB3LYP9x) − (ΔGB3LYP9j).
9a 2-Pyrrolyl H 16.8 −8.5
9b 2-Furyl H 21.0 −4.3
9c 4-Pyridyl-N-oxide H 21.1 −4.2
9d 2-Thiophenyl H 21.3 −3.8
9e 2-N-Boc-pyrroyl H 21.4 −3.9
9f 5-Benzo[d][1,3]dioxolyl H 23.0 −2.3
9g H2C[double bond, length as m-dash]CH2 (vinyl) H 24.5 −0.8
9h 3-Thiophenyl H 24.6 −0.7
9i 2-Thiazolyl H 25.3 0.0
9j Ph H 25.3 0.0
9k Ph Ph 25.8 +0.5
9l HC[triple bond, length as m-dash]C (alkynyl) H 25.9 +0.6
9m 4-Pyridyl H 26.4 +1.1
9n 2-Pyridyl H 28.2 +2.9
9o 2,6-Dimethylphenyl H 28.6 +3.3
9p CN H 29.3 +4.0
9q Me Me 31.8 +6.5
9r –CH2(CH2)3CH2 32.7 +7.4
9s Me H 33.9 +8.6
9t C6H11 H 34.0 +8.7
9u H H 37.2 +11.9



image file: c6sc01289b-f3.tif
Fig. 3 Difference in free energies of activation (ΔGB3LYP) between cyclopropane-substituted difluoro-VCP and reference 9j. Green = lower ΔG (<25.3 kcal mol−1), orange = 0–5 kcal mol−1 higher ΔG and red = >5 kcal mol−1 increase in ΔG.

We observed a dramatic rise in calculated activation energy when there was no additional substitution on the difluoro-VCP (9u, +11.9 kcal mol−1). Compounds with no aromatic functionality but one (9s and 9t) or two (9q and 9r) alkyl substituents were also found to have higher barriers for rearrangement (ranging from +6.5 to +8.6 kcal mol−1). We had observed previously that temperatures of 180 °C facilitated a [3,3]-rearrangement pathway with a calculated activation energy of 32.9 kcal mol−1GB3LYP, UB3LYP/6-31G*, gas phase, 298 K, Gaussian'09),3 suggesting that low temperature rearrangement with alkyl substituents would be unlikely. Higher reaction temperatures also have the potential to activate a competitive [1,5]-hydride shift pathway22 so these substitution patterns were ruled out of our synthetic study.

We observed a higher activation energy for 2,6-dimethylphenyl 9o (+3.3 kcal mol−1) due to steric interactions between the methyl protons and the proton (2.26 Å) and fluorine atoms (2.23 and 2.25 Å) attached to the cyclopropane ring in the transition state. These small steric interactions are tolerated in the transition state to ensure the benzylic radical remains coplanar with the aromatic ring. Unlike the bis-alkyl substituents, no barrier-lowering was observed for bis-phenyl 9k (+0.5 kcal mol−1). Steric repulsion between the ortho-protons in TS1k induces a slight rotation of each ring, forcing stabilising aryl groups out of the plane of the benzylic radical (Fig. 4). Sustmann and co-workers investigated this twist angle effect in more detail,23 but no further investigations into bis-arylated difluoro-VCP systems were undertaken because the ΔGB3LYP values were so similar to those for 9j.


image file: c6sc01289b-f4.tif
Fig. 4 Aromatic ring planarity effects on benzyl radical stabilisation in TS1j and TS1k (dihedral angles represented by ϕ, ΔG values (blue) calculated using UB3LYP/6-31G*, gas phase, 298 K, Spartan'10, all energy values are in kcal mol−1).

A wider range of effects was observed when R1 was an heteroaryl group; 2- and 4-pyridyl groups raised the predicted barrier. 2-Furyl, 2-pyrrolyl and 2-thiophenyl groups lowered it significantly, consistent with their known effects of free radical stabilisation as quantified by Creary.24 A piperonyl group also activated the system (see the ESI for a wider discussion of the systems screened and the substituent effects).

Effect of alkene substitution

We used 10 as a template for investigating modifications to the alkene fragment of the precursors, focusing on the effects of alkene configuration (R1versus R2) and radical stabilising substituents (Table 2 and Fig. 5).
Table 2 Effects of alkene substitution on VCP 10 (UB3LYP/6-31G*, gas phase, 298 K, Spartan'10)

image file: c6sc01289b-u2.tif

Entry VCP R1 R2 ΔGB3LYP ΔΔGB3LYPa
a ΔΔGB3LYP = (ΔGB3LYP10x) − (ΔGB3LYP9j).
1 10a CN H 23.2 −2.1
2 10b CON(Me)OMe H 24.4 −0.9
3 10c CH2OH H 25.1 −0.2
4 9j CO 2 Me H 25.3 0.0
5 10d Me H 25.5 0.2
6 10e H H 26.3 1.0
7 10f H CH2OH 28.5 3.2
8 10g H Me 29.7 4.4
9 10h H CN 29.7 4.4
10 10i CO2Me Me 28.3 3.0



image file: c6sc01289b-f5.tif
Fig. 5 Difference in free energy (ΔGB3LYP) between alkene substituted difluorovinylcyclopropanes and reference 9j. Green = lower ΔG (<25.3 kcal mol−1) and orange = 0–5 kcal mol−1 greater ΔG.

The dramatic reactivity difference between the E- and Z-alkene isomers of 9j previously reported3 was maintained when a wider range of substituents was examined. Calculated activation energies for E-substituted precursors (Table 2, entries 1–5) were lower than for the unsubstituted precursor 10e (Table 2, entry 6), whilst the barriers for the corresponding Z-diastereoisomers were higher (Table 2, entries 7–9). The narrow range of free energies of activation observed for the E-diastereoisomers (23.2 to 25.5 kcal mol−1) suggests that all of these precursors should rearrange at temperatures close to or below 100 °C. The narrow spread of barrier heights (ΔΔGB3LYP = 6.5 kcal mol−1) as the alkene substituents vary is half that observed when the similar changes were made on the cyclopropane (see Table 1). This suggests that a wider range of substituents could be tolerated on the alkene fragment since the radical is already stabilised through allylic resonance.

The transition states for the Z-diastereoisomers of Weinreb amide 10b and methyl ester 10i failed to optimise, but the alcohol (10f, +3.2 kcal mol−1), methyl (10g, +4.4 kcal mol−1) and cyanide (10h, +4.4 kcal mol−1) species all optimised with higher activation energies than the corresponding E-series. Disubstituted alkene 10i (Table 2, entry 10) had a higher free energy of activation than ester 9j (+3.0 kcal mol−1) but the introduction of the ester functionality lowered the activation energy for the rearrangement of Z-methyl 10g (difference of 1.4 kcal mol−1 between 10i and 10g).

Synthetic investigations

After the systematic examination of the functional group tolerance of the VCPR using electronic structure calculations, the synthesis of a selection of compounds was undertaken in order to test the computational predictions. While the literature describing the synthesis of difluorocyclopropanes by difluorocarbene transfer is extensive,17h,25 reaction substrates often seem to have been selected for electron-richness and low levels of functionality. We have challenged the published methods extensively in securing this set of compounds; our studies are described fully in the ESI. The approach described below is a pragmatic one, arrived at after extensive optimisation for this focussed set of compounds.

Synthesis of substituted difluoro-VCP

The efficient synthesis of phenyl-VCP 1a (73% over three steps) previously reported, relied on the successful difluorocyclopropanation of commercially available cinnamyl acetate with methyl 2,2-difluoro-2-(fluorosulfonyl)acetate (MDFA, 11). It was believed that the more electron-rich heteroarene substituents would help the cyclopropanation reaction by raising the nucleophilicity of the alkene. However, only decomposition was observed when 13a was exposed to the MDFA conditions (Scheme 3).
image file: c6sc01289b-s3.tif
Scheme 3 Synthesis and failed difluorocyclopropanation of 2-furyl allyl acetate 14a. Conditions: (i) (carbethoxymethylene)-triphenylphosphorane (1.1 eq.), DCM, r.t., 17 h (ii) DIBAL (3 eq.), toluene, −78 °C to r.t., 8 h (iii) Ac2O (1.2 eq.), pyridine (1.2 eq.), DCM, 45 °C, 5 h (iv) MDFA (2.5 eq.), TMSCl (2.5 eq.), KI (2.8 eq.), diglyme (1.17 eq.), 120 °C, 24 h.

A 2nd generation synthesis of precursor 16 from alkenoate 13 was proposed (Scheme 4); alkenoates 13 are easily accessible on a gram scale from the corresponding commercial or easy-to-prepare aldehydes and are more stable olefins for the high temperature reactions involving electron-rich aromatic substituents than allyl acetate 14.


image file: c6sc01289b-s4.tif
Scheme 4 1st and 2nd generation synthetic routes for accessing difluoro-VCP 16.

Wittig reactions of aldehydes 12a–h (see the ESI) as described previously afforded the desired alkenoates 13a–h in good to excellent yields (45–98%, Table 3, entries 1–10).

Table 3 Synthesis of difluorocyclopropyl allyl alcohols 17a–h from aldehydes 12a–ha

image file: c6sc01289b-u3.tif

Entry R1 x 13 15 17
Yieldb (%) Conv.c (%) Yieldb (%) Yieldb (%)
a Conditions: (i) (carbethoxymethylene)triphenylphosphorane (1.1–1.3 eq.), DCM, r.t., 6–20 h (ii) MDFA (2.5 eq.), TMSCl (2.5 eq.), KI (2.8 eq.), diglyme (1.17 eq.), 120 °C, 4 h (iii) DIBAL (3 eq.), toluene or DCM, −78 °C to r.t., 8 h. b Isolated yields. c Determined by 1H NMR. d Reaction time of 24 h. e Starting from commercial ethyl cinnamate 13b. f 13b and 15b were inseparable via column chromatography or distillation. g Aldehyde 12g was synthesised from thiazole and used crude in the olefination reaction (see ESI for details). h Calculated over two steps from thiazole, 4[thin space (1/6-em)]:[thin space (1/6-em)]1 mixture of E[thin space (1/6-em)]:[thin space (1/6-em)]Z-isomers. i Aldehyde synthesised in situ from the oxidation of 2-hydroxymethyl-3-methyl furan (see ESI for details). j Calculated over two steps from methyl 3-methyl-2-furoate.
1 2-Furyl a 82 66 40 94
2 Ph b 28d
3 29d,e n.d.f
4 73 37 n.d.f
5 2-Thiophenyl c 98 77 71 50
6 5-Benzo[d][1,3]-dioxole d 94 50 43 75
7 2-Pyridyl e 87 0
8 2-N-Boc-pyrrolyl f 60 92 54 6
9 2-Thiazolyl g 74h 0
10 3-Me-2-furyl h 45j 80 45 84


Furyl alkenoate 13a proved to be more stable under MDFA-mediated difluorocyclopropanation conditions than acetate 14a but only when a shorter reaction time of 4 hours was used (Table 3, entry 1); prolonged reaction times of 24 h resulted in a decreased conversion to ester 15a (66% compared with 50%, respectively). The conversion to 15a could be increased to 87% by using sodium chlorodifluoroacetate (Na-CDFA, 10 eq.) conditions,25e,26 but a drop in isolated yield was observed (22%) and attributed to product decomposition at the higher reaction temperature of 180 °C over a longer reaction time of 25 hours. Reactions with ethyl cinnamate 13b (Table 3, entries 2–4) showed that this 2nd generation route is less efficient than our previously published entry from cinnamyl acetate, due to the alkene reactivity rather than the conditions used to access alkenoates.

A wider range of heteroarene functionalised difluorocyclopropyl esters could be isolated in moderate to good yields (40–71%) using our shorter difluorocyclopropanation conditions from the corresponding alkenoates (Table 3). 2-Pyridyl 13e and 2-thiazoyl 13g failed to show any signs of reaction. Successful difluorocyclopropanation of a more reactive vinyl pyridine was reported in the patent literature using higher temperature decomposition of sodium chlorodifluoroacetate (Na-CDFA);27 these conditions also failed to cause alkenoate 13e to react. The lack of reactivity was attributed to side reactions between the substituents and difluorocarbene; a full discussion can be found the ESI.

DIBAL-mediated reduction of the isolated difluorocyclopropyl esters to the corresponding alcohols gave mixed results with moderate to excellent yields (50–94%) observed for furyl 17a, thiophenyl 17c, piperonyl 17d and 3′-methyl-furyl 17g analogues (Table 3). Unfortunately, N-Boc pyrrolyl 15f was unstable under reduction conditions and gave a poor 6% yield of 17f.

Vatele's room temperature tandem oxidation/olefination conditions were utilised as the final step in the precursor synthesis, minimising rearrangements during alkenoate formation.28 The one-pot method previously found to be successful with ethyl ester functionalisation only worked with a selection of commercially available phosphoranes (Scheme 5, Method A). However, aldehyde 18b could be isolated in good yield using the same room temperature oxidation chemistry, allowing direct Wittig reactions to proceed smoothly (Scheme 5, Method B).


image file: c6sc01289b-s5.tif
Scheme 5 Oxidation and olefination chemistry used to access a range of alkene functionality. Conditions: (i) BAIB (1.1 eq.), TEMPO (0.1 eq.), DCM, r.t., 5–6 h, (ii) Ph3P[double bond, length as m-dash]C(R1)R2, DCM, r.t., 15–16 h.

Three more alkenes were prepared in good yields from phenyl 17b; these were E-Weinreb amide 19a, α-methyl ester 20a and cyanide 21 (Table 4, entry 1–3, respectively). Furthermore, heteroarene-based building blocks all underwent functionalisation successfully but only piperonyl precursor 22a could be isolated (50%, Table 4, entry 4); the 2-furyl, 2-thiophenyl and 3-methyl-2-furyl congeners all rearranged at ambient temperature before isolation (Table 4, entry 5–10) could be completed.

Table 4 Complete synthesis of difluoro-VCP precursors 19–27a
Entry (b) R1 R2 R3 VCP Yieldc (%)
E-Alkene (a) Z-Alkene (b)
a Compounds represented by numbers and suffix a and b correspond to E- and Z-isomers, respectively. b Synthetic methodology based on Scheme 5. c Isolated yields unless otherwise stated. d 21a and 21b could not be separated by column chromatography and instead were isolated as a 3[thin space (1/6-em)]:[thin space (1/6-em)]2 mixture, respectively (mixture determined by 1H NMR). e 22b formed during the reaction but could not be separated from a mixture with 22a (21% isolated yield of 22a/22b mixture). f All precursors were successfully formed but reactions resulted in complex mixtures due to competing low temperature rearrangements.
1 (A) Ph CON(OMe)Me H 19 74 9
2 (A) Ph CO2Et Me 20 74 0
3 (B) Ph H/CN H/CN 21 84d
4 (A) 2-Piperonyl CO2Et H 22 50 n.d.e
5 (A) 2-Furyl CO2Et H 23 Full conversionf
6 (A) 2-Thiophenyl CO2Et H 24 Full conversionf
7 (A) 2-Furyl CO2Et Me 25 Full conversionf
8 (A) 2-Thiophenyl CO2Et Me 26 Full conversionf
9 (A) 3-Me-2-furyl CO2Et H 27 Full conversionf


Although 2-furyl 9b and 2-thiophenyl 9d had some of the lowest predicted VCPR activation energies (21.0 and 21.3 kcal mol−1, respectively), the observed low temperature rearrangements were still surprising; all reactions were conducted at room temperature and only warmed briefly during the evaporation of solvent (maximum temperature 40 °C). A repeated synthesis of VCP-23 (furan) evaporated the reaction solvent under a stream of nitrogen at ambient temperature and resulted in a similar mixture of products. These results suggest that the rearrangements occurred at room temperature (<20 °C) and not during work up. To our knowledge these are the first examples of VCP precursors undergoing low temperature thermolysis without additional additives; previously only transition metal mediated12a,29 or charge accelerated rearrangements11,30 occurred at temperatures close to ambient.

Deconvolution of crude oxidation/Wittig reaction mixtures for furyl and thiophenyl substrates was difficult due to the presence of triphenylphosphine oxide side product. All of these reaction mixtures were purified twice, first to remove reaction impurities in order to try to reveal the reaction outcome (see ESI for crude reaction spectra), then a final separation to isolate products free from phosphine oxide (Table 5).

Table 5 Reaction outcomes from furyl and thiophenyl precursors 23–26

image file: c6sc01289b-u4.tif

Entry VCP X R1 Crude observationsa (% conversion) Productb (%)
28/29 30/31
a Percentage conversion determined by 1H NMR after 1st purification. b Isolated yield. c Compound decomposed during purification attempts (see ESI for further information). d 50% purity (determined by 1H NMR) containing 29a.
1 23 O H 28a (>95), 30a (trace) n.a.c
2 24 S H 29a major, 31a minor, evidence of 24 29a (12) 31a (3)d
3 25 O Me 28b (100) 28b (48) 30b (0)
4 26 S Me 29b (100) 29b (55) 31b (0)


Precursor 23 (furyl) was absent after the reaction; instead mixed fractions from the first purification confirmed that mono-fluorinated cycloheptadiene 28a was the favoured rearrangement product over difluorocyclopentene 30a (Table 5, entry 1). Attempts at isolating rearrangement product failed due to product decomposition but distinctive 19F NMR signals which were consistent with isolated thiophene product were used for assigning rearrangement outcomes.

The synthesis of thiophenyl 24a/24b resulted in more complex fraction mixtures from the first purification, with NMR evidence for VCP 24a and 24b, as well as rearrangement products 29a and 31a (Table 5, entry 2). Further thermolysis reactions of mixtures containing VCP precursors showed that E-isomer 24a rearranged at 40 °C whereas the corresponding Z-isomer 24b required the higher temperature of 50 °C (see ESI for full discussion). The major mono-fluorinated cycloheptadiene 29a could be isolated in 12% yield and minor difluorocyclopentene 31a in a lower 3% yield (but of only 50% purity because of the presence of cycloheptadiene 29a).

Previous computational investigations predicted that an α-methyl substituent on the alkene for VCP precursor 10i would increase the activation barrier for VCPR rearrangement (vide supra). It was proposed that similar substitution would enable some temperature control of the rearrangement of highly reactive heterocyclic precursors. However, despite an increase in calculated activation energies for furyl 25aGB3LYP = 23.5 kcal mol−1, +2.5 cf.9b) and thiophenyl 26aGB3LYP = 23.7 kcal mol−1, +2.4 cf.9d), experimental results still gave room temperature rearrangements (Table 5, entries 3–4). Interestingly, both precursors reacted exclusively via the [3,3]-pathway, leading to moderate yields of both 28b (48%) and 29b (55%).

We also attempted the synthesis of 3-methyl furyl precursor 27 expecting that the methyl group would cause unfavourable steric interactions in the [3,3]-transition state and instead favour VCPR (Scheme 6a). However, like all other heteroarene substituted precursors, rearrangement was observed before isolation, but 19F NMR reaction monitoring of the tandem oxidation/olefination of alcohol 17h suggested that rearrangement had occurred from aldehyde 18h (see ESI for further discussion).


image file: c6sc01289b-s6.tif
Scheme 6 (a) Oxidation/olefination of alcohol 17h resulting in unexpected rearrangement of aldehyde 18h. (b) Electronic structure calculations for proposed rearrangement (UB3LYP/6-31G*, Spartan'10, gas phase, kcal mol−1). Conditions: (i) BAIB (1.1 eq.), TEMPO (10 mol%), DCM (3 mL), (ii) Ph3P[double bond, length as m-dash]CHCO2Et (1.3 eq.) (one pot).

The two major products were tentatively assigned as dihydrofuran 33 and benzooxepine 34 due to the strong similarities between 19F NMR chemical shifts reported by Hammond31 and ourselves,3 respectively. Electronic structure calculations for the [3,3]-rearrangement viaTS3 supported room temperature rearrangement with a low ΔGUB3LYP value of 20.8 kcal mol−1 (Scheme 6b). Analysis of the VCPR of aldehyde 18h through TS4 gave a higher calculated barrier for rearrangement (ΔGUB3LYP = 25.0 kcal mol−1), but a low spin operator value (S2 = 0.291) suggested that the rearrangement is more likely to be concerted or through a donor–acceptor ring opening/ring closing mechanism.32 Purification of the resulting crude reaction mixture failed, due to either decomposition or volatility of the proposed products. Further synthetic and computational investigations into these appealing fluorinated products are underway but they do not contribute to the development of the computational model in this study.

The observed experimental results with heterocyclic precursors link well to the calculated low activation energies for VCPR. However the lack of control, and in some cases dominance of the [3,3]-pathway, was surprising because temporary dearomatisation is required. Experimental activation energies were required from more controlled rearrangements in order to screen for the best computational methods for assessing this competing pathway.

Thermal rearrangement of isolated difluoro-VCP

The rearrangement temperatures for precursors which could be isolated were optimised to give full consumption of VCP after 17 hours (±5 °C) (Scheme 7). Normalising the reaction temperature against a fixed reaction time allows for the role that different substituents play on rearrangement rates (Table 6) to be appreciated more readily.
image file: c6sc01289b-s7.tif
Scheme 7 Rearrangement conditions for isolated-VCP.
Table 6 Thermal rearrangement of isolated VCP
Entry VCP precursor Temp. (°C) VCPR product
# X R1 R2 R3 # Conv.a (%) Yieldb (%)
a Conversion to product (determined by 1H or 19F NMR). b Isolated yields unless otherwise stated. c Full conversion of VCP precursor was observed. d Decomposition was observed. e Clean product could not be isolated from crude reaction mixture. f 3[thin space (1/6-em)]:[thin space (1/6-em)]2 mixture of 21a and 21b, respectively. g Crude mixture also contains 26% 21b and 3% cis-42 (determined by 19F NMR). h Using crude reaction mixture from entry 4. i 6[thin space (1/6-em)]:[thin space (1/6-em)]1 ratio of difluorocyclopentene 37 and alkene isomer 43 (by 1H NMR). j 22% of cis-cyclopentene 38b was also isolated. k Crude reaction mixture also contains 58% of [3,3]-product (44, by 19F NMR). l 11% of 44 was also isolated.
1 19a F Ph CON(Me)OMe H 95 35 100 97
2 19b F Ph H CON(Me)OMe 160 0c d
3 20a F Ph CO2Et Me 155 36 0c e
4 21a/bf F Ph H/CN H/CN 90 37 71g
5h 21b F Ph H CN 160 100 48i
6 47 H Ph CO2Et H 220 38 100 40j
7 19a F Piperonyl CO2Et H 70 39 42k 18l


Weinreb amide 19a rearranged at 95 °C, 5 °C lower than the corresponding ethyl ester 1a, to afford a near quantitative yield of difluorocyclopentene 35 (Table 6, entry 1), consistent with the slightly lower calculated activation energy (ΔG = −0.9 kcal mol−1). Minor stereoisomer 19b required higher temperatures to induce rearrangement, favouring what seemed to be a [3,3]-pathway over VCPR, consistent with previously investigated Z-alkenoates. Diene 40 was proposed as one of the major products (Fig. 6) but these harsher conditions resulted in some decomposition and poor recovery of observed products (Table 6, entry 2). Despite the slightly lower temperature of 155 °C required for methylated alkenoate 20a (Table 6, entry 3), decomposition was again observed and only firm 19F NMR identification of diene 41 could be achieved (see ESI for further discussion). The higher temperature required for rearrangement of 20a compared with ester 1a was entirely consistent with an increased calculated barrier height (+3.0 kcal mol−1, vide supra).


image file: c6sc01289b-f6.tif
Fig. 6 Potential side products from the rearrangement of difluorinated VCP discussed in Table 6.

It was predicted that the E-isomer (21a) would rearrange more rapidly in the isolated mixture of nitriles-21a and 21b since the calculated activation energy for 21a was 6.5 kcal mol−1 lower than that of the Z-isomer (21b). Thermolysis of the mixture at 90 °C over 17 hours gave full conversion of 21a to difluorocyclopentene 37; 21b only showed cyclopropane stereoisomerisation to cis-42 at this temperature (Table 6, entry 4). Only when the resulting mixture was re-heated to 160 °C was full conversion of the Z-isomers observed (Table 6, entry 5), affording difluorocyclopentene 37 in a 48% yield in a 6[thin space (1/6-em)]:[thin space (1/6-em)]1 ratio with alkene regioisomer 43.

The computational triage was not limited to fluorinated precursors; VCP 47 had a calculated activation barrier of 33.4 kcal mol−1 for rearrangement through TS5 (Scheme 8). Simmons–Smith cyclopropanation33 of commercially available cinnamyl alcohol 46, followed by tandem oxidation/olefination, afforded the desired precursor 47 in an unoptimised 35% yield over two steps. A much higher temperature of 220 °C was required for full conversion of 47 to cyclopentene 38 over 17 hours (Table 6, entry 6), 120 °C higher than required for the corresponding difluorinated precursor 1a. A mixture of trans-38a and cis-38b diastereoisomers was observed and chromatographic separation gave 40% and 22% yields of the two products, respectively. Electronic structure calculations showed that transition states representing the VCPR from both trans-47 (TS5a, ΔGB3LYP = 33.4 kcal mol−1) and the corresponding cis-isomer (TS5b, ΔGB3LYP = 34.5 kcal mol−1) had similar energies and could therefore be competitive at high temperatures.


image file: c6sc01289b-s8.tif
Scheme 8 Two step synthesis of non-fluorinated VCP-47; activation energies to TS5G value (in blue) were calculated using (U)B3LYP/6-31G* from intermediate 47, gas phase, 298 K, Spartan'10, using conformationally simpler Me ester, units are kcal mol−1). Conditions: (i) ZnEt2 (1 M in hexane, 5 eq.), CH2I2 (10 eq.), 0 °C to r.t., 2.5 h; (ii) TEMPO (0.1 eq.), BAIB (1.1 eq.), DCM, r.t., 3.5 h; (iii) Ph3P[double bond, length as m-dash]C(H)CO2Et (1.3 eq.), 20 h.

This result conclusively shows the accelerative effect of gem-difluorination, and justifies the decision not to invest time in the synthesis of precursors predicted to have activation energies greater than 30 kcal mol−1. In fact, experimental results from compounds 20a and 21b suggest that the maximum temperature for synthetically useful VCPR with fluorinated precursors could be much lower than expected.

Piperonyl species 22a had a much lower predicted activation energy (ΔGB3LYP = 23.0 kcal mol−1) and subsequently rearranged at much lower optimised temperature of 70 °C (Table 6, entry 7). However, like compounds containing the highly activating heteroarene substituents, both difluorocyclopentene 39 and heptadiene 44 were observed and could be isolated after preparative HPLC in low yields (18% and 11% respectively). Diene 44 was the exclusive product of the [3,3]-rearrangement, despite the possibility of the formation of regioisomer 45via reaction at aromatic carbon centre C4 (Scheme 9). Electronic structure calculations were consistent with experimental results, showing that the observed pathway through TS6a had a lower activation energy than that through TS6b (ΔΔG = 5.3 kcal mol−1).


image file: c6sc01289b-s9.tif
Scheme 9 Electronic structure calculations used to predict the favoured [3,3]-rearrangement pathway from 48G values are relative to 48 (blue) and calculated using UB3LYP/6-31G*, gas phase, 298 K, Spartan'10, units are in kcal mol−1).

The greater thermal control observed from the rearrangement of 22a, allowed the competition between VCPR and [3,3]-pathways to be examined fully and for the first time, using VT 19F NMR spectroscopy (344–373 K in [D8]toluene (see the ESI for experimental details and a full discussion on kinetic modelling)). Unlike the VCPR of phenyl 1a monitored previously, no evidence of the cis-VCP was observed during thermolysis because the [3,3]-pathway was now competitive with VCPR. At this higher rearrangement temperature (100 °C), full consumption of 22a was observed after 30 minutes, contrasting with the 10 hours required for full conversion of 1a and providing further experimental evidence that the activation energy for the latter was higher, as predicted. An Arrhenius determination of activation parameters was carried out; the value for VCPR Ea from piperonyl 22a of (23.4 ± 0.2) kcal mol−1 was very close to the calculated ΔGB3LYP (23.0 kcal mol−1). A slightly higher Ea value of (24.9 ± 0.3) kcal mol−1 was calculated for the [3,3]-pathway; these values were used to screen for the best computational method for treating this manifold of reactions.

Predictive electronic structure calculations

We carried out an extensive set of calculations to determine the most appropriate methods for the treatment of VCPR and competing [3,3]-rearrangement. The study is reported fully in the ESI and only the main findings will be summarised here. The UM05-2X/6-31+G* method gave the closest agreement between experimental values for the VCPR of 1a and 22a, while the M06-2X/6-31G* method was the best for the [3,3]-rearrangement of 48; these Minnesota functionals were selected as the most accurate methods for assessing the two pathways. Because of its consistency of performance, the (U)B3LYP/6-31G* was also retained as a low cost method able to handle both rearrangements at a useful level of accuracy. The calculated ΔG values for VCPR using each of these methods were corrected by the average error (ΔG) from 9f and 9j; these values are +2.7 kcal mol−1 for UB3LYP and −1.2 kcal mol−1 for UM05-2X. Since the experimental activation energy for the [3,3]-rearrangement could only be determined for 48, the correcting values are based solely on this compound's ΔG values, which are +2.4 kcal mol−1 for B3LYP and +0.5 kcal mol−1 for M06-2X.

The differences in corrected ΔGVCPR and corrected ΔG[3,3] values ((ΔGVCPRΔ) − (ΔG[3,3])) were then used to predict which rearrangement pathway would be favoured by either B3LYP or Minnesota methods; values greater than +1.0 kcal mol−1 would favour VCPR and values less than −1.0 kcal mol−1 would favour sigmatropic rearrangement. Small differences between these values could result from computational error (±0.5 kcal mol−1) and would represent the limit of our ability to predict the composition of rearrangement product mixtures or the identity of the dominant pathway (Table 7).

Table 7 Predicting thermal rearrangement pathways (VCPR or [3,3]) from VCPs using corrected free energies of activation (ΔG)
a Calculated ΔG + 2.7 kcal mol−1 (UB3LYP/6-31G*). b Calculated ΔG + 2.4 kcal mol−1 (B3LYP/6-31G*). c ΔG = (corrected VCPR ΔGUB3LYP) − (corrected [3,3] ΔGB3LYP). d Calculated ΔG − 1.2 kcal mol−1 (UM05-2X/6-31+G*). e Calculated ΔG + 0.5 kcal mol−1 (M06-2X/6-31G*). f ΔG = (corrected VCPR ΔGUM052X) − (corrected [3,3] ΔGM062X). All calculations were performed in gas phase at 298 K using Gaussian'09, units are in kcal mol−1. Blue values predict VCPR pathway, orange values predict a mixture of pathways and green values predict [3,3]-rearrangement.
image file: c6sc01289b-u5.tif


Since the correction factors were derived from their experimental values, it was no surprise that phenyl 1a was correctly predicted to undergo VCPR and that piperonyl 22a was predicted to give a mixture of products (Table 7, entry 1 and 2, respectively). More pleasingly, the major rearrangement pathways for the nine out of the ten novel difluoro-VCP systems examined were all predicted correctly using B3LYP/6-31G* calculations; for phenyl 1a, the VCPR pathway was correctly predicted as the major one but is within computational error (ΔG = +0.9 kcal mol−1, Fig. 7).


image file: c6sc01289b-f7.tif
Fig. 7 Prediction plot of the difference between corrected ΔGVCPR (UB3LYP/6-31G* or UM05-2X/6-31+G*) and corrected ΔG[3,3] (B3LYP/6-31G* or M06-2X/6-31G*, respectively) for synthesised VCP. Colour used to represent the correct (green) or incorrect (red) predictions compared with experimental observations.

Both methods failed to deal with the non-fluorinated system, predicting a mixture of rearrangement products instead of only the VCPR observed experimentally. The Minnesota functionals also failed to deal with more sterically hindered heteroarene VCPs 25 and 26, but it is unknown whether the errors arise from the VCPR or [3,3]-rearrangement calculations, or from a combination of both. These results strongly suggest that the lowest cost method is comparable and in some case better than the more expensive Minnesota methods (we refer to implementation in Spartan), consistent with studies carried out by Simón and Goodman.34 Since all experimental rearrangements were optimised to give full conversion of VCP after 17 hours, a strong trend was apparent between the corrected ΔG values for VCPR with reaction temperatures using either M05-2X/6-31+G* (Fig. 8) or B3LYP/6-31G* methods (see ESI).


image file: c6sc01289b-f8.tif
Fig. 8 Correlation between corrected ΔGUM052X values (kcal mol−1) for VCPR and the optimised reaction temperatures (K) which gave 100% conversion of VCP. Error in ΔGUM052X values = ±0.5 kcal mol−1. Error in temperature = ±5 K.

As more compounds are synthesised, the error associated with these models may be reduced. However, from the small set of varied difluorinated vinylcyclopropanes that were examined, the most effective computational models look reliable enough to be used with confidence in the design and assessment of new precursors before any synthetic commitments are required.

A set of VCP from the literature were selected to test the predictive capability of the developed computational models and the rearrangement pathway for all four compounds was successfully predicted using the lower cost UB3LYP/6-31G* method (Fig. 9). The indole-vinylcyclopropane rearrangement of 508b and divinylcyclopropane rearrangement of 498c were both favoured over VCPR, whereas Sustmann and co-workers bis-aryl VCP 51 was correctly predicted to undergo VCPR.35 The calculated ΔGVCPR value of 31.0 ± 0.5 kcal mol−1 for 51 is within error of the reported experimental activation energy of 32.8 ± 1.6 kcal mol−1 (corrected to 298 K). Our temperature prediction suggests that the rearrangement could give full conversion after 17 hours at 130 °C, 70 °C lower than the reported conditions (reaction time of 1.5 hours). Difluoro-VCP 52 could only undergo VCPR and a calculated ΔGVCPR value of 28.7 kcal mol−1 was very similar to phenyl-VCP 1a and predicted to rearrange at the same temperature of 100 °C. These results strongly suggest that the Ni-catalyst present during the reaction of 52 at 140 °C does not facilitate VCPR, a factor that was not obvious from experimental observations.17h


image file: c6sc01289b-f9.tif
Fig. 9 Testing the predictive capability of electronic structure models against compounds which undergo selective VCPR or [3,3]-rearrangement (simplified computational models were used for 50 (Ts replaced with Ms) and 52 (TBS replaced with TMS)). Free energies of activation (ΔG) calculated on Gaussian'09 using UB3LYP/6-31G* (gas phase, 298 K) and quoted in kcal mol−1. Predicted temperature derived from the straight line equation y = 0.0721x + 2.0341 from a plot of corrected ΔGB3LYP against VCPR rearrangement temperatures (see ESI).

Conclusion

A low cost computational assessment (UB3LYP/6-31G*) of substituent effects on the VCPR of difluorinated vinylcyclopropanes was used to guide the synthesis of a test set of novel difluorocyclopropanes. The VCPR system examined is most unusual in that reactions involving open shell singlets and triplet species are concurrent with a more classically concerted pericyclic rearrangement, posing a significant challenge to available computational methods.

Radical stabilising groups, specifically heteroarenes, were found to lower calculated free energies of activation more when bound directly to the cyclopropane instead of to the alkene, consistent with the open shell singlet mechanism for VCPR, and bringing about rearrangements at room temperature in some cases. The most reactive systems based on heteroarenes, underwent rearrangements at unexpectedly low temperatures, through competing VCPR and aromatic-vinylcyclopropane rearrangements to give access to both novel difluorocyclopentenes and fluorinated benzocycloheptadienes, respectively. Optimised rearrangement temperatures for isolated precursors showed a good trend with calculated activation energies, allowing estimates of rearrangement temperatures to be made before synthesis. Comparison of predictions made by electronic structure calculations with experimental activation energies for piperonyl 22a and literature examples showed that the (U)M05-2X/6-31+G* method remained the most accurate for assessing VCPR, but M06-2X/6-31G* calculations were better for the aromatic-vinylcyclopropane rearrangement. No single method stood out overall but the consistency in error observed with (U)B3LYP/6-31G* calculations for both pathways, meant that it came closest to a universal method for dealing with the reaction manifold. There was no simple relationship between the amount of HF-exchange and the accuracy of the predictions. The selectively for rearrangement pathways could be predicted accurately using electronic structure calculations, either with the Minnesota functionals or lower cost DFT methods ((U)B3LYP/6-31G*). The computational design model developed was tested against literature compounds and was found to predict observed experimental results correctly.

The ability to determine whether a VCP molecule will rearrange thermally at a synthetically useful temperature, and to predict which pathway it will take, is an extremely powerful development, and shows that effective triage of synthetic chemistry programmes is not only effective, but also practicable by non-specialists.

Acknowledgements

We thank Dr Vipulkumar Patel (GSK, Flexible Discovery Unit) and the University of Strathclyde (studentship to D.O.), for funding, the EPRSC National Mass Spectrometry Service Centre, Swansea for accurate mass measurements, Craig Irving for help with variable temperature NMR kinetic work and Dr Tell Tuttle (University of Strathclyde) and Dr Colin Edge (GSK, Department of Platform Technology and Science) for helpful suggestions.

References

  1. (a) S. M. Bachrach, Computational Organic Chemistry, John Wiley & Sons, Inc., Hoboken, New Jersey, 2nd edn, 2014 Search PubMed; (b) C. J. Cramer, Essentials of Computational Chemistry: Theories and Models, John Wiley & Sons, Ltd., Chichester, West Sussex, 2nd edn, 2004 Search PubMed.
  2. (a) Y. Liang, J. L. Mackey, S. A. Lopez, F. Liu and K. N. Houk, J. Am. Chem. Soc., 2012, 134, 17904–17907 CrossRef CAS PubMed; (b) G. Jiménez-Osés, A. J. Brockway, J. T. Shaw and K. N. Houk, J. Am. Chem. Soc., 2013, 135, 6633–6642 CrossRef PubMed; (c) K. N. Houk and P. H.-Y. Cheong, Nature, 2008, 455, 309–313 CrossRef CAS PubMed; (d) Q. N. N. Nguyen and D. J. Tantillo, Chem.–Asian J., 2014, 9, 674–680 CrossRef CAS PubMed; (e) A. S. K. Tsang, I. A. Sanhueza and F. Schoenebeck, Chem. –Eur. J., 2014, 20, 16432–16441 CrossRef CAS PubMed; (f) R. N. Straker, Q. Peng, A. Mekareeya, R. S. Paton and E. A. Anderson, Nat. Commun., 2016, 7, 10109 CrossRef CAS PubMed.
  3. D. Orr, J. M. Percy, T. Tuttle, A. R. Kennedy and Z. A. Harrison, Chem. –Eur. J., 2014, 20, 14305–14316 CrossRef CAS PubMed.
  4. T. Hudlicky and J. W. Reed, Angew. Chem., Int. Ed., 2010, 49, 4864–4876 CrossRef CAS PubMed.
  5. Y. Tu-Hsin and L. A. Paquette, Tetrahedron Lett., 1982, 23, 3227–3230 CrossRef.
  6. C. K. Murray, D. C. Yang and W. D. Wulff, J. Am. Chem. Soc., 1990, 112, 5660–5662 CrossRef CAS.
  7. E. A. Ilardi, C. E. Stivala and A. Zakarian, Chem. Soc. Rev., 2009, 38, 3133–3148 RSC.
  8. (a) T. Hudlicky, R. Fan, J. W. Reed and K. G. Gadamasetti, Org. React., 2004, 41, 1–133 Search PubMed; (b) S. Krüger and T. Gaich, Beilstein J. Org. Chem., 2014, 10, 163–193 CrossRef PubMed; (c) P. Erbes and W. Boland, Helv. Chim. Acta, 1992, 75, 766–772 CrossRef CAS.
  9. (a) J. P. Olson and H. M. L. Davies, Org. Lett., 2008, 10, 573–576 CrossRef CAS PubMed; (b) J. P. Olson and H. M. L. Davies, Org. Lett., 2010, 12, 1144 CrossRef CAS.
  10. P. J. Gritsch, E. Stempel and T. Gaich, Org. Lett., 2013, 15, 5472–5475 CrossRef CAS PubMed.
  11. R. L. Danheiser, C. Martinez-Davila and J. M. Morin, J. Org. Chem., 1980, 45, 1340–1341 CrossRef CAS.
  12. (a) G. Zuo and J. Louie, Angew. Chem., Int. Ed., 2004, 43, 2277–2279 CrossRef CAS PubMed; (b) S. C. Wang and D. J. Tantillo, J. Organomet. Chem., 2006, 691, 4386–4392 CrossRef CAS.
  13. (a) W. R. Dolbier, B. H. Al-Sader, S. F. Sellers and H. Koroniak, J. Am. Chem. Soc., 1981, 103, 2138–2139 CrossRef CAS; (b) W. R. Roth, W. Kirmse, W. Hoffmann and H.-W. Lennartz, Chem. Ber., 1982, 115, 2508–2515 CrossRef CAS.
  14. (a) S. W. Benson, A. N. Bose and P. Nangia, J. Am. Chem. Soc., 1963, 85, 1388–1389 CrossRef CAS; (b) D. N. Zeiger and J. F. Liebman, J. Mol. Struct., 2000, 556, 83–94 CrossRef CAS.
  15. W. R. Dolbier, Acc. Chem. Res., 1981, 14, 195–200 CrossRef CAS.
  16. (a) D. K. Lewis, D. J. Charney, B. L. Kalra, A.-M. Plate, M. H. Woodard, S. J. Cianciosi and J. E. Baldwin, J. Phys. Chem. A, 1997, 101, 4097–4102 CrossRef CAS; (b) W. R. Dolbier and S. F. Sellers, J. Am. Chem. Soc., 1982, 104, 2494–2497 CrossRef CAS.
  17. (a) L. A. Buttle and W. B. Motherwell, Tetrahedron Lett., 1994, 35, 3995–3998 CrossRef CAS; (b) F. Barth and C. O. Yang, Tetrahedron Lett., 1991, 32, 5873–5876 CrossRef CAS; (c) T. Taguchi, M. Kurishita, A. Shibuya and K. Aso, Tetrahedron, 1997, 53, 9497–9508 CrossRef CAS; (d) T. Mase, I. N. Houpis, A. Akao, I. Dorziotis, K. Emerson, T. Hoang, T. Iida, T. Itoh, K. Kamei, S. Kato, Y. Kato, M. Kawasaki, F. Lang, J. Lee, J. Lynch, P. Maligres, A. Molina, T. Nemoto, S. Okada, R. Reamer, J. Z. Song, D. Tschaen, T. Wada, D. Zewge, R. P. Volante, P. J. Reider and K. Tomimoto, J. Org. Chem., 2001, 66, 6775–6786 CrossRef CAS PubMed; (e) P. E. Harrington, L. Li and M. A. Tius, J. Org. Chem., 1999, 64, 4025–4029 CrossRef CAS; (f) C. Qing-Yun, J. Fluorine Chem., 1995, 72, 241–246 CrossRef; (g) D. Munemori, K. Narita, T. Nokami and T. Itoh, Org. Lett., 2014, 16, 2638–2641 CrossRef CAS PubMed; (h) T. Aono, H. Sasagawa, K. Fuchibe and J. Ichikawa, Org. Lett., 2015, 17, 5736–5739 CrossRef CAS PubMed.
  18. (a) G. Theodoridis, in Advances in Fluorine Science, ed. T. Alain, Elsevier Science, 1st edn, 2006, vol. 2, pp. 121–175 Search PubMed; (b) S. Swallow, in Fluorine in Pharmaceutical and Medicinal Chemistry From Bioaspects to Clinical Applications, ed. V. Gouverneur and K. Müller, Imperial College Press, London, 2012, pp. 141–174 Search PubMed; (c) M. Bremer, P. Kirsch, M. Klasen-Memmer and K. Tarumi, Angew. Chem., Int. Ed., 2013, 52, 8880–8896 CrossRef CAS PubMed; (d) J. Wang, M. Sánchez-Roselló, J. L. Aceña, C. del Pozo, A. E. Sorochinsky, S. Fustero, V. A. Soloshonok and H. Liu, Chem. Rev., 2013, 114, 2432–2506 CrossRef PubMed.
  19. J. M. Percy, Top. Curr. Chem., 1997, 193, 131–195 CrossRef CAS.
  20. Reviewers have pointed out that spin projection methods or more elaborate calculations (CAS-SCF) can be used to deal with the high degrees of spin contamination which arose in our earlier publication. We wished to treat full-sized experimental systems and therefore evaluated the effectiveness of DFT calculations. Our approach is precedented by Cremer et al.. See J. Grafenstein, E. Kraka, M. Filatov and D. Cremer, Int. J. Mol. Sci., 2002, 3, 360–394. CrossRef Cremer states “However, if the overlap between the open-shell orbitals is small (the single electrons are located at different atomic centers) errors become small and, then, the principal failure of UDFT in these cases is not apparent”.
  21. (a) A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS; (b) C. Lee, W. Yang and R. G. Parr, Phys. Rev. B: Condens. Matter Mater. Phys., 1988, 37, 785–789 CrossRef CAS; (c) S. H. Vosko, L. Wilk and M. Nusair, Can. J. Phys., 1980, 58, 1200–1211 CrossRef CAS; (d) P. J. Stephens, F. J. Devlin, C. F. Chabalowski and M. J. Frisch, J. Phys. Chem., 1994, 98, 11623–11627 CrossRef CAS; (e) the seminal work was carried out by Houk and co-workers and who state that “major pathway involves a single diradicaloid transition state” associated with 〈S2〉 = 0.75 (characteristic of a doublet state). We use the term “open shell singlet” to cohere with the sense of Houk's original expression (thanking a reviewer for the recommendation); see K. N. Houk, M. Nendel, O. Wiest and J. W. Storer, J. Am. Chem. Soc., 1997, 119, 10545–10546 CrossRef CAS; (f) for a discussion of spin operator values for radical systems, see K. C. Mondal, S. Roy and H. W. Roesky, Chem. Soc. Rev., 2016, 45, 1080–1111 RSC.
  22. R. J. Ellis and H. M. Frey, J. Chem. Soc., 1964, 5578–5583 RSC.
  23. W. Sicking, R. Sustmann, J. Mulzer and R. Huisgen, Helv. Chim. Acta, 2011, 94, 1389–1405 CrossRef CAS.
  24. X. Creary, Acc. Chem. Res., 2006, 39, 761–771 CrossRef CAS PubMed.
  25. (a) W. R. Dolbier Jr, F. Tian, J.-X. Duan, A.-R. Li, S. Ait-Mohand, O. Bautista, S. Buathong, J. Marshall Baker, J. Crawford, P. Anselme, X. H. Cai, A. Modzelewska, H. Koroniak, M. A. Battiste and Q.-Y. Chen, J. Fluorine Chem., 2004, 125, 459–469 CrossRef; (b) S. Eusterwiemann, H. Martinez and W. R. Dolbier, J. Org. Chem., 2012, 77, 5461–5464 CrossRef CAS PubMed; (c) F. Wang, T. Luo, J. Hu, Y. Wang, H. S. Krishnan, P. V. Jog, S. K. Ganesh, G. K. S. Prakash and G. A. Olah, Angew. Chem., Int. Ed., 2011, 50, 7153–7157 CrossRef CAS PubMed; (d) F. Wang, W. Zhang, J. Zhu, H. Li, K.-W. Huang and J. Hu, Chem. Commun., 2011, 47, 2411–2413 RSC; (e) Y. Fujioka and H. Amii, Org. Lett., 2008, 10, 769–772 CrossRef CAS PubMed; (f) K. Oshiro, Y. Morimoto and H. Amii, Synthesis, 2010, 12, 2080–2084 Search PubMed.
  26. (a) R. Csuk and L. Eversmann, Tetrahedron, 1998, 54, 6445–6456 CrossRef CAS; (b) A. Shibuya, A. Sato and T. Taguchi, Bioorg. Med. Chem. Lett., 1998, 8, 1979–1984 CrossRef CAS PubMed.
  27. E. J. Brnardic, A. Converso, M. E. Fraley, R. M. Garbaccio and S. Y. Huang, PCT Int. Appl., WO 2009140163 A1, 2009.
  28. J.-M. Vatèle, Tetrahedron Lett., 2006, 47, 715–718 CrossRef.
  29. V. Aris, J. M. Brown, J. A. Conneely, B. T. Golding and D. H. Williamson, J. Chem. Soc., Perkin Trans. 2, 1975, 4–10 RSC.
  30. (a) S. D. Larsen, J. Am. Chem. Soc., 1988, 110, 5932–5934 CrossRef CAS; (b) T. Hudlicky, N. E. Heard and A. Fleming, J. Org. Chem., 1990, 55, 2570–2572 CrossRef CAS.
  31. S. Arimitsu and G. B. Hammond, J. Org. Chem., 2007, 72, 8559–8561 CrossRef CAS PubMed.
  32. T. F. Schneider, J. Kaschel and D. B. Werz, Angew. Chem., Int. Ed., 2014, 53, 5504–5523 CrossRef CAS PubMed.
  33. S. E. Denmark and J. P. Edwards, J. Org. Chem., 1991, 56, 6974–6981 CrossRef CAS.
  34. L. Simón and J. M. Goodman, Org. Biomol. Chem., 2011, 9, 689–700 Search PubMed.
  35. J. Mulzer, R. Huisgen, V. Arion and R. Sustmann, Helv. Chim. Acta, 2011, 94, 1359–1388 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available: Computational methodology, Cartesian coordinates for intermediates and transition states, synthetic procedures and characterisation spectra for novel compounds and further discussion on points highlighted in the text. See DOI: 10.1039/c6sc01289b

This journal is © The Royal Society of Chemistry 2016
Click here to see how this site uses Cookies. View our privacy policy here.