1
|
Gómez S, Flórez E, Acelas N, Cappelli C, Hadad C, Restrepo A. Encapsulation of charged halogens by the 5 12 water cage. Phys Chem Chem Phys 2024; 26:15426-15436. [PMID: 38747303 DOI: 10.1039/d4cp01340a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Elizabeth Flórez
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Nancy Acelas
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | - Cacier Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| |
Collapse
|
2
|
Gómez S, Gómez S, Rojas-Valencia N, Hernández JG, Ardila-Fierro KJ, Gómez T, Cárdenas C, Hadad C, Cappelli C, Restrepo A. Interactions and reactivity in crystalline intermediates of mechanochemical cyclorhodation reactions. Phys Chem Chem Phys 2024; 26:2228-2241. [PMID: 38165158 DOI: 10.1039/d3cp04201d] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
There is experimental evidence that solid mixtures of the rhodium dimer [Cp*RhCl2]2 and benzo[h] quinoline (BHQ) produce two different polymorphic molecular cocrystals called 4α and 4β under ball milling conditions. The addition of NaOAc to the mixture leads to the formation of the rhodacycle [Cp*Rh-(BHQ)Cl], where the central Rh atom retains its tetracoordinate character. Isolate 4β reacts with NaOAc leading to the same rhodacycle while isolate 4α does not under the same conditions. We show that the puzzling difference in reactivity between the two cocrystals can be traced back to fundamental aspects of the intermolecular interactions between the BHQ and [Cp*RhCl2]2 fragments in the crystalline environment. To support this view, we report a number of descriptors of the nature and strength of chemical bonds and intermolecular interactions in the extended solids and in a cluster model. We calculate formal quantum mechanical descriptors based on electronic structure, electron density, and binding and interaction energies including an energy decomposition analysis. Without exception, all descriptors point to 4β being a transient structure higher in energy than 4α with larger local and global electrophilic and nucleophilic powers, a more favorable spatial and energetic distribution of the frontier orbitals, and a more fragile crystal structure.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Santiago Gómez
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Natalia Rojas-Valencia
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - José G Hernández
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Karen J Ardila-Fierro
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Avenida Pedro de Valdivia 425, Santiago, Chile
| | - Carlos Cárdenas
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile
- Centro para el desarrollo de las Nanociencias y Nanotecnología, CEDENNA, Av. Ecuador 3493, Santiago, Chile
| | - Cacier Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| |
Collapse
|
3
|
Rojas-Valencia N, Gómez S, Giovannini T, Cappelli C, Restrepo A, Núñez Zarur F. Water Maintains the UV-Vis Spectral Features During the Insertion of Anionic Naproxen and Ibuprofen into Model Cell Membranes. J Phys Chem B 2023; 127:2146-2155. [PMID: 36877579 DOI: 10.1021/acs.jpcb.2c08332] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
UV-vis spectra of anionic ibuprofen and naproxen in a model lipid bilayer of the cell membrane are investigated using computational techniques in combination with a comparative analysis of drug spectra in purely aqueous environments. The simulations aim at elucidating the intricacies behind the negligible changes in the maximum absorption wavelength in the experimental spectra. A set of configurations of the systems constituted by lipid, water, and drugs or just water and drugs are obtained from classical Molecular Dynamics simulations. UV-vis spectra are computed in the framework of atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) approaches together with Time-Dependent Density Functional Theory (TD-DFT). Our results suggest that the molecular orbitals involved in the electronic transitions are the same, regardless of the chemical environment. A thorough analysis of the contacts between the drug and water molecules reveals that no significant changes in UV-vis spectra are a consequence of ibuprofen and naproxen molecules being permanently microsolvated by water molecules, despite the presence of lipid molecules. Water molecules microsolvate the charged carboxylate group as expected but also microsolvate the aromatic regions of the drugs.
Collapse
Affiliation(s)
- Natalia Rojas-Valencia
- Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026, Medellín, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia, UdeA, Calle 70 No. 52-21 050010, Medellín, Colombia
| | - Francisco Núñez Zarur
- Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026, Medellín, Colombia
| |
Collapse
|
4
|
Gómez S, Gómez S, David J, Guerra D, Cappelli C, Restrepo A. Dissecting Bonding Interactions in Cysteine Dimers. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27248665. [PMID: 36557799 PMCID: PMC9786917 DOI: 10.3390/molecules27248665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Revised: 12/01/2022] [Accepted: 12/02/2022] [Indexed: 12/13/2022]
Abstract
Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.
Collapse
Affiliation(s)
- Santiago Gómez
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín 050010, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Jorge David
- Escuela de Ciencias y Humanidades, Departamento de Ciencias Básicas, Universidad Eafit, AA 3300, Medellín 050022, Colombia
| | - Doris Guerra
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín 050010, Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Correspondence: (C.C.); (A.R.)
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín 050010, Colombia
- Correspondence: (C.C.); (A.R.)
| |
Collapse
|
5
|
Flórez E, Gómez S, Acelas N, Hadad C, Restrepo A. Microsolvation versus Encapsulation in Mono, Di, and Trivalent Cations. Chemphyschem 2022; 23:e202200456. [PMID: 35962558 DOI: 10.1002/cphc.202200456] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 08/04/2022] [Indexed: 01/05/2023]
Abstract
The effects of the formal charge in the stability and bonding of water cavities when solvating a cation are studied here using [X(H2 O)20 ]q+ clusters starting with the well known 512 isomer of (water)20 , placing a single mono, di, or trivalent Xq+ cation at the interior, and then optimizing and characterizing the resulting clusters. Highly correlated interaction and deformation energies are calculated using the CCSD(T)-DLPNO formalism. Bonding interactions are characterized using the tools provided by the quantum theory of atoms in molecules, natural bond orbitals, and non-covalent surfaces. Our results indicate that water to water hydrogen bonds are sensibly strengthened resulting in strong cooperative effects, which amount to ≈ 2 ${ \approx 2}$ kcal/mol per hydrogen bond in the bare cavity and to larger values for the systems including the cations. Approximate encapsulation, that is, surrounding the cation by a network of hydrogen bonds akin to the well known methane clathrate seems to be preferred by cations with smaller charge densities while microsolvation, that is, cluster structures having explicit X⋯O contacts seem to be preferred by cations with larger charge densities which severely deform the cavity.
Collapse
Affiliation(s)
- Elizabeth Flórez
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, Medellín, 050026, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Nancy Acelas
- Grupo de Materiales con Impacto, Mat&mpac. Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, Medellín, 050026, Colombia
| | - Cacier Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| |
Collapse
|
6
|
Kenouche S, Martínez-Araya JI. A combined QTAIM/IRI topological analysis of the effect of axial/equatorial positions of NH 2 and CN substituents in the [(PY 5Me 2)MoO] + complex. J Mol Graph Model 2022; 116:108273. [PMID: 35930821 DOI: 10.1016/j.jmgm.2022.108273] [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: 12/24/2021] [Revised: 05/09/2022] [Accepted: 07/08/2022] [Indexed: 12/15/2022]
Abstract
By means of the Interaction Region Indicator (IRI) and Quantum Theory of Atoms in Molecules (QTAIM), the influence exerted by NH2 (amino) and CN (cyano) as electron donor and electron acceptor substituent groups, respectively, located at para-positions of axial and equatorial pyridine rings of derivatized complexes coming from the [(PY5Me2)MoO]+ complex during the hydrogen molecular release in the gas phase was analyzed. In any case, a H-H covalent bond is forming at the transition state, with a strengthening of the electron density of 5.5% when the substituent group involved is NH2 at the para-position of the axial pyridine ring. However, there was no difference between NH2 and CN when these substituent groups are located at the para-positions of the equatorial pyridine rings. The topological properties of electron densities from the QTAIM are not perturbed by the electron donor and electron acceptor nature of the substituents, even when these substituent groups are located at the axial or equatorial pyridine rings of the Mo-based complex.
Collapse
Affiliation(s)
- Samir Kenouche
- Group of Modeling of Chemical Systems Using Quantum Calculations, Applied Chemistry Laboratory (LCA), University M. Khider of Biskra, 07000 Biskra, Algeria
| | - Jorge I Martínez-Araya
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andres Bello (UNAB), Av. República 275, 8370146, Santiago, Chile.
| |
Collapse
|
7
|
Moreno N, Hadad CZ, Restrepo A. Microsolvation of electrons by a handful of ammonia molecules. J Chem Phys 2022; 157:134301. [PMID: 36209021 DOI: 10.1063/5.0107245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Microsolvation of electrons in ammonia is studied here via anionic NH3 n - clusters with n = 2-6. Intensive samplings of the corresponding configurational spaces using second-order perturbation theory with extended basis sets uncover rich and complex energy landscapes, heavily populated by many local minima in tight energy windows as calculated from highly correlated coupled cluster methods. There is a marked energetical preference for structures that place the excess electron external to the molecular frame, effectively coordinating it with the three protons from a single ammonia molecule. Overall, as the clusters grow in size, the lowest energy dimer serves as the basic motif over which additional ammonia molecules are attached via unusually strong charge-assisted hydrogen bonds. This is a priori quite unexpected because, on electrostatic grounds, the excess electron would be expected to be in contact with as many protons as possible. Accordingly, a full quantum mechanical treatment of the bonding interactions under the tools provided by the quantum theory of atoms in molecules is carried out in order to dissect and understand the nature of intermolecular contacts. Vertical detachment energies reveal bound electrons even for n = 2.
Collapse
Affiliation(s)
- Norberto Moreno
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Cacier Z Hadad
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| |
Collapse
|
8
|
Arias A, Gómez S, Rojas-Valencia N, Núñez-Zarur F, Cappelli C, Murillo-López JA, Restrepo A. Formation and evolution of C-C, C-O, C[double bond, length as m-dash]O and C-N bonds in chemical reactions of prebiotic interest. RSC Adv 2022; 12:28804-28817. [PMID: 36320504 PMCID: PMC9549586 DOI: 10.1039/d2ra06000k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 09/27/2022] [Indexed: 11/18/2022] Open
Abstract
A series of prebiotic chemical reactions yielding the precursor building blocks of amino acids, proteins and carbohydrates, starting solely from HCN and water is studied here. We closely follow the formation and evolution of the pivotal C-C, C-O, C[double bond, length as m-dash]O, and C-N bonds, which dictate the chemistry of the molecules of life. In many cases, formation of these bonds is set in motion by proton transfers in which individual water molecules act as catalysts so that water atoms end up in the products. Our results indicate that the prebiotic formation of carbon dioxide, formaldehyde, formic acid, formaldimine, glycolaldehyde, glycine, glycolonitrile, and oxazole derivatives, among others, are best described as highly nonsynchronous concerted single step processes. Nonetheless, for all reactions involving double proton transfer, the formation and breaking of O-H bonds around a particular O atom occur in a synchronous fashion, apparently independently from other primitive processes. For the most part, the first process to initiate seems to be the double proton transfer in the reactions where they are present, then bond breaking/formation around the reactive carbon in the carbonyl group and finally rupture of the C-N bonds in the appropriate cases, which are the most reluctant to break. Remarkably, within the limitations of our non-dynamical computational model, the wide ranges of temperature and pressure in which these reactions occur, downplay the problematic determination of the exact constraints on the early Earth.
Collapse
Affiliation(s)
- Alejandro Arias
- Instituto de Química, Universidad de Antioquia UdeACalle 70 No. 52-21MedellínColombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di ScienzePiazza dei Cavalieri 7Pisa56126Italy
| | - Natalia Rojas-Valencia
- Instituto de Química, Universidad de Antioquia UdeACalle 70 No. 52-21MedellínColombia,Escuela de Ciencias y Humanidades, Departamento de Ciencias Básicas, Universidad EafitMedellínAA 3300Colombia
| | - Francisco Núñez-Zarur
- Facultad de Ciencias Básicas, Universidad de MedellínCarrera 87 No. 30-65Medellín050026Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di ScienzePiazza dei Cavalieri 7Pisa56126Italy
| | - Juliana A. Murillo-López
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andres BelloAutopista, Concepción-TalcahuanoTalcahuano 7100Chile
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeACalle 70 No. 52-21MedellínColombia
| |
Collapse
|
9
|
Correa E, Montaño D, Restrepo A. Cation ⋯anion bonding interactions in 1–Ethyl–3–Methylimidazolium based ionic liquids. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111648] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
10
|
Analysis of Conformational Preferences in Caffeine. Molecules 2022; 27:molecules27061937. [PMID: 35335301 PMCID: PMC8949453 DOI: 10.3390/molecules27061937] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 12/13/2021] [Accepted: 12/13/2021] [Indexed: 12/10/2022] Open
Abstract
High level DLPNO−CCSD(T) electronic structure calculations with extended basis sets over B3LYP−D3 optimized geometries indicate that the three methyl groups in caffeine overcome steric hindrance to adopt uncommon conformations, each one placing a C−H bond on the same plane of the aromatic system, leading to the C−H bonds eclipsing one carbonyl group, one heavily delocalized C−N bond constituent of the fused double ring aromatic system, and one C−H bond from the imidazole ring. Deletion of indiscriminate and selective non-Lewis orbitals unequivocally show that hyperconjugation in the form of a bidirectional −CH3 ⇆ aromatic system charge transfer is responsible for these puzzling conformations. The structural preferences in caffeine are exclusively determined by orbital interactions, ruling out electrostatics, induction, bond critical points, and density redistribution because the steric effect, the allylic effect, the Quantum Theory of Atoms in Molecules (QTAIM), and the non-covalent interactions (NCI), all predict wrong energetic orderings. Tiny rotational barriers, not exceeding 1.3 kcal/mol suggest that at room conditions, each methyl group either acts as a free rotor or adopts fluxional behavior, thus preventing accurate determination of their conformations. In this context, our results supersede current experimental ambiguity in the assignation of methyl conformation in caffeine and, more generally, in methylated xanthines and their derivatives.
Collapse
|
11
|
Porto CM, Santana LC, Morgon NH. Theoretical investigation of the cooperative effect of solvent: a case study. Phys Chem Chem Phys 2022; 24:14603-14615. [DOI: 10.1039/d1cp05388d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effect of solvent was investigated at the DFT level, M06-2X/6-31++G(d,p), for the implicit, namely the universal solvent model based on solute electron density (SMD), and hybrid solvation models, and...
Collapse
|
12
|
David J, Gómez S, Guerra D, Guerra D, Restrepo A. A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers. Chemphyschem 2021; 22:2401-2412. [PMID: 34554628 DOI: 10.1002/cphc.202100585] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Revised: 09/22/2021] [Indexed: 12/14/2022]
Abstract
High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.
Collapse
Affiliation(s)
- Jorge David
- Escuela de Ciencias, Departamento de Ciencias Físicas, Universidad Eafit, AA 3300, Medellín, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Doris Guerra
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Dario Guerra
- Departamento de Educación y Ciencias Básicas, Instituto Tecnológico Metropolitano, Calle 73 No. 76 A-354, Medellín, Colombia
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| |
Collapse
|
13
|
Florez E, Acelas N, Gomez S, Hadad C, Restrepo A. To be or not to be? that is the entropic, enthalpic, and molecular interaction dilemma in the formation of (water)20 clusters and methane clathrate. Chemphyschem 2021; 23:e202100716. [PMID: 34761856 DOI: 10.1002/cphc.202100716] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Revised: 11/08/2021] [Indexed: 11/09/2022]
Abstract
A detailed analysis under a comprehensive set of theoretical and computational tools of the thermodynamical factors and of the intermolecular interactions behind the stabilization of a well known set of (water)20 cavities and of the methane clathrate is offered in this work. Beyond the available reports of experimental characterization at extreme conditions of most of the systems studied here, all clusters should be amenable to experimental detection at 1 atm and moderate temperatures since 280 K marks the boundary at which, ignoring reaction paths, formation of all clusters is no longer spontaneous from the 20H2O → (H2O)20 and CH4 + 20H2O → CH4@512 processes. As a function of temperature, a complex interplay leading to the free energy of formation occurs between the destabilizing entropic contributions, mostly due to cluster vibrations, and the stabilizing enthalpic contributions, due to intermolecular interactions and the PV term, is best illustrated by the highly symmetric 512 cage consistently showing signs of stronger intermolecular bonding despite having smaller binding energy than the other clusters. A fluxional wall of attractive non-covalent interactions, arising because of the cumulative effect of a large number of tiny individual charge transfers to the interstitial region, plays a pivotal role stabilizing the CH4@512 clathrate.
Collapse
Affiliation(s)
- Elizabeth Florez
- Universidad de Medellín: Universidad de Medellin, Ciencias basicas, COLOMBIA
| | - Nancy Acelas
- Universidad de Medellín: Universidad de Medellin, Ciencias Basicas, COLOMBIA
| | - Sara Gomez
- Scuola Normale Superiore Classe di Scienze, Chemistry, ITALY
| | - Cacier Hadad
- Universidad de Antioquía: Universidad de Antioquia, Chemistry, COLOMBIA
| | - Albeiro Restrepo
- Universidad de Antioquia, Chemistry, AA 1226, 00000, Medellin, COLOMBIA
| |
Collapse
|
14
|
Uribe L, Gómez S, Giovannini T, Egidi F, Restrepo A. An efficient and robust procedure to calculate absorption spectra of aqueous charged species applied to NO 2. Phys Chem Chem Phys 2021; 23:14857-14872. [PMID: 34223573 DOI: 10.1039/d1cp00652e] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Accurate calculation of absorption spectra of aqueous NO2- requires rigorously sampling the quantum potential energy surfaces for microsolvation of NO2- with at least five explicit water molecules and embedding the resulting clusters in a continuum solvent accounting for the statistical weighted contributions of individual isomers. This method, which we address as ASCEC + PCM, introduces several desired features when compared against MD simulations derived QM/MM spectra: comparatively fewer explicit solvent molecules to be treated with expensive QM methods, the identification of equilibrium structures in the quantum PES to be used in further vibrational spectroscopy, and the unequivocal identification of cluster orbitals undergoing electronic transitions and charge transfer that originate the spectral bands.
Collapse
Affiliation(s)
- Lina Uribe
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Franco Egidi
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy.
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia.
| |
Collapse
|
15
|
Gómez S, Rojas-Valencia N, Gómez SA, Cappelli C, Merino G, Restrepo A. A molecular twist on hydrophobicity. Chem Sci 2021; 12:9233-9245. [PMID: 34276953 PMCID: PMC8261874 DOI: 10.1039/d1sc02673a] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2021] [Accepted: 05/19/2021] [Indexed: 12/22/2022] Open
Abstract
A thorough exploration of the molecular basis for hydrophobicity with a comprehensive set of theoretical tools and an extensive set of organic solvent S/water binary systems is discussed in this work. Without a single exception, regardless of the nature or structure of S, all quantum descriptors of bonding yield stabilizing S⋯water interactions, therefore, there is no evidence of repulsion and thus no reason for etymological hydrophobicity at the molecular level. Our results provide molecular insight behind the exclusion of S molecules by water, customarily invoked to explain phase separation and the formation of interfaces, in favor of a complex interplay between entropic, enthalpic, and dynamic factors. S⋯water interfaces are not just thin films separating the two phases; instead, they are non-isotropic regions with density gradients for each component whose macroscopic stability is provided by a large number of very weak dihydrogen contacts. We offer a definition of interface as the region in which the density of the components in the A/B binary system is not constant. At a fundamental level, our results contribute to better current understanding of hydrophobicity. Notwithstanding the very weak nature of individual contacts, it is the cumulative effect of a large number of interactions (green NCI surfaces) which provides macroscopic stability to the interfaces.![]()
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze Piazza dei Cavalieri 7 56126 Pisa Italy
| | - Natalia Rojas-Valencia
- Instituto de Química, Universidad de Antioquia UdeA Calle 70 No. 52-21 Medellín Colombia .,Escuela de Ciencias y Humanidades, Departamento de Ciencias Básicas, Universidad Eafit AA 3300 Medellín Colombia
| | - Santiago A Gómez
- Instituto de Química, Universidad de Antioquia UdeA Calle 70 No. 52-21 Medellín Colombia
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze Piazza dei Cavalieri 7 56126 Pisa Italy
| | - Gabriel Merino
- Departamento de Física Aplicada, Centro de Investigación y de Estudios Avanzados Unidad Mérida. Km 6 Antigua Carretera a Progreso. Apdo. Postal 73, Cordemex 97310 Mérida Yucatan Mexico
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA Calle 70 No. 52-21 Medellín Colombia
| |
Collapse
|
16
|
Kaviani S, Shahab S, Sheikhi M, Khaleghian M, Al Saud S. Characterization of the binding affinity between some anti-Parkinson agents and Mn2+, Fe3+ and Zn2+ metal ions: A DFT insight. INORG CHEM COMMUN 2021. [DOI: 10.1016/j.inoche.2021.108582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
17
|
|
18
|
Tobón P, Gómez S, Restrepo A, Núñez-Zarur F. Role of Substrate Substituents in Alkene Metathesis Mediated by a Ru Alkylidene Catalyst. Organometallics 2021. [DOI: 10.1021/acs.organomet.0c00482] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Paola Tobón
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín 50010, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín 50010, Colombia
| | - Francisco Núñez-Zarur
- Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026 Medellín, Colombia
| |
Collapse
|