Quantum-Chemical Modeling of The
Total Page:16
File Type:pdf, Size:1020Kb
Article Cite This: J. Phys. Chem. A 2018, 122, 1643−1657 pubs.acs.org/JPCA Quantum-Chemical Modeling of the First Steps of the Strecker Synthesis: From the Gas-Phase to Water Solvation Published as part of The Journal of Physical Chemistry virtual special issue “Manuel Yań̃ez and Otilia Mó Festschrift”. V. Riffet, G. Frison, and G. Bouchoux* LCM, CNRS, Ecole Polytechnique, UniversitéParis-Saclay, 91128 Palaiseau, France *S Supporting Information ABSTRACT: The two first steps of the Strecker synthesis of glycine, namely (a) addition of ammonia to formaldehyde to give aminomethanol and (b) its dehydration leading to methaneimine, is studied using high level quantum chemistry computations (G3B3). Water solvation is modeled by considering the effect of adding a discrete number of active or passive molecules of water (up to four) and by immersing the identified water-solute complexes in a conductor-like polarizable continuum solvent model. Activation of the reactants by protonation is also examined. Exhaustive search of microhydrated neutral and protonated aminomethanol has been performed using a combination of hierarchical and genealogical approaches. Critical energies associated with all the elementary reaction steps were estimated using the accurate G3B3 composite method thus providing benchmarks to discuss the possible occurrence of Strecker synthesis in prebiotic chemistry. 1. INTRODUCTION byproducts of this synthesis route) in carbonaceous chondrites The Strecker reaction1 is probably the most straightforward has been considered as an indication that the Strecker reaction has been involved in the formation of amino acids in the parent route from simple compounds such as ketone, aldehyde, 9,18,22 ammonia, and hydrogen cyanide to α-amino acids. For this bodies of this family of meteorites. Recently, laboratory − reason, it has been extensively studied2 and reviewed.3 5 experiments on interstellar ice analogues show that amino- Briefly, the Strecker synthesis proceeds through several methanol, HOCH2NH2, (in parallel with hydroxyacetonitrile, ∼ successive steps as depicted in Scheme 1: (a) addition of HOCH2CN) may be produced at low temperature (i.e., ammonia to an aldehyde RCHO giving rise to the 40K) from ices containing formaldehyde and ammonia.23,24 corresponding aminoalcohol RCH(OH)NH2, (b) dehydration Previous mechanistic studies of the Strecker reaction using of the later producing the imine RCHNH, then, (c) reaction of quantum chemical calculations have been devoted to step (a) in − the imine with hydrogen cyanide leading to formation of an α- the case R = H.25 31 Depending upon the theoretical level aminonitrile RCH(NH2)CN. Additionally, in aqueous solution, used, the computed critical energy for the formation of α −1 -aminonitrile formation is accompanied by a series of other H2NCH2OH was situated between 115 kJ mol and 200 kJ processes constituting a complex network of equilibrated mol−1 with respect to the reactants CH O+NH. The second 2 2 3 reactions from which emerge two noticeable irreversible step (b) has been predicted by theory to be more endothermic α − pathways: (d) the hydrolytic transformation leading to - than step (a) by ca. 40−60 kJ mol 1.27,29,30 Consideration of a amino acid, and (e) the reaction with CO2 giving hydantoins ̈ − 6 water solvent environment has been made by means of the through the Bucherer Berg reaction. isodensity surface polarizable model, IPCM, but without Occurrence of the abiotic formation of α-amino acids on the fi optimization of the geometry at this level, and inclusion of early Earth through Strecker synthesis was rst advocated by 26,28 7 one or two discrete water molecules or a 6 water ice Miller. This hypothesis has been retained by a number of 30 31 researchers who suggested that this pathway may be efficient in model and clusters including 4, 9, or 12 water molecules. In − many other environments in the Universe.8 19 Indeed, since its parallel, the reactivity of methaneimine with hydrogen cyanide (step (c), Scheme 1) has been studied both in vacuum and in reactants (hydrogen cyanide, ammonia, and aldehydes such as 32,33 30 formaldehyde and acetaldehyde) have been identified in aqueous solution or in ice environment. interstellar or circumstellar media,20,21 the occurrence of Strecker reaction can be conceivable as a synthetic pathway Received: October 24, 2017 in interstellar media, even in the gas-phase or on the surface of Revised: January 22, 2018 ice grains. Moreover, the presence of imino-diacids (usual Published: January 22, 2018 © 2018 American Chemical Society 1643 DOI: 10.1021/acs.jpca.7b10534 J. Phys. Chem. A 2018, 122, 1643−1657 The Journal of Physical Chemistry A Article Scheme 1. A Simplified View of the Strecker Synthesis of Amino Acids After consideration of these previous results, it appears that times do not bring any new minima. Finally, this procedure led (i) only a limited number of structures have been investigated to several hundred local minima. More precisely, 82 and 306 fi so far and that (ii) the level of theory used was not necessarily minima have been identi ed for (HOCH2NH2)(H2O)n and ffi + su cient to accurately describe the energetic aspects of these 102 and 216 minima for (HOCH2NH3 )(H2O)n with n = 3 and reactions. As a consequence, we propose here a large 4, respectively. All these minima, with up to 4 explicit water investigation at a uniform and high level of theory, including molecules, were obtained from MD simulations that do not optimization of the geometries in the solvent environment, in include any bulk water solvation effect. order to offer a realistic view of the considered systems. For this All structures B obtained from calculations detailed above purpose, we examine here the first steps (a) and (b) (A → B→ were optimized at the DFT level B3LYP using the double-ζ C, Scheme 1) of the Strecker synthesis of glycine (R = H, quality basis set 6-31G(d) and vibrational frequencies were Scheme 1) by quantum chemical calculations, using the calculated at the same DFT level to obtain enthalpies at 0 and 34 accurate composite methods G3B3. Starting from reactions 298 K using the Gaussian09 suite of programs.39 Among these −1 between isolated species CH2 O and NH3 in the gas phase, local minima, those situated in the first ca. 20 kJ mol potential we progressively examine the modifications brought by energy range (B3LYP/6-31G(d) level) were retained in order explicitly including several molecules of water acting as active to construct the potential energy profiles. At this stage, a catalysis and, further, a polarizable medium to mimic the bulk Darwinian tree is built on all minima identified in order to solvent or ice effects. Extensive explorations of the conforma- control that each generation of conformations involving n water tional spaces were done using methods that combine stationary − 50,35 molecule, Gn, is related to the previous generation Gn 1. and dynamic approaches developed in our group. Starting from the most stable structures B, all reaction routes Activation of the reactants by protonation both in the gas- toward reactant A and products C were investigated and the phase and in aqueous media is also considered. Conclusions as corresponding transition and final structures identified. The to whether the Strecker synthesis is a possible pathway leading search for the transition structures B → A and B → C was also to the primitive amino acids are drawn. done at the B3LYP/6-31G(d) level, and their connectivity to reactant and product was systematically controlled by IRC 2. COMPUTATIONAL METHODS calculations. Finally, energy computations were accurately 34 In order to construct the various potential energy profiles performed using the composite G3B3 method. presented in this work, geometries of the central structures B, It may be recalled that, in the G3B3 procedure, the + geometries of stable and transition structures are fully i.e., (HOCH2NH2)(H2O)n and (HOCH2NH3 )(H2O)n (n =0, 1, 2, 3, and 4), were first generated. For low n values, i.e., from 0 optimized with the hybrid functional B3LYP using the 6- to 2, structures B have been built manually, by identifying 31G(d) basis sets. The same level of theory is used (i) to fi calculate harmonic frequencies in order to obtain zero point hydrogen bond acceptor and donor sites able to xH2O molecules. If chemical intuition may be sufficient to explore energies and correction to enthalpies and (ii) to calculate exhaustively the conformational landscapes of structures B,itis entropies in the rigid-rotor harmonic oscillator approximation. not the case for the higher n values where molecular dynamic Improved energies are then obtained by making use of several (MD) simulation is necessary. For these last cases only (n > 2), successive more elaborated computations taking into account we used the methodology combining the hierarchical approach larger basis sets (6-31+G(d), 6-31G(2df,p), and G3large) and and the Darwin’s family tree building.51 Indeed, this method- larger electron correlation levels (MP2, MP4, and QCISD(T)). ology has already been used successfully to study the The final energy is calculated assuming additivity of the various microhydration of protonated Nα-acetylhistidine.35 For the contributions and including a higher-order correction (HLC) hierarchical approach, we used the Tinker molecular modeling to remove the remaining deficiencies in the energy calculations. package36 associated with the AMOEBA37,38 polarizable force The G3B3 method has been assessed on a test set of 299 fi eld. The parameters of the neutral CH2NH2 and CH2OH energies (enthalpies of formation, ionization energies, electron termination of ethylamine and ethanol allowed us to adjust the affinities, and proton affinities).34 The average absolute parameters of HOCH2NH2, respectively. In parallel, the deviations from experiments of these 299 energies are less + −1 parameters of HOCH2NH3 were adjusted on the basis of than 5 kJ mol .