Proton transfers in the Strecker reaction revealed by DFT calculations.

The Strecker reaction of acetaldehyde, NH3, and HCN to afford alanine was studied by DFT calculations for the first time, which involves two reaction stages. In the first reaction stage, the aminonitrile was formed. The rate-determining step is the deprotonation of the NH3 (+) group in MeCH(OH)-NH3 (+) to form 1-aminoethanol, which occurs with an activation energy barrier (ΔE (≠)) of 9.6 kcal/mol. The stereochemistry (R or S) of the aminonitrile product is determined at the NH3 addition step to the carbonyl carbon of the aldehyde. While the addition of CN(-) to the carbon atom of the protonated imine 7 appears to scramble the stereochemistry, the water cluster above the imine plane reinforces the CN(-) to attack the imine group below the plane. The enforcement hinders the scrambling. In the second stage, the aminonitrile transforms to alanine, where an amide Me-CH(NH2)-C(=O)-NH2 is the key intermediate. The rate-determining step is the hydrolysis of the cyano group of N(amino)-protonated aminonitrile which occurs with an ΔE (≠) value of 34.7 kcal/mol. In the Strecker reaction, the proton transfer along the hydrogen bonds plays a crucial role.

Introduction

In 1850, Adolph Strecker reported a reaction that affords alanine from n class="Chemical">acetaldehyde, ammonia and hydrogen cyanide [1]. The original form of Strecker amino acid synthesis is shown in Scheme 1(a). In this reaction, the aldehyde reacts with hydrogen cyanide to form an aminonitrile, which undergoes hydrolysis to afford alanine in the acidic solution. The traditional Strecker reaction gave racemic α-aminonitriles (mixtures of equal amounts of R and S forms), where an imine RCH=NH was considered to be the key intermediate [2]. Three typical reactions are presented in Scheme 1(b) [3].

Scheme 1

The general form of the Strecker reaction. The reaction (b) is taken from [2].

In 1963, Harada reported the first asymmetric Strecker reaction, in which an (S)-α-phenylethylamine was employed as the chiral auxiliary [4]. In this reaction, he obtained a chiral n class="Chemical">alanine with 95% optically activity; see Scheme 2. In 1996, Lipton et al. succeeded in a series of asymmetric Strecker reactions by employing a chiral catalyst, a cyclic dipeptide [5]. In these reactions, N-substituted imines react with HCN to yield (S)-α-aminonitriles with remarkably high enantiomeric excess (ee). One example is shown in Scheme 3.

Scheme 2

The first asymmetric Strecker reaction [4].

Scheme 3

The first asymmetric synthesis of α-aminonitirles via a chiral catalyst [5].

However, when benzaldehyde and n class="Chemical">NH3 instead of the N-substituted imine were employed as the substrates, the reaction afforded an initial product Ph-CH(NH2)-CN of configurational instability[5]. In the following, Sigman and Jacobsen used a parallel combinatorial library synthesis for the discovery and optimization of a chiral catalyst for the reaction of imines and HCN [6]. From then on, various catalytic asymmetric Strecker reactions have been reported to gain high enantioselectivity of the hydrocyanation reaction of imines [7-12]. However, the origin of the enantioselectivity in the asymmetric Strecker reactions has not been clarified.

To our knowledge, the elementary processes of the whole Strecker reaction have not been elucidated. As shown in Scheme 1, the Strecker reaction includes two reaction stages. The first reaction stage is the condensation of aldehydes with n class="Chemical">ammonia and hydrogen cyanide leading to α-aminonitriles . The second reaction stage is the hydrolysis of the nitrile group. In these reactions, K+ (or Na+) and Cl− ions are not involved, as shown in Scheme 1(a). Therefore, it is suitable to the theoretical investigation of the reaction mechanism, because the effect of counter ions does not need to be considered.

Actually, several theoretical studies were reported of the last step of the first reaction stage of the Strecker reaction [10-16], i.e. the hydrocyanation of imines (or protonated n class="Chemical">imines + CN−) to aminonitriles. In those works, how the nucleophile CN− is generated has not been examined. Because HCN is a very weak acid with a dissociation constant of Ka = 1.3 × 10−9 mol/L (in water, 18 °C), direct dissociation reaction HCN → H+ + CN− is difficult to occur.

In the second reaction stage, i.e., the acid-catalyzed hydrolysis of the cyano group, protonation of the group appears to cause the addition of OH2 to the cyano n class="Chemical">carbon:

R-CH(NH2)-CN + H+ → R-n class="Chemical">CH(NH2)-CNH+

R-CH(NH2)-CNH+ + n class="Chemical">OH2 → R-CH(NH2)-C(OH)=NH + H+

However, the proton affinity (PA) of the nitrile is much smaller than that of the amino group, for example, the n class="Chemical">PAs of the cyano and amino groups of 2-amino-propanonitrile (Me-CH(NH2)-CN) are 190.7 and 199.6 kcal/mol, respectively. Thus, in the acidic solution (2H2O + 2HCl), Me-CH(NH3+)-CN should be afforded; see Scheme 1(b). The reaction mechanism of this hydrolysis is also unclear.

To address the above issues, we performed DFT calculations of the Strecker reaction shown in Scheme 1(a). Here, ten specific water molecules were considered.

Methods of calculation

Geometry optimizations were performed by density functional theory (DFT) with the B3LYP [17-18] functional. The basis set 6-311+G(d,p) was employed for all the atoms in the calculations. The solution (water) effect was considered by the Polarizable Continuum Model (PCM) [19-21]. Vibrational analyses were carried out to make sure whether a stationary point is an equilibrium structure or a transition state (TS). From TSs, reaction paths were traced by the intrinsic reaction n class="Chemical">coordinate (IRC) method [22-23] to obtain the energy-minimum geometries. All the calculations were carried out using the GAUSSIAN 09 [24] program package. Throughout this paper, the discussion was presented based on the potential energy changes with zero-point vibrational energy (ZPE) correction unless otherwise noted.

Results and Discussion

Formation reaction of aminonitrile (the first stage)

The reaction model is shown in Scheme 4. In the model, lone-pair electrons of the oxygen and n class="Chemical">nitrogen atoms participate in hydrogen bonds. In calculating each TS, ten water molecules were placed so that the large hydrogen-bond stabilization is gained. After geometries of TSs were determined, those of energy minima were obtained by IRC and the subsequent optimizations. By the use of the similar water cluster models, TS geometries and activation energies (Ea's) in the base promoted ester hydrolyses were calculated [25]. In Ph-COOEt + OH−(H2O)→Ph-COO− + HO-Et + (H2O), Ea = +14.7 kcal/mol (n = 5), +16.3 (n = 8), +16.3 (n = 12) and +15.6 (n = 16) were obtained, where the experimental Ea is + 14.6 kcal/mol. Also, in the Bamberger rearrangement Ph-NH(OH) + (H3O+)2(H2O)13 → para-HO-C6H4-NH3+ + H3O+(H2O)14 [26], Ea = +26.3 kcal/mol was calculated, where the experimental Ea is +24.8 kcal/mol.

Scheme 4

A reaction model composed of Me-CH=O, HCN, NH3 and (H2O)10 for geometry optimizations to trace elementary processes. Broken lines stand for hydrogen bonds.

A reaction model composed of Me-CH=O, HCN, n class="Chemical">NH3 and (H2O)10 for geometry optimizations to trace elementary processes. Broken lines stand for hydrogen bonds.

Our proposed reaction pathways are shown in Scheme 5. Geometries of TSs are shown in Figure 1 and those of precursor 1, intermediates and product 8 are provided in Supporting Information File 1 Figure S1.

Scheme 5

Possible pathways for the formation of aminonitrile from acetaldehyde.

Figure 1

Geometries of transition states along the reaction from acetaldehyde (1) to the aminonitrile 8. Distances are in Å. TS means, for instance, a transition state for the step 1 → 2.

Possible pathways for the formation of aminonitrile from n class="Chemical">acetaldehyde.

Geometries of transition states along the reaction from acetaldehyde (1) to the n class="Chemical">aminonitrile 8. Distances are in Å. TS means, for instance, a transition state for the step 1 → 2.

As shown in Figure S1, Supporting Information File 1, MeCH=O, n class="Chemical">NH3 and HCN are separated by the H2O cluster in the precursor complex 1. The reaction begins with the addition of NH3 to the carbonyl carbon of acetaldehyde to form a Mulliken charge-transfer complex 2. This complex was firstly proposed here. In 2, the C–O bond is elongated to1.345 Å, which has an alkoxide character and the complex is not stable in the gas phase. However, it is more stable than the precursor complex by 3.5 kcal/mol when ten specific water molecules are considered; see Figure 2. This result indicates that the consideration of water molecules in the reaction is necessary to describe the step 1 → 2. Then, the alkoxide oxygen atom captures a proton from a surrounding water molecule to form Me(H)C(OH)-NH3+ and a remaining OH− anion in the surrounding. They form an ion pair 3 [Me(H)C(OH)-NH3+ and OH−]. Next, the hydroxide ion catches a proton from HCN through the transition state TS to afford a more stable ion-pair intermediate 4 [Me(H)C(OH)-NH3+ and CN−]. This step is exothermic by 5.5 kcal/mol with a small activation energy barrier (ΔE≠) of 1.4 kcal/mol. Starting from 4, a proton of the NH3+ group migrates to one water molecule to form 5 [Me(H)C(OH)-NH2(H3O+) and CN−] via a transition state TS. After that, the proton migrates from H3O+ to the hydroxy group of Me(H)C(OH)-NH2 via TS to yield Me(H)C(OH2+)-NH2 6. From 6, H2O is easily eliminated through TS to afford the protonated imine 7 with a ΔE≠ value of 0.3 kcal/mol. At last, CN− nucleophilically attacks the carbon atom of MeCH=NH2+ through TS to afford a 2-aminopropanonitrile 8.

Figure 2

Energy changes along elementary processes from acetaldehyde to aminonitrile. Bold numbers are defined in Scheme 5.

Energy changes along elementary processes from acetaldehyde to n class="Chemical">aminonitrile. Bold numbers are defined in Scheme 5.

Starting from 4, a concerted SN2-type pathway was also examined, which directly leads to the n class="Chemical">nitrile compound 8; see Scheme 6. In this pathway, the proton is transferred from the NH3+ group to the hydroxy group via a two-water-molecule bridge. At the same time, the H2O elimination and the approach of CN− to 4 concomitantly take place with a Walden inversion. The transition state TS was successfully located; see Figure 1. However, this pathway needs a large energy barrier of 58.8 = [+51.3 − (−7.5)] kcal/mol, indicating that it is difficult to occur.

Scheme 6

A short-cut path by the nucleophilic displacement and the concomitant proton transfer. “The first bypass” in Scheme 5.

Starting from the ion-pair intermediate 3, we also investigated the possibility that the hydroxide ion captures a proton from the n class="Chemical">NH3+ group to form 9, Me-C(H)(OH)-NH3+ + OH− → Me-C(H)(OH)-NH2 + OH2. This reaction step occurs through a transition state TS with a ΔE≠ value of 1.9 kcal/mol, which is comparable to that (1.4 kcal/mol) of the proton transfer step from HCN to OH−. However, 9 is less stable than 4 by 1.8 kcal/mol, indicating that the OH− prefers to capture a proton from HCN rather than from the NH3+ group.

As shown in Figure 2, the rate-determining step of this reaction stage is the proton migration from the NH3+ group to the n class="Chemical">water cluster (from 4 to 5), where the energy barrier is 17.1= [+9.6 − (−7.5)] kcal/mol. Other proton transfer steps facilely occur. These results are in consistent with the room-temperature experimental condition in Scheme 1. For the rate-determining step, we also checked an extended model "TS–ext", where ten water molecules are added (the molecular formula, C3H48N2O21). The geometry of TS–ext is shown in Supporting Information File 1 Figure S2. The geometrical parameters of the proton-transfer region of TS–ext are similar to that of TS. Also, the energy difference between TS–ext and 4–ext (17.3 kcal/mol) is very close to that (17.1 kcal/mol) between TS and 4.

As shown in Figure 1, the product, aminonitrile 8, is in an S form. However, racemic α-n class="Chemical">aminonitriles are obtained experimentally. This stereochemical scrambling is explicable on the basis of the computational results, as follows. In TS, the NH3 molecule may add to MeCH=O from both upper and lower directions equivalently, which leads to the racemic products. However, in TS, the nucleophile CN− is obligated to attack MeC(H)=NH2+ at the plane opposite to the OH2 dissociating side (see Scheme 7). The addition model of TS was examined by the use of the amine in Scheme 2. The activation energy of the less crowded TS–A is 1.8 kcal/mol smaller than the more crowded TS–B; see Figure 3. This calculation result is consistent with Harada's work that a chiral product [4] was obtained in the Strecker reaction.

Scheme 7

A contrast of the nucleophilic addition.

Figure 3

Two transition states (A and B) of the nucleophilic addition of (S)-α-phenylethylamine to acetaldehyde. (H2O)10 is also included, and the molecular formula of the reaction system is C11H36N2O11.

Two transition states (A and B) of the nucleophilic addition of (S)-α-phenylethylamine to n class="Chemical">acetaldehyde. (H2O)10 is also included, and the molecular formula of the reaction system is C11H36N2O11.

Hydrolysis of amino nitrile to amino acid (the second stage)

In the acidic hydrolysis of amino nitrile, we take n class="Chemical">2-amino-propanonitrile +H3O+(H2O)10 (8) as a precursor complex; see Figure S3 in Supporting Information File 1 for the geometry of 8. Our proposed pathways are shown in Scheme 8.

Scheme 8

Elementary processes of the acid-catalyzed hydrolysis of 2-amino-propanonitrile.

There are two competitive pathways from the precursor complex: One (I) is the protonation of the amino group to form a N(on amino)-protonated n class="Chemical">aminonitrile 10. This step occurs through TS with the ΔE≠ and ΔE values of 4.1 and −6.8 kcal/mol, respectively; see Figure 4.

Figure 4

Energy changes along elementary processes from 2-amino nitrile 8 to 2-amino acid 16. Brown-color lines stand for the most favorable route.

Energy changes along elementary processes from 2-amino nitrile 8 to n class="Chemical">2-amino acid 16. Brown-color lines stand for the most favorable route.

The other pathway (II) is the hydrolysis of the C≡N group to form a compound 18 through TS, where the OH group is added to the carbon atom and the n class="Chemical">hydrogen atom attachs to the nitrogen atom. This step needs a considerably large energy barrier of 32.0 kcal/mol with an endothermicity of 0.2 kcal/mol. Obviously, the protonation of the amino group (I) is much more favorable than the hydrolysis of the cyano group (II). As a result, compound 10 is the starting point for the following reactions. From 10, a water trimer reacts with the cyano group through TS to afford the N(on amino)-protonated 2-amino-1-hydroxypropanimine 11. At TS, the hydroxy group adds to the carbon atom of the cyano carbon nucleophilically. Simultaneously, the proton migrates to the nitrogen atom through a two-water-bridge; see Figure 5 for the geometry of TS. The C≡N group in 10 convers to a C(OH)=NH group in 11. This reaction step is endothermic by 4.6 = [−2.2 − (−6.8)] kcal/mol with a large ΔE≠ value of 34.7 = [+27.9 − (−6.8)] kcal/mol. Although this ΔE≠ value is apparently larger than that of TS, TS lies lower than TS by 4.1 kcal/mol when taking the energy of the precursor complex as a reference; see Figure 4. We examined the role of the NH3+ group in the reaction by investigating the hydrolysis of the cyano group of a methyl-substituted model 10(Me). In this model, we replaced the NH3+ group in 10 with a methyl group. The ΔE≠ value of the hydrolysis step increases to 34.6 kcal/mol. It indicates that the NH3+ group enhances the electrophilicity of the cyano carbon, which is favorable for the OH2 addition. The TS geometry of the methyl substituted model, TS(Me), is shown in Supporting Information File 1 Figure S6.

Figure 5

Geometries of transition states along the most favorable route from 2-aminonitrile 8 to 2-amino acid 16.

Geometries of transition states along the most favorable route from 2-aminonitrile 8 to n class="Chemical">2-amino acid 16.

In the following, a proton on the NH3+ group in 11 is transferred to the n class="Chemical">imine nitrogen to afford 12. This proton transfer step is facilitated by a water molecule bridge. The ΔE≠ and ΔE values are 5.5 = [+3.3 − (−2.2)] and −3.8 = [−6.0 −(2.2)] kcal/mol, respectively. Starting from 12, there are two possible pathways, paths A and B, to form the final alanine product; see Scheme 8. In path A, the deprotonation of the amino group occurs to produce an amide intermediate 13 with an exothermicity of −8.3 = [−14.3 − (−6.0)] kcal/mol; see Figure 4. In the following, the protonation of the amide nitrogen atom occurs to produce a cationic species MeC(NH2)H-C(=O)-NH3+ 14. The amide carbon in 14 is subject to the OH2 addition to afford a zwitterion compound MeC(NH2)H-C(OH)(O−)-NH3+ 15. The ΔE≠ and ΔE values of the H2O addition step relative to the energy of 13 are 30.5 and 25.8 kcal/mol, respectively. After that, the NH3 moiety is ready to dissociate from 15 to produce the product (R)-alanine 16 with a ΔE≠ value of 0.7 kcal/mol of TS. In path B, the second H2O molecule is added to the C=N double bond in 12 to form an intermediate 17 through TS. The ΔE≠ and ΔE values of this step are 25.5 and 8.6 kcal/mol, respectively. However, TS lies higher than TS by 3.1 kcal/mol and 17 in path B is much more unstable than the intermediate 13 in path A by 17.0 kcal/mol. These energy differences suggest that the deprotonation of the amino group occurs more favorably. Then, from 17 the elimination of the NH4+ group takes place to afford the alanine 16 with a ΔE≠ value of 2.5 kcal/mol. According to the above discussion, the path-B is less favorable than the path A.

The most favorable pathway for the second stage of the Strecker reaction was shown in brown color in Figure 4. The rate-determining step is the OH2 addition to the cyano group (TS) with the activation energy of 27.9 kcal/mol. The competitive TS has an energy barrier of 32.0 kcal/mol. Geometries along the unfavorable routes pathway II and pathway IB, are shown in Figures S4 and S5, Supporting Information File 1, respectively. The relative stability of these two transition states were checked with extended models TS-ext and TS-ext, which have a molecular formula of C3H49N2O21+. Their geometries are shown in Figure S7. TS-ext lies higher than TS-ext by 3.4 kcal/mol, which is consistent with the small model system.

Conclusion

In this work, we theoretically investigated the whole Strecker reaction shown in Scheme 1(a), which includes two reaction stages. The most favorable pathways are summarized in Scheme 9.

Scheme 9

Summary of the present computational work expressed by minimal models.

As shown in the upper half of Scheme 9, the first reaction stage, acetaldehyde + n class="Chemical">NH3 + HCN + (H2O)10 (1) → 2-aminopropanonitrile (H2O)11(8), is composed of seven elementary processes. The rate-determining step is the deprotonation of the NH3+ group in MeCH(OH)-NH3+ to form 1-aminoethanol, which occurs with an activation energy barrier of 9.6 kcal/mol. The stereochemistry (R or S) of the product aminonitrile is determined by equal addition of NH3 to the carbonyl carbon of the aldehyde in both sides. While the addition of CNˉ to the carbon atom of the protonated imine 7 appears to give the scrambling of the stereochemistry, the water cluster above the imine plane reinforces the CNˉ to attack the carbon atom below the plane; see Scheme 7. While HCN is a very weak acid, CN− may be generated by the proton transfer, HCN + OH− → CN− + H2O (3 → TS → 4) in this reaction stage. As shown in the lower half of Scheme 9, the second reaction stage, aminonitrile + H3O+(H2O)10 → alanine + NH4+(H2O)9, is also composed of seven elementary processes. In this reaction stage, the protonation to the amino nitrogen occurs first, which enhances the subsequent hydrolysis of the cyano group to form an imine C(OH)=NH moiety. Its rate-determining step is the hydrolysis of the cyano group of N(amino)-protonated aminonitrile to afford a N(amino)-protonated 1-hydroxypropanimine 11. The ΔE≠ value of this step is 34.7 kcal/mol and the large value corresponds to the high temperature conditions in Scheme 1(b).

Supporting Information

Suppoting Information File 1:

File Format: PDF.

Cartesian coordinates of optimized geometries in Figures 1, 3, and 6 and Figures S1–S7.