1
|
Martínez-Haya B, Avilés-Moreno JR, Gámez F, Martens J, Oomens J, Berden G. Correlated proton dynamics in hydrogen bonding networks: the benchmark case of 3-hydroxyglutaric acid. Phys Chem Chem Phys 2023; 26:198-208. [PMID: 38053486 DOI: 10.1039/d3cp04514e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
Proton and hydrogen-bonded networks sustain a broad range of structural and charge transfer processes in supramolecular materials. The modelling of proton dynamics is however challenging and demands insights from prototypical benchmark systems. The intramolecular H-bonding networks induced by either protonation or deprotonation of 3-hydroxyglutaric acid provide intriguing case studies of correlated proton dynamics. The vibrational signatures associated with the fluxional proton bonding and its coupling with the hydroxyglutaric backbone are investigated here with infrared action ion spectroscopy experiments and Born-Oppenheimer molecular dynamics (BOMD) computations. Despite the formally similar symmetry of protonated and deprotonated hydroxyglutaric acid, the relative proton affinities of the oxygen centers of the carboxylic and carboxylate groups with respect to that of the central hydroxyl group lead to distinct proton dynamics. In the protonated acid, a tautomeric arrangement of the type HOCO·[HOH]+·OCOH is preferred with the proton binding tighter to the central hydroxyl moiety and the electronic density being shared between the two nearly symmetric H-bonds with the carboxylic end groups. In the deprotonated acid, the asymmetric [OCO]-·HO·HOCO configuration is more stable, with a stronger H-bonding on the bare carboxylate end. Both systems display active backbone dynamics and concerted Grothuss-like proton motions, leading to diffuse band structures in their vibrational spectra. These features are accurately reproduced by the BOMD computations.
Collapse
Affiliation(s)
- Bruno Martínez-Haya
- Center for Nanoscience and Sustainable Technologies (CNATS), Universidad Pablo de Olavide, E-41013 Seville, Spain.
| | - Juan Ramón Avilés-Moreno
- Department of Applied Physical Chemistry, Universidad Autónoma de Madrid, E-28049, Madrid, Spain
| | - Francisco Gámez
- Departamento de Química Física, Universidad Complutense, E-28040 Madrid, Spain
| | - Jonathan Martens
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525ED Nijmegen, The Netherlands
| | - Jos Oomens
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525ED Nijmegen, The Netherlands
| | - Giel Berden
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525ED Nijmegen, The Netherlands
| |
Collapse
|
2
|
Molina FL, Jara-Toro RA, Noble JA, Dedonder-Lardeux C, Jouvet C, Pino GA. Photodetachment of Deprotonated R-Mandelic Acid: The Role of Proton Delocalization on the Radical Stability. Chemphyschem 2023; 24:e202200324. [PMID: 36000956 DOI: 10.1002/cphc.202200324] [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: 05/10/2022] [Revised: 08/23/2022] [Indexed: 01/05/2023]
Abstract
The photodetachment and stability of R-Mandelate, the deprotonated form of the R-Mandelic acid, was investigated by observing the neutral species issued from either simple photodetachment or dissociative photodetachment in a cold anions set-up. R-Mandalate has the possibility to form an intramolecular ionic hydrogen-bond between adjacent hydroxyl and carboxylate groups. The potential energy surface along the proton transfer (PT) coordinate between both groups (O- …H+ …- OCO) features a single local minima, with the proton localized on the O- group (OH…- OCO). However, the structure with the proton localized on the - OCO group (O- …HOCO) is also observed because it falls within the extremity of the vibrational wavefunction of the OH…- OCO isomer along the PT coordinate. The stability of the corresponding radicals, produced upon photodetachment, is strongly dependent on the position of the proton in the anion: the radicals produced from the OH…- OCO isomer decarboxylate without barrier, while the radicals produced from the O- …HOCO isomer are stable.
Collapse
Affiliation(s)
- Franco L Molina
- INFIQC: Instituto de Investigaciones en Fisicoquímica de Córdoba (CONICET - UNC) -, Haya de la Torre s/n, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Departamento de Fisicoquímica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba -, Haya de la Torre y Medina Allende, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Centro Láser de Ciencias Moleculares -, Universidad Nacional de Córdoba, Haya de la Torre s/n, Pabellón Argentina, Ciudad Universitaria, X5000HUA, Córdoba, Argentina
| | - Rafael A Jara-Toro
- INFIQC: Instituto de Investigaciones en Fisicoquímica de Córdoba (CONICET - UNC) -, Haya de la Torre s/n, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Departamento de Fisicoquímica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba -, Haya de la Torre y Medina Allende, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Centro Láser de Ciencias Moleculares -, Universidad Nacional de Córdoba, Haya de la Torre s/n, Pabellón Argentina, Ciudad Universitaria, X5000HUA, Córdoba, Argentina
| | - Jennifer A Noble
- Physique des Interactions Ioniques et Moléculaires (PIIM), CNRS, Aix-Marseille Université, 13397, Marseille, France
| | - Claude Dedonder-Lardeux
- Physique des Interactions Ioniques et Moléculaires (PIIM), CNRS, Aix-Marseille Université, 13397, Marseille, France
| | - Christophe Jouvet
- Physique des Interactions Ioniques et Moléculaires (PIIM), CNRS, Aix-Marseille Université, 13397, Marseille, France
| | - Gustavo A Pino
- INFIQC: Instituto de Investigaciones en Fisicoquímica de Córdoba (CONICET - UNC) -, Haya de la Torre s/n, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Departamento de Fisicoquímica, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba -, Haya de la Torre y Medina Allende, Ciudad Universitaria, X5000HUA, Córdoba, Argentina.,Centro Láser de Ciencias Moleculares -, Universidad Nacional de Córdoba, Haya de la Torre s/n, Pabellón Argentina, Ciudad Universitaria, X5000HUA, Córdoba, Argentina
| |
Collapse
|
3
|
Clavaguéra C, Thaunay F, Ohanessian G. Manifolds of low energy structures for a magic number of hydrated sulfate: SO 42-(H 2O) 24. Phys Chem Chem Phys 2021; 23:24428-24438. [PMID: 34693943 DOI: 10.1039/d1cp03123f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Low energy structures of SO42-(H2O)24 have been obtained using a combination of classical molecular dynamics simulations and refinement of structures and energies by quantum chemical calculations. Extensive exploration of the potential energy surface led to a number of low-energy structures, confirmed by accurate calibration calculations. An overall analysis of this large set was made after devising appropriate structural descriptors such as the numbers of cycles and their combinations. Low energy structures bear common motifs, the most prominent being fused cycles involving alternatively four and six water molecules. The latter adopt specific conformations which ensure the appropriate surface curvature to form a closed cage without dangling O-H bonds and at the same time provide 12-coordination of the sulfate ion. A prominent feature to take into account is isomerism via inversion of hydrogen bond orientations along cycles. This generates large families of ca. 100 isomers for this cluster size, spanning energy windows of 10-30 kJ mol-1. This relatively ignored isomerism must be taken into account to identify reliably the lowest energy minima. The overall picture is that the magic number cluster SO42-(H2O)24 does not correspond to formation of a single, remarkable structure, but rather to a manifold of structural families with similar stabilities. Extensive calculations on isomerization mechanisms within a family indicate that large barriers are associated to direct inversion of hydrogen bond networks. Possible implications of these results for magic number clusters of other anions are discussed.
Collapse
Affiliation(s)
- Carine Clavaguéra
- Institut de Chimie Physique, Université Paris-Saclay - CNRS, UMR 8000, 91405 Orsay, France.
| | - Florian Thaunay
- Laboratoire de Chimie Moléculaire (LCM), CNRS, École Polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France.
| | - Gilles Ohanessian
- Laboratoire de Chimie Moléculaire (LCM), CNRS, École Polytechnique, Institut Polytechnique de Paris, 91120 Palaiseau, France.
| |
Collapse
|
4
|
Thomas DA, Taccone M, Ober K, Mucha E, Meijer G, von Helden G. Helium Nanodroplet Infrared Action Spectroscopy of the Proton-Bound Dimer of Hydrogen Sulfate and Formate: Examining Nuclear Quantum Effects. J Phys Chem A 2021; 125:9279-9287. [PMID: 34652165 PMCID: PMC8558860 DOI: 10.1021/acs.jpca.1c05705] [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] [Indexed: 01/29/2023]
Abstract
The proton-bound dimer of hydrogen sulfate and formate is an archetypal structure for ionic hydrogen-bonding complexes that contribute to biogenic aerosol nucleation. Of central importance for the structure and properties of this complex is the location of the bridging proton connecting the two conjugate base moieties. The potential energy surface for bridging proton translocation features two local minima, with the proton localized at either the formate or hydrogen sulfate moiety. However, electronic structure methods reveal a shallow potential energy surface governing proton translocation, with a barrier on the order of the zero-point energy. This shallow potential complicates structural assignment and necessitates a consideration of nuclear quantum effects. In this work, we probe the structure of this complex and its isotopologues, utilizing infrared (IR) action spectroscopy of ions captured in helium nanodroplets. The IR spectra indicate a structure in which a proton is shared between the hydrogen sulfate and formate moieties, HSO4-···H+···-OOCH. However, because of the nuclear quantum effects and vibrational anharmonicities associated with the shallow potential for proton translocation, the extent of proton displacement from the formate moiety remains unclear, requiring further experiments or more advanced theoretical treatments for additional insight.
Collapse
Affiliation(s)
- Daniel A Thomas
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Martín Taccone
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Katja Ober
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Eike Mucha
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gerard Meijer
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
5
|
Inakollu VS, Geerke DP, Rowley CN, Yu H. Polarisable force fields: what do they add in biomolecular simulations? Curr Opin Struct Biol 2020; 61:182-190. [PMID: 32044671 DOI: 10.1016/j.sbi.2019.12.012] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 12/18/2019] [Accepted: 12/19/2019] [Indexed: 12/11/2022]
Abstract
The quality of biomolecular simulations critically depends on the accuracy of the force field used to calculate the potential energy of the molecular configurations. Currently, most simulations employ non-polarisable force fields, which describe electrostatic interactions as the sum of Coulombic interactions between fixed atomic charges. Polarisation of these charge distributions is incorporated only in a mean-field manner. In the past decade, extensive efforts have been devoted to developing simple, efficient, and yet generally applicable polarisable force fields for biomolecular simulations. In this review, we summarise the latest developments in accounting for key biomolecular interactions with polarisable force fields and applications to address challenging biological questions. In the end, we provide an outlook for future development in polarisable force fields.
Collapse
Affiliation(s)
- Vs Sandeep Inakollu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia
| | - Daan P Geerke
- AIMMS Division of Molecular and Computational Toxicology, Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, the Netherlands.
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador, Canada.
| | - Haibo Yu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia.
| |
Collapse
|