Gas Phase Synthesis of Protonated Glycine by Chemical Dynamics Simulations

Computational studies on the possible formation of glycine via open shell gas-phase chemistry in the interstellar medium

Glycine is the simplest proteinogenic amino acid. It has significant astrobiological implications owing to the ongoing investigation for its detection in the interstellar medium (ISM). Hence, a suitable mechanistic elucidation for its formation in the ISM is of current research interest. In the present work, by employing electronic structure calculations [UCCSD(T) and density functional theory (DFT)], various plausible chemical pathways in the gas phase have been examined for the formation of glycine (whose existence has been indirectly proposed in the ISM) and other simple amino acids (yet to be detected in the ISM) from some simpler molecules present in the ISM. This work suggests that step 1: HO-CO (radical) + CH₂NH → NHCH₂COOH (radical) and step 2a: NHCH₂COOH (radical) + H₂ → glycine + H (radical) have very small barriers of 0.14 kcal mol^-1 and ∼3 kcal mol^-1, respectively (easily surmountable at a temperature of ∼50 K under putative interstellar conditions). Hence this should likely be feasible in interstellar gas-phase chemistry. Therefore, HO-CO (radical), CH₂NH, and H₂ could be the possible precursors for the formation of glycine (subject to the presence of the HO-CO radical). The energetic information related to the interstellar reactions, and how this work takes the putative interstellar conditions into account are presented. This paper also highlights how a reaction found to be unsuitable for interstellar molecular evolution via surface chemistry could nonetheless occur via gas-phase chemistry. Based on our results, this work also recommends detecting three new open-shell molecules, HO-CO radical, NHCH₂COOH radical, and NH₂CHCOOH radical, in the ISM.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations

Glycine (Gly), NH2CH2COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

Study of the Reactivity of CH3COOH+• and COOH+ Ions with CH3NH2: Evidence of the Formation of New Peptide-like C(O)-N Bonds

Acetamide, a small organic compound containing a peptide bond, was observed in the interstellar medium, but reaction pathways leading to the formation of this prebiotic molecule remain uncertain. We investigated the possible formation of a peptide-like bond from the reaction between acetic acid (CH₃-COOH) and methylamine (CH₃-NH₂) that were identified in the interstellar medium. From an experimental point of view, a quadrupole/octopole/quadrupole mass spectrometer was used in combination with synchrotron radiation as a tunable source of VUV photons for monitoring the reactivity of selected ions. Acetic acid was photoionized, and the reactivity of CH₃COOH^+• as well as COOH⁺ (produced from either acetic acid or formic acid) ions with neutral CH₃NH₂ was further studied. With no surprise, charge transfer, proton transfer, and concomitant dissociation processes were found to largely dominate the reactivity. However, a C(O)-N bond formation process between the two reactants was also evidenced, with a weak cross section reaction. From a theoretical point of view, results concerning reactivity and barrier heights were obtained using density functional theory, with the LC-ωPBE range-separated functional in combination with the 6-311++G(d,p) Pople basis set and are in perfect agreement with the experimental data.

Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories

In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.