tions were performed with the Gaussian03 package.22 Long-
range dispersion interactions were also included through an
energy refinement step with the double hybrid B2PLYP
functional23 in conjunction with a basis set of triple-ú quality
with polarization functions (TZVP)24 for all of the atoms.25
The latter single-point energy refinement was performed with
the Orca package.26
does take place, but the reaction mechanism detours con-
siderably from the direct attack of the carbamate nitrogen
to C2 (see Scheme 1). The electrophilic selenium activates
the indole ring in a first step that occurs with neither facial
nor positional selectivity, in clear contrast with the textbook
hypothesis summoned initially (vide supra). For the very
unsymmetrical indole C2-C3 bond, the activation occurs
by MeSe+ attacking either the C2 or the C3 positions, with
formation of a benzyl carbenium ion or an iminium ion,
respectively (see A1-A2 and B1-B2 species in Scheme 1).
These four species are dynamically exchanging in fast
equilibrium processes. Actually, it is remarkable that only
the benzyl carbenium ions are active regarding the seleno-
cyclization reaction, whereas A2 and B2 are unproductive
species. Since the carbamate nitrogen captures the benzyl
carbenium ions, the attack occurs on C3 instead of C2, thus
forming a spirocyclic intermediate (Int-1A and Int-1B).
Subsequently, a concerted double rearrangement of the N
and the Se atoms with ring expansion of the azetidine to the
final fused pyrrol unit was located.27
In keeping with the textbook explanation, the first step in
this mechanism is seemingly the activation of the C2-C3
bond of R by the electrophilic species, with formation of a
charged three-membered ring that is subsequently opened
by the nucleophile. For the case at hand, the activation step
is likely not diastereoselective since the steric cluttering upon
selenonium ion formation is similar on both faces of the
indole side chain due to its considerable flexibility. The good
diastereomeric ratios (dr) of the selenocyclization reactions
(Table 1) argue against the activation step being rate-limiting
since low dr should be expected in that case. The mechanistic
proposal by Danishefsky involves the attack of the carbamate
N lone pair to the C2 position of the electrophile-activated
intermediate. All attempts to locate this transition structure
and its associated cyclic intermediate were unsuccessful.
These calculations, in turn, suggest that the carbamate lone
pair is very compromised in its conjugation with the carbonyl
group and lacks the nucleophilicity required for an effective
attack. Since these cyclization reactions are performed under
acid catalysis (pTSA or PPTS), we instead considered likely
the involvement of the carbamate keto-enol equilibrium,
thus leaving the N lone pair in the enol form isolated from
the carbonyl group and enhancing its nucleophilic character.
In fact, when this tautomer is proposed as the reactive species
of the tryptophan model system R, the nucleophilic attack
The formation of the azetidine is the rate-limiting step with
an activation barrier of ca. 23 kcal/mol (see Table 2). There
Table 2. Relative Free Energies in Solution (CH2Cl2) for the
Seleno- and Bromocyclization Mechanismsa
selenocyclization
bromocyclization
exo
A1
A2
TS-1A
Int-1A
TS-2A
P-1A
0.00
1.61
23.14
11.56
19.01
-5.83
0.00
-2.30
15.69
-0.41
18.29
-15.11
(21) (a) Tomasi, J.; Persico, M. Chem. ReV. 1994, 94, 2027. (b) Mineva,
T.; Russo, N.; Sicilia, E. J. Comput. Chem. 1998, 19, 290.
endo
B1
B2
TS-1B
Int-1B
TS-2B
P-1B
-3.81
-4.12
24.16
11.86
20.33
-6.37
-10.13
-3.10
16.17
0.22
19.76
-15.74
(22) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb,
M. A.; Cheeseman, J. R.; Montgomery, J. A., Jr.; Vreven, T.; Kudin, K.
N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.;
Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.;
Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.;
Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li,
X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.;
Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.;
Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.;
Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich,
S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A.
D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A.
G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.;
Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham,
M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.;
Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian
03, revision C.02; Gaussian, Inc.: Wallingford, CT, 2004.
(23) Grimme and co-workers have demonstrated that intramolecular long-
range dispersion is very important to obtain an accurate electronic description
of molecules in general and oligopeptides in particular (Schwabe, T.;
Grimme, S. Phys. Chem. Chem. Phys. 2007, 9, 3397). They introduced
long-range dispersion by perturbative methods (Grimme, S. J. Chem. Phys.
2006, 124, 034108) and developed a semiempirical scheme to further
improve the latter functionals (Grimme, S. J. Comput. Chem. 2004, 25,
1463). In particular, the double hybrid B2PLYP functional with dispersion
corrections included semiempirically has proved to yield errors considerably
smaller than those of the venerable B3LYP for the G3 set.
a Relative free energies in kcal/mol.
is about 1 kcal/mol difference in activation energy favoring
TS-1A over TS-1B, which we surmise is the source of the
final exo/endo selectivity. Danishefsky’s finding that the exo
adduct is the kinetic product in the rearrangement of 3 (Table
1) is in good agreement with our calculations.
After thorough inspection of these two transition structures,
we concluded that there are no obvious structural features
responsible for the difference in activation energies. Both
transition structures exhibit similar geometries (bond lengths,
angles, and dihedrals); in fact, the energy difference between
both transition structures is only ∼0.5 kcal/mol in the gas
phase, which means that half of the energy difference
(24) Schaefer, A.; Horn, H.; Ahlrichs, R. J. Chem. Phys. 1992, 97, 2571.
(25) To make the perturbative part more affordable given the large basis
functions employed, the resolution of the identity (RI) was used as
implemented in the Orca program: (a) Weigend, F.; Ha¨ser, M. Theor. Chem.
Acc. 1997, 97, 331. (b) Weigend, F.; Ko¨hn, A.; Ha¨ttig, C. J. Chem. Phys.
2002, 116, 3175.
(27) Similar ring expansions on indole motifs can be found in recent
literature: (a) Li, C.; Chan, C.; Heimann, A. C.; Danishefsky, S. J. Angew.
Chem., Int. Ed. 2007, 46, 1444. (b) Baran, P. S.; Shenvi, R. A. J. Am.
Chem. Soc. 2006, 128, 14028.
Org. Lett., Vol. 10, No. 1, 2008
79