1
|
Lee AS, Elliott S, Harb H, Ward L, Foster I, Curtiss L, Assary RS. Emin: A First-Principles Thermochemical Descriptor for Predicting Molecular Synthesizability. J Chem Inf Model 2024; 64:1277-1289. [PMID: 38359461 DOI: 10.1021/acs.jcim.3c01583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
Predicting the synthesizability of a new molecule remains an unsolved challenge that chemists have long tackled with heuristic approaches. Here, we report a new method for predicting synthesizability using a simple yet accurate thermochemical descriptor. We introduce Emin, the energy difference between a molecule and its lowest energy constitutional isomer, as a synthesizability predictor that is accurate, physically meaningful, and first-principles based. We apply Emin to 134,000 molecules in the QM9 data set and find that Emin is accurate when used alone and reduces incorrect predictions of "synthesizable" by up to 52% when used to augment commonly used prediction methods. Our work illustrates how first-principles thermochemistry and heuristic approximations for molecular stability are complementary, opening a new direction for synthesizability prediction methods.
Collapse
Affiliation(s)
- Andrew S Lee
- Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Sarah Elliott
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Hassan Harb
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Logan Ward
- Data Science and Learning Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Ian Foster
- Data Science and Learning Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Larry Curtiss
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Rajeev S Assary
- Materials Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
2
|
Petelski AN, Marquez J, Pamies SC, Sosa GL, Peruchena NM. Understanding the Chloride Affinity of Barbiturates for Anion Receptor Design. Chemphyschem 2021; 22:665-674. [PMID: 33538090 DOI: 10.1002/cphc.202100008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Revised: 02/02/2021] [Indexed: 11/09/2022]
Abstract
Due to their potential binding sites, barbituric acid (BA) and its derivatives have been used in metal coordination chemistry. Yet their abilities to recognize anions remain unexplored. In this work, we were able to identify four structural features of barbiturates that are responsible for a certain anion affinity. The set of coordination interactions can be finely tuned with covalent decorations at the methylene group. DFT-D computations at the BLYP-D3(BJ)/aug-cc-pVDZ level of theory show that the C-H bond is as effective as the N-H bond to coordinate chloride. An analysis of the electron charge density at the C-H⋅⋅⋅Cl- and N-H⋅⋅⋅Cl- bond critical points elucidates their similarities in covalent character. Our results reveal that the special acidity of the C-H bond shows up when the methylene group moves out of the ring plane and it is mainly governed by the orbital interaction energy. The amide and carboxyl groups are the best choices to coordinate the ion when they act together with the C-H bond. We finally show how can we use this information to rationally improve the recognition capability of a small cage-like complex that is able to coordinate NaCl.
Collapse
Affiliation(s)
- Andre Nicolai Petelski
- Departamento de Ingeniería Química, Grupo de Investigación en Química Teórica y Experimental (QUITEX), Universidad Tecnológica Nacional, Facultad Regional Resistencia, French 414, H3500CHJ), Resistencia, Chaco, Argentina.,Instituto de Química Básica y Aplicada del Nordeste Argentino, IQUIBA-NEA, UNNE-CONICET, Avenida Libertad 5460, 3400, Corrientes, Argentina
| | - Josefina Marquez
- Departamento de Ingeniería Química, Grupo de Investigación en Química Teórica y Experimental (QUITEX), Universidad Tecnológica Nacional, Facultad Regional Resistencia, French 414, H3500CHJ), Resistencia, Chaco, Argentina
| | - Silvana Carina Pamies
- Departamento de Ingeniería Química, Grupo de Investigación en Química Teórica y Experimental (QUITEX), Universidad Tecnológica Nacional, Facultad Regional Resistencia, French 414, H3500CHJ), Resistencia, Chaco, Argentina
| | - Gladis Laura Sosa
- Departamento de Ingeniería Química, Grupo de Investigación en Química Teórica y Experimental (QUITEX), Universidad Tecnológica Nacional, Facultad Regional Resistencia, French 414, H3500CHJ), Resistencia, Chaco, Argentina.,Instituto de Química Básica y Aplicada del Nordeste Argentino, IQUIBA-NEA, UNNE-CONICET, Avenida Libertad 5460, 3400, Corrientes, Argentina
| | - Nélida María Peruchena
- Instituto de Química Básica y Aplicada del Nordeste Argentino, IQUIBA-NEA, UNNE-CONICET, Avenida Libertad 5460, 3400, Corrientes, Argentina.,Área de Química Física - Departamento de Química, Laboratorio de Estructura Molecular y Propiedades (LEMyP), Facultad de Ciencias Exactas y Naturales y Agrimensura, Universidad Nacional del Nordeste, Avenida Libertad 5460, 3400, Corrientes, Argentina
| |
Collapse
|
3
|
Kuznetsov ML. Strength of the [Z-I···Hal] - and [Z-Hal···I] - Halogen Bonds: Electron Density Properties and Halogen Bond Length as Estimators of Interaction Energy. Molecules 2021; 26:2083. [PMID: 33916483 PMCID: PMC8038634 DOI: 10.3390/molecules26072083] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Revised: 03/29/2021] [Accepted: 03/31/2021] [Indexed: 11/29/2022] Open
Abstract
Bond energy is the main characteristic of chemical bonds in general and of non-covalent interactions in particular. Simple methods of express estimates of the interaction energy, Eint, using relationships between Eint and a property which is easily accessible from experiment is of great importance for the characterization of non-covalent interactions. In this work, practically important relationships between Eint and electron density, its Laplacian, curvature, potential, kinetic, and total energy densities at the bond critical point as well as bond length were derived for the structures of the [Z-I···Hal]- and [Z-Hal···I]- types bearing halogen bonds and involving iodine as interacting atom(s) (totally 412 structures). The mean absolute deviations for the correlations found were 2.06-4.76 kcal/mol.
Collapse
Affiliation(s)
- Maxim L. Kuznetsov
- Centro de Química Estrutural, Instituto Superior Técnico, Universidade de Lisboa, Avenida Rovisco Pais, 1049-001 Lisbon, Portugal; ; Tel.: +351-218-419-236
- Institute of Chemistry, Saint Petersburg State University, Universitetskaya Nab. 7/9, 199034 Saint Petersburg, Russia
| |
Collapse
|
4
|
Jiménez‐Grávalos F, Gallegos M, Martín Pendás Á, Novikov AS. Challenging the electrostatic
σ
‐hole picture of halogen bonding using minimal models and the interacting quantum atoms approach. J Comput Chem 2021; 42:676-687. [DOI: 10.1002/jcc.26488] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Revised: 01/10/2021] [Accepted: 01/18/2021] [Indexed: 12/14/2022]
Affiliation(s)
| | - Miguel Gallegos
- Department of Analytical and Physical Chemistry University of Oviedo Oviedo Spain
| | - Ángel Martín Pendás
- Department of Analytical and Physical Chemistry University of Oviedo Oviedo Spain
| | | |
Collapse
|
5
|
Zarechnaya OM, Anisimov AA, Belov EY, Burakov NI, Kanibolotsky AL, Mikhailov VA. Polycentric binding in complexes of trimethylamine- N-oxide with dihalogens. RSC Adv 2021; 11:6131-6145. [PMID: 35423161 PMCID: PMC8694807 DOI: 10.1039/d0ra08165e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Accepted: 01/13/2021] [Indexed: 02/01/2023] Open
Abstract
Dihalogens readily interact with trimethylamine-N-oxide under ambient conditions. Accordingly, herein, stable 1 : 1 adducts were obtained in the case of iodine chloride and iodine bromide. The crystal and molecular structure of the trimethylamine-N-oxide-iodine chloride adduct was solved. Furthermore, the geometry and electronic structure of the trimethylamine-N-oxide-dihalogen complexes were studied computationally. Only molecular ensembles were found in the global minimum for the 1 : 1 stoichiometry. The O⋯X-Y halogen bond is the main factor for the thermodynamic stability of these complexes. Arguments for electrostatic interactions as the driving force for this noncovalent interaction were discussed. Also, the equilibrium structures are additionally stabilised by weak C-H⋯X hydrogen bonds. Consequently, formally monodentate ligands are bound in a polycentric manner.
Collapse
Affiliation(s)
- Olga M Zarechnaya
- L.M. Litvinenko Institute of Physical Organic and Coal Chemistry R. Luxemburg St., 70 Donetsk Ukraine
| | - Aleksei A Anisimov
- A.N. Nesmeyanov Institute of Organoelement Compounds, Russian Academy of Sciences 28 Vavilov St. 119991 Moscow Russia
- D.I. Mendeleev Russian Chemical Technological University 9 Miusskaya Sq. 125047 Moscow Russia
| | - Eugenii Yu Belov
- L.M. Litvinenko Institute of Physical Organic and Coal Chemistry R. Luxemburg St., 70 Donetsk Ukraine
| | - Nikolai I Burakov
- L.M. Litvinenko Institute of Physical Organic and Coal Chemistry R. Luxemburg St., 70 Donetsk Ukraine
| | | | - Vasilii A Mikhailov
- L.M. Litvinenko Institute of Physical Organic and Coal Chemistry R. Luxemburg St., 70 Donetsk Ukraine
| |
Collapse
|
6
|
Guevara-Vela JM, Francisco E, Rocha-Rinza T, Martín Pendás Á. Interacting Quantum Atoms-A Review. Molecules 2020; 25:E4028. [PMID: 32899346 PMCID: PMC7504790 DOI: 10.3390/molecules25174028] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 08/24/2020] [Accepted: 08/26/2020] [Indexed: 12/16/2022] Open
Abstract
The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.
Collapse
Affiliation(s)
- José Manuel Guevara-Vela
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., Mexico City 04510, Mexico; (J.M.G.-V.); (T.R.-R.)
| | - Evelio Francisco
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006 Oviedo, Spain;
| | - Tomás Rocha-Rinza
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P., Mexico City 04510, Mexico; (J.M.G.-V.); (T.R.-R.)
| | - Ángel Martín Pendás
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006 Oviedo, Spain;
| |
Collapse
|
7
|
Miranda MO, Duarte DJR, Alkorta I. Anion‐Anion Complexes Established between Aspartate Dimers. Chemphyschem 2020; 21:1052-1059. [DOI: 10.1002/cphc.201901200] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Revised: 02/28/2020] [Indexed: 01/20/2023]
Affiliation(s)
- Matias O. Miranda
- Laboratorio de Estructura Molecular y Propiedades Departamento de Química Facultad de Ciencias Exactas y Naturales y AgrimensuraUniversidad Nacional del Nordeste Avenida Libertad 5460 3400 Corrientes Argentina
- IQUIBA-NEA (UNNE-CONICET) Avenida Libertad 5460 3400 Corrientes Argentina
| | - Darío J. R. Duarte
- Laboratorio de Estructura Molecular y Propiedades Departamento de Química Facultad de Ciencias Exactas y Naturales y AgrimensuraUniversidad Nacional del Nordeste Avenida Libertad 5460 3400 Corrientes Argentina
- IQUIBA-NEA (UNNE-CONICET) Avenida Libertad 5460 3400 Corrientes Argentina
| | - Ibon Alkorta
- Instituto de Química Médica, CSIC Juan de la Cierva, 3 28006 Madrid Spain
| |
Collapse
|
8
|
Souza GK, Kischkel B, Freitas CF, Negri M, Back D, Johann G, Hioka N, Schuquel ITA, Santin SMO, Pomini AM. Antiproliferative activity and energy calculations of a new triterpene isolated from the palm tree Acrocomia totai. Nat Prod Res 2019; 35:4225-4234. [PMID: 31773984 DOI: 10.1080/14786419.2019.1696331] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Acrocomia totai Mart (Arecaceae) is a palm tree native to South America, widely studied for biodiesel production. The aim of this work was to perform the first phytochemical study of A. totai leaves, as well as to do biological assays against human cancer cell lines. A new triterpene of the hopane class named totaiol (1), three known triterpenes (2-4), and two phytosteroids (5-6) were identified. The new natural product was characterized using 1 D and 2 D NMR, single crystal X-ray diffraction analises, and high resolution mass spectrometry. The intercontacts in the crystal packing were also analised. Complete stereochemical characterization of compound 1 revealed an unusual positioning pattern for methyl and isopropenyl groups in the polycyclic skeleton. Compounds 1-5 were evaluated for the first time in antiproliferative assays against Ca Ski, MCF-7 and MCF-10 cells. The new natural product was active against Ca Ski cells with IC50 ≤ 6.25 µg mL-1.
Collapse
Affiliation(s)
- Gredson Keiff Souza
- Department of Chemistry, State University of Maringá, Maringá, Brazil.,Faculty of Administration and Economic Sciences, Cianorte, Brazil
| | - Brenda Kischkel
- Clinical Analysis Department, State University of Maringá, Maringá, Brazil
| | | | - Melyssa Negri
- Clinical Analysis Department, State University of Maringá, Maringá, Brazil
| | - Davi Back
- Department of Chemistry, Federal University of Santa Maria, Santa Maria, Brazil
| | - Gracielle Johann
- Federal Technological University of Paraná, Dois Vizinhos, Brazil
| | - Noboru Hioka
- Department of Chemistry, State University of Maringá, Maringá, Brazil
| | | | | | | |
Collapse
|
9
|
Del Bene JE, Alkorta I, Elguero J. What Types of Noncovalent Bonds Stabilize Dimers (XCP)2, for X = CN, Cl, F, and H? J Phys Chem A 2019; 123:10086-10094. [DOI: 10.1021/acs.jpca.9b08141] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Janet E. Del Bene
- Department of Chemistry, Youngstown State University, Youngstown, Ohio 44555, United States
| | - Ibon Alkorta
- Instituto de Química Médica (CSIC), Juan de la Cierva, 3, E-28006 Madrid, Spain
| | - José Elguero
- Instituto de Química Médica (CSIC), Juan de la Cierva, 3, E-28006 Madrid, Spain
| |
Collapse
|
10
|
|
11
|
Can Combined Electrostatic and Polarization Effects Alone Explain the F···F Negative-Negative Bonding in Simple Fluoro-Substituted Benzene Derivatives? A First-Principles Perspective. COMPUTATION 2018. [DOI: 10.3390/computation6040051] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The divergence of fluorine-based systems and significance of their nascent non-covalent chemistry in molecular assemblies are presented in a brief review of the field. Emphasis has been placed to show that type-I and -II halogen-centered F···F long-ranged intermolecular distances viable between the entirely negative fluorine atoms in some fluoro-substituted dimers of C6H6 can be regarded as the consequence of significant non-covalent attractive interactions. Such attractive interactions observed in the solid-state structures of C6F6 and other similar fluorine-substituted aromatic compounds have frequently been underappreciated. While these are often ascribed to crystal packing effects, we show using first-principles level calculations that these are much more fundamental in nature. The stability and reliability of these interactions are supported by their negative binding energies that emerge from a supermolecular procedure using MP2 (second-order Møller-Plesset perturbation theory), and from the Symmetry Adapted Perturbation Theory, in which the latter does not determine the interaction energy by computing the total energy of the monomers or dimer. Quantum Theory of Atoms in Molecules and Reduced Density Gradient Non-Covalent Index charge-density-based approaches confirm the F···F contacts are a consequence of attraction by their unified bond path (and bond critical point) and isosurface charge density topologies, respectively. These interactions can be explained neither by the so-called molecular electrostatic surface potential (MESP) model approach that often demonstrates attraction between sites of opposite electrostatic surface potential by means of Coulomb’s law of electrostatics, nor purely by the effect of electrostatic polarization. We provide evidence against the standalone use of this approach and the overlooking of other approaches, as the former does not allow for the calculation of the electrostatic potential on the surfaces of the overlapping atoms on the monomers as in the equilibrium geometry of a complex. This study thus provides unequivocal evidence of the limitation of the MESP approach for its use in gaining insight into the nature of reactivity of overlapped interacting atoms and the intermolecular interactions involved.
Collapse
|
12
|
Petković M, Nakarada Đ, Etinski M. When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms. J Comput Chem 2018; 39:1868-1877. [PMID: 29799128 DOI: 10.1002/jcc.25359] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2018] [Revised: 04/16/2018] [Accepted: 04/27/2018] [Indexed: 01/01/2023]
Abstract
Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra-atomic and inter-atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra-atomic and inter-atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange-correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange-correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| | - Đura Nakarada
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, Belgrade, 11 158, Serbia
| |
Collapse
|