1
|
Gimferrer M, Salvador P. Exact decompositions of the total KS-DFT exchange-correlation energy into one- and two-center terms. J Chem Phys 2023; 158:234105. [PMID: 37326158 DOI: 10.1063/5.0142778] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 05/25/2023] [Indexed: 06/17/2023] Open
Abstract
In the so-called Interacting Quantum Atoms (IQA) approach, the molecular energy is numerically decomposed as a sum of atomic and diatomic contributions. While proper formulations have been put forward for both Hartree-Fock and post-Hartree-Fock wavefunctions, this is not the case for the Kohn-Sham density functional theory (KS-DFT). In this work, we critically analyze the performance of two fully additive approaches for the IQA decomposition of the KS-DFT energy, namely, the one from Francisco et al., which uses atomic scaling factors, and that from Salvador and Mayer based upon the bond order density (SM-IQA). Atomic and diatomic exchange-correlation (xc) energy components are obtained for a molecular test set comprising different bond types and multiplicities and along the reaction coordinate of a Diels-Alder reaction. Both methodologies behave similarly for all systems considered. In general, the SM-IQA diatomic xc components are less negative than the Hartree-Fock ones, which is in good agreement with the known effect of electron correlation upon (most) covalent bonds. In addition, a new general scheme to minimize the numerical error of the sum of two-electron energy contributions (i.e., Coulomb and exact exchange) in the framework of overlapping atoms is described in detail.
Collapse
Affiliation(s)
- Martí Gimferrer
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| | - Pedro Salvador
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| |
Collapse
|
2
|
The IQA Energy Partition in a Drug Design Setting: A Hepatitis C Virus RNA-Dependent RNA Polymerase (NS5B) Case Study. Pharmaceuticals (Basel) 2022; 15:ph15101237. [PMID: 36297349 PMCID: PMC9609620 DOI: 10.3390/ph15101237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Revised: 09/30/2022] [Accepted: 10/01/2022] [Indexed: 11/05/2022] Open
Abstract
The interaction of the thumb site II of the NS5B protein of hepatitis C virus and a pair of drug candidates was studied using a topological energy decomposition method called interacting quantum atoms (IQA). The atomic energies were then processed by the relative energy gradient (REG) method, which extracts chemical insight by computation based on minimal assumptions. REG reveals the most important IQA energy contributions, by atom and energy type (electrostatics, sterics, and exchange–correlation), that are responsible for the behaviour of the whole system, systematically from a short-range ligand–pocket interaction until a distance of approximately 22 Å. The degree of covalency in various key interatomic interactions can be quantified. No exchange–correlation contribution is responsible for the changes in the energy profile of both pocket–ligand systems investigated in the ligand–pocket distances equal to or greater than that of the global minimum. Regarding the hydrogen bonds in the system, a “neighbour effect” was observed thanks to the REG method, which states that a carbon atom would rather not have its covalent neighbour oxygen form a hydrogen bond. The combination of IQA and REG enables the automatic identification of the pharmacophore in the ligands. The coarser Interacting Quantum Fragments (IQF) enables the determination of which amino acids of the pocket contribute most to the binding and the type of energy of said binding. This work is an example of the contribution topological energy decomposition methods can make to fragment-based drug design.
Collapse
|
3
|
Zapata-Acevedo CA, Guevara-Vela JM, Popelier PLA, Rocha Rinza T. Binding Energy Partition of Promising IRAK-4 Inhibitor (Zimlovisertib) for the Treatment of COVID-19 Pneumonia. Chemphyschem 2022; 23:e202200455. [PMID: 36044560 PMCID: PMC9538207 DOI: 10.1002/cphc.202200455] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Revised: 08/19/2022] [Indexed: 01/05/2023]
Abstract
The technique of Fragment-Based Drug Design (FBDD) considers the interactions of different moieties of molecules with biological targets for the rational construction of potential drugs. One basic assumption of FBDD is that the different functional groups of a ligand interact with a biological target in an approximately additive, that is, independent manner. We investigated the interactions of different fragments of ligands and Interleukin-1 Receptor-Associated Kinase 4 (IRAK-4) throughout the FBDD design of Zimlovisertib, a promising anti-inflammatory, currently in trials to be used for the treatment of COVID-19 pneumonia. We utilised state-of-the-art methods of wave function analyses mainly the Interacting Quantum Atoms (IQA) energy partition for this purpose. By means of IQA, we assessed the suitability of every change to the ligand in the five stages of FBDD which led to Zimlovisertib on a quantitative basis. We determined the energetics of the interaction of different functional groups in the ligands with the IRAK-4 protein target and thereby demonstrated the adequacy (or lack thereof) of the changes made across the design of this drug. This analysis permits to verify whether a given alteration of a prospective drug leads to the intended tuning of non-covalent interactions with its protein objective. Overall, we expect that the methods exploited in this paper will prove valuable in the understanding and control of chemical modifications across FBDD processes.
Collapse
Affiliation(s)
- César Arturo Zapata-Acevedo
- Tecnologico de Monterrey: Instituto Tecnologico y de Estudios Superiores de MonterreyChemistryAv. Carlos Lazo 100Santa Fe, La Loma01389Álvaro ObregónMEXICO
| | | | - Paul L. A. Popelier
- UoM: The University of ManchesterChemistryOxford RoadM13 9PLManchesterUNITED KINGDOM
| | - Tomás Rocha Rinza
- Institute Of Chemistry, National Autonomous University of MexicoDepartment of Physical ChemistryCircuito Exterior, Ciudad Universitaria04510Mexico CityMEXICO
| |
Collapse
|
4
|
Jara-Cortés J, Leal-Sánchez E, Francisco E, Pérez-Pimienta JA, Martín Pendás Á, Hernández-Trujillo J. Implementation of the interacting quantum atom energy decomposition using the CASPT2 method. Phys Chem Chem Phys 2021; 23:27508-27519. [PMID: 34874377 DOI: 10.1039/d1cp02837e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Edith Leal-Sánchez
- Departamento de Física y Química Teórica, Facultad de Química, UNAM, México City 04510, Mexico
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | - José A Pérez-Pimienta
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | | |
Collapse
|
5
|
Polestshuk P. TWOE Code: An Efficient Tool for Explicit Partition of Coupled Cluster and Configuration Interaction Energies into Atomic and Diatomic Contributions. J Phys Chem A 2021; 125:3198-3207. [PMID: 33848147 DOI: 10.1021/acs.jpca.1c00675] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The efficient implementation of the TWOE program for evaluating the atomic and interatomic energy components at post-HF level was developed. The systematic convergence of these terms up to a near full-CI limit was performed for the first time for a series of coupled cluster methods: CCSD → CCSDT → CCSDTQ → CCSDTQP. A comparison with corresponding CI approaches (up to fifth excitation level) is additionally discussed. For a set of diatomic systems, it was demonstrated that, along with a full molecular energy convergence, all its components are also converged but with different patterns. It was found that not all components are decreased in their values at increasing computational rank. For instance, atomic energy parts are decreased while interatomic (interaction) energies are increased as the limiting level is approached. Two schemes were employed for atomic partition of molecules: the Baders approach and planes dissection. Influence of dynamical correlation effects on atomic energy components was analyzed in detail. Current TWOE implementation allows one, in principle, to work with any ab initio method providing the two-particle density matrix. It is believed that the developed program will be a useful tool for a real space energy decomposition that helps to reveal the most peculiar points in the structure of the total and correlation energies of a molecule.
Collapse
|
6
|
Fernández-Alarcón A, Guevara-Vela JM, Casals-Sainz JL, Costales A, Francisco E, Martín Pendás Á, Rocha Rinza T. Photochemistry in Real Space: Batho- and Hypsochromism in the Water Dimer. Chemistry 2020; 26:17035-17045. [PMID: 32822523 DOI: 10.1002/chem.202002854] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Indexed: 11/09/2022]
Abstract
The development of chemical intuition in photochemistry faces several difficulties that result from the inadequacy of the one-particle picture, the Born-Oppenheimer approximation, and other basic ideas used to build models. It is shown herein how real-space approaches can be efficiently used to gain valuable insights in photochemistry through a simple example of red and blue shift effects: the double hypso- and bathochromic shifts in the low-lying valence excited states of (H2 O)2 . It is demonstrated that 1) the use of these techniques allows the perturbative language used in the theory of intermolecular interactions, even in the strongly interacting short-range regime, to be maintained; 2) one and only one molecule is photoexcited in each of the addressed excited states and 3) the electrostatic interaction between the in-the-cluster molecular dipoles provides a fairly intuitive rationalisation of the observed batho- and hypsochromism. The methods exploited and illustrated herein are able to maintain the individuality and properties of the interacting entities in a molecular aggregate, and thereby they allow chemical intuition in general states, at any geometry and using a broad variety of electronic structure methods to be kept and built.
Collapse
Affiliation(s)
- Alberto Fernández-Alarcón
- Institute of Chemistry, National Autonomous University of Mexico, 04510, Mexico City, Mexico.,Department of Analytical and Physical Chemistry, University of Oviedo, 33006, Oviedo, Spain
| | | | - José Luis Casals-Sainz
- Department of Analytical and Physical Chemistry, University of Oviedo, 33006, Oviedo, Spain
| | - Aurora Costales
- Department of Analytical and Physical Chemistry, University of Oviedo, 33006, Oviedo, Spain
| | - Evelio Francisco
- Department of Analytical and Physical Chemistry, University of Oviedo, 33006, Oviedo, Spain
| | - Ángel Martín Pendás
- Department of Analytical and Physical Chemistry, University of Oviedo, 33006, Oviedo, Spain
| | - Tomás Rocha Rinza
- Institute of Chemistry, National Autonomous University of Mexico, 04510, Mexico City, Mexico
| |
Collapse
|
7
|
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: 55] [Impact Index Per Article: 13.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
|
8
|
Vincent MA, Silva AF, Popelier PLA. A Comparison of the Interacting Quantum Atoms (IQA) Analysis of the Two‐Particle Density‐Matrices of MP4SDQ and CCSD. Z Anorg Allg Chem 2020. [DOI: 10.1002/zaac.202000169] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Mark A. Vincent
- Manchester Institute of Biotechnology The University of Manchester M1 7DN Manchester UK
- Department of Chemistry The University of Manchester M13 9PL Manchester UK
| | - Arnaldo F. Silva
- Manchester Institute of Biotechnology The University of Manchester M1 7DN Manchester UK
- Department of Chemistry The University of Manchester M13 9PL Manchester UK
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology The University of Manchester M1 7DN Manchester UK
- Department of Chemistry The University of Manchester M13 9PL Manchester UK
| |
Collapse
|
9
|
An Interacting Quantum Atoms (IQA) and Relative Energy Gradient (REG) Study of the Halogen Bond with Explicit Analysis of Electron Correlation. Molecules 2020; 25:molecules25112674. [PMID: 32526931 PMCID: PMC7321288 DOI: 10.3390/molecules25112674] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 06/03/2020] [Accepted: 06/03/2020] [Indexed: 01/18/2023] Open
Abstract
Energy profiles of seven halogen-bonded complexes were analysed with the topological energy partitioning called Interacting Quantum Atoms (IQA) at MP4(SDQ)/6-31+G(2d,2p) level of theory. Explicit interatomic electron correlation energies are included in the analysis. Four complexes combine X2 (X = Cl or F) with HCN or NH3, while the remaining three combine ClF with HCN, NH3 or N2. Each complex was systematically deformed by translating the constituent molecules along its central axis linking X and N, and reoptimising its remaining geometry. The Relative Energy Gradient (REG) method (Theor. Chem. Acc. 2017, 136, 86) then computes which IQA energies most correlate with the total energy during the process of complex formation and further compression beyond the respective equilibrium geometries. It turns out that the covalent energy (i.e., exchange) of the halogen bond, X…N, itself drives the complex formation. When the complexes are compressed from their equilibrium to shorter X…N distance then the intra-atomic energy of N is in charge. When the REG analysis is restricted to electron correlation then the interatomic correlation energy between X and N again drives the complex formation, and the complex compression is best described by the destabilisation of the through-space correlation energy between N and the "outer" halogen.
Collapse
|
10
|
Casals‐Sainz JL, Guevara‐Vela JM, Francisco E, Rocha‐Rinza T, Martín Pendás Á. Efficient implementation of the interacting quantum atoms energy partition of the second‐order Møller–Plesset energy. J Comput Chem 2020; 41:1234-1241. [DOI: 10.1002/jcc.26169] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Revised: 01/30/2020] [Accepted: 01/30/2020] [Indexed: 12/19/2022]
Affiliation(s)
| | - José Manuel Guevara‐Vela
- Institute of Chemistry, National Autonomous University of Mexico, Circuito ExteriorCiudad Universitaria Mexico City Mexico
| | - Evelio Francisco
- Department of Analytical and Physical ChemistryUniversity of Oviedo Oviedo Spain
| | - Tomás Rocha‐Rinza
- Institute of Chemistry, National Autonomous University of Mexico, Circuito ExteriorCiudad Universitaria Mexico City Mexico
| | - Ángel Martín Pendás
- Department of Analytical and Physical ChemistryUniversity of Oviedo Oviedo Spain
| |
Collapse
|
11
|
Vincent MA, Silva AF, Popelier PLA. Atomic Partitioning of the MPn (n = 2, 3, 4) Dynamic Electron Correlation Energy by the Interacting Quantum Atoms Method: A Fast and Accurate Electrostatic Potential Integral Approach. J Comput Chem 2019; 40:2793-2800. [PMID: 31373709 PMCID: PMC6900022 DOI: 10.1002/jcc.26037] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2019] [Revised: 07/10/2019] [Accepted: 07/11/2019] [Indexed: 11/13/2022]
Abstract
Recently, the quantum topological energy partitioning method called interacting quantum atoms (IQA) has been extended to MPn (n = 2, 3, 4) wave functions. This enables the extraction of chemical insight related to dynamic electron correlation. The large computational expense of the IQA-MPn approach is compensated by the advantages that IQA offers compared to older nontopological energy decomposition schemes. This expense is problematic in the construction of a machine learning training set to create kriging models for topological atoms. However, the algorithm presented here markedly accelerates the calculation of atomically partitioned electron correlation energies. Then again, the algorithm cannot calculate pairwise interatomic energies because it applies analytical integrals over whole space (rather than over atomic volumes). However, these pairwise energies are not needed in the quantum topological force field FFLUX, which only uses the energy of an atom interacting with all remaining atoms of the system that it is part of. Thus, it is now feasible to generate accurate and sizeable training sets at MPn level of theory. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Mark A. Vincent
- Manchester Institute of BiotechnologyThe University of ManchesterManchesterM1 7DNUK
- School of ChemistryThe University of ManchesterManchesterM13 9PLUK
| | - Arnaldo F. Silva
- Manchester Institute of BiotechnologyThe University of ManchesterManchesterM1 7DNUK
- School of ChemistryThe University of ManchesterManchesterM13 9PLUK
| | - Paul L. A. Popelier
- Manchester Institute of BiotechnologyThe University of ManchesterManchesterM1 7DNUK
- School of ChemistryThe University of ManchesterManchesterM13 9PLUK
| |
Collapse
|
12
|
Maulén B, Echeverri A, Gómez T, Fuentealba P, Cárdenas C. Electron Localization Function in Excited States: The Case of the Ultrafast Proton Transfer of the Salicylidene Methylamine. J Chem Theory Comput 2019; 15:5532-5542. [DOI: 10.1021/acs.jctc.9b00691] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Boris Maulén
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile
| | - Andrea Echeverri
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile
| | - Tatiana Gómez
- Instituto de Ciencias Químicas Aplicadas, Theoretical and Computational Chemistry Center, Facultad de Ingeniería, Universidad Autónoma de Chile, Av. El Llano Subercaseaux 2801, San Miguel, Santiago, Chile
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile
- Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Avda. Ecuador 3493, Santiago 9170124, Chile
| | - Carlos Cárdenas
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile
- Centro para el Desarrollo de la Nanociencia y la Nanotecnología (CEDENNA), Avda. Ecuador 3493, Santiago 9170124, Chile
| |
Collapse
|
13
|
Tetrel Interactions from an Interacting Quantum Atoms Perspective. Molecules 2019; 24:molecules24122204. [PMID: 31212835 PMCID: PMC6632095 DOI: 10.3390/molecules24122204] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Revised: 06/05/2019] [Accepted: 06/07/2019] [Indexed: 02/06/2023] Open
Abstract
Tetrel bonds, the purportedly non-covalent interaction between a molecule that contains an atom of group 14 and an anion or (more generally) an atom or molecule with lone electron pairs, are under intense scrutiny. In this work, we perform an interacting quantum atoms (IQA) analysis of several simple complexes formed between an electrophilic fragment (A) (CH3F, CH4, CO2, CS2, SiO2, SiH3F, SiH4, GeH3F, GeO2, and GeH4) and an electron-pair-rich system (B) (NCH, NCO-, OCN-, F-, Br-, CN-, CO, CS, Kr, NC-, NH3, OC, OH2, SH-, and N3-) at the aug-cc-pvtz coupled cluster singles and doubles (CCSD) level of calculation. The binding energy ( E bind AB ) is separated into intrafragment and inter-fragment components, and the latter in turn split into classical and covalent contributions. It is shown that the three terms are important in determining E bind AB , with absolute values that increase in passing from electrophilic fragments containing C, Ge, and Si. The degree of covalency between A and B is measured through the real space bond order known as the delocalization index ( δ AB ). Finally, a good linear correlation is found between δ AB and E xc AB , the exchange correlation (xc) or covalent contribution to E bind AB .
Collapse
|
14
|
Cukrowski I. Reliability of HF/IQA, B3LYP/IQA, and MP2/IQA data in interpreting the nature and strength of interactions. Phys Chem Chem Phys 2019; 21:10244-10260. [DOI: 10.1039/c9cp00243j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Using CCSD/BBC1/IQA as a reference to describe interactions in glycol conformers, including O4⋯H6 and O3⋯O4 interactions, B3LYP/IQA data produced an exact qualitative picture and exceptionally reliable IQA energy changes, resulting in a perfectly comparable FAMSEC-based quantitative interpretation.
Collapse
Affiliation(s)
- Ignacy Cukrowski
- Department of Chemistry
- Faculty of Natural and Agricultural Sciences
- University of Pretoria
- Pretoria
- South Africa
| |
Collapse
|
15
|
Fernández-Alarcón A, Casals-Sainz JL, Guevara-Vela JM, Costales A, Francisco E, Martín Pendás Á, Rocha-Rinza T. Partition of electronic excitation energies: the IQA/EOM-CCSD method. Phys Chem Chem Phys 2019; 21:13428-13439. [PMID: 30942218 DOI: 10.1039/c9cp00530g] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Different developments in chemistry and emerging technologies have generated a renewed interest in the properties of molecular excited states. We present herein the partition of black-box, size-consistent equation-of-motion coupled cluster singles and doubles (EOM-CCSD) excitation energies within the framework of the interacting quantum atoms (IQA) formalism. We denote this method as IQA/EOM-CCSD. We illustrate this approach by considering small molecules used often in the study of excited states. This investigation shows how the combination of IQA and EOM-CCSD may provide valuable insights into the molecular changes induced by electron excitation via the real space distribution of the energy of an absorbed photon in a molecular system. Our results reveal (i) the most energetically deformed atomic basins and (ii) the most affected covalent and non-covalent interactions within a molecule due to a given photoexcitation. In other words, this kind of analysis provides insights into the spatial energetic redistribution accompanying an electronic excitation, with interesting foreseeable applications in the rational design of photoexcitations with tailored chemical effects. Altogether, we expect that the IQA/EOM-CCSD excitation energy partition will prove useful in the understanding of systems and processes of interest in photophysics and photochemistry.
Collapse
Affiliation(s)
- Alberto Fernández-Alarcón
- Universidad Nacional Autónoma de México, Instituto de Química, Ciudad Universitaria, Ciudad de México, 04510, Mexico.
| | | | | | | | | | | | | |
Collapse
|
16
|
Tognetti V, Silva AF, Vincent MA, Joubert L, Popelier PLA. Decomposition of Møller–Plesset Energies within the Quantum Theory of Atoms-in-Molecules. J Phys Chem A 2018; 122:7748-7756. [DOI: 10.1021/acs.jpca.8b05357] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Vincent Tognetti
- Normandy University, COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesniére, 76821 Mont St Aignan, Cedex, France
| | - Arnaldo F. Silva
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Mark A. Vincent
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Laurent Joubert
- Normandy University, COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesniére, 76821 Mont St Aignan, Cedex, France
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
17
|
Thacker JCR, Vincent MA, Popelier PLA. Using the Relative Energy Gradient Method with Interacting Quantum Atoms to Determine the Reaction Mechanism and Catalytic Effects in the Peptide Hydrolysis in HIV-1 Protease. Chemistry 2018; 24:11200-11210. [PMID: 29802794 PMCID: PMC6099506 DOI: 10.1002/chem.201802035] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2018] [Indexed: 11/10/2022]
Abstract
The reaction mechanism in an active site is of the utmost importance when trying to understand the role that an enzyme plays in biological processes. In a recently published paper [Theor. Chem. Acc. 2017, 136, 86], we formalised the Relative Energy Gradient (REG) method for automating an Interacting Quantum Atoms (IQA) analysis. Here, the REG method is utilised to determine the mechanism of peptide hydrolysis in the aspartic active site of the enzyme HIV-1 Protease. Using the REG method along with the IQA approach we determine the mechanism of peptide hydrolysis without employing any arbitrary parameters and with remarkable ease (albeit at large computational cost: the system contains 133 atoms, which means that there are 17 689 individual IQA terms to be calculated). When REG and IQA work together it is possible to determine a reaction mechanism at atomistic resolution from data directly derived from quantum calculations, without arbitrary parameters. Moreover, the mechanism determined by this novel method gives concrete insight into how the active site residues catalyse peptide hydrolysis.
Collapse
Affiliation(s)
- Joseph C. R. Thacker
- Manchester Institute of Biotechnology (MIB)131 Princess StreetManchesterM1 7DNUK
- School of ChemistryUniversity of ManchesterOxford RoadManchesterM13 9PLUK
| | - Mark A. Vincent
- Manchester Institute of Biotechnology (MIB)131 Princess StreetManchesterM1 7DNUK
- School of ChemistryUniversity of ManchesterOxford RoadManchesterM13 9PLUK
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB)131 Princess StreetManchesterM1 7DNUK
- School of ChemistryUniversity of ManchesterOxford RoadManchesterM13 9PLUK
| |
Collapse
|
18
|
Pastorczak E, Corminboeuf C. Perspective: Found in translation: Quantum chemical tools for grasping non-covalent interactions. J Chem Phys 2018; 146:120901. [PMID: 28388098 DOI: 10.1063/1.4978951] [Citation(s) in RCA: 70] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Today's quantum chemistry methods are extremely powerful but rely upon complex quantities such as the massively multidimensional wavefunction or even the simpler electron density. Consequently, chemical insight and a chemist's intuition are often lost in this complexity leaving the results obtained difficult to rationalize. To handle this overabundance of information, computational chemists have developed tools and methodologies that assist in composing a more intuitive picture that permits better understanding of the intricacies of chemical behavior. In particular, the fundamental comprehension of phenomena governed by non-covalent interactions is not easily achieved in terms of either the total wavefunction or the total electron density, but can be accomplished using more informative quantities. This perspective provides an overview of these tools and methods that have been specifically developed or used to analyze, identify, quantify, and visualize non-covalent interactions. These include the quantitative energy decomposition analysis schemes and the more qualitative class of approaches such as the Non-covalent Interaction index, the Density Overlap Region Indicator, or quantum theory of atoms in molecules. Aside from the enhanced knowledge gained from these schemes, their strengths, limitations, as well as a roadmap for expanding their capabilities are emphasized.
Collapse
Affiliation(s)
- Ewa Pastorczak
- Laboratory for Computational Molecular Design, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
19
|
Rodríguez-Mayorga M, Via-Nadal M, Solà M, Ugalde JM, Lopez X, Matito E. Electron-Pair Distribution in Chemical Bond Formation. J Phys Chem A 2018; 122:1916-1923. [PMID: 29381071 DOI: 10.1021/acs.jpca.7b12556] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The chemical formation process has been studied from relaxation holes, Δh(u), resulting from the difference between the radial intracule density and the nonrelaxed counterpart, which is obtained from atomic radial intracule densities and the pair density constructed from the overlap of the atomic densities. Δh(u) plots show that the internal reorganization of electron pairs prior to bond formation and the covalent bond formation from electrons in separate atoms are completely recognizable processes from the shape of the relaxation hole, Δh(u). The magnitude of Δh(u), the shape of Δh(u) ∀ u < Req, and the distance between the minimum and the maximum in Δh(u) provide further information about the nature of the chemical bond formed. A computational affordable approach to calculate the radial intracule density from approximate pair densities has been also suggested, paving the way to study electron-pair distributions in larger systems.
Collapse
Affiliation(s)
- M Rodríguez-Mayorga
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC). P.K. 1072 , 20080 Donostia, Euskadi, Spain.,Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, University of Girona , C/ Maria Aurèlia Capmany, 69, 17003 Girona, Catalonia, Spain
| | - M Via-Nadal
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC). P.K. 1072 , 20080 Donostia, Euskadi, Spain
| | - M Solà
- Institut de Química Computacional i Catàlisi (IQCC) and Departament de Química, University of Girona , C/ Maria Aurèlia Capmany, 69, 17003 Girona, Catalonia, Spain
| | - J M Ugalde
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC). P.K. 1072 , 20080 Donostia, Euskadi, Spain
| | - X Lopez
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC). P.K. 1072 , 20080 Donostia, Euskadi, Spain
| | - E Matito
- Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC). P.K. 1072 , 20080 Donostia, Euskadi, Spain.,IKERBASQUE, Basque Foundation for Science , 48013 Bilbao, Euskadi, Spain
| |
Collapse
|
20
|
Romero-Montalvo E, Guevara-Vela JM, Costales A, Pendás ÁM, Rocha-Rinza T. Cooperative and anticooperative effects in resonance assisted hydrogen bonds in merged structures of malondialdehyde. Phys Chem Chem Phys 2018; 19:97-107. [PMID: 27722301 DOI: 10.1039/c6cp04877c] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We analyzed non-additive effects in resonance assisted hydrogen bonds (RAHBs) in different β-enolones, which are archetypal compounds of these types of interactions. For this purpose, we used (i) potential energy curves to compute the formation energy, ΔE, of the RAHBs of interest in different circumstances along with (ii) tools offered by quantum chemical topology, namely, the Quantum Theory of Atoms In Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) electronic energy partition. We established the effect that a given H-bond exerts over ΔE associated with another RAHB, determining in this way the cooperativity or the anticooperativity of these interactions. The mesomeric structures and the QTAIM delocalisation indices are consistent with the determined cooperative or anticooperative character of two given RAHBs. The HB cooperativity and anticooperativity studied herein are directly reflected in the IQA interaction energy E, but they are modulated by the surrounding hydrocarbon chain. The IQA decomposition of ΔEcoop, a measure of the cooperativity between a pair of interacting RAHBs, indicates that the analyzed H-bond cooperative/anticooperative effects are associated with greater/smaller (i) strengthening of the pseudo-bicyclic structure of the compounds of interest and (ii) electron localisations with its corresponding changes in the intra and intermolecular exchange-correlation contributions to ΔE. Overall, we expect that this investigation will provide valuable insights into the interplay among hydrogen bonded atoms and the π system in RAHBs contributing in this way to the understanding of the general features of H-bonds.
Collapse
Affiliation(s)
- Eduardo Romero-Montalvo
- Institute of Chemistry, National Autonomous University of Mexico, Circuito Exterior, Ciudad Universitaria, Delegación Coyoacán C.P. 04510, Mexico City, Mexico.
| | | | - Aurora Costales
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006, Oviedo, Spain
| | - Ángel Martín Pendás
- 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. 04510, Mexico City, Mexico.
| |
Collapse
|
21
|
McDonagh JL, Silva AF, Vincent MA, Popelier PLA. Machine Learning of Dynamic Electron Correlation Energies from Topological Atoms. J Chem Theory Comput 2017; 14:216-224. [PMID: 29211469 DOI: 10.1021/acs.jctc.7b01157] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
We present an innovative method for predicting the dynamic electron correlation energy of an atom or a bond in a molecule utilizing topological atoms. Our approach uses the machine learning method Kriging (Gaussian Process Regression with a non-zero mean function) to predict these dynamic electron correlation energy contributions. The true energy values are calculated by partitioning the MP2 two-particle density-matrix via the Interacting Quantum Atoms (IQA) procedure. To our knowledge, this is the first time such energies have been predicted by a machine learning technique. We present here three important proof-of-concept cases: the water monomer, the water dimer, and the van der Waals complex H2···He. These cases represent the final step toward the design of a full IQA potential for molecular simulation. This final piece will enable us to consider situations in which dispersion is the dominant intermolecular interaction. The results from these examples suggest a new method by which dispersion potentials for molecular simulation can be generated.
Collapse
Affiliation(s)
- James L McDonagh
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain
| | - Arnaldo F Silva
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain
| | - Mark A Vincent
- School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, Great Britain.,School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
22
|
Casalz-Sainz JL, Guevara-Vela JM, Francisco E, Rocha-Rinza T, Martín Pendás Á. Where Does Electron Correlation Lie? Some Answers from a Real Space Partition. Chemphyschem 2017; 18:3553-3561. [DOI: 10.1002/cphc.201700940] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2017] [Revised: 10/04/2017] [Indexed: 11/09/2022]
Affiliation(s)
- José Luis Casalz-Sainz
- Departament of Analytical and Physical Chemistry; University of Oviedo; E-33006 Oviedo Spain
| | | | - Evelio Francisco
- Departament 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.; 04510 Mexico City Mexico
| | - Ángel Martín Pendás
- Departament of Analytical and Physical Chemistry; University of Oviedo; E-33006 Oviedo Spain
| |
Collapse
|
23
|
Maxwell PI, Popelier PLA. Unfavorable regions in the ramachandran plot: Is it really steric hindrance? The interacting quantum atoms perspective. J Comput Chem 2017; 38:2459-2474. [PMID: 28841241 PMCID: PMC5659141 DOI: 10.1002/jcc.24904] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2017] [Revised: 07/10/2017] [Accepted: 07/18/2017] [Indexed: 01/06/2023]
Abstract
Accurate description of the intrinsic preferences of amino acids is important to consider when developing a biomolecular force field. In this study, we use a modern energy partitioning approach called Interacting Quantum Atoms to inspect the cause of the φ and ψ torsional preferences of three dipeptides (Gly, Val, and Ile). Repeating energy trends at each of the molecular, functional group, and atomic levels are observed across both (1) the three amino acids and (2) the φ/ψ scans in Ramachandran plots. At the molecular level, it is surprisingly electrostatic destabilization that causes the high-energy regions in the Ramachandran plot, not molecular steric hindrance (related to the intra-atomic energy). At the functional group and atomic levels, the importance of key peptide atoms (Oi-1 , Ci , Ni , Ni+1 ) and some sidechain hydrogen atoms (Hγ ) are identified as responsible for the destabilization seen in the energetically disfavored Ramachandran regions. Consistently, the Oi-1 atoms are particularly important for the explanation of dipeptide intrinsic behavior, where electrostatic and steric destabilization unusually complement one another. The findings suggest that, at least for these dipeptides, it is the peptide group atoms that dominate the intrinsic behavior, more so than the sidechain atoms. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Peter I. Maxwell
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford RoadManchesterGreat BritainM13 9PL
| | - Paul L. A. Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain and School of Chemistry, University of Manchester, Oxford RoadManchesterGreat BritainM13 9PL
| |
Collapse
|
24
|
McDonagh JL, Vincent MA, Popelier PL. Partitioning dynamic electron correlation energy: Viewing Møller-Plesset correlation energies through Interacting Quantum Atom (IQA) energy partitioning. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.09.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
25
|
Ruiz I, Matito E, Holguín-Gallego FJ, Francisco E, Martín Pendás Á, Rocha-Rinza T. Fermi and Coulomb correlation effects upon the interacting quantum atoms energy partition. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1957-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|