1
|
Kwan V, Consta S, Malek SMA. Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit. J Phys Chem B 2024; 128:193-207. [PMID: 38127582 DOI: 10.1021/acs.jpcb.3c05701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.
Collapse
Affiliation(s)
- Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, ON, Canada N6A 5B7
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, ON, Canada N6A 5B7
| | - Shahrazad M A Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's, NL, Canada A1B 3X7
| |
Collapse
|
2
|
Oliveira A, Filipe HAL, Ramalho JP, Salvador A, Geraldes CFGC, Moreno MJ, Loura LMS. Modeling Gd 3+ Complexes for Molecular Dynamics Simulations: Toward a Rational Optimization of MRI Contrast Agents. Inorg Chem 2022; 61:11837-11858. [PMID: 35849762 PMCID: PMC9775472 DOI: 10.1021/acs.inorgchem.2c01597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The correct parametrization of lanthanide complexes is of the utmost importance for their characterization using computational tools such as molecular dynamics simulations. This allows the optimization of their properties for a wide range of applications, including medical imaging. Here we present a systematic study to establish the best strategies for the correct parametrization of lanthanide complexes using [Gd(DOTA)]- as a reference, which is used as a contrast agent in MRI. We chose the bonded model to parametrize the lanthanide complexes, which is especially important when considering the study of the complex as a whole (e.g., for the study of the dynamics of its interaction with proteins or membranes). We followed two strategies: a so-called heuristic approach employing strategies already published by other authors and another based on the more recent MCPB.py tool. Adjustment of the Lennard-Jones parameters of the metal was required. The final topologies obtained with both strategies were able to reproduce the experimental ion to oxygen distance, vibrational frequencies, and other structural properties. We report a new strategy to adjust the Lennard-Jones parameters of the metal ion in order to capture dynamic properties such as the residence time of the capping water (τm). For the first time, the correct assessment of the τm value for Gd-based complexes was possible by recording the dissociative events over up to 10 μs all-atom simulations. The MCPB.py tool allowed the accurate parametrization of [Gd(DOTA)]- in a simpler procedure, and in this case, the dynamics of the water molecules in the outer hydration sphere was also characterized. This sphere was divided into the first hydration layer, an intermediate region, and an outer hydration layer, with a residence time of 18, 10 and 19 ps, respectively, independent of the nonbonded parameters chosen for Gd3+. The Lennard-Jones parameters of Gd3+ obtained here for [Gd(DOTA)]- may be used with similarly structured gadolinium MRI contrast agents. This allows the use of molecular dynamics simulations to characterize and optimize the contrast agent properties. The characterization of their interaction with membranes and proteins will permit the design of new targeted contrast agents with improved pharmacokinetics.
Collapse
Affiliation(s)
- Alexandre
C. Oliveira
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal
| | - Hugo A. L. Filipe
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,CPIRN-IPG-Center
of Potential and Innovation of Natural Resources, Polytechnic Institute of Guarda, 6300-559 Guarda, Portugal
| | - João P.
Prates Ramalho
- Hercules
Laboratory, LAQV, REQUIMTE, Department of Chemistry, School of Science
and Technology, University of Évora, 7000-671 Évora, Portugal
| | - Armindo Salvador
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,CNC−Center
for Neuroscience and Cell Biology, University
of Coimbra, P-3004-517 Coimbra, Portugal,Institute
for Interdisciplinary Research - University of Coimbra, Casa Costa Alemão- Polo II, Rua D. Francisco de Lemos, 3030-789 Coimbra, Portugal
| | - Carlos F. G. C. Geraldes
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,Department
of Life Sciences, University of Coimbra, Calçada Martim de Freitas, 3000-393 Coimbra, Portugal,CIBIT/ICNAS
- Instituto de Ciências Nucleares Aplicadas à Saúde, Pólo das Ciências
da Saúde, Azinhaga de Santa Comba, 3000-548 Coimbra, Portugal
| | - Maria João Moreno
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,Department
of Chemistry, University of Coimbra, 3004-535 Coimbra, Portugal,
| | - Luís M. S. Loura
- Coimbra
Chemistry Center - Institute of Molecular Sciences (CQC-IMS), University of Coimbra, 3004-535 Coimbra, Portugal,Faculty
of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal,
| |
Collapse
|
3
|
Phan LX, Lynch CI, Crain J, Sansom MS, Tucker SJ. Influence of effective polarization on ion and water interactions within a biomimetic nanopore. Biophys J 2022; 121:2014-2026. [DOI: 10.1016/j.bpj.2022.05.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 03/25/2022] [Accepted: 05/04/2022] [Indexed: 11/16/2022] Open
|
4
|
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
|
5
|
Malek SMA, Kwan V, Saika-Voivod I, Consta S. Low Density Interior in Supercooled Aqueous Nanodroplets Expels Ions to the Subsurface. J Am Chem Soc 2021; 143:13113-13123. [PMID: 34375522 DOI: 10.1021/jacs.1c04142] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interaction between water and ions within droplets plays a key role in the chemical reactivity of atmospheric and man-made aerosols. Here we report direct computational evidence that in supercooled aqueous nanodroplets a lower density core of tetrahedrally coordinated water expels the cosmotropic ions to the denser and more disordered subsurface. In contrast, at room temperature, depending on the nature of the ion, the radial distribution in the droplet core is nearly uniform or elevated toward the center. We analyze the spatial distribution of a single ion in terms of a reference electrostatic model. The energy of the system in the analytical model is expressed as the sum of the electrostatic and surface energy of a deformable droplet. The model predicts that the ion is subject to a harmonic potential centered at the droplet's center of mass. We name this effect "electrostatic confinement". The model's predictions are consistent with the simulation findings for a single ion at room temperature but not at supercooling. We anticipate this study to be the starting point for investigating the structure of supercooled (electro)sprayed droplets that are used to preserve the conformations of macromolecules originating from the bulk solution.
Collapse
Affiliation(s)
- Shahrazad M A Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada
| | - Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Ivan Saika-Voivod
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's A1B 3X7, Canada.,Department of Applied Mathematics, Western University, London, Ontario N6A 3K7, Canada
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
6
|
Das B, Mondal S, Chandra A. Two-Dimensional Infrared Spectroscopy of Aqueous Solutions of Metal Nitrates: Slowdown of Spectral Diffusion in the Presence of Divalent Cations. J Phys Chem B 2020; 124:7391-7404. [PMID: 32790404 DOI: 10.1021/acs.jpcb.0c03471] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The hydrogen-bonded network of water can be affected both structurally and dynamically by the presence of ions. In the present study, we have considered three aqueous solutions of metal nitrates to investigate the effects of divalent cations (Mg2+ and Ca2+), compared to that of monovalent Na+ ions, on hydrogen-bond fluctuations and vibrational spectral diffusion through calculations of linear and two-dimensional infrared spectra of these solutions at room temperature. We have employed the methods of molecular dynamics simulations using effective polarizable models of ions combined with quantum mechanical calculations of transition variables and statistical mechanical calculations of spectral response functions of vibrational spectroscopy. Divalent cations are found to have much stronger and longer-ranged effects on the structure and dynamics of the hydrogen-bonded network than that induced by the monovalent sodium ions. The blue shifts in the calculated linear spectra are found to follow the Hofmeister trend for the cations. The 2D-IR spectral lineshape and intensity corresponding to three-pulse echo peak shift (3PEPS) experiments are calculated. The timescales of these nonlinear spectral responses and also frequency-time correlations show significant slowing down of spectral diffusion for solutions containing divalent Mg2+ and Ca2+ ions compared to the corresponding dynamics of the solution containing Na+ ions. Unlike NaNO3 solution, the relaxation of frequency and dipole orientational fluctuations of anion-bound water in Mg(NO3)2 and Ca(NO3)2 solutions are found to be somewhat slower than bulk water, which can be attributed to the presence of divalent cations whose effects go beyond their first solvation shells. This is also seen in the dynamics of bulk water in these solutions which is found to be notably slower for the solutions containing divalent cations than that in the NaNO3 solution. Unlike Mg2+ and Ca2+ ions, no specific cationic effect is observed for the Na+ ions.
Collapse
Affiliation(s)
- Banshi Das
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Subhadip Mondal
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
7
|
Le Barbu-Debus K, Bowles J, Jähnigen S, Clavaguéra C, Calvo F, Vuilleumier R, Zehnacker A. Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches. Phys Chem Chem Phys 2020; 22:26047-26068. [DOI: 10.1039/d0cp03869e] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.
Collapse
Affiliation(s)
- Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| | - Jessica Bowles
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Sascha Jähnigen
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Carine Clavaguéra
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Florent Calvo
- Université Grenoble Alpes
- CNRS
- LiPhy
- F-38000 Grenoble
- France
| | - Rodolphe Vuilleumier
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| |
Collapse
|
8
|
Kwan V, Malevanets A, Consta S. Where Do the Ions Reside in a Highly Charged Droplet? J Phys Chem A 2019; 123:9298-9310. [DOI: 10.1021/acs.jpca.9b03368] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Victor Kwan
- Department of Chemistry, The University of Western Ontario, London N6A 5B7, Ontario, Canada
| | - Anatoly Malevanets
- Department of Electrical and Computer Engineering, The University of Western Ontario, London N6A 5B9, Ontario, Canada
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London N6A 5B7, Ontario, Canada
| |
Collapse
|
9
|
Thaunay F, Jana C, Clavaguéra C, Ohanessian G. Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops. J Phys Chem A 2018; 122:832-842. [PMID: 29266957 DOI: 10.1021/acs.jpca.7b10554] [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
Hydrated ions are ubiquitous in environmental and biological media. Understanding the perturbation exerted by an ion on the water hydrogen bond network is possible in the nanodrop regime by recording vibrational spectra in the O-H bond stretching region. This has been achieved experimentally in recent years by forming gaseous ions containing tens to hundreds of water molecules and recording their infrared photodissociation spectra. In this paper, we demonstrate the capabilities of a modeling strategy based on an extension of the AMOEBA polarizable force field to implement water atomic charge fluctuations along with those of intramolecular structure along the dynamics. This supplementary flexibility of nonbonded interactions improves the description of the hydrogen bond network and, therefore, the spectroscopic response. Finite temperature IR spectra are obtained from molecular dynamics simulations by computing the Fourier transform of the dipole moment autocorrelation function. Simulations of 1-2 ns are required for extensive sampling in order to reproduce the experimental spectra. Furthermore, bands are assigned with the driven molecular dynamics approach. This method package is shown to compare successfully with experimental spectra for 11 ions in water drops containing 36-100 water molecules. In particular, band frequency shifts of the free O-H stretching modes at the cluster surface are well reproduced as a function of both ion charge and drop size.
Collapse
Affiliation(s)
- Florian Thaunay
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| | - Chandramohan Jana
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| | - Carine Clavaguéra
- Laboratoire de Chimie Physique, Université Paris Sud - CNRS, Université Paris Saclay , 15, avenue Jean Perrin, 91405 Orsay Cedex, France
| | - Gilles Ohanessian
- LCM, CNRS, Ecole Polytechnique, Université Paris Saclay , 91128 Palaiseau, France
| |
Collapse
|