1
|
Magalhães IOM, Cabral BJC, Martins JBL. Ab Initio Approach to the Structure, Vibrational Properties, and Electron Binding Energies of H 2S∙∙∙SO 2. Molecules 2023; 28:6656. [PMID: 37764431 PMCID: PMC10535889 DOI: 10.3390/molecules28186656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Revised: 09/07/2023] [Accepted: 09/14/2023] [Indexed: 09/29/2023] Open
Abstract
The present study employs high-level ab initio calculations to investigate the structure, vibrational frequencies, and electronic properties of H2S∙∙∙SO2. The analysis of vibrational frequencies reveals an intramolecular vibrational energy transfer phenomenon, where energy from the stretching modes of H2S is transferred to the ν1s mode of SO2. At the CCSD(T)/aug-cc-pVQZ level, the interaction energy between H2S and SO2 is predicted to be 2.78 kcal/mol. Electron propagator theory calculations yield a HOMO-LUMO gap of 8.24 eV for H2S∙∙∙SO2. Furthermore, by utilizing ab initio results for the adiabatic ionization energy and electron affinity, the electrophilicity of H2S∙∙∙SO2 is estimated to be 2.01 eV. This value is similar to the electrophilicity of SO2, suggesting comparable reactivity and chemical behavior. The non-covalent interaction (NCI) analysis of the H2S∙∙∙SO2 complex emphasizes the significant contribution of non-covalent van der Waals interactions in its energetic stabilization.
Collapse
Affiliation(s)
- Isaac O. M. Magalhães
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia 70910900, DF, Brazil; (I.O.M.M.); (B.J.C.C.)
| | - Benedito J. C. Cabral
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia 70910900, DF, Brazil; (I.O.M.M.); (B.J.C.C.)
- Biosystems and Integrative Sciences Institute, BioISI, Faculdade de Ciências de Lisboa, 1749-016 Lisboa, Portugal
| | - João B. L. Martins
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia 70910900, DF, Brazil; (I.O.M.M.); (B.J.C.C.)
| |
Collapse
|
2
|
Nascimento LA, Nascimento ÉCM, Martins JBL. In silico study of tacrine and acetylcholine binding profile with human acetylcholinesterase: docking and electronic structure. J Mol Model 2022; 28:252. [PMID: 35947248 DOI: 10.1007/s00894-022-05252-2] [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: 03/14/2022] [Accepted: 07/25/2022] [Indexed: 11/30/2022]
Abstract
Alzheimer disease (AD) is a neurodegenerative process, one of the most common and incident dementia in the population over 60 years. AD manifests the presence of complex biochemical processes involved in neuronal degeneration, such as the formation of senile plaques containing amyloid-β peptides, the development of intracellular neurofibrillary tangles, and the suppression of the acetylcholine neurotransmitter. In this way, we performed a set of theoretical tests of tacrine ligand and acetylcholine neurotransmitter against the human acetylcholinesterase enzyme. Molecular docking was used to understand the most important interactions of these molecules with the enzyme. Computational chemistry calculation was carried out using MP2, DFT, and semi-empirical methods, starting from molecular docking structures. We have also performed studies regarding the non-covalent interactions, electron localization function, molecular electrostatic potential and explicit water molecule influence. For Trp86 residue, we show two main interactions in accordance to the results of the literature for TcAChE. First, intermolecular interactions of the cation-π and sigma-π type were found. Second, close stacking interactions were stablished between THA+ and Trp86 residue on one side and with Tyr337 residue on the other side.
Collapse
Affiliation(s)
- Letícia A Nascimento
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia, DF, 70910-900, Brazil
| | - Érica C M Nascimento
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia, DF, 70910-900, Brazil
| | - João B L Martins
- Computational Chemistry Laboratory, Institute of Chemistry, University of Brasilia, Brasilia, DF, 70910-900, Brazil.
| |
Collapse
|
3
|
Martins JBL, Quintino RP, Politi JRDS, Sethio D, Gargano R, Kraka E. Computational analysis of vibrational frequencies and rovibrational spectroscopic constants of hydrogen sulfide dimer using MP2 and CCSD(T). SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 239:118540. [PMID: 32502813 DOI: 10.1016/j.saa.2020.118540] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2020] [Revised: 05/18/2020] [Accepted: 05/24/2020] [Indexed: 06/11/2023]
Abstract
Previous studies have shown that the weakly bonded H2S dimer demands high level quantum chemical calculations to reproduce experimental values. We investigated the hydrogen bonding of H2S dimer using MP2 and CCSD(T) levels of theory in combination with aug-cc-pV(D,T,Q)Z basis sets. More precisely, the binding energies, potential energy curves, rovibrational spectroscopic constants, decomposition lifetime, and normal vibrational frequencies were calculated. In addition, we introduced the local mode analysis of Konkoli-Cremer to quantify the hydrogen bonding in the H2S dimer as well as providing for the first time the comprehensive decomposition of normal vibrational modes into local modes contributions, and a decomposition lifetime based on rate constant. The local mode force constant of the H2S dimer hydrogen bond is smaller than that of the water dimer, in accordance with the weaker hydrogen bonding in the H2S dimer.
Collapse
Affiliation(s)
- João B L Martins
- Institute of Chemistry, University of Brasília, Brasília, DF 70910-900, Brazil.
| | - Rabeshe P Quintino
- Institute of Chemistry, University of Brasília, Brasília, DF 70910-900, Brazil
| | - José R Dos S Politi
- Institute of Chemistry, University of Brasília, Brasília, DF 70910-900, Brazil
| | - Daniel Sethio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, United States
| | - Ricardo Gargano
- Institute of Physics, University of Brasília, Brasília, DF 70910-900, Brazil
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, United States
| |
Collapse
|
4
|
Jouypazadeh H, Farrokhpour H, Solimannejad M. Evaluation of one-dimensional potential energy surfaces for prediction of spectroscopic properties of hydrogen bonds in linear bonded complexes. J Mol Model 2017; 23:157. [PMID: 28397087 DOI: 10.1007/s00894-017-3336-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2017] [Accepted: 03/27/2017] [Indexed: 10/19/2022]
Abstract
This work evaluated the reliability of the one-dimensional potential energy surface for calculating the spectroscopic properties (rovibrational constants and rotational line energies) of hydrogen bonds in linear bonded complexes by comparing theoretical results with the corresponding experimental results. For this purpose, two hydrogen bonded complexes were selected: the HCN···HCN homodimer and the HCN···HF heterodimer. The one-dimensional potential energy surfaces related to the hydrogen bonds in these complexes were calculated using different computational methods and basis sets. The calculated potential curve of each complex was fitted to an analytical one-dimensional potential function to obtain the potential parameters. The obtained analytical potential function of each complex was used in a two-particle Schrödinger equation to obtain the rovibrational energy levels of the hydrogen bond. Using the calculated rovibrational levels, the rovibrational spectra and constants of each complex were calculated and compared with experimental data available from the literature. Compared with experimental data, the calculated one-dimensional potential energy surface at the QCISD/aug-cc-pVDZ level of theory was found to predict the spectroscopic properties of hydrogen bonds better than the potential curves obtained using other computational methods, especially for the HCN···HCN homodimer complex. Generally, the results obtained for the HCN···HCN homodimer complex were closer to experimental data than those obtained for the HCN···HF heterodimer complex. The investigation performed in this work showed that the one-dimensional potential curve related to the hydrogen bond between two linear molecules can be used to predict the spectroscopic constants of hydrogen bonds. Graphical abstract Potential energy curves of HCN···HCN and HCN···HF complexes calculated at the different computational levels.
Collapse
Affiliation(s)
- Hamidreza Jouypazadeh
- Department of Chemistry, Isfahan University of Technology, Isfahan, 84156-83111, Iran.
| | - Hossein Farrokhpour
- Department of Chemistry, Isfahan University of Technology, Isfahan, 84156-83111, Iran.
| | - Mohammad Solimannejad
- Department of Chemistry, Faculty of Sciences, Arak University, Arak, 38156-8-8349, Iran
| |
Collapse
|
5
|
Bonfim VS, Borges NM, Martins JBL, Gargano R, Politi JRDS. Quantum Monte Carlo with density matrix: potential energy curve derived properties. J Mol Model 2017; 23:104. [DOI: 10.1007/s00894-017-3272-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Accepted: 02/06/2017] [Indexed: 11/30/2022]
|
6
|
Solimannejad M, Jouypazadeh H, Farrokhpour H. Rovibrational energy and spectroscopic constant calculations of complexes pairing via dihydrogen bonds. J Mol Model 2015; 21:119. [PMID: 25877527 DOI: 10.1007/s00894-015-2662-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2015] [Accepted: 03/23/2015] [Indexed: 11/25/2022]
Abstract
In the present investigation, we performed a thorough study of potential energy curves, rovibrational spectra, and spectroscopic constants for complexes pairing via dihydrogen bonds. In particular, we dealt with LiH···HX (X = F, CN, CCH, CCF, CCCl) complexes by employing accurate electronic energy calculations at the MP2/aug-cc-pVDZ level of theory. Following this, the Numerov method was applied to solve the nuclear Schrödinger equation, thus obtaining spectroscopic constants and rovibrational spectra. Good linear correlation between the magnitudes of the interaction energies for interaction of HX with LiH, and the most positive electrostatic potentials of hydrogen in HX, was established.
Collapse
Affiliation(s)
- Mohammad Solimannejad
- Department of Chemistry, Faculty of Sciences, Arak University, Arak, 38156-8-8349, Iran,
| | | | | |
Collapse
|