Figure 2. Plot of the coefficients.
Table 1. Factors and ranges
In order to understand the response surface better, we plot-
ted the predicted yields as a function of two variables in a
series of contour plots as seen in Figure 3. In the four plots,
yield is plotted as a function of aniline and catalyst stoichi-
ometry. Four combinations of temperature and time are em-
ployed in the four different plots. An analysis of the plots
shows that the strongest dependency of yield is on tempera-
ture, with higher temperatures giving higher yields. Similarly,
higher stoichiometries of both aniline and catalyst give higher
yields. As per the model, in order to improve our yield be-
yond 65% (the highest yield we obtained from the initial set
of experiments), we needed to increase the stoichiometry of
aniline, catalyst, temperature, or various combinations thereof.
Figure 4 shows a contour plot of aniline versus catalyst at
130 °C with an incubation time of 24 h. It can be seen that
the model predicts yields in the 90% range with high equiva-
lents of both aniline and catalyst. This prompted us to modify
the initial ranges that we set for the factors aniline (1-12
equiv) and catalyst (0-9 equiv) while we maintained the
temperature at the boiling point of the solvent (130 °C). The
results of these experiments are shown in Table 3. As seen
from Table 3, at 9 equiv each of aniline and catalyst we
obtained ∼70% yield as predicted by the model. This was
an acceptable yield for our purposes, but such high equiva-
lents of aniline and catalyst were not too appealing since it
would involve additional work to remove excess reagents at
the end of the reaction, particularly during scale-up. We were
then left with temperature and incubation time as the only
remaining factors. Reaction times of >24 h are inconvenient,
thus leaving temperature as the only remaining factor to vary.
The maximum temperature used in the design is indeed
the boiling point of the solvent, 2-ethoxyethanol. To achieve
higher temperatures under conventional heating, we would
need to use a different, higher-boiling solvent. Switching the
solvent at this stage of optimization might cause unforesee-
able variations and was deemed undesirable. Taking these
factors into consideration, we decided to try the reaction
under microwave irradiation conditions. Use of microwave
irradiation in organic synthesis has gained popularity in
recent years, particularly for parallel synthesis. The higher
reaction temperatures that can be attained using microwave
irradiationswell above the boiling point of the solventswas
especially attractive in our situation. Table 4 shows results
factor
conc. units
range
aniline
equiv
equiv
°C
1-4
0-5
50-130
8-24
catalyst
temperature
time
h
conditions to maximize the yield. This approach has been
extensively used in several disciplines3 and is now commonly
being adopted by synthetic chemists as well.4
Results and Discussion
Experimental design, model building, analysis, and op-
timization of reaction conditions were all performed using
the MODDE5 software application. Four different factors6
were considered in our studysstoichiometry of aniline 7,
stoichiometry of pyridine hydrochloride (henceforth referred
to as catalyst), temperature, and reaction time. The initial
ranges of these factors are shown in Table 1. A two-level
full factorial design was chosen to screen the factor space.
A set of 16 experiments plus three replicates of the center
point was carried out to develop a model of the response
surface. The design and the yields from these experiments
are shown in Table 2. The three center point experiments
had a mean value of 26.1 with a standard deviation of 3.5.
The data was used to build a predictive model using multiple
linear regression. We obtained a good fit, with R2 ) 0.95
and Q2 ) 0.79. Figure 2 shows a plot of the coefficients of
the regression model. It can be seen that temperature,
stoichiometry of aniline, stoichiometry of catalyst, and two
cross terms are significant.
(3) (a) Morris, N. D.; Mallouk, T. E. J. Am. Chem. Soc. 2002, 124, 11114-
111121. (b) Le Mapihan, K.; Vial, J.; Jardy, A. J. Chromatogr., A 2004,
1061, 149-158. (c) Arauzo-Bravo, M. J.; Shimizu, K. J. Biotechnol. 2003,
105, 117-133. (d) Halliwell, C. M.; Cass, A. E. Anal. Chem. 2001, 73,
2476-2483.
(4) (a) Stazi, F.; Palmisano, G.; Turconi, M.; Clini, S.; Marco, S. J. Org. Chem.
2004, 69, 1097-1103. (b) Chen, J. J.; Nugent, T. C.; Lu, C. V.; Kondapally,
S.; Giannousis, P.; Wang, Y.; Wilmot, J. T. Org. Process Res. DeV. 2003,
7, 313-317. (c) Gooding, O. W. Curr. Opin. Chem. Biol. 2004, 8, 297-
304.
(6) Preliminary experiments carried out in parallel using the conditions detailed
in method A, but varying the concentration (0.025 M, 0.05 M, 0.1 M and
0.2 M) did not show any change. As a result, a detailed study was not
carried out using concentration as a factor.
Vol. 11, No. 3, 2007 / Organic Process Research & Development
•
451