Preliminary Oligomerization in a Glycolic Acid–Glycine Mixture: A Free Energy Map

Early Steps of Glycolonitrile Oligomerization: A Free-Energy Map

Building on our previous free-energy map (J. Phys. Chem. A 2018, 122, 6769-6779) examining the reactions of monomeric glycolonitrile, we explore the formation of its dimers and trimers in aqueous solution under neutral conditions. While 5-membered rings are kinetically favored, open-chain oligomers with ester or amide linkages are thermodynamically favored. Accessing the 5-membered rings also provides a potential route to glyoxal that bypasses preforming glycolamide, the thermodynamic sink for monomers. However, finding a kinetically accessible route to glycine starting from glycolonitrile in the absence of added ammonia at room temperature proved challenging; the best case involved an intramolecular nucleophilic substitution reaction in a dimer containing neighboring imine and amide groups. Our free-energy map also examines routes to experimentally proposed moieties, explaining why some are observed in low yield or not at all.

Role of Acid in the Co-oligomerization of Formaldehyde and Pyrrole

Building on previous work (J. Phys. Chem. A 2017, 121, 8154-8166) under neutral conditions, we examined the co-oligomerization of CH₂O and pyrrole to form porphryinogen under acidic conditions using density functional theory (B3LYP//6-311G**). Thermodynamically, we found that azafulvene intermediates were significantly stabilized under highly acidic conditions. Kinetically, energy barriers were lowered for C-C bond formation, discriminating in favor of reactions that lead to porphyrinogen. However, it was challenging to satisfactorily combine our thermodynamic and kinetic profiles into a unified free-energy profile because of difficulties in optimizing transition states of cationic species involving proton hops. Instead, we used neutral carboxylic acids as a proxy to study how energy barriers changed. By combining data from both neutral and acidic conditions, we estimate a free-energy profile for the initial steps of oligomerization under milder acidic conditions more relevant to prebiotic chemistry.

Exploring Free Energy Profiles of Uracil and Cytosine Reactions with Formaldehyde

Simple polymers can be potentially formed by the co-oligomerization of pyrimidine nucleobases, uracil and cytosine, with the small molecule formaldehyde. Using density functional calculations, we have constructed a free energy map outlining the thermodynamics and kinetics for (1) the addition of formaldehyde to uracil and cytosine to form hydroxymethylated uracil (HMU) and hydroxymethylated cytosine (HMC), (2) the deamination of cytosine and HMC to uracil and HMU, respectively, and (3) the initial oligomerization of 5-HMU. For the initial formation of monomeric HMU, addition of formaldehyde to the C5 and C6 positions is thermodynamically favored over N1 and N3, but faces higher kinetic barriers, and explains why 5-HMU is the main product observed experimentally. Oligomerization of 5-HMU is thermodynamically favorable although decreasingly so at the tetramer stage. In addition, decreasing concentrations of initial monomer shifts the equilibrium disfavoring oligomer formation.

Reactions of Glycolonitrile with Ammonia and Water: A Free Energy Map

Glycolonitrile, the product of combining CH₂O and HCN, is an intermediate in the Strecker reaction leading to the synthesis of the amino acid glycine. However, besides glycine, a plethora of other compounds are also generated when CH₂O and HCN react in the presence of ammonia and water. As a starting point to analyze the possible components of this complex mixture, we have employed density functional theory to construct a free energy map of all two-carbon (C2) species that may be present when glycolonitrile participates in addition or elimination reactions with ammonia and water. By identifying thermodynamic sinks and kinetic barriers, we find that the myriad C2 species can be grouped into three broad regions across the free energy landscape. This allows us to trace possible routes to glycine and other molecules of interest in the reaction mixture. The present map also extends our previous work on one-carbon (C1) species. We had previously found one issue with our computational protocol in the C1 map; however, our present C2 map provides a larger data set that supports using an empirical correction to our original protocol for imidic acid to amide transformations, without increasing the computational cost, while retaining the original protocol for other classes of reactions.

Studies of intramolecular H-bond interactions and solvent effects in the conformers of glycolic acid — A quantum chemical study

In this work, density functional theory is applied to understand the conformational stability and solvent effects on glycolic acid conformers in different solvents. In addition, the role of intramolecular hydrogen bond (H-bond) interactions in the stability of conformers are investigated. The molecular geometries of selected conformers are optimized using B3LYP and PBE0 functionals with 6-311[Formula: see text]G(d,p) basis set. The effects of solvent on the geometrical parameters, relative stability, dipole moment, chemical hardness, chemical potential, etc. are studied for the conformers of glycolic acid. Our calculations show that the order of stability of the SSC and AAT conformers does not change in liquid phase. However, the energy of SSC and AAT conformers is very close to each other in water media. In water media, strong intramolecular H-bond interaction is present in AAT conformer which causes the energy of AAT conformer to be very close to that of SSC conformer. This may be due to the influence of water media.