1
|
Stylianakis I, Zervos N, Lii JH, Pantazis DA, Kolocouris A. Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory. J Comput Aided Mol Des 2023; 37:607-656. [PMID: 37597063 PMCID: PMC10618395 DOI: 10.1007/s10822-023-00513-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Accepted: 06/03/2023] [Indexed: 08/21/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ioannis Stylianakis
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Nikolaos Zervos
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece
| | - Jenn-Huei Lii
- Department of Chemistry, National Changhua University of Education, Changhua City, Taiwan
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| | - Antonios Kolocouris
- Department of Medicinal Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimioupolis Zografou, 15771, Athens, Greece.
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece.
| |
Collapse
|
2
|
Neier R, Baak M, Jaun BM, Belaj F. Conformations of 4‐tert.‐butyloxy‐, 4‐trimethyl silyloxy‐ and 4‐trimethyl stannyloxy‐6‐methyl‐2‐pyrones in the crystalline state and in solution. Helv Chim Acta 2022. [DOI: 10.1002/hlca.202200029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Reinhard Neier
- University of Neuchatel Department of Chemistry Av. Bellevaux 51 CH-2000 Neuchatel SWITZERLAND
| | - Marcel Baak
- BFH TI: Berner Fachhochschule Technik und Informatik Abt MNG Allgemeinbildung Quellgasse 21 2501 Biel SWITZERLAND
| | - Bernhard Mathias Jaun
- ETH-Zürich LOC: Eidgenossische Technische Hochschule Zurich Laboratorium fur Organische Chemie LOC: Laboratorium für organische Chemie Vladimir-Prelog-Weg 1-5/10HCI D314 8093 Zürich SWITZERLAND
| | - Ferdinand Belaj
- Karl-Franzens-Universität Graz: Karl-Franzens-Universitat Graz Institut für Chemie Schubertstraße 1/III 8010 Graz AUSTRIA
| |
Collapse
|
3
|
Schleif T, Mieres-Perez J, Henkel S, Mendez-Vega E, Inui H, McMahon RJ, Sander W. Conformer-Specific Heavy-Atom Tunneling in the Rearrangement of Benzazirines to Ketenimines. J Org Chem 2019; 84:16013-16018. [PMID: 31730349 DOI: 10.1021/acs.joc.9b02482] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
5-Methoxy-2H-benzazirine was prepared via irradiation of the corresponding phenyl azide, isolated in an argon matrix at cryogenic temperatures. It undergoes ring expansion to the corresponding ketenimine in the dark at T < 30 K despite a calculated activation barrier of 4.9 kcal mol-1 [B3LYP/6-311++G(d,p)]. Since this rearrangement proceeds with a rate constant in the order of 10-4 s-1, exhibiting only a shallow temperature dependence, the results are interpreted in terms of heavy-atom tunneling. Of the four isomeric benzazirines resulting from the initial photolysis, only one can be observed to rearrange; this conformer specificity is explained by the other potentially observable rearrangements being either too fast or too slow to be detected due to the differences in heights and widths of their respective activation barriers.
Collapse
Affiliation(s)
- Tim Schleif
- Lehrstuhl für Organische Chemie II , Ruhr-Universität Bochum , 44801 Bochum , Germany
| | - Joel Mieres-Perez
- Lehrstuhl für Organische Chemie II , Ruhr-Universität Bochum , 44801 Bochum , Germany
| | - Stefan Henkel
- Lehrstuhl für Organische Chemie II , Ruhr-Universität Bochum , 44801 Bochum , Germany
| | - Enrique Mendez-Vega
- Lehrstuhl für Organische Chemie II , Ruhr-Universität Bochum , 44801 Bochum , Germany
| | - Hiroshi Inui
- Department of Chemistry, School of Science , Kitasato University , 1-15-1 Kitasato, Minami-ku , Sagamihara , Kanagawa 252-0373 , Japan
| | - Robert J McMahon
- Department of Chemistry , University of Wisconsin-Madison , 1101 University Avenue , Madison , Wisconsin 53706-1322 , United States
| | - Wolfram Sander
- Lehrstuhl für Organische Chemie II , Ruhr-Universität Bochum , 44801 Bochum , Germany
| |
Collapse
|
4
|
Janesko BG, Villegas H. Attractive Nonbonded Interactions Help Stabilize the Z Form of Alkenyl Anions. J Org Chem 2018; 83:8208-8213. [PMID: 30003784 DOI: 10.1021/acs.joc.8b00960] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Alkenyl anions ((-)H2C-CH═CH-Y; Y = aliphatic, aromatic, amine, ether) are often thermodynamically stable in the sterically congested Z form. This preference for Z structures is used to control regiochemistry in organometallic and Grignard reactions, allyl amine rearrangements to enamines, and allyl ether rearrangements to enol ethers. Explanations for Z stability in alkenes (Y = CH2R) typically invoke through-space attraction (Coulomb or charge transfer) between the formally anionic carbon C1 and the Y = CH2R hydrogens. However, this explanation is difficult to generalize to amines and ethers. We suggest that the orbital-driven so-called "attractive nonbonded interactions" suggested to stabilize the Z forms of 1,2-difluoroethylene and 1-substituted propenes also help stabilize Z alkenyl anions. We present electronic structure calculations and surveys of the experimental literature to show these effects' relevance for alkenyl anions. Our results suggest new approaches for regiocontrol in reactions with alkenyl anion intermediates and motivate revisiting older orbital-based theories of "attractive nonbonded interactions".
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry and Biochemistry , Texas Christian University , Fort Worth , Texas 76129 , United States
| | - Hector Villegas
- Department of Chemistry and Biochemistry , Texas Christian University , Fort Worth , Texas 76129 , United States
| |
Collapse
|
5
|
Matta CF. How dependent are molecular and atomic properties on the electronic structure method? Comparison of Hartree-Fock, DFT, and MP2 on a biologically relevant set of molecules. J Comput Chem 2010; 31:1297-311. [PMID: 19882732 DOI: 10.1002/jcc.21417] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This article compares molecular properties and atomic properties defined by the quantum theory of atoms in molecules (QTAIM) obtained from three underlying levels of theory: MP2(full), density functional theory (DFT) (B3LYP), and Hartree-Fock (H-F). The same basis set (6-311++G(d,p)) has been used throughout the study. The calculations and comparisons were applied to a set of 30 small molecules representing common fragments of biological molecules. The molecular properties investigated are the energies and the electrostatic moments (up to and including the quadrupoles), and the atomic properties include electron populations (and atomic charge), atomic dipolar and quadrupolar polarizations, atomic volumes, and corrected and raw atomic energies. The Cartesian distance between dipole vectors and the Frobenius distance between the quadrupole tensors calculated at the three levels of theory provide a measure of their correlation (or lack thereof). With the exception of energies (atomic and molecular), it is found that both DFT and H-F are in excellent agreement with MP2, especially with regards to the electrostatic mutipoles up to the quadrupoles, but DFT and MP2 agree better in almost all studied properties (with the exception of molecular geometries). QTAIM properties whether obtained from H-F, DFT(B3LYP), or MP2 calculations when used in the construction of empirical correlations with experiment such as quantitative structure-activity-(or property)-relationships (QSAR/QSPR) are equivalent (because the properties calculated at the three levels are very highly correlated among themselves with r(2) typically >0.95, and therefore preserving trends). These results suggest that the massive volume of results that were published in the older literature at the H-F level is valid especially when used to study trends or in QSAR or QSPR studies, and, as long as our test set of molecules is representative, there is no pressing need to re-evaluate them at other levels of theory except when inadequate basis sets were used by today's standards. Extensive tabulation of molecular and atomic properties at the three theoretical levels is available in the Supporting Information, including optimized geometries, molecular energies, virial ratios, molecular electrostatic moments up to and including hexadecapoles, atomic populations, atomic volumes, atomic electrostatic moments up to and including the quadrupoles, and atomic energies.
Collapse
Affiliation(s)
- Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada B3M 2J6.
| |
Collapse
|
6
|
|
7
|
Wolstenholme DJ, Cameron TS. A comparison of the energetic and topological properties of weak interactions in molecular organic crystals. CAN J CHEM 2007. [DOI: 10.1139/v07-090] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Several molecular organic crystals [tetraphenylphosphonium squarate, bianthrone, and bis(benzophenone) azine] have been analyzed topologically following a multipole refinement of each compound. This has resulted in the complete characterization of all the weak interactions present in these systems (C-H···O, C-H···C, and C-Hδ+···δ+H-C). Correlations have been found to exist between the energetic and topological properties for each type of interaction. The coexistence of the three types of interactions provides a unique opportunity to examine their energetic differences and similarities. This examination can also be used to help understand the complex nature of C-Hδ+···δ+H-C type interactions, in particular those characterized as H–H bonds. All this leads to a better understanding of the role H–H bonds play in molecular organic crystals compared with other weak interactions such as hydrogen bonds and van der Waals interactions.Key words: energetics, multipole refinement, H–H bonds, hydrogen bonds.
Collapse
|
8
|
Wolstenholme DJ, Matta CF, Cameron TS. Experimental and Theoretical Electron Density Study of a Highly Twisted Polycyclic Aromatic Hydrocarbon: 4-Methyl-[4]helicene. J Phys Chem A 2007; 111:8803-13. [PMID: 17705352 DOI: 10.1021/jp071002+] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Helicenes are molecules of considerable interest in view of their aromaticity which persists despite a marked departure from planarity and because of the extreme potency of some of their metabolites as tumor and mutation promoters. In this study, the electron density of 4-methyl-[4]helicene (or 4-methylbenzo[c]phenanthrene) is studied topologically with an emphasis on the fjord region since this region is where metabolic activation is initiated. The molecule consists of four fused aromatic rings that assume a twisted geometry. This geometry brings two hydrogen atoms into close proximity in the fjord region of the molecule accompanied by the appearance of an intramolecular C-Hdelta+...delta+H-C bond path (an interaction termed hydrogen-hydrogen or H- H bonding to distinguish it from dihydrogen bonding from which it is qualitatively distinct). In addition to the intramolecular H-H interaction, a number of intermolecular interactions are shown to be involved in the packing of this molecule in the crystalline state. The effect of the nonplanarity of the molecule on the local aromaticity of each ring is also discussed.
Collapse
Affiliation(s)
- David J Wolstenholme
- Department of Chemistry, Dalhousie University, Halifax, Nova Scotia B3H 4J3, Canada
| | | | | |
Collapse
|
9
|
Ignatyev IS, Montejo M, Sundius T, Partal Ureña F, López González JJ. Structure and vibrational spectra of vinyl ether conformers. The comparison of B3LYP and MP2 predictions. Chem Phys 2007. [DOI: 10.1016/j.chemphys.2007.01.013] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
10
|
da Silva G, Kim CH, Bozzelli JW. Thermodynamic Properties (Enthalpy, Bond Energy, Entropy, and Heat Capacity) and Internal Rotor Potentials of Vinyl Alcohol, Methyl Vinyl Ether, and Their Corresponding Radicals. J Phys Chem A 2006; 110:7925-34. [PMID: 16789782 DOI: 10.1021/jp0602878] [Citation(s) in RCA: 80] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Vinyl alcohols (enols) have been discovered as important intermediates and products in the oxidation and combustion of hydrocarbons, while methyl vinyl ethers are also thought to occur as important combustion intermediates. Vinyl alcohol has been detected in interstellar media, while poly(vinyl alcohol) and poly(methyl vinyl ether) are common polymers. The thermochemical property data on these vinyl alcohols and methyl vinyl ethers is important for understanding their stability, reaction paths, and kinetics in atmospheric and thermal hydrocarbon-oxygen systems. Enthalpies , entropies , and heat capacities (C(p)()(T)) are determined for CH(2)=CHOH, C(*)H=CHOH, CH(2)=C(*)OH, CH(2)=CHOCH(3), C(*)H=CHOCH(3), CH(2)=C(*)OCH(3), and CH(2)=CHOC(*)H(2). Molecular structures, vibrational frequencies, , and C(p)(T) are calculated at the B3LYP/6-31G(d,p) density functional calculation level. Enthalpies are also determined using the composite CBS-Q, CBS-APNO, and G3 methods using isodesmic work reactions to minimize calculation errors. Potential barriers for internal rotors are calculated at the B3LYP/6-31G(d,p) level and used to determine the hindered internal rotational contributions to entropy and heat capacity. The recommended ideal gas phase values calculated in this study are the following (in kcal mol(-1)): -30.0, -28.9 (syn, anti) for CH(2)=CHOH; -25.6, -23.9 for CH(2)=CHOCH(3); 31.3, 33.5 for C(*)H=CHOH; 27.1 for anti-CH(2)=C(*)OH; 35.6, 39.3 for C(*)H=CHOCH(3); 33.5, 32.2 for CH(2)=C(*)OCH(3); 21.3, 22.0 for CH(2)=CHOC(*)H(2). Bond dissociation energies (BDEs) and group additivity contributions are also determined. The BDEs reveal that the O-H, O-CH(3), C-OH, and C-OCH(3) bonds in vinyl alcohol and methyl vinyl ether are similar in energy to those in the aromatic molecules phenol and methyl phenyl ether, being on average around 3 kcal mol(-1) weaker in the vinyl systems. The keto-enol tautomerization enthalpy for the interconversion of vinyl alcohol to acetaldehyde is determined to be -9.7 kcal mol(-1), while the activation energy for this reaction is calculated as 55.9 kcal mol(-1); this is the simplest keto-enol tautomerization and is thought to be important in the reactions of vinyl alcohol. Formation of the formyl methyl radical (vinoxy radical/vinyloxy radical) from both vinyl alcohol and methyl vinyl ether is also shown to be important, and its reactions are discussed briefly.
Collapse
Affiliation(s)
- Gabriel da Silva
- New Jersey Institute of Technology, Department of Chemistry and Environmental Science, Newark, 07102, USA
| | | | | |
Collapse
|