Cantillo and Kappe
JOCArticle
methods to compute the free energy barriers for the reaction,
is interesting but only compares the energetics of the direct
proton transfer pathways (via a four-membered transition
structure) and the proton 1,3-shift assisted by water, omit-
ting the possibility of the assistance by a second molecule of
aldehyde or other protic species.
SCHEME 2. The DABCO-Catalyzed Morita-Baylis-
Hillman (MBH) Reaction of Aryl Aldehydes (4) with Methyl
Acrylate (2)
In the present work we present a detailed computational
and experimental reinvestigation on the amine-catalyzed
MBH reaction of benzaldehyde with methyl acrylate. In
particular, issues relating to the hitherto neglected tempera-
ture dependence of the thermodynamic properties;which
prevents the use of elevated temperatures as an activation
method in the otherwise sluggish MBH reaction;are ad-
dressed. In addition, we computationally reinvestigate, em-
ploying the recently introduced M06-2X method14 that pro-
vides accurate thermodynamics for carbon-carbon bond
forming reactions,15 all the previously suggested mechanistic
proposals for this transformation, emphasizing the most
controversial step: the proton migration from the R-carbon
to the oxygen derived from the aldehyde counterpart (Figure 1).
To fully address the mechanistic conundrum surrounding the
MBH reaction, we also take into consideration the reagents/
additives reported in the literature to activate this proton
migration step: a second molecule of aldehyde,7 alcohols,16
water,16,17 phenol,18 and the observed autocatalysis.8 Grati-
fyingly, these new calculations can rationalize all existing
experimental evidence, and most importantly, can settle the
argument about the character of the hydrogen migration
step, in particular the Aggarwal versus McQuade proposals
(Figure 1). On the basis of the present theoretical study and a
reinterpretation of the available kinetic data it can be con-
cluded that both mechanisms are competing reactions, and
depending on the reaction progress and conditions either of
the two pathways is favored.
stationary points. Ground and transition states were character-
ized by none and one imaginary frequencies, respectively.
All the presented relative energies are free energies with
respect to the reactants. As in previous studies,7,8,10,13 the
reaction of benzaldehyde with methyl acrylate and DABCO
was chosen as the model for the MBH reaction. In addition,
4-nitrobenzaldehyde and 4-methoxybenzaldehyde were also
employed in order to observe substituent effects in the tempera-
ture dependence studies. Since along the mechanistic pathways
several diastereomeric transition structures can be formed, all
the possibilities have been investigated and are presented in the
Supporting Information.
Results and Discussion
Temperature Dependence of Thermodynamic Properties.
For the present study we have chosen the reaction between
benzaldehyde (4a) and methyl acrylate (2) catalyzed by
diazabicyclo[2.2.2]octane (DABCO) (1) in methanol as the
model system, since this reaction can be considered the
archetypal MBH reaction and is one of the most experimen-
tally studied examples (Scheme 2).1-10 The most relevant
experimental and theoretical mechanistic studies on the
MBH reaction by McQuade7 and Aggarwal8,10 were also
carried out with this system. Additional experimental data
with 4-nitrobenzaldehyde (4b) instead of benzaldehyde have
also been reported.7
Computational Details
The available computational data on this transformation,
with similar energies for the reactants and the product, point
to a reversible character for the MBH reaction.10 However,
in a previous study7 attempts to experimentally observe the
reversibility of this model reaction (the reaction of 4a was
chosen in this study) were unsuccessful. To clarify this issue
for our model system, a purified sample of MBH adduct 6a
was resubjected to the reaction conditions. After dissolving
adduct 6a in methanol in the presence of DABCO (2 equiv)
and heating the mixture at 120 °C for 1 h the formation of
aldehyde substrate 4a in considerable amounts (∼68% by
GC-FID) was observed. Importantly, by letting the reaction
mixture stand at room temperature and monitoring the
progress continuously by GC-FID the equilibrium again
shifted to the product side, with 72% MBH adduct 6a being
observable after 24 h. These simple experimental studies not
only verify the reversible character of the MBH reaction, but
also point to a strong temperature dependence of the equi-
librium constant, in the sense that the adduct (6a) is favored
at low temperatures, while the reactants 4a þ 2 are the major
species at elevated temperatures.
All geometries and energies, as well as frequency calculations,
were computed with use of the Gaussian09 package.19 The
MP220 method along with the 6-311þG(d,p) basis set was
employed for the temperature dependence studies. In the case
of the pathway simulations, the M06-2X14 density-functional
method with the 6-311G(d,p) basis was used.21 The geometries
were optimized including the solvation effect. For this purpose
the SMD22 solvation method was employed, using methanol or
tetrahydrofuran as solvent. The frequency calculations on all
the stationary points were carried out at the same level of theory
as the geometry optimizations to ascertain the nature of the
(14) Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215–241.
(15) Pieniazek, S. N.; Clemente, F. R.; Houk, K. N. Angew. Chem., Int.
Ed. 2008, 47, 7746–7749.
(16) Aggarwal, V. K.; Emme, I. S.; Fulford, Y. J. Org. Chem. 2003, 68,
692–700.
(17) Cai, J.; Zhou, Z.; Zhao, G.; Tang, C. Org. Lett. 2002, 4, 4723–4725.
(b) Aggarwal, V. K.; Dean, D. K.; Mereu, A.; Williams, R. J. Org. Chem.
2002, 67, 510–514.
(18) Faltin, C.; Fleming, E. M.; Connon, S. J. J. Org. Chem. 2004, 69,
6496–6499.
(19) Gaussian 09, Revision A.1; Frisch, M. J. et al. ; Gaussian, Inc.,
Wallingford, CT, 2009.
(20) Head-Gordon, M.; Pople, J. A.; Frisch, M. J. Chem. Phys. Lett. 1988,
153, 503–506.
(21) Single points calculations for key transition structures with the
6-311þG(d,p) basis set were also performed to ensure that the addition of
diffuse functions does not affect the drawn conclusions.
(22) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. J. Phys. Chem. B
2009, 113, 6378–6396.
This temperature dependence can in principle be observed
in computational studies on the MBH reaction. In the
previous theoretical study by Aggarwal and Harvey10 no
frequency analysis of the stationary points is presented, thus
avoiding the chance to determine any temperature influence.
J. Org. Chem. Vol. 75, No. 24, 2010 8617