Jordi Vila1, Miquel Solà1, Anna Pla-Quintana1, Anna Roglans1. 1. Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, Facultat de Ciències, Universitat de Girona (UdG), C/ Maria Aurèlia Capmany, 69, Girona 17003, Catalonia, Spain.
Abstract
The synthesis of spiro compounds featuring seven- and six-membered rings in the spirobicyclic motif is successfully achieved through a cascade process encompassing a rhodium(I)-catalyzed cycloisomerization followed by a highly selective Diels-Alder homodimerization. The scope of the reaction is analyzed based on a series of synthetic substrates, and control experiments and DFT calculations led us to justify the exquisite degree of selectivity observed.
The synthesis of spiro compounds featuring seven- and six-membered rings in the spirobicyclic motif is successfully achieved through a cascade process encompassing a rhodium(I)-catalyzed cycloisomerization followed by a highly selective Diels-Alder homodimerization. The scope of the reaction is analyzed based on a series of synthetic substrates, and control experiments and DFT calculations led us to justify the exquisite degree of selectivity observed.
The synthesis of spiro
compounds, structures containing two rings
connected through a single carbon atom (referred to as spiro atoms),
has received significant attention. This can be explained by their
broad spectrum of biological and pharmaceutical activity[1,2] as well as their applications in organic optoelectronics[3] and catalysis.[4] The
presence of a spirobicyclic scaffold in a molecule imposes an inherent
rigidity that greatly influences its chemical and physical properties.
In drug discovery, for example, the rigidity of the spirobicyclic
motif has been used to design the 3D orientation of the pharmacophore
more efficiently. This maximizes intermolecular interactions and thus
improves the recognition process.[5] However,
the presence in a molecule of a quaternary spiro atom, which can present
central or axial chirality, poses major challenges for its synthesis.[6]Plethoras of natural products contain spirobicyclic
moieties, but
spiro compounds containing a seven-membered ring in the spirobicyclic
motif are not particularly common. However, those that contain a nitrogen
atom in the seven-membered ring present very interesting pharmaceutical
properties (Figure ). Regarding derivatives of the spirooxindole family,[6c] spirodiazepineoxindole a has been
evaluated for its antimicrobial and antianxiety activity, providing
notable results,[7] and spiroazepineindole b shows antiplasmodial activity toward Plasmodium
falciparum, the most relevant malaria parasite.[8] A peptidomimetic containing the spiroazepinoxindole
moiety (c, Figure ) has also shown antiproliferative activity on human prostatic
carcinoma cell lines.[9] On the other hand,
if we focus on spirobicyclic motifs containing six and seven cyclic
rings, then galanthamine (d) is a natural product that
is used to slow down the process of neurological degeneration in Alzheimer’s
disease.[10] Furthermore, spirolide derivatives[11] such as pinnatoxins A–D and pteriatoxins
A–C (e), isolated from Pinna muricata(12) and Pteria penguin(13) shellfish, respectively, are potent
marine toxins containing a spirocyclic seven-membered imine that is
key for its toxic activity.[14,15]
Figure 1
Bioactive molecules that
contain seven-membered azaspiro scaffolds.
Bioactive molecules that
contain seven-membered azaspiro scaffolds.A straightforward strategy to construct a spirocyclic moiety containing
a six-membered ring is based on the Diels–Alder cycloaddition
between an exocyclic dienophile and a diene. In an ongoing project
of our research group aimed at developing rhodium-catalyzed cyclization
reactions involving allenes,[16] we recently
described[17] a rhodium-catalyzed cycloisomerization/Diels–Alder
cascade reaction of bisallenes and alkenes that afforded dihydroazepine-
and dihydrooxepine-fused ring systems in good yields. In a control
experiment excluding the alkene (that was designed to obtain mechanistic
information), we observed that cycloheptatriene (cHT),
obtained as a non-isolable intermediate in the reaction, dimerized
through a Diels–Alder reaction to afford spirocyclic derivatives
in a highly selective manner (Scheme ). This dimerization is remarkably analogous to the
postulated biosynthesis of xylopidimers A and B through [4 + 2] Diels–Alder
cycloaddition of two guaiane moieties (Scheme ), from which different orientations explain
the formation of the various regioisomers that are isolated.[18] Due to the interest both in the process and
the properties of the products obtained, here we report the preparation
of seven-membered azaspiro compounds and a full analysis of the reasons
behind the exquisite degree of selectivity in their formation.
Scheme 1
Synthesis of Spiro Derivatives through Diels–Alder Reactions
Results and Discussion
We started
our study with N-tosyl-tethered bisallene 1a (Table ). Based
on our earlier study,[17] we carried
out the dimerization of 1a using the previously optimized
reaction conditions (entry 1, Table ). Employing 10 mol % of cationic rhodium complex [Rh(cod)2]BF4 with DTBM-Segphos in THF/CH2Cl2 (4:1), two different compounds were obtained: compound 2a (47% yield) resulting from a dimerization reaction of intermediate cHT, isolated as a single regioisomer, and 3a (14% yield) derived from a [2 + 2] cyclization of the two internal
double bonds of the two allene moieties of 1a.[19] Since the in situ-generated
cycloheptatriene cHT has three potential double bonds
to be involved in the Diels–Alder reaction as dienophiles,
we first proceeded to fully characterize the compound obtained. Figure shows all the possible
cycloadducts that can be formed in the homo-Diels–Alder cycloaddition
of cHT ordered according to the double bond that acts
as a dienophile (dp).
Table 1
Optimization of the Rhodium(I)-Catalyzed
Dimerization of Bisallene 1a
The solvent mixture THF/DCM was
non-anhydrous and degassed.
Figure 2
Possible cycloadducts that can be formed
in the homo-Diels–Alder
of cycloheptatriene cHT.
Possible cycloadducts that can be formed
in the homo-Diels–Alder
of cycloheptatriene cHT.L1 = (R)-DTBM-Segphos/L2
= BINAP/L3 = Tol-BINAP/L4 = BIPHEP.The solvent mixture THF/DCM was
non-anhydrous and degassed.The molecular formula of dimeric compound 2a was determined
by HRMS, showing a peak at m/z =
573.1836 corresponding to the [M + Na]+ adduct. 1D and
2D NMR spectroscopic experiments clearly showed the formation of only
one cycloadduct, and its complete analysis was carried out to ascertain
the structure of the product formed. The 1H NMR chemical
shift of the six olefinic protons and their multiplicity were first
analyzed. Two pairs of doublets, at δ = 4.71 and 6.27 ppm and
δ = 4.84 and 6.64 ppm, are characteristic of a structure that
has two cis endocyclic double bonds. In addition, only two singlets
at δ = 4.52 and 4.63 ppm are observed, corresponding to the
two geminal protons of a single exocyclic double bond. This allowed
us to discard cycloadducts 1A and 1B. In
order to distinguish between cycloadducts formed by the reaction of dp_2 and dp_3, 2D NMR experiments were conducted.
The HMBC spectrum displays a three-bond correlation between the spiro
carbon atom and the proton closest to the nitrogen of one of the cis
endocyclic olefins (C7-H5, in Figure ). This is, in fact, a clear demonstration that it
is the exocyclic double bond conjugated to the endocyclic double bond
that is involved in the Diels–Alder reaction (dp_2), and thus cycloadducts 3A and 3B can
be discarded. Finally, the analysis of the HMBC, COSY, and NOESY experiments
allowed us to distinguish between regioisomers 2A and 2B (Figure ). NOE contacts between H9 (H8 and H9 identified as the two contiguous
methylenic groups by COSY) and H11 protons led us to assign 2A as the single cycloadduct that formed. The HMBC crosspeak
between H11 and C9 also supports the formation of the cycloadduct 2A. Selected signals in the 1H and 13C NMR spectra of 2a are shown in Figure .
Figure 3
COSY and HMBC crosspeaks and NOE contacts observed
in 2a confirming the formation of regioisomer 2A and selected 1H and 13C (in italics) NMR shifts
of 2a.
COSY and HMBC crosspeaks and NOE contacts observed
in 2a confirming the formation of regioisomer 2A and selected 1H and 13C (in italics) NMR shifts
of 2a.Since a highly chemo-
and regioselective cascade process was found,
we searched for the best reaction conditions to allow the efficient
dimerization of bisallene 1a while avoiding the formation
of byproduct 3a (Table ).First, the effect of the ligand was examined.
BINAP, Tol-BINAP,
and BIPHEP bisphosphines were tested, but spiro derivative 2a was obtained in low yields (entries 2–4, Table ). DPEphos, xantphos, and Xphos
did not promote the reaction. Therefore, the bulky phosphine, DTBM-Segphos,
was established as the ligand of choice. The reaction was then tested
at both a lower and higher concentration of 1a (entries
5 and 6, Table ),
and the yield rose up to 53% when the concentration was increased
to 36 mM. The effect of an additive using the dimeric rhodium complex
[Rh(cod)Cl]2 was next examined (entries 7 and 8, Table ). The mixture of
[Rh(cod)Cl]2/NaBArF maintained the yield of 2a at 46% but decreased the yield of byproduct 3a (6%)
considerably. We then evaluated the effect of water on the reaction
mixture. Using non-anhydrous solvents taken directly from the bottles
but degassed, an 89% yield of 2a was obtained (entry
9, Table ), resulting
in our set of optimized conditions. Since a chiral ligand was used,
we checked the optical purity of 2a; however, no enantioinduction
was observed. This result points to a mechanism involving a rhodium-catalyzed
cycloisomerization coupled to a Diels–Alder cycloaddition in
which the chiral rhodium complex does not participate (vide
infra), in line with our previous study.[17]The scope of the reaction was then evaluated (Figure ). The nature of
the substituents
at the phenyl ring of the sulfonamide tether in bisallene 1a was explored. The reaction proceeded efficiently with both electron-donating
(2b) and electron-withdrawing groups (2c), and the substitution at the ortho position of the phenyl ring
(2c) did not hamper the reaction. A bisallene bearing
the 5-methyl-2-pyridinesulfonyl group provided 2d in
75% yield, indicating that the presence of a potentially coordinating
nitrogen atom did not poison the catalyst. Sulfonamide tethers with
aliphatic substitution (tert-butyl and trimethylsilylethyl)
were also efficient, delivering spiro derivatives 2e and 2f in 68 and 55% yields, respectively. Tethers other than
sulfonamides were also tested. Carbon-tethered bisallene 1g and N-Boc bisallene 1h also participate
in the cascade process, affording 2g and 2h with 80 and 42% yields, respectively.
Figure 4
Scope of the cascade
process. Byproducts 3a–3f (5–13%
yields) and 3g (20% yield) were
observed in the 1H NMR of the reaction crude.
Scope of the cascade
process. Byproducts 3a–3f (5–13%
yields) and 3g (20% yield) were
observed in the 1H NMR of the reaction crude.In addition, 1 mmol-scale reactions using NTs-tethered bisallene 1a and N-Boc-tethered bisallene 1h were performed, affording 84 and 42% yields of 2a and 2h, respectively. In all the reactions, byproducts 3 were formed in low yields. Unfortunately, when we attempted heterodimerization
reactions, a complex mixture was observed by NMR.To gain further
understanding on the chemo- and regioselectivity
on the formation of cycloadducts 2, we completed our
study by performing DFT calculations of the DA reaction of the cHT-1g (for the formation of cycloheptatriene cHT, see ref (17)). The Gibbs energy profile computed at 313.15 K and 1 atm with the
(U)B3LYP-D3/cc-pVTZ/SMD(76% THF and 24% CH2Cl2)[20]//(U)B3LYP-D3/cc-pVDZ method is depicted
in Figure ; the Gibbs
energy barriers and the molecular structures of all TSs can be found
in Figures S7 and S8, respectively, and
the closed-shell and open-shell singlet/triplet energies are in Table S1.
Figure 5
(a) Gibbs energy profile in kcal·mol–1 for
the transformation of cHT (from 1g) into 2g. (b) Molecular structures of TS and TS. (c) Values of the condensed Fukui functions f and f (units are electrons)
on the HOMO and LUMO orbitals.
(a) Gibbs energy profile in kcal·mol–1 for
the transformation of cHT (from 1g) into 2g. (b) Molecular structures of TS and TS. (c) Values of the condensed Fukui functions f and f (units are electrons)
on the HOMO and LUMO orbitals.We first evaluated the formation of all the Diels–Alder
cycloadducts shown in Figure , considering both the endo and exo approximations for each
of them (see Figure S7 in the SI) using
closed-shell (cs) calculations. We found that the
lowest energy transition states are those that lead to the formation
of our experimentally observed product, TS and TS, having barriers of 23.9 and 26.6 kcal·mol–1, respectively (see Figure S7 for detailed
Gibbs energy barriers of all reaction paths). The geometries of these
transition states present a high asynchronous character as shown in Figure b, suggesting the
possibility of a two-step process. However, these transition states
lead to product 2g releasing 10.8 kcal·mol–1, and an intermediate could not be localized with the closed-shell
calculations. In contrast to this, when performing open-shell (os) calculations, the reaction was found to occur in a two-step
manner through a biradical intermediate. The transition states for
the 2A and 2A approximations are slightly lower
in energy than their closed-shell counterparts, with TS (⟨S⟩2 = 0.05) being the
lowest one by 0.4 kcal·mol–1 and preferred
over TS (⟨S⟩2 =
0.19) by 2.1 kcal·mol–1. This step generates int(2A in an endergonic process
by 13.1 kcal·mol–1. High delocalization of
the unpaired electrons α and β of int(2A over the conjugated system has been
observed for these intermediates (Figure S9). From this point, the collapse of the biradical intermediate int(2A leads to 2g through TS, surpassing a Gibbs energy barrier of
7.3 kcal·mol–1 and releasing 23.9 kcal·mol–1. Although the two-step biradical pathway through
the endo approximation is the lowest in energy, the difference of
0.4 kcal·mol–1 as compared to the concerted
asynchronous mechanism for the endo approach is not sufficient to
be able to distinguish between the two pathways and they are probably
competing to generate the final product 2g. However,
it should be noted that the biradical pathway clearly helps in rationalizing
the regio- and chemoselectivity observed, which is clearly governed
by the stability of the biradical intermediate and has its two radicals
delocalized in a double allylic position.The condensed Fukui
functions[21] calculated
from the natural population analysis (NPA) charges (f4 = 0.266; f4 = 0.287, Figure c) support the formation of the first bond between
the unsubstituted terminus of the doubly conjugated exocyclic double
bond of the two cHT units (C4–C4, atom labels
in Figure a), leading
to the selective formation of cycloadduct 2g (Figure ). The same conclusion
can be extracted by looking at the HOMO and LUMO orbitals, as the
best molecular orbital overlap is the one that comes from the coupling
of the same methylene positions.To obtain experimental evidence
of the proposed two-step biradical
process, we performed the reaction in the presence of 1.5 equiv of
TEMPO as a radical trapping agent.[22] The
reaction yield toward the formation of 2a was drastically
reduced from 89 to 44%, although the yield of 3a was
not significantly reduced. Furthermore, a fraction could be isolated
from the reaction mixture that, when analyzed by ESI-MS, showed the
incorporation into int(2A) of either two TEMPO moieties
(giving rise to an adduct at m/z = 863.5 for [int(2A) + 2TEMPO + H]+ that
fragmented to a species at m/z =
706.3 by the loss of neutral TEMPO-H as confirmed by MS/MS) or a hydrogen
radical and a TEMPO (giving rise to an adduct at m/z = 708.3 for [int(2A) + TEMPO+H·
+ H]+) (Scheme ). This experiment gives experimental evidence of the intermediacy
of biradical intermediates in the reaction under study.
Scheme 2
Species
Detected by ESI(+)-MS when the Reaction Is Run in the Presence
of TEMPO
Conclusions
In
summary, a method that enables rapid access to seven- and six-membered
spirocyclic compounds in a complete chemo- and regioselective manner
from 1,5-bisallenes has been developed. The whole process involves
an initial rhodium-catalyzed cycloisomerization of 1,5-bisallenes
leading to the non-isolable cycloheptatrienes followed by Diels–Alder
homodimerization. A set of aromatic and aliphatic sulfonamide-tethered
1,5-bisallenes, as well as a carbon-tethered 1,5-bisallene, gave the
final spirocyclic product with complete chemo- and regioselectivity.
The DFT calculations have demonstrated the selectivity observed, arising
from the highly favored homocoupling of the unsubstituted terminus
of the doubly conjugated double bond. Additionally, although experimental
evidence has been found for the formation of the proposed biradical
intermediate, the computational study suggests that both mechanisms,
the two-step biradical and the concerted asynchronous, compete to
generate the spirocyclic products.
Experimental
Section
General Information
Unless otherwise noted, materials
were obtained from commercial suppliers and used without further purification.
Bisallenes 1a (X = NTs), 1b (X = p-MeO-PhSO2N), 1c (X = o-CF3-PhSO2N), 1d (X = 5-Me-2-PySO2N), 1e (X = Bu-SO2N), 1f (X = TMS-CH2CH2-SO2N), 1g (X = C(COOEt)2), and 1h (X = N-Boc) were prepared from the corresponding
bisalkynes using Crabbé homologation reaction. Experimental
procedures and full characterization have been described previously
by us.[16f]CH2Cl2 and THF were dried under nitrogen by passing through solvent purification
columns (MBraun, SPS-800). Reaction progress during the preparation
of all compounds was monitored using thin-layer chromatography on
Macherey-Nagel Xtra SIL G/UV254 silica gel plates. Solvents were removed
under reduced pressure with a rotary evaporator. Reaction mixtures
were chromatographed on silica gel using an automated purification
instrument Interchim PuriFlash XS 520 Plus equipped with a quaternary
gradient pump (up to 300 mL/min, 20 bar) and a UV–vis 200–800
nm diode array detector. All 1H and 13C NMR
spectra were recorded on a Bruker ASCEND 400 spectrometer equipped
with a 5 mm BBFO probe using CDCl3 as a deuterated solvent.
Chemical shifts for 1H and 13C NMR are reported
in ppm (δ) relative to residual solvent signals (7.26 ppm for
1H and 77.16 ppm for 13C). Coupling constants are given
in hertz (Hz). 1H and 13C NMR signals were assigned
based on 2D-NMR HSQC, HMBC, COSY, and NOESY experiments. Electrospray
mass spectrometry analyses were recorded on an Esquire 6000 ion trap
mass spectrometer (Bruker) equipped with an electrospray ion source.
Electrospray ionization high-resolution mass spectrometry was performed
using a Bruker microTOF-Q II instrument. Both mass instruments were
operated in the positive ESI(+) ion mode. IR spectra were recorded
on an Agilent Cary 630 FT-IR spectrometer equipped with an ATR sampling
accessory. Melting points were measured in an SMP10 apparatus from
Stuart without any correction.
Computational Details
Geometries of all stationary
points were optimized without symmetry constraint with the Gaussian
16 program[23] using the DFT B3LYP hybrid
exchange-correlation functional[24] employing
the all-electron cc-pVDZ basis set,[25] and TSexo) and TSendo) were optimized
again using UB3LYP/cc-pVDZ. The electronic energies were improved
(singlet and triplet states) by performing single-point energy calculations
with the cc-pVTZ basis set and the (U)B3LYP hybrid exchange-correlation
functional including solvent effect corrections computed with the
solvent model based on density (SMD) continuum solvation.[25] To mimic the experimental solvent mixture with
a molar fraction ratio of 76:24 of THF:CH2Cl2, the values of the solvent descriptors used in the SMD solvation
model were redefined on the basis of a linear behavior with the molar
fraction. Using the “Solvent = (Generic,Read)” options
of the Gaussian16 SCRF keyword, the solvent mixture was defined employing
the following solvent descriptors: dynamic dielectric constant = 1.410;
static dielectric constant = 1.4085; Abraham’s hydrogen bond
acidity = 0.02424; Abraham’s hydrogen bond basicity = 0.3758;
surface tension = 39.3697; carbon aromaticity = 0; electronegativity
halogenicity = 0.1617. The D3 Grimme energy corrections for dispersion[26] with the original damping function were added
in all (U)B3LYP/cc-pVDZ and (U)B3LYP/cc-pVTZ calculations. Analytical
Hessians were computed to determine the nature of stationary points
(one and zero imaginary frequencies for TSs and minima, respectively)
and to calculate unscaled zero-point energies as well as thermal corrections
and entropy effects using the standard statistical-mechanics relationships
for an ideal gas.[27] These two latter terms
were computed at 313.15 K and 1 atm to provide the reported relative
Gibbs energies. As a summary, the reported Gibbs energies contain
electronic energies including solvent effects calculated at the (U)B3LYP-D3/cc-pVTZ//(U)B3LYP-D3/cc-pVDZ
level together with gas-phase thermal and entropic contributions computed
at 313.15 K and 1 atm with the (U)B3LYP-D3/cc-pVDZ method. All stationary
points were unambiguously confirmed by IRC calculations. For the condensed
Fukui functions, the natural charges were obtained by performing the
natural population analysis of the neutral, cationic, and anionic cHT at the (U)B3LYP-D3/cc-pVTZ//(U)B3LYP-D3/cc-pVDZ theory
level.
General Procedure for the Synthesis of Spiro Derivatives 2
In a 10 mL capped vial, a mixture of [Rh(cod)Cl]2 (2.2 mg, 0.004 mmol, 0.05 equiv), (R)-DTBM-Segphos
(11.3 mg, 0.01 mmol, 0.10 equiv), and NaBArF (8.5 mg, 0.01 mmol, 0.10
equiv) was purged with nitrogen and dissolved in anhydrous CH2Cl2 (4 mL). Hydrogen gas was bubbled into the catalyst
solution, and the mixture was stirred for 30 min. The mixture was
then concentrated to dryness under a stream of hydrogen, dissolved
again in degassed CH2Cl2 (0.5 mL), and transferred
via a syringe into a solution of bisallene 1a–g (0.09 mmol, 1 equiv) in degassed THF (2 mL) under an inert
atmosphere at 40 °C (aluminum heating block). The resulting mixture
was stirred for 16 h at 40 °C. The solvent was removed under
reduced pressure, and the crude reaction mixture was purified by column
chromatography on silica gel using mixtures of hexane/EtOAc as the
eluent (90:10 to 60:40 v/v).
Authors: Bryan K S Yeung; Bin Zou; Matthias Rottmann; Suresh B Lakshminarayana; Shi Hua Ang; Seh Yong Leong; Jocelyn Tan; Josephine Wong; Sonja Keller-Maerki; Christoph Fischli; Anne Goh; Esther K Schmitt; Philipp Krastel; Eric Francotte; Kelli Kuhen; David Plouffe; Kerstin Henson; Trixie Wagner; Elizabeth A Winzeler; Frank Petersen; Reto Brun; Veronique Dartois; Thierry T Diagana; Thomas H Keller Journal: J Med Chem Date: 2010-07-22 Impact factor: 7.446
Authors: Jacob S Bestwick; David J Jones; Helen E Jones; Panagiotis G Kalomenopoulos; Rafal Szabla; Andrew L Lawrence Journal: Angew Chem Int Ed Engl Date: 2022-06-28 Impact factor: 16.823