1
|
Thorat A, Chauhan R, Sartape R, Singh MR, Shah JK. Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol. J Phys Chem B 2024; 128:3707-3719. [PMID: 38572661 PMCID: PMC11033864 DOI: 10.1021/acs.jpcb.3c08480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 03/08/2024] [Accepted: 03/18/2024] [Indexed: 04/05/2024]
Abstract
Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion-ion and ion-solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K+ ion. Diffusion-based Nernst-Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.
Collapse
Affiliation(s)
- Amey Thorat
- School
of Chemical Engineering, Oklahoma State
University, Stillwater, Oklahoma 74078, United States
| | - Rohit Chauhan
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Rohan Sartape
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Meenesh R. Singh
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Jindal K. Shah
- School
of Chemical Engineering, Oklahoma State
University, Stillwater, Oklahoma 74078, United States
| |
Collapse
|
2
|
Mamatkulov S, Polák J, Razzokov J, Tomaník L, Slavíček P, Dzubiella J, Kanduč M, Heyda J. Unveiling the Borohydride Ion through Force-Field Development. J Chem Theory Comput 2024; 20:1263-1273. [PMID: 38227434 PMCID: PMC10867804 DOI: 10.1021/acs.jctc.3c01020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Revised: 12/15/2023] [Accepted: 12/18/2023] [Indexed: 01/17/2024]
Abstract
The borohydride ion, BH4-, is an essential reducing agent in many technological processes, yet its full understanding has been elusive, because of at least two significant challenges. One challenge arises from its marginal stability in aqueous solutions outside of basic pH conditions, which considerably limits the experimental thermodynamic data. The other challenge comes from its unique and atypical hydration shell, stemming from the negative excess charge on its hydrogen atoms, which complicates the accurate modeling in classical atomistic simulations. In this study, we combine experimental and computer simulation techniques to devise a classical force field for NaBH4 and deepen our understanding of its characteristics. We report the first measurement of the ion's activity coefficient and extrapolate it to neutral pH conditions. Given the difficulties in directly measuring its solvation free energies, owing to its instability, we resort to quantum chemistry calculations. This combined strategy allows us to derive a set of nonpolarizable force-field parameters for the borohydride ion for classical molecular dynamics simulations. The derived force field simultaneously captures the solvation free energy, the hydration structure, as well as the activity coefficient of NaBH4 salt across a broad concentration range. The obtained insights into the hydration shell of the BH4- ion are crucial for accurately modeling and understanding its interactions with other molecules, ions, materials, and interfaces.
Collapse
Affiliation(s)
- Shavkat Mamatkulov
- Institute
of Material Science of AS, Ch.Aytmatov str.2B, 100084 Tashkent, Uzbekistan
| | - Jakub Polák
- Department
of Physical Chemistry, University of Chemistry
and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Jamoliddin Razzokov
- Institute
of Fundamental and Applied Research, National
Research University TIIAME, Kori Niyoziy 39, 100000 Tashkent, Uzbekistan
- School
of Engineering, Akfa University, Milliy Bog Street 264, 111221 Tashkent, Uzbekistan
| | - Lukáš Tomaník
- Department
of Physical Chemistry, University of Chemistry
and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Petr Slavíček
- Department
of Physical Chemistry, University of Chemistry
and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Joachim Dzubiella
- Applied
Theoretical Physics–Computational Physics, Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder-Str. 3, D-79104 Freiburg, Germany
| | - Matej Kanduč
- Jožef
Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
| | - Jan Heyda
- Department
of Physical Chemistry, University of Chemistry
and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| |
Collapse
|
3
|
Blazquez S, Bourg IC, Vega C. Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2. J Chem Phys 2024; 160:046101. [PMID: 38270237 DOI: 10.1063/5.0186233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Accepted: 12/27/2023] [Indexed: 01/26/2024] Open
Abstract
In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Ian C Bourg
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, USA
- High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, USA
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
4
|
Tesei G, Hsiao YW, Dabkowska A, Grönberg G, Yanez Arteta M, Ulkoski D, Bray DJ, Trulsson M, Ulander J, Lund M, Lindfors L. Lipid shape and packing are key for optimal design of pH-sensitive mRNA lipid nanoparticles. Proc Natl Acad Sci U S A 2024; 121:e2311700120. [PMID: 38175863 PMCID: PMC10786277 DOI: 10.1073/pnas.2311700120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 11/27/2023] [Indexed: 01/06/2024] Open
Abstract
The ionizable-lipid component of RNA-containing nanoparticles controls the pH-dependent behavior necessary for an efficient delivery of the cargo-the so-called endosomal escape. However, it is still an empirical exercise to identify optimally performing lipids. Here, we study two well-known ionizable lipids, DLin-MC3-DMA and DLin-DMA using a combination of experiments, multiscale computer simulations, and electrostatic theory. All-atom molecular dynamics simulations, and experimentally measured polar headgroup pKa values, are used to develop a coarse-grained representation of the lipids, which enables the investigation of the pH-dependent behavior of lipid nanoparticles (LNPs) through Monte Carlo simulations, in the absence and presence of RNA molecules. Our results show that the charge state of the lipids is determined by the interplay between lipid shape and headgroup chemistry, providing an explanation for the similar pH-dependent ionization state observed for lipids with headgroup pKa values about one-pH-unit apart. The pH dependence of lipid ionization is significantly influenced by the presence of RNA, whereby charge neutrality is achieved by imparting a finite and constant charge per lipid at intermediate pH values. The simulation results are experimentally supported by measurements of α-carbon 13C-NMR chemical shifts for eGFP mRNA LNPs of both DLin-MC3-DMA and DLin-DMA at various pH conditions. Further, we evaluate the applicability of a mean-field Poisson-Boltzmann theory to capture these phenomena.
Collapse
Affiliation(s)
- Giulio Tesei
- Structural Biology and NMR Laboratory & The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, CopenhagenDK-2200, Denmark
- Department of Chemistry, Division of Computational Chemistry, Lund University, LundSE-221 00, Sweden
| | - Ya-Wen Hsiao
- The Hartree Centre, Science and Technology Facilities Council (STFC) Daresbury Laboratory, WarringtonWA4 4AD, United Kingdom
| | - Aleksandra Dabkowska
- Advanced Drug Delivery, Pharmaceutical Sciences, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| | - Gunnar Grönberg
- Medicinal Chemistry, Early Respiratory & Immunology, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| | - Marianna Yanez Arteta
- Advanced Drug Delivery, Pharmaceutical Sciences, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| | - David Ulkoski
- Advanced Drug Delivery, Pharmaceutical Sciences, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| | - David J. Bray
- The Hartree Centre, Science and Technology Facilities Council (STFC) Daresbury Laboratory, WarringtonWA4 4AD, United Kingdom
| | - Martin Trulsson
- Department of Chemistry, Division of Computational Chemistry, Lund University, LundSE-221 00, Sweden
| | - Johan Ulander
- Data Science and Modelling, Pharmaceutical Sciences, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| | - Mikael Lund
- Department of Chemistry, Division of Computational Chemistry, Lund University, LundSE-221 00, Sweden
| | - Lennart Lindfors
- Advanced Drug Delivery, Pharmaceutical Sciences, BioPharmaceuticals Research and Development, AstraZeneca, Gothenburg, Mölndal431 83, Sweden
| |
Collapse
|
5
|
Sun Y, Zhang D, Bashir A, Li C, Fan Z. Scaling Solute-Solvent Distances to Improve Solubility and Ion Paring Predictions in Rigid Ion Models. J Phys Chem B 2023; 127:9575-9586. [PMID: 37906589 DOI: 10.1021/acs.jpcb.3c05993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
Force fields based on the rigid ion model (RIM) have been developed to accurately predict the various physical and chemical properties of salts and water. However, the combined use of these models often fails to accurately predict the solubility of salts in water. To address this issue, several approaches, such as charge scaling or reparameterization, have been proposed. Nevertheless, these methods require laborious reparameterization of nonbonded force field parameters. In this article, we propose a scaling solute-solvent distance (SSSD) method to improve force fields in predicting salt solubility without changing the solute-solute and solvent-solvent interactions in the original force fields. This method can also tune the ion pairing of salt in water. One main advantage of the SSSD method is that reparameterization of the crystal and water models is not needed. We use two RIMs for the NaCl-water system (JC-SPC/E and SD-SPC/E) and the CHARMM force field for the KCl-water system to demonstrate the improved accuracy in predicting solubility by the SSSD method. Furthermore, we use the RDG-SPC/Fw force field to show that the SSSD method can also be used to tune the ion pairing of CaCO3 in water. Limitations of this method are also discussed.
Collapse
Affiliation(s)
- Ying Sun
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Dan Zhang
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- Department of Chemistry, Shanghai University, Shanghai 200444, P. R. China
| | - Ayesha Bashir
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Chuncheng Li
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
| | - Zhaochuan Fan
- Suzhou Institute of Nano-Tech and Nano-Bionics, Chinese Academy of Sciences, Suzhou 215123, P. R. China
- School of Nano-Tech and Nano-Bionics, University of Science and Technology of China, Hefei 230026, P. R. China
| |
Collapse
|
6
|
Gao Y, Li M, Zhan C, Zhang H, Yin M, Lu W, Xu B. A Nanoconfined Water-Ion Coordination Network for Flexible Energy-Dissipation Devices. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2303759. [PMID: 37410996 DOI: 10.1002/adma.202303759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2023] [Revised: 07/03/2023] [Accepted: 07/05/2023] [Indexed: 07/08/2023]
Abstract
Water-ion interaction in a nanoconfined environment that deeply constrains spatial freedoms of local atomistic motion with unconventional coupling mechanisms beyond that in a free, bulk state is essential to spark designs of a broad spectrum of nanofluidic devices with unique properties and functionalities. Here, it is reported that the interaction between ions and water molecules in a hydrophobic nanopore forms a coordination network with an interaction density that is nearly fourfold that of the bulk counterpart. Such strong interaction facilitates the connectivity of the water-ion network and is uncovered by corroborating the formation of ion clusters and the reduction of particle dynamics. A liquid-nanopore energy-dissipation system is designed and demonstrated in both molecular simulations and experiments that the formed coordination network controls the outflow of confined electrolytes along with a pressure reduction, capable of providing flexible protection for personnel and devices and instrumentations against external mechanical impact and attack.
Collapse
Affiliation(s)
- Yuan Gao
- State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan, 430074, China
- Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA, 22904, USA
| | - Mingzhe Li
- Department of Civil and Environmental Engineering, Michigan State University, East Lansing, MI, 48824, USA
| | - Chi Zhan
- Department of Civil and Environmental Engineering, Michigan State University, East Lansing, MI, 48824, USA
| | - Haozhe Zhang
- Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA, 22904, USA
| | - Mengtian Yin
- Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA, 22904, USA
| | - Weiyi Lu
- Department of Civil and Environmental Engineering, Michigan State University, East Lansing, MI, 48824, USA
- Department of Chemical Engineering and Materials Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Baoxing Xu
- Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA, 22904, USA
| |
Collapse
|
7
|
Blazquez S, Abascal JLF, Lagerweij J, Habibi P, Dey P, Vlugt TJH, Moultos OA, Vega C. Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property. J Chem Theory Comput 2023; 19:5380-5393. [PMID: 37506381 PMCID: PMC10448725 DOI: 10.1021/acs.jctc.3c00562] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Indexed: 07/30/2023]
Abstract
In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.
Collapse
Affiliation(s)
- Samuel Blazquez
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L. F. Abascal
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jelle Lagerweij
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Parsa Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Carlos Vega
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
8
|
Bendová M, Heyda J, Wagner Z, Feder-Kubis J, Polák J, Tankam T, Sýkorová A. Aqueous solutions of chiral ionic liquids based on (–)-menthol: An experimental and computational study of volumetric and transport properties. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2023]
|
9
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
10
|
Effect of hydrated shell layers on surface tension of electrolyte solutions: Insights from interpretable machine learning. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120887] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
11
|
McIvor JAP, Larsen DS, Mercadante D. Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM. J Phys Chem B 2022; 126:7833-7846. [PMID: 36125334 DOI: 10.1021/acs.jpcb.2c03837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the αR region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.
Collapse
Affiliation(s)
- Jordan A P McIvor
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Danaé S Larsen
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, 23 Symonds Street, Auckland 1010, New Zealand
| |
Collapse
|
12
|
Pethes I. Towards the correct microscopic structure of aqueous CsCl solutions with a comparison of classical interatomic potential models. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
13
|
Borkowski AK, Thompson WH. Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions. J Phys Chem B 2022; 126:6700-6712. [PMID: 36004804 DOI: 10.1021/acs.jpcb.2c03957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The Hofmeister series is frequently used to rank ions based on their behavior from chaotropes ("structure breakers"), which weaken the surrounding hydrogen-bond network, to kosmotropes ("structure makers"), which enhance it. Here, we use fluctuation theory to investigate the energetic and entropic driving forces underlying the Hofmeister series for aqueous alkali-halide solutions. Specifically, we exploit the OH stretch infrared (IR) spectrum in isotopically dilute HOD/D2O solutions as a probe of the effect of the salt on the water properties for different concentrations and choice of halide anion. Fluctuation theory is used to calculate the temperature derivative of these IR spectra, including decomposition of the derivative into different energetic contributions. These contributions are used to determine the thermodynamic driving forces in terms of effective internal energy and entropic contributions. This analysis implicates entropic contributions as the key factor in the Hofmeister series behavior of the OH stretch IR spectra, while the effective internal energy is nearly ion-independent.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
14
|
Fan K, Zhang Y, Qiu Y, Zhang H. Impacts of targeting different hydration free energy references on the development of ion potentials. Phys Chem Chem Phys 2022; 24:16244-16262. [PMID: 35758314 DOI: 10.1039/d2cp01237e] [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
Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H+). The choice of ΔG°(H+) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K+, Mg2+, Al3+, and Cl- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H+) for the cation and an increase in ΔG°(H+) for the anion. The increase in ΔG°(H+) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H+) value of -259.8 kcal mol-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H+) for Mg2+ model is better for a reasonable description of the metal binding network.
Collapse
Affiliation(s)
- Kun Fan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| |
Collapse
|
15
|
Singh O, Das BK, Chakraborty D. Influence of ion specificity and concentration on the conformational transition of intrinsically disordered sheep prion peptide. Chemphyschem 2022; 23:e202200211. [PMID: 35621322 DOI: 10.1002/cphc.202200211] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Revised: 05/26/2022] [Indexed: 11/12/2022]
Abstract
The structural sensitivity of the IDPs with the ions has been observed experimentally; however, it is still unclear how the presence of different metal ions affects structural stability. We performed atomistic molecular dynamics simulation of sheep prion peptide (142-167) in presence of different monovalent, divalent ions at various concentrations to find out the effect of the size, charge, and ionic concentration on the structure of the peptide. It is found that Li + ions have a higher survival probability compared to Na + , K + and Mg2 + affecting the solvation structure of the protein leading to the alpha-helix structure. At high concentration, due to the increase in the ion-solvent and ion-counter interactions, the effect of the ions is screened on the surface of the protein and hence no ion specificity is observed. This study demonstrates how ions can be used to regulate the protein structure and function that can help in designing drugs.
Collapse
Affiliation(s)
- Omkar Singh
- National Institute of Technology Karnataka, Chemistry, INDIA
| | | | - Debashree Chakraborty
- National Institute of Technology Karnataka, Chemistry, Department of Chemistry, Science Block, NIT K Surathkal, 575025, Mangalore, INDIA
| |
Collapse
|
16
|
Venetsanos F, Anogiannakis SD, Theodorou DN. Mixing Thermodynamics and Flory–Huggins Interaction Parameter of Polyethylene Oxide/Polyethylene Oligomeric Blends from Kirkwood–Buff Theory and Molecular Simulations. Macromolecules 2022. [DOI: 10.1021/acs.macromol.2c00642] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fotis Venetsanos
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Stefanos D. Anogiannakis
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| | - Doros N. Theodorou
- School of Chemical Engineering, National Technical University of Athens, 9 Heroon Polytechniou Street, Athens 15780, Greece
| |
Collapse
|
17
|
Grotz KK, Schwierz N. Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E. J Chem Phys 2022; 156:114501. [DOI: 10.1063/5.0087292] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.
Collapse
Affiliation(s)
- Kara K. Grotz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| | - Nadine Schwierz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| |
Collapse
|
18
|
Blazquez S, Conde MM, Abascal JLF, Vega C. The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: Extension to the ions F−, Br−, I−, Rb+, and Cs+. J Chem Phys 2022; 156:044505. [DOI: 10.1063/5.0077716] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- S. Blazquez
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M. M. Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - J. L. F. Abascal
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
19
|
Sevilla M, Cortes-Huerto R. Connecting density fluctuations and Kirkwood–Buff integrals for finite-size systems. J Chem Phys 2022; 156:044502. [DOI: 10.1063/5.0076744] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Mauricio Sevilla
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | | |
Collapse
|
20
|
Grotz KK, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D. J Chem Theory Comput 2022; 18:526-537. [PMID: 34881568 PMCID: PMC8757469 DOI: 10.1021/acs.jctc.1c00791] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Indexed: 12/14/2022]
Abstract
Magnesium is essential in many vital processes. To correctly describe Mg2+ in physiological processes by molecular dynamics simulations, accurate force fields are fundamental. Despite the importance, force fields based on the commonly used 12-6 Lennard-Jones potential showed significant shortcomings. Recently progress was made by an optimization procedure that implicitly accounts for polarizability. The resulting microMg and nanoMg force fields (J. Chem. Theory Comput. 2021, 17, 2530-2540) accurately reproduce a broad range of experimental solution properties and the binding affinity to nucleic acids in TIP3P water. Since countless simulation studies rely on available water models and ion force fields, we here extend the optimization and provide Mg2+ parameters in combination with the SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D water models. For each water model, the Mg2+ force fields reproduce the solvation free energy, the distance to oxygens in the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, and the binding affinity and distance to the phosphate oxygens on nucleic acids. We present two parameter sets: MicroMg yields water exchange on the microsecond time scale and matches the experimental exchange rate. Depending on the water model, nanoMg yields accelerated water exchange in the range of 106 to 108 exchanges per second. The nanoMg parameters can be used to enhance the sampling of binding events, to obtain converged distributions of Mg2+, or to predict ion binding sites in biomolecular simulations. The parameter files are freely available at https://github.com/bio-phys/optimizedMgFFs.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical
Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|
21
|
Rubinson KA. Measurement of Enthalpies and Entropies of Activation as a Function of Pairwise Distance for the Pairwise Relative Diffusion of SrI 2 in Water over Lengthscales from 6 Å to 40 Å. NEW J CHEM 2022. [DOI: 10.1039/d2nj00712f] [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
Quasi-elastic neutron scattering (QENS) has been used to provide insight into the motions within solutions for many decades, and coherent QENS (CQENS) determines the motions of one ion relative to...
Collapse
|
22
|
Hantal G, Sega M, Horvai G, Jedlovszky P. Contribution of Different Molecules and Moieties to the Surface Tension in Aqueous Surfactant Solutions. II: Role of the Size and Charge Sign of the Counterions. J Phys Chem B 2021; 125:9005-9018. [PMID: 34319728 DOI: 10.1021/acs.jpcb.1c04216] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the role of the counterion species in surfactant solutions is a complicated task, made harder by the fact that, experimentally, it is not possible to vary independently bulk and surface quantities. Here, we perform molecular dynamics simulations at constant surface coverage of the liquid/vapor interface of lithium, sodium, potassium, rubidium, and cesium dodecyl sulfate aqueous solutions. We investigate the effect of counterion type and charge sign on the surface tension of the solution, analyzing the contribution of different species and moieties to the lateral pressure profile. The observed trends are qualitatively compatible with the Hofmeister series, with the notable exception of sodium. We point out a possible shortcoming of what is at the moment, in our experience, the most realistic nonpolarizable force field (CHARMM36) that includes the parametrization for the whole series of alkali counterions. In the artificial system where the counterion and surfactant charges are inverted in sign, the counterions become considerably harder. This charge inversion changes considerably the surface tension contributions of the counterions, surfactant headgroups, and water molecules, stressing the key role of the hardness of the counterions in this respect. However, the hydration free energy gain of the counterions, occurring upon charge inversion, is compensated by the concomitant free energy loss of the headgroups and water molecules, leading to a negligible change in the surface tension of the entire system.
Collapse
Affiliation(s)
- György Hantal
- Institute of Physics and Materials Science, University of Natural Resources and Life Sciences, Peter Jordan Straße 82, A-1190 Vienna, Austria.,Department of Chemistry, Eszterházy Károly University, Leányka utca 6, H-3300 Eger, Hungary
| | - Marcello Sega
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11),Fürther Straße 248, D-90429 Nürnberg, Germany
| | - George Horvai
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H-1111 Budapest, Hungary
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka utca 6, H-3300 Eger, Hungary
| |
Collapse
|
23
|
Saravi SH, Panagiotopoulos AZ. Individual Ion Activity Coefficients in Aqueous Electrolytes from Explicit-Water Molecular Dynamics Simulations. J Phys Chem B 2021; 125:8511-8521. [PMID: 34319101 DOI: 10.1021/acs.jpcb.1c04019] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
We compute individual ion activity coefficients (IIACs) in aqueous NaCl, KCl, NaF, and KF solutions from explicit-water molecular dynamics simulations. Free energy changes are obtained from insertion of single ions-accompanied by uniform neutralizing backgrounds-into solution by gradually turning on first Lennard-Jones interactions, followed by Coulombic interactions using Ewald electrostatics. Simulations are performed at multiple system sizes, and all results are extrapolated to the thermodynamic limit, thus eliminating any possible artifacts from the neutralizing backgrounds. Because of controversies associated with measurements of IIACs from electrochemical cells with ion-selective electrodes, the reported experimental data are not widely accepted; thus there remains a knowledge gap with respect to the contributions of individual ions to solution nonidealities. Our results are in good qualitative agreement with these reported measurements, though significantly larger in magnitude. In particular, the relative positioning for the activity coefficients of anions and cations matches the experimental ordering for all four systems. This work establishes a robust thermodynamic framework, without a need to invoke extra hypotheses, that sheds light on the behavior of individual ions and their contributions to nonidealities of aqueous electrolyte solutions.
Collapse
Affiliation(s)
- Sina Hassanjani Saravi
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | | |
Collapse
|
24
|
Proteins maintain hydration at high [KCl] concentration regardless of content in acidic amino acids. Biophys J 2021; 120:2746-2762. [PMID: 34087206 PMCID: PMC8390907 DOI: 10.1016/j.bpj.2021.05.015] [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] [Received: 01/06/2021] [Revised: 04/22/2021] [Accepted: 05/19/2021] [Indexed: 11/23/2022] Open
Abstract
Proteins of halophilic organisms, which accumulate molar concentrations of KCl in their cytoplasm, have a much higher content in acidic amino acids than proteins of mesophilic organisms. It has been proposed that this excess is necessary to maintain proteins hydrated in an environment with low water activity, either via direct interactions between water and the carboxylate groups of acidic amino acids or via cooperative interactions between acidic amino acids and hydrated cations. Our simulation study of five halophilic proteins and five mesophilic counterparts does not support either possibility. The simulations use the AMBER ff14SB force field with newly optimized Lennard-Jones parameters for the interactions between carboxylate groups and potassium ions. We find that proteins with a larger fraction of acidic amino acids indeed have higher hydration levels, as measured by the concentration of water in their hydration shell and the number of water/protein hydrogen bonds. However, the hydration level of each protein is identical at low (bKCl = 0.15 mol/kg) and high (bKCl = 2 mol/kg) KCl concentrations; excess acidic amino acids are clearly not necessary to maintain proteins hydrated at high salt concentration. It has also been proposed that cooperative interactions between acidic amino acids in halophilic proteins and hydrated cations stabilize the folded protein structure and would lead to slower dynamics of the solvation shell. We find that the translational dynamics of the solvation shell is barely distinguishable between halophilic and mesophilic proteins; if such a cooperative effect exists, it does not have that entropic signature.
Collapse
|
25
|
Klein T, Piszko M, Kankanamge CJ, Kasapis G, Fröba AP. Fick Diffusion Coefficient in Binary Mixtures of [HMIM][NTf 2] and Carbon Dioxide by Dynamic Light Scattering and Molecular Dynamics Simulations. J Phys Chem B 2021; 125:5100-5113. [PMID: 33960781 DOI: 10.1021/acs.jpcb.1c01616] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Dynamic light scattering (DLS) experiments and equilibrium molecular dynamics (EMD) simulations were performed in the saturated liquid phase of the binary mixture of 1-hexyl-3-methylimidazolium bis(trifluormethylsulfonyl)imide ([HMIM][NTf2]) and carbon dioxide (CO2) to access the Fick diffusion coefficient (D11). The investigations were performed within or close to saturation conditions at temperatures between (298.15 and 348.15) K and CO2 mole fractions (xCO2) up to 0.81. The DLS experiments were combined with polarization-difference Raman spectroscopy (PDRS) to simultaneously access the composition of the liquid phase. For the first time in an electrolyte-based system, D11 was directly calculated from EMD simulations by accessing the Maxwell-Stefan (MS) diffusion coefficient and the thermodynamic factor. Agreement within combined uncertainties was found between D11 from DLS and EMD simulations for CO2 mole fractions up to 0.5. In general, an increasing D11 with increasing xCO2 could be observed, with a local maximum present at a CO2 mole fraction of about 0.75. The local maximum could be explained by an increasing MS diffusion coefficient with increasing xCO2 over the entire studied composition range and a decreasing thermodynamic factor at xCO2 above 0.7. Finally, PDRS and EMD simulations were combined to investigate the influence of the fluid structure on the diffusive process.
Collapse
Affiliation(s)
- Tobias Klein
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Maximilian Piszko
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Chathura J Kankanamge
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| | - Georgios Kasapis
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany.,Institute for Multiscale Thermofluids, School of Engineering, University of Edinburgh, Robert Stevenson Road, Edinburgh EH9 3FB, U.K
| | - Andreas P Fröba
- Institute of Advanced Optical Technologies - Thermophysical Properties (AOT-TP), Department of Chemical and Biological Engineering (CBI) and Erlangen Graduate School in Advanced Optical Technologies (SAOT), Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Paul-Gordan-Straße 8, 91052 Erlangen, Germany
| |
Collapse
|
26
|
Cruz-León S, Grotz KK, Schwierz N. Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations. J Chem Phys 2021; 154:171102. [PMID: 34241062 DOI: 10.1063/5.0048113] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca2+ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg2+ improves the structure of two specifically bound Mg2+ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Kara K Grotz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
27
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
28
|
Grotz KK, Cruz-León S, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties. J Chem Theory Comput 2021; 17:2530-2540. [PMID: 33720710 PMCID: PMC8047801 DOI: 10.1021/acs.jctc.0c01281] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Indexed: 12/31/2022]
Abstract
Magnesium ions play an essential role in many vital processes. To correctly describe their interactions in molecular dynamics simulations, an accurate parametrization is crucial. Despite the importance and considerable scientific effort, current force fields based on the commonly used 12-6 Lennard-Jones interaction potential fail to reproduce a variety of experimental solution properties. In particular, no parametrization exists so far that simultaneously reproduces the solvation free energy and the distance to the water oxygens in the first hydration shell. Moreover, current Mg2+ force fields significantly underestimate the rate of water exchange leading to unrealistically slow exchange kinetics. In order to make progress in the development of improved models, we systematically optimize the Mg2+ parameters in combination with the TIP3P water model in a much larger parameter space than previously done. The results show that a long-ranged interaction potential and modified Lorentz-Berthelot combination rules allow us to accurately reproduce multiple experimental properties including the solvation free energy, the distances to the oxygens of the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, the self-diffusion coefficient, and the binding affinity to the phosphate oxygen of RNA. Matching this broad range of thermodynamic properties, we present two sets of optimal parameters: MicroMg yields water exchange on the microsecond timescale in agreement with experiments. NanoMg yields water exchange on the nanosecond timescale facilitating the direct observation of ion-binding events. As shown for the example of the add A-riboswitch, the optimized parameters correctly reproduce the structure of specifically bound ions and permit the de novo prediction of Mg2+-binding sites in biomolecular simulations.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Sergio Cruz-León
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|
29
|
Duboué-Dijon E, Javanainen M, Delcroix P, Jungwirth P, Martinez-Seara H. A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization. J Chem Phys 2021; 153:050901. [PMID: 32770904 DOI: 10.1063/5.0017775] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.
Collapse
Affiliation(s)
- E Duboué-Dijon
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - M Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - P Delcroix
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - P Jungwirth
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - H Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| |
Collapse
|
30
|
Primorac T, Požar M, Sokolić F, Zoranić L. The influence of binary mixtures' structuring on the calculation of Kirkwood-Buff integrals: A molecular dynamics study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114773] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
31
|
Hantal G, Horváth RA, Kolafa J, Sega M, Jedlovszky P. Surface Affinity of Alkali and Halide Ions in Their Aqueous Solution: Insight from Intrinsic Density Analysis. J Phys Chem B 2020; 124:9884-9897. [PMID: 33084342 DOI: 10.1021/acs.jpcb.0c05547] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The surface tension of all aqueous alkali halide solutions is higher than that of pure water. According to the Gibbs adsorption equation, this indicates a net depletion of these ions in the interfacial region. However, simulations and experiments show that large, soft ions, such as I-, can accumulate at the liquid/vapor interface. The presence of a loose hydration shell is usually considered to be the reason for this behavior. In this work, we perform computer simulations to characterize the liquid-vapor interface of aqueous alkali chloride and sodium halide solutions systematically, considering all ions from Li+ to Cs+ and from F- to I-. Using computational methods for the removal of surface fluctuations, we analyze the structure of the interface at a dramatically enhanced resolution, showing that the positive excess originates in the very first molecular layer and that the next 3-4 layers account for the net negative excess. With the help of a fictitious system with charge-inverted ion pairs, we also show that it is not possible to rationalize the surface affinity of ions in solutions in terms of the properties of anions and cations separately. Moreover, the surface excess is generally dominated by the smaller of the two ions.
Collapse
Affiliation(s)
- György Hantal
- Institute of Physics and Materials Science, University of Natural Resources and Life Sciences, Peter Jordan Straße 82, A-1190 Vienna, Austria
| | - Réka A Horváth
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H-1111 Budapest, Hungary
| | - Jiří Kolafa
- Department of Physical Chemistry, University of Chemistry and Technology, 166 28 Prague 6, Czech Republic
| | - Marcello Sega
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Straße 248, D-90429 Nürnberg, Germany
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka utca 6, H-3300 Eger, Hungary
| |
Collapse
|
32
|
Smirnov PR. Structure of the Nearest Environment of
Na+, K+,
Rb+, and Cs+ Ions in
Oxygen-Containing Solvents. RUSS J GEN CHEM+ 2020. [DOI: 10.1134/s1070363220090169] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
33
|
Panagiotopoulos AZ. Simulations of activities, solubilities, transport properties, and nucleation rates for aqueous electrolyte solutions. J Chem Phys 2020; 153:010903. [DOI: 10.1063/5.0012102] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
34
|
Rubinson KA, Mountain RD. Ion and water transport reasonably involves rotation and pseudorotation: measurement and modeling the temperature dependence of small-angle neutron scattering from aqueous SrI 2. Phys Chem Chem Phys 2020; 22:13479-13488. [PMID: 32525150 DOI: 10.1039/d0cp02088e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
X-ray and neutron scattering have provided insight into the short range (<8 Å) structures of ionic solutions for over a century. For longer distances, single scattering bands have, however, been seen. For the non-hydrolyzing salt SrI2 in aqueous (D2O) solution, a structure sufficient to scatter slow neutrons has been seen to persist down to a concentration of 0.1 mol L-1 where the measured average spacing between scatterers is over 20 Å. Theoretical studies of such long distance solution structures are difficult, and these difficulties are discussed. The width of the distribution in distances between the scatterers (ions, ion pairs, etc.) remains less than 10 Å, which approximates the average size of the ions and their first hydration shell. Here, we measure the temperature dependence from 10 °C to 90 °C of the small angle neutron scattering (SANS) by a 0.5 molar SrI2 solution in D2O and find that this surprisingly narrow distribution of the distances remains constant within experimental uncertainty. This structure of the ions in the solution appears to endure because changes in interion distances along any single spatial dimension require displacements near the size of a water molecule. Together, the experimental measurements support a rotatory mechanism for simultaneous ion transport and water countertransport. Since rotation minimizes displacement of the solution framework, it is suggested that water transport alone also involves rotation of multimolecular structures, and that the interpretation of single-molecule water rotation is confounded by pseudorotation that results from paired picosecond proton exchanges. It is pointed out that NMR-determined millisecond to microsecond proton exchange times of chelated-metal-ion bound waters and the much faster chelate rotational correlation times around 10 picoseconds, both of which require making and breaking of hydrogen bonds, are difficult to impossible to reconcile.
Collapse
Affiliation(s)
- Kenneth A Rubinson
- Department of Biochemistry and Molecular Biology, Wright State University, Dayton, OH 45435, USA. and NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.
| | - Raymond D Mountain
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.
| |
Collapse
|
35
|
Stenqvist B, Aspelin V, Lund M. Generalized Moment Correction for Long-Ranged Electrostatics. J Chem Theory Comput 2020; 16:3737-3745. [PMID: 32315176 PMCID: PMC7588037 DOI: 10.1021/acs.jctc.9b01003] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
![]()
Describing
long-ranged electrostatics using short-ranged pair potentials
is appealing because the computational complexity scales linearly
with the number of particles. The foundation of the approach presented
here is to mimic the long-ranged medium response by cancelling electric
multipoles within a small cutoff sphere. We propose a rigorous and
formally exact new method that cancels up to infinitely many multipole moments and is free of operational damping parameters
often required in existing theories. Using molecular dynamics simulations
of water with and without added salt, we discuss radial distribution
functions, Kirkwood–Buff integrals, dielectrics, diffusion
coefficients, and angular correlations in relation to existing electrostatic
models. We find that the proposed method is an efficient and accurate
alternative for handling long-ranged electrostatics as compared to
Ewald summation schemes. The methodology and proposed parameterization
are applicable also for dipole–dipole interactions.
Collapse
Affiliation(s)
- Björn Stenqvist
- Division of Physical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden.,Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Vidar Aspelin
- Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Mikael Lund
- Division of Theoretical Chemistry, Department of Chemistry, Lund University, 221 00 Lund, Sweden.,LINXS-Lund Institute of Advanced Neutron and X-ray Science, Scheelevägen 19, SE-223 70 Lund, Sweden
| |
Collapse
|
36
|
Saric D, Kohns M, Vrabec J. Dielectric constant and density of aqueous alkali halide solutions by molecular dynamics: A force field assessment. J Chem Phys 2020; 152:164502. [PMID: 32357782 DOI: 10.1063/1.5144991] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The concentration dependence of the dielectric constant and the density of 11 aqueous alkali halide solutions (LiCl, NaCl, KCl, RbCl, CsCl, LiI, NaI, KI, CsI, KF, and CsF) is investigated by molecular simulation. Predictions using eight non-polarizable ion force fields combined with the TIP4P/ε water model are compared to experimental data. The influence of the water model and the temperature on the results for the NaCl brine are also addressed. The TIP4P/ε water model improves the accuracy of dielectric constant predictions compared to the SPC/E water model. The solution density is predicted well by most ion models. Almost all ion force fields qualitatively capture the decline of the dielectric constant with the increase of concentration for all solutions and with the increase of temperature for NaCl brine. However, the sampled dielectric constant is mostly in poor quantitative agreement with experimental data. These results are related to the microscopic solution structure, ion pairing, and ultimately the force field parameters. Ion force fields with excessive contact ion pairing and precipitation below the experimental solubility limit generally yield higher dielectric constant values. An adequate reproduction of the experimental solubility limit should therefore be a prerequisite for further investigations of the dielectric constant of aqueous electrolyte solutions by molecular simulation.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics and Energy Technology, University of Paderborn, 33098 Paderborn, Germany
| | - Maximilian Kohns
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67633 Kaiserslautern, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, Technical University Berlin, 10587 Berlin, Germany
| |
Collapse
|
37
|
Role of the Counterions in the Surface Tension of Aqueous Surfactant Solutions. A Computer Simulation Study of Alkali Dodecyl Sulfate Systems. COLLOIDS AND INTERFACES 2020. [DOI: 10.3390/colloids4020015] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
We have investigated the surface tension contributions of the counterions, surfactant headgroups and tails, and water molecules in aqueous alkali dodecyl sulfate (DS) solutions close to the saturated surface concentration by analyzing the lateral pressure profile contribution of these components using molecular dynamics simulations. For this purpose, we have used the combination of two popular force fields, namely KBFF for the counterions and GROMOS96 for the surfactant, which are both parameterized for the SPC/E water model. Except for the system containing Na+ counterions, the surface tension of the surfactant solutions has turned out to be larger rather than smaller than that of neat water, showing a severe shortcoming of the combination of the two force fields. We have traced back this failure of the potential model combination to the unphysically strong attraction of the KBFF counterions, except for Na+, to the anionic head of the surfactants. Despite this failure of the model, we have observed a clear relation between the soft/hard character (in the sense of the Hofmeister series) and the surface tension contribution of the counterions, which, given the above limitations of the model, can only be regarded as an indicative result. We emphasize that the obtained results, although in a twisted way, clearly stress the crucial role the counterions of ionic surfactants play in determining the surface tension of the aqueous surfactant solutions.
Collapse
|
38
|
Yagasaki T, Matsumoto M, Tanaka H. Lennard-Jones Parameters Determined to Reproduce the Solubility of NaCl and KCl in SPC/E, TIP3P, and TIP4P/2005 Water. J Chem Theory Comput 2020; 16:2460-2473. [PMID: 32207974 DOI: 10.1021/acs.jctc.9b00941] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Most classical nonpolarizable ion potential models underestimate the solubility values of NaCl and KCl in water significantly. We determine Lennard-Jones parameters of Na+, K+, and Cl- that reproduce the solubility as well as the hydration free energy in dilute aqueous solutions for three water potential models, SPC/E, TIP3P, and TIP4P/2005. The ion-oxygen distance in the solution and the cation-anion distance in salt are also considered in the parametrization. In addition to the target properties, the hydration enthalpy, hydration entropy, self-diffusion coefficient, coordination number, lattice energy, enthalpy of solution, density, viscosity, and number of contact ion pairs are calculated for comparison with 17 frequently used or recently developed ion potential models. The overall performance of each ion model is represented by a global score using a scheme that was originally developed for comparison of water potential models. The global score is better for our models than for the other 17 models not only because of the quite good prediction for the solubility but also because of the relatively small deviation from the experimental value for many of the other properties.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Masakazu Matsumoto
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| | - Hideki Tanaka
- Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan.,Department of Chemistry, Faculty of Science, Okayama University, Okayama 700-8530, Japan
| |
Collapse
|
39
|
Polák J, Ondo D, Heyda J. Thermodynamics of N-Isopropylacrylamide in Water: Insight from Experiments, Simulations, and Kirkwood-Buff Analysis Teamwork. J Phys Chem B 2020; 124:2495-2504. [PMID: 32118440 DOI: 10.1021/acs.jpcb.0c00413] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The behavior of thermoresponsive polymer poly(N-isopropylacrylamide) (PNiPAM), an essential building block in the design of smart soft materials, in aqueous solutions has attracted much interest, which contrasts with our knowledge of N-isopropylacrylamide (NiPAM) monomer. Strikingly, the physicochemical properties of aqueous NiPAM are similarly rich, and their understanding is far from being complete. This stems from the lack of accurate thermodynamic data and quantitative model for atomistic simulations. In this joint study, we have probed the thermodynamic behavior of aqueous NiPAM by experimental methods, molecular dynamics (MD) simulations, and Kirkwood-Buff (KB) analysis at ambient conditions. From the partial molar volumes and simultaneously correlated osmotic coefficients, with excess partial molar enthalpies of NiPAM in water, the concentration and temperature dependence of KB integrals was determined. For the purpose of this work, we have developed and employed a novel NiPAM force field, which not only reproduces KB integrals (Gij) and adequately captures macroscopic thermodynamic quantities but also provides more accurate structural insight than the original force fields. We revealed in the vicinity of NiPAM the competing effect of amide hydration with interaction between nonpolar regions. This microscopic picture is reflected in the experimentally observed NiPAM-NiPAM association, which is present from highly dilute conditions up to the solubility limit and is evidenced by G22. From intermediate concentrations, it is accompanied by the existence of apparent dense-water regions, as indicated by positive G11 values. The here-employed KB-based framework provided a mutually consistent thermodynamic and microscopic insight into the NiPAM solution and may be further extended for ion-specific effects. Moreover, our findings contribute to the understanding of thermodynamic grounds behind PNiPAM collapse transition.
Collapse
Affiliation(s)
- Jakub Polák
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Daniel Ondo
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| | - Jan Heyda
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague 6, Czech Republic
| |
Collapse
|
40
|
|
41
|
Lbadaoui-Darvas M, Takahama S. Water Activity from Equilibrium Molecular Dynamics Simulations and Kirkwood-Buff Theory. J Phys Chem B 2019; 123:10757-10768. [DOI: 10.1021/acs.jpcb.9b06735] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mária Lbadaoui-Darvas
- ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Satoshi Takahama
- ENAC/IIE Swiss Federal Institute of Technology Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
42
|
Zeron IM, Abascal JLF, Vega C. A force field of Li +, Na +, K +, Mg 2+, Ca 2+, Cl -, and SO 4 2- in aqueous solution based on the TIP4P/2005 water model and scaled charges for the ions. J Chem Phys 2019; 151:134504. [PMID: 31594349 DOI: 10.1063/1.5121392] [Citation(s) in RCA: 120] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li+, Na+, K+, Mg2+, and Ca2+ and the anions Cl- and SO4 2-. These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt.
Collapse
Affiliation(s)
- I M Zeron
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J L F Abascal
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
43
|
Yue S, Panagiotopoulos AZ. Dynamic properties of aqueous electrolyte solutions from non-polarisable, polarisable, and scaled-charge models. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1645901] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ, USA
| | | |
Collapse
|
44
|
Saravi SH, Ravichandran A, Khare R, Chen CC. Bridging two-liquid theory with molecular simulations for electrolytes: An investigation of aqueous NaCl solution. AIChE J 2019. [DOI: 10.1002/aic.16521] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Sina H. Saravi
- Dept. of Chemical Engineering; Texas Tech University; Lubbock TX 79409
| | | | - Rajesh Khare
- Dept. of Chemical Engineering; Texas Tech University; Lubbock TX 79409
| | - Chau-Chyun Chen
- Dept. of Chemical Engineering; Texas Tech University; Lubbock TX 79409
| |
Collapse
|
45
|
Petris PC, Anogiannakis SD, Tzounis PN, Theodorou DN. Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory. J Phys Chem B 2019; 123:247-257. [PMID: 30516991 DOI: 10.1021/acs.jpcb.8b10425] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δmix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.
Collapse
Affiliation(s)
- Panagiotis C Petris
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | - Stefanos D Anogiannakis
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| | | | - Doros N Theodorou
- School of Chemical Engineering , National Technical University of Athens , GR 15780 Athens , Greece
| |
Collapse
|
46
|
Zhang F, Yu Z, Rondinone AJ, Huang J, Sumpter BG, Qiao R. Adsorption of Molecular Nitrogen in Electrical Double Layers near Planar and Atomically Sharp Electrodes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:14552-14561. [PMID: 30411900 DOI: 10.1021/acs.langmuir.8b03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The adsorption of gas molecules at electrode-electrolyte interfaces is an important step in electrochemical reactions. Using molecular dynamics simulations, we investigate the adsorption of dissolved N2 in the electrical double layers (EDLs) of an aqueous electrolyte near planar and 1 nm radius spherical carbon electrodes. The adsorption of N2 is found to be overall enriched near neutral electrodes regardless of their surface curvature, although it can be locally enriched or depleted depending on the distance from the electrode surface. In comparison, the adsorption of N2 in the EDL near negatively charged electrodes is found to increase under a moderate surface charge density, but decrease under a high surface charge density, especially near a planar electrode. By analyzing the potential of mean force for dissolved N2, the solvent-induced effects are found to play important roles in influencing the adsorption of N2 in the EDLs. The adsorption behavior of N2 molecules, especially their dependence on the surface charge and curvature of electrodes, is further rationalized by examining the structure of interfacial water molecules, their interference with the hydration shell of N2, and their modification by the electrification of electrodes.
Collapse
Affiliation(s)
- Fei Zhang
- Department of Mechanical Engineering , Virginia Tech , Blacksburg , Virginia 24061 , United States
| | - Zhou Yu
- Department of Mechanical Engineering , Virginia Tech , Blacksburg , Virginia 24061 , United States
| | | | | | | | - Rui Qiao
- Department of Mechanical Engineering , Virginia Tech , Blacksburg , Virginia 24061 , United States
| |
Collapse
|
47
|
Properties of Ion Complexes and Their Impact on Charge Transport in Organic Solvent-Based Electrolyte Solutions for Lithium Batteries: Insights from a Theoretical Perspective. BATTERIES-BASEL 2018. [DOI: 10.3390/batteries4040062] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Electrolyte formulations in standard lithium ion and lithium metal batteries are complex mixtures of various components. In this article, we review molecular key principles of ion complexes in multicomponent electrolyte solutions in regards of their influence on charge transport mechanisms. We outline basic concepts for the description of ion–solvent and ion–ion interactions, which can be used to rationalize recent experimental and numerical findings concerning modern electrolyte formulations. Furthermore, we discuss benefits and drawbacks of empirical concepts in comparison to molecular theories of solution for a more refined understanding of ion behavior in organic solvents. The outcomes of our discussion provide a rational for beneficial properties of ions, solvent, co-solvent and additive molecules, and highlight possible routes for further improvement of novel electrolyte solutions.
Collapse
|
48
|
Michalowsky J, Zeman J, Holm C, Smiatek J. A polarizable MARTINI model for monovalent ions in aqueous solution. J Chem Phys 2018; 149:163319. [PMID: 30384758 DOI: 10.1063/1.5028354] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a new polarizable coarse-grained martini force field for monovalent ions, called refIon, which is developed mainly for the accurate reproduction of electrostatic properties in aqueous electrolyte solutions. The ion model relies on full long-range Coulomb interactions and introduces satellite charges around the central interaction site in order to model molecular polarization effects. All force field parameters are matched to reproduce the mass density and the static dielectric permittivity of aqueous NaCl solutions, such that experimental values are well-reproduced up to moderate salt concentrations of 2 m o l / l . In addition, an improved agreement with experimentally measured ionic conductivities is observed. Our model is validated with regard to analytic solutions for the ion distribution around highly charged rod-like polyelectrolytes in combination with atomistic simulations and experimental results concerning structural properties of lipid bilayers in the presence of distinct salt concentrations. Further results regarding the coordination numbers of counterions around dilute poly(styrene sulfonate) and poly(diallyldimethylammonium) polyelectrolyte chains also highlight the applicability of our approach. The introduction of our force field allows us to eliminate heuristic scaling factors, as reported for previous martini ion models in terms of effective salt concentrations, and in consequence provides a better agreement between simulation and experimental results. The presented approach is specifically useful for recent martini attempts that focus on highly charged systems-such as models of DNA, polyelectrolytes or polyelectrolyte complexes-where precise studies of electrostatic effects and charge transport processes are essential.
Collapse
Affiliation(s)
- Julian Michalowsky
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| | - Johannes Zeman
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| | - Christian Holm
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| | - Jens Smiatek
- Institut für Computerphysik, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany
| |
Collapse
|
49
|
Pylaeva S, Brehm M, Sebastiani D. Salt Bridge in Aqueous Solution: Strong Structural Motifs but Weak Enthalpic Effect. Sci Rep 2018; 8:13626. [PMID: 30206276 PMCID: PMC6133928 DOI: 10.1038/s41598-018-31935-z] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2018] [Accepted: 08/29/2018] [Indexed: 12/03/2022] Open
Abstract
Salt bridges are elementary motifs of protein secondary and tertiary structure and are commonly associated with structural driving force that increases stability. Often found on the interface to the solvent, they are highly susceptible to solvent–solute interactions, primarily with water but also with other cosolvents (especially ions). We have investigated the interplay of an Arginine–Aspartic acid salt bridge with simple salt ions in aqueous solution by means of molecular dynamics simulations. Besides structural and dynamical features at equilibrium, we have computed the mean force along the dissociation pathway of the salt bridge. We demonstrate that solvated ions influence the behavior of the salt bridge in a very specific and local way, namely the formation of tight ionic pairs Li+/Na+–Asp−. Moreover, our findings show that the enthalpic relevance of the salt bridge is minor, regardless of the presence of solvated ions.
Collapse
Affiliation(s)
- Svetlana Pylaeva
- Martin Luther University Halle-Wittenberg, Institute of Physical Chemistry, 06120, Halle, Saale, Germany
| | - Martin Brehm
- Martin Luther University Halle-Wittenberg, Institute of Physical Chemistry, 06120, Halle, Saale, Germany
| | - Daniel Sebastiani
- Martin Luther University Halle-Wittenberg, Institute of Physical Chemistry, 06120, Halle, Saale, Germany.
| |
Collapse
|
50
|
Pethes I. The structure of aqueous lithium chloride solutions at high concentrations as revealed by a comparison of classical interatomic potential models. J Mol Liq 2018. [DOI: 10.1016/j.molliq.2018.05.044] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
|