Received date: August 18, 2015; Accepted date: September 22, 2015; Published date: September 24, 2015
Citation: Ramkissoon C, Szórád JJ, Rágyanszki A, et al. Screening Constitutional Isomers of C3H2O2 for Thermodynamic Properties Obtained from Systematic Electron Structure Calculations. Chem Inform. 2015, 1:2.
Propiolic acid is widely employed in the synthesis and preparation of many reagents, like bis-alkynes, which has revealed cytotoxicity against B16 melanoma cancer cells. In addition, propiolic acid has antifungal and antimicrobial properties. These properties have inspired us to explore the thermodynamic properties of the isomers of propiolic acid, C3H2O2, using electron structure calculations. In all 35 isomers were studied. The thermodynamic data can be useful in the pharmaceutical and medicinal industries.
Bis-alkynes; Molecular isomers; Geometries; Transition state
Propiolic acid is an acetylenic compound with derivatives that may be used in different organic synthesis . Acetylenic compounds can be isolated from plants, fungi and insects . Some have antimicrobial and antifungal properties, which are used by organisms such as the solider beetle, Chaudiognathus leconteri, as a defensive mechanism against its predators . There is ongoing research that focuses on the potential of acetylenic fatty acids as the forefronts in generating new topical treatments for fungal infections .
Propiolic acid is also used in the preparation of bis-alkynes which has revealed noticeable cytotoxicity against certain types of cancer cells such as the B16 melanoma cancer cells . Computational techniques have been used to predict molecular properties of propiolic acid  and in the elucidation of reaction mechanism for decarboxylation of various carboxylic acids . Relevant properties of propiolic acid in industrial context are found in ref. [8,9].
Traditional research in organic chemistry started with the synthesis of a given molecule. Subsequently its various properties, like thermodynamic stability, were determined experimentally. With the arrival of high-performance computers and with the associated software development one can generate all possible isomeric structures computationally via combinatorial chemistry. This allows one to construct virtual libraries of molecular isomers. Quantum chemical molecular computations make it possible to determine the thermodynamic stability of each of these isomers even before one would attempt their synthesis in the laboratory. Various areas of non-accessible regions of chemical research would benefit from this new approach, such as drug-design and the study of molecular evolution to mention the two most obvious aspects of "molecular design". In other words this new approach can help us to study large number of molecular isomers in a variety of areas of organic chemistry.
The present study starts with a very small molecular composition of C3H2O2 that has only a relatively few isomers. In this case it is easier to demonstrate the various principles of this method of research.
All of the isomers for C3H2O2 with a singlet electronic state were generated using valence restrictions and geometry parameters from the MM2 molecular mechanics program. The hydrogen, oxygen and carbon atoms were considered to form bond with 1, 2 and 4 valences, respectively. The C3H2O2 composition could occur either as a single molecule (29 cases) or as a molecular complex formed from two of its fragments (6 cases). Complexes with more than two fragments were not found. Based on graph theory Molgen 5.0 program  was employed to generate all possible singlet state initial structures using the above mentioned restrictions. These structures were used as the primary structures to determine the local minima of the potential energy surface of C3H2O2. Geometries were optimized at the B3LYP/6-31G(d) level of theory. The maximum optimization step size was restricted to 0.01 Bohr and the maximum number of the optimization steps were set to 300. During optimization analytical first and second derivatives were calculated in each step. The normal mode analysis was performed to secure only local minima of the multidimensional potential energy surface. After optimization, the G3MP2B3 composite method  was used to derive standard thermochemical properties. Subsequently, the enthalpy of formation, ΔfΗ°, was found using the atomization scheme. The relative Gibbs free energy values were then used to rank the thermodynamic stability of the thirty-five isomers.
29 out of the 35 calculated molecular structures have conventional chemical bonding while the remaining 6 are bimolecular complexes consisting of two molecular fragments. The fragments were small molecules such as CO, CO2, H2O2, O2, C2H2, C2H2O, CH2O, C3.
Although, the structural generator software (Molgen 5.0) is expected to produce only tetravalent carbons, occasionally divalent carbene type carbon atoms have also appeared in some of the computed structures. This happened both in the open chained compounds as well as among the cyclic compounds. The structures of all of the 35 isomers are shown in Table 1 along with their thermodynamic data. It appears that that there are many ring systems among the isomers. There are five three membered rings, eight four membered rungs and six five membered rings. The clustering of the isomers in unimolecular structures and bimolecular complexes is shown in the entropy funnel representation in Figure 1, which illustrates the molecular network of stable isomers associated with the potential energy hyper-surface (PEHS) that describes the C3H2O2 empirical formula. The left-hand side of the double funnel collects all the unimolecular structures while the right hand side collects those consisting of two fragments. Only those isomers of the network can be joined by a line that may be directly inter-converted by a chemical reaction. For example structure #4 in Table 1 (H-CCCOOH) may undergo decarboxylation to form structure #1 ([CO2, H2C2])
Data for the conformer where the OH group is oriented towards the CCH moiety
|28||3-ethenylidenedioxirane or 3-vinylidenedioxirane
|31||Complex of dioxiren-1-ium-3-ide
compound with acetylene
Table 1: List of the investigated 35 isomers of C3H2O2. ΔrelG° (kJ/mol) is the standard free energy relative to that of the global minimum of [CO2, H2C2]. S° (J/(mol*K) is the standard entropy and ΔfH° (kJ/mol) is
the standard enthalpy of formation. ΔrelG° and S° were calculated at the G3MP2B3 level while ΔfH° was obtained using the atomization scheme.
H-CC-COOH → [CO2, H2C2]
In that case structures #4 and #1, in Figure 1, may be interconnected by a linear, or curvy-linear, path that has a transition state somewhere along this line.
Since the propiolic acid is a well-studied compound one might expect that it was the most stable among the isomers, but it is not. The isomer with the lowest free energy (the global minimum) turns out to be a planar complex consisting of H2C2 and CO2. The complex has C2v symmetry. For convenience the stability of all isomers are listed relative to the global minimum as ΔrelGº. According to the G3MP2B3 calculations, propiolic acid exists in two conformations, differing in the orientation of the OH group relative to the CCH fragment. The conformer where the OH group is oriented away from the CCH fragment (conformer A) has a free energy which is 15 kJ/mol lower than the other isomer (conformer B). In Table 1, we list all the 35 isomers ranked according to ΔrelGº. Table 1 includes also ΔfHº and Sº along with an indication of the structure.
Transition states for interconversion
Topological properties of the potential energy hyper-surface: The C3H2O2 structure has seven atoms (n=7) implying the number of internal coordinates, 3n-6, is 15. This means that the PEHS has 15 independent variables:
The 35 minima are located on the PEHS, however, without explicit knowledge of the topology of the hyper-surface it is impossible to know which pair of minima is adjacent and which ones are further separated from each other. Those that are adjacent have a transition state in between to enhance the conversion of one minimum to another. Other non-adjacent minima can only be converted through multiple conversions involving more than one transition state.
case study for interconversion: The transition from the global minimum, [CO2, C2H2], to the propionic acid isomer has been determined using the QST2 method implemented in Gausian09 . The transition state was determined for both conformers of propiolic acid. It turns out that the two transition states are identical. Investigation of the transition state using the internal reaction coordinate method  shows that it connects the global minimum to the conformer, where the OH group is oriented towards the CCH moiety of propionic acid (conformer B). The structure of the transition state is depicted in Figure 2. The free energy barrier for this transition at standard condition is 321 kJ/mol at the G3MP2B3 level. The barrier is so high that the reaction of CO2, C2H2], to the propionic acid isomer has been determined usin and C2H2 to form propionic acid will not happen - in line with observations. Calculations of the transition state using other levels of theory lead to similar estimates for the barrier height as reported in Table 2. Conformer B of propionic acid is connected to conformer A (which has lower free energy). A QST2 calculation leads to an estimate for the barrier height of 32 kJ/ mol (G3MP2B3). Thus the transition from the global ground state to the most stable isomer of propionic acid is a two state process.
Figure 2: Transition state connecting the global minimum, [CO2, C2H2], among the C3H2O2 isomers with the propiolic acid isomer, HCC-COOH, in a conformation where the OH group is oriented towards the CCH moiety. The transition state is calculated at the G3MP2B3 level of theory. The angle O2-C1-O3 is 144°. The distances 7H-2O and 7H-5C are 1.32 Å and 1.33 Å, respectively. The structure of the transition state is planar.
|Level of theory||E1>||E2|
Table 2: Barrier heights (kJ/mol) for the transition from the global minimum, [CO2, C2H2], to propiolic acid calculated at various levels of theory. E1 is the barrier from the global minimum to conformer B of propionic acid and E2 is the barrier from conformer B to conformer A.
It was shown that new molecular structures can be predicted by molecular computations and the relative stabilities can be assessed by computing thermodynamic functions. Even though this is a small step on a long road, it is hoped that such research would be useful in designing new drugs and other desirable molecules.
All Published work is licensed under a Creative Commons Attribution 4.0 International License
Copyright © 2019 All rights reserved. iMedPub LTD Last revised : April 21, 2019