Structures and conformational dynamics of monomethylated derivatives of acrolein: A quantum-chemical study

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Theoretical Spectroscopic Study of Two Ketones of Atmospheric Interest: Methyl Glyoxal (CH3COCHO) and Methyl Vinyl Ketone (CH3COCH═CH2)

Two ketones of atmospheric interest, methyl glyoxal and methyl vinyl ketone, are studied using explicitly correlated coupled cluster theory and core-valence correlation-consistent basis sets. The work focuses on the far-infrared region. At the employed level of theory, the rotational constants can be determined to within a few megahertz of the experimental data. Both molecules present two conformers, trans/cis and antiperiplanar (A_p)/synperiplanar (S_p), respectively. trans-Methyl glyoxal and A_p-methyl vinyl ketone are the preferred structures. cis-Methyl glyoxal is a secondary minimum of very low stability, which justifies the unavailability of experimental data in this form. In methyl vinyl ketone, the two conformers are almost isoenergetic, but the interconversion implies a relatively high torsional barrier of 1798 cm^-1. A very low methyl torsional barrier was estimated for trans-methyl glyoxal (V₃ = 273.6 cm^-1). Barriers of 429.6 and 380.7 cm^-1 were computed for A_p- and S_p-methyl vinyl ketone. Vibrational second-order perturbation theory was applied to determine the rovibrational parameters. The far-infrared region was explored using a variational procedure of reduced dimensionality. For trans-methyl glyoxal, the ground vibrational state was estimated to split by 0.067 cm^-1, and the two low excited energy levels (1 0) and (0 1) were found to lie at 89.588 cm^-1/88.683 cm^-1 (A₂/E) and 124.636 cm^-1/123.785 cm^-1 (A₂/E). For A_p- and S_p-methyl vinyl ketone, the ground vibrational state splittings were estimated to be 0.008 and 0.017 cm^-1, respectively.

Conformer-specific VUV-MATI spectroscopy of methyl vinyl ketone: stabilities and cationic structures of the s-trans and s-cis conformers

Methyl vinyl ketone (MVK), a volatile compound with photochemical activity, has received considerable attention in the fields of environmental chemistry and atmospheric chemistry. We explored the conformational stabilities of MVK in the neutral S0 and the cationic D0 states using conformer-specific vacuum ultraviolet mass-analyzed threshold ionization (VUV-MATI) spectroscopy, which provided identifiable vibrational spectra for cationic MVK conformers. Based on the origin bands of the two individual conformers of MVK identified in the MATI spectra under different supersonic expansion conditions, the accurate adiabatic ionization energies of the s-trans and the s-cis conformers were determined to be 77 867 ± 4 (9.6543 ± 0.0005 eV) and 78 222 ± 4 cm-1 (9.6983 ± 0.0005 eV), respectively. The identifiable vibrational spectra of the two cationic conformers were further confirmed using vibrational assignments based on the Franck-Condon fit. Accordingly, precise cationic structures of the MVK conformers could be determined. The structural changes of the two conformers upon ionization could be attributed to the removal of an electron from the highest occupied molecular orbital of each conformer, which consists of nonbonding molecular orbitals on the oxygen atom in the carbonyl group interacting with the σ orbitals in the molecular plane. Consequently, the s-trans conformer was preferred by 48 ± 18 and 403 ± 18 cm-1 in the neutral ground S0 and the cationic D0 states, respectively, which was supported by density-corrected density functional theory calculations and natural bond orbital analysis.

Rate coefficients and product branching ratios for (E)-2-butenal + H reactions

Thermal rate constants for the hydrogen abstraction reactions of (E)-2-butenal by hydrogen atoms were calculated, for the first time, using the multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). After a torsional potential energy surface exploration, ten conformations of the transition states (including the mirror images) were found and separated into four conformational reaction channels (CRCs). Individual energy paths of each CRC were built, recrossing and quantum tunneling effects estimated, and the thermal rate constants obtained. Due to the hindered rotors, the torsional anharmonicity was incorporated in the rate coefficient through the calculations of the rovibrational partition functions using the extended two-dimensional torsional method (E2DT). For comparison, the one-well (1W-CVT/SCT) and harmonic multipath (MP-CVT/SCT) thermal rate constants were also estimated. In addition, kinetic Monte Carlo (KMC) simulations were performed to predict the product branching ratios. For all kinetic approaches, the formation of products of (R1) is predominant. Compared to the harmonic multipath estimation, the percentage of reaction (R4) increases by approximately 9% when the torsional anharmonicity is taken into account. For the reactions (R2) and (R3), the product branching ratio is slightly decreased when compared with the harmonic simulation.

Identification and composition of conformational isomers and their cations in crotonaldehyde by VUV-MATI spectroscopy

Composition of the crotonaldehyde sample at ambient temperature was determined as 93.0%, 3.4%, 3.4%, and 0.2% for the tt-/tc-/ct-/cc-conformers, respectively, for the first time.

Tautomerism in large databases

We have used the Chemical Structure DataBase (CSDB) of the NCI CADD Group, an aggregated collection of over 150 small-molecule databases totaling 103.5 million structure records, to conduct tautomerism analyses on one of the largest currently existing sets of real (i.e. not computer-generated) compounds. This analysis was carried out using calculable chemical structure identifiers developed by the NCI CADD Group, based on hash codes available in the chemoinformatics toolkit CACTVS and a newly developed scoring scheme to define a canonical tautomer for any encountered structure. CACTVS's tautomerism definition, a set of 21 transform rules expressed in SMIRKS line notation, was used, which takes a comprehensive stance as to the possible types of tautomeric interconversion included. Tautomerism was found to be possible for more than 2/3 of the unique structures in the CSDB. A total of 680 million tautomers were calculated from, and including, the original structure records. Tautomerism overlap within the same individual database (i.e. at least one other entry was present that was really only a different tautomeric representation of the same compound) was found at an average rate of 0.3% of the original structure records, with values as high as nearly 2% for some of the databases in CSDB. Projected onto the set of unique structures (by FICuS identifier), this still occurred in about 1.5% of the cases. Tautomeric overlap across all constituent databases in CSDB was found for nearly 10% of the records in the collection.