Accounting for Electronic Polarization Effects in Aqueous Sodium Chloride via Molecular Dynamics Aided by Neutron Scattering

Effective Inclusion of Electronic Polarization Improves the Description of Electrostatic Interactions: The prosECCo75 Biomolecular Force Field

prosECCo75 is an optimized force field effectively incorporating electronic polarization via charge scaling. It aims to enhance the accuracy of nominally nonpolarizable molecular dynamics simulations for interactions in biologically relevant systems involving water, ions, proteins, lipids, and saccharides. Recognizing the inherent limitations of nonpolarizable force fields in precisely modeling electrostatic interactions essential for various biological processes, we mitigate these shortcomings by accounting for electronic polarizability in a physically rigorous mean-field way that does not add to computational costs. With this scaling of (both integer and partial) charges within the CHARMM36 framework, prosECCo75 addresses overbinding artifacts. This improves agreement with experimental ion binding data across a broad spectrum of systems─lipid membranes, proteins (including peptides and amino acids), and saccharides─without compromising their biomolecular structures. prosECCo75 thus emerges as a computationally efficient tool providing enhanced accuracy and broader applicability in simulating the complex interplay of interactions between ions and biomolecules, pivotal for improving our understanding of many biological processes.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

Anomalous Concentration Dependence of Viscosity: Hidden Role of Cross-Correlations in Aqueous Electrolyte Solutions

The viscosity of aqueous electrolyte solutions exhibits well-known composition-dependent anomalies that show certain definitive trends and universal features. The viscosity of LiCl and NaCl solutions increases with concentration in a monotonic fashion, while solutions of KCl, RbCl, and CsCl exhibit a more complex behavior. Here, the viscosity first decreases and then increases with increasing concentration, with a rather broad minimum at intermediate concentrations (ca. 1-3 m). To unearth the origin of such puzzling behavior, we carried out detailed molecular-level analyses by interrogating the exact Green-Kubo expression of viscosity in terms of the stress-stress time correlation function (SS-TCF). The total SS-TCF can be decomposed into a collection of three self- and three cross-SS-TCFs arising from the three constituent components (water, cations, and anions). Mode coupling theory (MCT) analysis for the friction on ions and the viscosity of the solution suggests the possible importance of two-particle static and time-dependent cross-correlations between water and the ions. We calculate the viscosity and other dynamical properties for all five electrolyte (LiCl, NaCl, KCl, RbCl, and CsCl) solutions over a range of concentrations, using two models of water (SPC/E and TIP4P/2005). The total viscosity derives non-negligible contributions from all of the terms. The cross-correlations are found to be surprisingly large and seen to play a hidden role in the concentration dependence. However, the importance of cross-correlations is often not discussed. Our study leads to a theoretical understanding of the microscopic origin of the observed anomalies in the composition dependence of viscosity across all five electrolytes.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.

Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials

The diffusivity of water in aqueous cesium iodide solutions is larger than that in neat liquid water and vice versa for sodium chloride solutions. Such peculiar ion-specific behavior, called anomalous diffusion, is not reproduced in typical force field based molecular dynamics (MD) simulations due to inadequate treatment of ion-water interactions. Herein, this hurdle is tackled by using machine learned atomic potentials (MLPs) trained on data from density functional theory calculations. MLP based atomistic MD simulations of aqueous salt solutions reproduce experimentally determined thermodynamic, structural, dynamical, and transport properties, including their varied trends in water diffusivities across salt concentration. This enables an examination of their intermolecular structure to unravel the microscopic underpinnings of the differences in their transport properties. While both ions in CsI solutions contribute to the faster diffusion of water molecules, the competition between the heavy retardation by Na ions and the slight acceleration by Cl ions in NaCl solutions reduces their water diffusivity.

Anions and Cations Affect Amino Acid Dissociation Equilibria via Distinct Mechanisms

Salts reduce the pKa of weak acids by a mechanism sensitive to ion identity and concentration via charge screening of the deprotonated state. In this study, we utilize constant pH molecular dynamics simulations to understand the molecular mechanism behind the salt-dependent dissociation of aspartic acid (Asp). We calculate the pKa of Asp in the presence of a monovalent salt and investigate Hofmeister ion effects by systematically varying the ionic radii. We observe that increasing the anion size leads to a monotonic decrease in Asp pKa. Conversely, the cation size affects the pKa nonmonotonically, interpretable in the context of the law of matching water affinity. The net effect of salt on Asp acidity is governed by an interplay of solvation and competing ion interactions. The proposed mechanism is rather general and can be applicable to several problems in Hofmeister ion chemistry, such as pH effects on protein stability and soft matter interfaces.

Evolution of the electrical double layer with electrolyte concentration probed by second harmonic scattering

Investigating the electrical double layer (EDL) structure has been a long-standing challenge and has seen the emergence of several sophisticated techniques able to probe selectively the few molecular layers of a solid/water interface. While a qualitative estimation of the thickness of the EDL can be obtained using simple theoretical models, following experimentally its evolution is not straightforward and can be even more complicated in nano- or microscale systems, particularly when changing the ionic concentration by several orders of magnitude. Here, we bring insight into the structure of the EDL of SiO2 nanoparticle suspensions and its evolution with increasing ionic concentration using angle-resolved second harmonic scattering (AR-SHS). Below millimolar salt concentrations, we can successively characterize inner-sphere adsorption, diffuse layer formation, and outer-sphere adsorption. Moreover, we show for the first time that, by appropriately selecting the nanoparticle size, it is possible to retrieve information also in the millimolar range. There, we observe a decrease in the magnitude of the surface potential corresponding to a compression in the EDL thickness, which agrees with the results of several other electroanalytical and optical techniques. Molecular dynamics simulations suggest that the EDL compression mainly results from the diffuse layer compression rather than outer-sphere ions (Stern plane) moving closer to the surface.

Nonaqueous Ion Pairing Exemplifies the Case for Including Electronic Polarization in Molecular Dynamics Simulations

The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

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.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

Common Cations Are Not Polarizable: Effects of Dispersion Correction on Hydration Structures from Ab Initio Molecular Dynamics

We employed density functional theory-based ab initio molecular dynamics simulations to examine the hydration structure of several common alkali and alkali earth metal cations. We found that the commonly used atom pairwise dispersion correction scheme D3, which assigns dispersion coefficients based on the neutral form of the atom rather than its actual oxidation state, leads to inaccuracies in the hydration structures of these cations. We evaluated this effect for lithium, sodium, potassium, and calcium and found that the inaccuracies are particularly pronounced for sodium and potassium compared to the experiment. To remedy this issue, we propose disabling the D3 correction specifically for all cation-including pairs, which leads to a much better agreement with experimental data.

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water

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.

Ion-Induced Transient Potential Fluctuations Facilitate Pore Formation and Cation Transport through Lipid Membranes

Unassisted ion transport through lipid membranes plays a crucial role in many cell functions without which life would not be possible, yet the precise mechanism behind the process remains unknown due to its molecular complexity. Here, we demonstrate a direct link between membrane potential fluctuations and divalent ion transport. High-throughput wide-field non-resonant second harmonic (SH) microscopy of membrane water shows that membrane potential fluctuations are universally found in lipid bilayer systems. Molecular dynamics simulations reveal that such variations in membrane potential reduce the free energy cost of transient pore formation and increase the ion flux across an open pore. These transient pores can act as conduits for ion transport, which we SH image for a series of divalent cations (Cu2+, Ca2+, Ba2+, Mg2+) passing through giant unilamellar vesicle (GUV) membranes. Combining the experimental and computational results, we show that permeation through pores formed via an ion-induced electrostatic field is a viable mechanism for unassisted ion transport.

Structure of POPC Lipid Bilayers in OPLS3e Force Field

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters S_CH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured S_CH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca²⁺. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca²⁺ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Hydration Structure of Na+ and K+ Ions in Solution Predicted by Data-Driven Many-Body Potentials

The hydration structure of Na⁺ and K⁺ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na⁺ and K⁺ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

Ionic Strength and Solution Composition Dictate the Adsorption of Cell-Penetrating Peptides onto Phosphatidylcholine Membranes

Adsorption of arginine-rich positively charged peptides onto neutral zwitterionic phosphocholine (PC) bilayers is a key step in the translocation of those potent cell-penetrating peptides into the cell interior. In the past, we have shown both theoretically and experimentally that polyarginines adsorb to the neutral PC-supported lipid bilayers in contrast to polylysines. However, comparing our results with previous studies showed that the results often do not match even at the qualitative level. The adsorption of arginine-rich peptides onto 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) may qualitatively depend on the actual experimental conditions where binding experiments have been performed. In this work, we systematically studied the adsorption of R₉ and K₉ peptides onto the POPC bilayer, aided by molecular dynamics (MD) simulations and fluorescence cross-correlation spectroscopy (FCCS) experiments. Using MD simulations, we tested a series of increasing peptide concentrations, in parallel with increasing Na⁺ and Ca²⁺ salt concentrations, showing that the apparent strength of adsorption of R₉ decreases upon the increase of peptide or salt concentration in the system. The key result from the simulations is that the salt concentrations used experimentally can alter the picture of peptide adsorption qualitatively. Using FCCS experiments with fluorescently labeled R₉ and K₉, we first demonstrated that the binding of R₉ to POPC is tighter by almost 2 orders of magnitude compared to that of K₉. Finally, upon the addition of an excess of either Na⁺ or Ca²⁺ ions with R₉, the total fluorescence correlation signal is lost, which implies the unbinding of R₉ from the PC bilayer, in agreement with our predictions from MD simulations.

Shining (Infrared) Light on the Hofmeister Series: Driving Forces for Changes in the Water Vibrational Spectra in Alkali-Halide Salt Solutions

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/D₂O 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.

The lung surfactant activity probed with molecular dynamics simulations

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

The "good," the "bad," and the "hidden" in neutron scattering and molecular dynamics of ionic aqueous solutions

We characterize a concentrated 7.3 m CaCl₂ solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca²⁺-Cl^- interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca²⁺-O_w and Ca²⁺-H_w contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models.

Where Lennard-Jones Potentials Fail: Iterative Optimization of Ion-Water Pair Potentials Based on Ab Initio Molecular Dynamics Data

The use of the Lennard-Jones (LJ) potential in computer simulations of aqueous electrolyte solutions is widespread. The standard approach is to parametrize LJ potential parameters against thermodynamic solution properties, but problems in representing the local structural and dynamic properties of ion hydration shells remain. The r^-12-term in the LJ potential is responsible for this as it leads to overly repulsive ion-water interactions at short range. As a result, the LJ potential predicts blue-shifted vibrational peaks of the cations' rattling mode and too large negative ion hydration entropies. We demonstrate that cation-water effective pair potentials derived from ab initio MD data have softer short-range repulsions and represent hydration shell properties significantly better. Our findings indicate that replacing the LJ potential with these effective pair potentials offers a promising route to represent thermodynamic solution properties and local interactions of specific ions with nonpolarizable force field models.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field

In this work, we studied the effect of Li⁺, Na⁺, K⁺, Mg²⁺, and Ca²⁺ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na⁺, Ca²⁺, and SO₄ ^2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior. In particular, the diffusion of water or the salt solubility are often not well reproduced by most of the existing force-fields. Recently, the Madrid-2019 model was proposed, and it was shown that this force-field, which uses the TIP4P/2005 model for water and non-integer charges for the ions, provides a good description of a large number of properties, including the solution densities, viscosities, and the diffusion of water. In this work, we assess the performance of this force-field on the evaluation of the freezing point depression. Although the freezing point depression is a colligative property that at low salt concentrations depends solely on properties of pure water, a good model for the electrolytes is needed to accurately predict the freezing point depression at moderate and high salt concentrations. The coexistence line between ice and several salt aqueous solutions (NaCl, KCl, LiCl, MgCl₂, and Li₂SO₄) up to the eutectic point is estimated from direct coexistence molecular dynamics simulations. Our results show that this force-field reproduces fairly well the experimentally measured freezing point depression with respect to pure water freezing for all the salts and at all the compositions considered.

Ab Initio Molecular Dynamics Study of Aqueous Solutions of Magnesium and Calcium Nitrates: Hydration Shell Structure, Dynamics and Vibrational Echo Spectroscopy

Ab initio molecular dynamics simulations are performed to study the hydration shell structure, dynamics, and vibrational echo spectroscopy of aqueous Mg(NO₃)₂ and Ca(NO₃)₂ solutions. The hydration shell structure is probed through calculations of various ion-ion and ion-water radial and spatial distribution functions. On the dynamical side, calculations have been made for the hydrogen bond dynamics of hydration shells and also residence dynamics and lifetimes of water in different solvation environments. Subsequently, we looked at the dynamics of frequency fluctuations of OD modes of heavy water in different hydration environments. Specifically, the temporal decay of spectral observables of two-dimensional infrared (2DIR) spectroscopy, three pulse echo peak shift (3PEPS) measurements and also of time correlations of frequency fluctuations are calculated to investigate the dynamics of vibrational spectral diffusion of water in different hydration environments in these solutions. The OD stretch frequencies of water molecules in the vicinity of both divalent cations are found to be red-shifted and also fluctuating at a slower rate than other water molecules present in the solutions. The Mg²⁺ ions are found to be strongly hydrated which can be linked to their lower tendency to form contact ion-pairs and essentially no water exchange between the cationic hydration shells and bulk during the time scale of the current simulations. The stronger hydration of Mg²⁺ ions make their hydration shells structurally and dynamically more rigid and make the dynamics of hydrogen bonds and vibrational spectral diffusion, as revealed through spectral observables of 2DIR and 3PEPS slower than that for the Ca²⁺ ions. The structural and spectral dynamics of water molecules outside the cationic solvation shells in the Mg(NO₃)₂ solution are also found to be relatively slower than that of the Ca(NO₃)₂ solution and pure water which show the effects of stronger electric fields of Mg²⁺ ions extending beyond their first hydration shells. Also, water molecules in the hydration shells of the NO₃^- ions are found to relax at a slower rate in the Mg(NO₃)₂ solution which manifests the effect countercations have on anionic hydration shells for divalent metal nitrate solutions.

Changes in the Local Conformational States Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Seven Force Fields with and without NBFIX and ECC Corrections

Electrostatic interactions have a determining role in the conformational and dynamic behavior of polyelectrolyte molecules. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in a water solution with the most commonly used K+ or Na+ counterions, were investigated using atomistic molecular dynamics (MD) simulations. We performed a comparison of seven popular force fields, namely AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and using two common corrections for overbinding of ions, the non-bonded fix (NBFIX), and electronic continuum corrections (ECC). These corrections were originally introduced to correct for the often-reported problem concerning the overbinding of ions to the charged groups of polyelectrolytes. In this work, a comparison of the simulation results with existing experimental data revealed several differences between the investigated force fields. The data from these simulations and comparisons with previous experimental data were then used to determine the limitations and strengths of these force fields in the context of the structural and dynamic properties of anionic polyamino acids. Physical properties, such as molecular sizes, local structure, and dynamics, were studied using two types of common counterions, namely potassium and sodium. The results show that, in some cases, both the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of the ions and charged side chain groups. The local structures and dynamics are more sensitive to dihedral angle parameterization, resulting in a preference for defined monomer conformations and the type of correction used. We also provide recommendations based on the results.

Salt-in-water and water-in-salt electrolytes: the effects of the asymmetry in cation and anion valence on their properties

We investigated the structural, dynamic, energetic, and electrostatic properties of electrolytes based on the ion pairs LiCl and Li₂SO₄. Atomistic molecular dynamics simulations were used to simulate these aqueous electrolytic solutions at two different concentrations 2 M (normal) and 21 M (superconcentrated, WiSE). The effects of the valence asymmetry of the Li₂SO₄ electrolyte were also discussed for both salt concentrations. Our results differ in the physical aspect of pure electrolytes, showing the drastic effect of high concentration, in particular on the viscosity, which is dramatically increased in WiSE. This is a consequence of their reduced ionic mobility and has a direct effect on ionic conductivity. Also, our results for graphene-based supercapacitors, as indicated by some experimental work, do not indicate any better performance of WiSEs over normal electrolytes. In fact, the differences in the total capacitance, due to the concentration of ions, presented by both electrolytes are negligible. The valence asymmetry can be clearly observed in some properties but for most of them its effects could not be quantified or isolated.

Calcium Ion Binding at the Lipid-Water Interface Alters the Ion Permeability of Phospholipid Bilayers

Calcium ions (Ca²⁺) play a fundamental role in membrane-associated physiological processes. Ca²⁺ can also significantly modulate the physicochemical properties of phospholipid bilayers, but whether this occurs at physiologically relevant concentrations is difficult to determine because of the uncertainty in the reported affinity of Ca²⁺ for phospholipid bilayers. In this article, we determine the apparent affinity of Ca²⁺ for zwitterionic phospholipid bilayers using tethered bilayer lipid membranes (tBLMs) used in conjunction with swept-frequency electrical impedance spectroscopy (EIS). We report that Ca²⁺ binds to phospholipid bilayers at physiologically relevant concentrations and modulates membrane permeability. We present direct experimental evidence that this effect is governed by specific interactions with select lipid headgroup moieties, which is supported by data from molecular dynamics (MD) simulations. This is the first reported use of tBLM/EIS to estimate cation-membrane affinity. Combined with MD simulations, this technique provides a novel methodology to elucidate the molecular details of cation-membrane interactions at the water-phospholipid interface.

Surfactant-induced wettability reversal on oil-wet calcite surfaces: Experimentation and molecular dynamics simulations with scaled-charges

HYPOTHESIS

Surfactant flooding is the leading approach for reversing the wettability of oil-wet carbonate reservoirs, which is critical for the recovery of the remaining oil. Combination of molecular dynamics (MD) simulations with experiments on simplified model systems can uncover the molecular mechanisms of wettability reversal and identify key molecular properties for systematic design of new, effective chemical formulations for the enhanced oil recovery.

EXPERIMENTS/SIMULATIONS

Wettability reversal by a series of surfactant solutions was studied experimentally using contact angle measurements on aged calcite chips, and a novel MD simulation methodology with scaled-charges that provides superior description of the ionic interactions in aqueous solutions.

FINDINGS

The MD simulation results were in excellent agreement with the experiments. Cationic surfactants were the most effective in reversing the calcite wettability, resulting in complete detachment of the oil from the surface. Some nonionic surfactants also altered the wettability, but to a lesser degree, while the amphoteric and anionic surfactants had no effect. From the tested cationic surfactants, the double-tailed one was the least effective, but the experiments were inconclusive due to its poor solubility. Contributions of specific interactions to the wettability reversal process and implications for the design and optimization of surfactants for the enhanced oil recovery are discussed.

Second Harmonic Scattering Reveals Ion-Specific Effects at the SiO2 and TiO2 Nanoparticle/Aqueous Interface

Ion-specific effects play a crucial role in controlling the stability of colloidal systems and regulating interfacial processes. Although mechanistic pictures have been developed to explain the electrostatic structure of solid/water colloidal interfaces, ion-specific effects remain poorly understood. Here we quantify the average interfacial water orientation and the electrostatic surface potential around 100 nm SiO2 and TiO2 colloidal particles in the presence of NaCl, RbCl, and CaCl2 using polarimetric angle-resolved second harmonic scattering. We show that these two parameters can be used to establish the ion adsorption mechanism in a low ionic strength regime (<1 mM added salt). The relative differences between salts as a function of the ionic strength demonstrate cation- and surface-specific preferences for inner- vs outer-sphere adsorption. Compared to monovalent Rb+ and Na+, Ca2+ is found to be preferentially adsorbed as outer-sphere on SiO2 surfaces, while a dominant inner-sphere adsorption is observed for Ca2+ on TiO2. Molecular dynamics simulations performed on crystalline SiO2 and TiO2 surfaces support the experimental conclusions. This work contributes to the understanding of the electrostatic environment around colloidal nanoparticles on a molecular level by providing insight into ion-specific effects with micromolar sensitivity.

Transferable Ion Force Fields in Water from a Simultaneous Optimization of Ion Solvation and Ion-Ion Interaction

The poor performance of many existing nonpolarizable ion force fields is typically blamed on either the lack of explicit polarizability, the absence of charge transfer, or the use of unreduced Coulomb interactions. However, this analysis disregards the large and mostly unexplored parameter range offered by the Lennard-Jones potential. We use a global optimization procedure to develop water-model-transferable force fields for the ions K⁺, Na⁺, Cl^–, and Br^– in the complete parameter space of all Lennard-Jones interactions using standard mixing rules. No extra-thermodynamic assumption is necessary for the simultaneous optimization of the four ion pairs. After an optimization with respect to the experimental solvation free energy and activity, the force fields reproduce the concentration-dependent density, ionic conductivity, and dielectric constant with high accuracy. The force field is fully transferable between simple point charge/extended and transferable intermolecular potential water models. Our results show that a thermodynamically consistent force field for these ions needs only Lennard-Jones and standard Coulomb interactions.

Ion Size Dependences of the Salting-Out Effect: Reversed Order of Sodium and Lithium Ions

A general trend of the salting-out effect on hydrophobic solutes in aqueous solution is that the smaller the size of a dissolved ion, the larger the effect of reducing the solubility of a hydrophobe. An exception is that Li⁺, the smallest in alkali metal ions, has a notably weaker effect than Na⁺. To understand the reversed order in the cation series, we performed molecular dynamics simulations of aqueous solutions of salt ions and calculated the Setschenow coefficient of methane with the ionic radius of either a cation or an anion varied in a wide range. It is confirmed that the Setschenow coefficient is correlated with the packing fraction of salt solution, as observed in earlier studies, and also correlated with the partial molar volume of an ion. Analyses of correlation function integrals, packing fractions of solvation spheres, and orientations of water molecules surrounding an ion reveal the key differences in microscopic properties between the cation and anion series, which give rise to the reversed order in the cation series of the partial molar volumes of ions and ultimately that of the Setschenow coefficients.

Expression of BMP2-Hydrophobin fusion protein in the tobacco plant and molecular dynamic evaluation of its simulated model

Plants are one of the ideal models for therapeutic protein production, however the recombinant protein purification problems in them must be overcome. Bone Morphogenetic Protein2 (BMP2) is employed for the restoration and construction of bone tissues. Hydrophobin is a fungal based protein with high hydrophobic characteristics. Due to this specificity, it is suitable for the purification of chimer protein from complex solutions when is fused to a protein utilizing an aqueous two-phase (A2P) technique. The plant optimized mature human BMP2 gene was designed and evaluated by in silico method. This process involves simulating molecular dynamics using the RMSD, RMSF and Gyration radius indexes. The synthesized Hyd-BMP2 gene was cloned into a pTRAkc-ERH plasmid and Transferred into Agrobacterium (Gv3101). The Nicotiana benthamiana plant leaves were co-agroinfiltrated with HA-Hyd-BMP2 and P19-pCambia1304 containing silencing suppressor. After purification of plant extract utilizing the A2P method, the sample was subjected to SDS-PAGE and Western-blot. By in silico study, the simulated fusion protein profitably shows reasonable protein compactness and the effect of amino acid substitution on protein-protein interaction is not remarkable. Western-blotting using anti HA tag has shown that the A2P technique partially purified the two 22 kDa and 44 kDa forms of Hydrophobin-BMP2. These results confirmed the presence of monomer and dimer forms of Hydrophobin-BMP2 proteins. Moreover, the expression level of the protein using P19 silencing suppressor increased six times and to 0.018% as shown by ELISA. This study presents a fast and easy technique for the purification of transient expressed pharmaceutical proteins from plants.

Efficient Parametrization of Force Field for the Quantitative Prediction of the Physical Properties of Ionic Liquid Electrolytes

The prediction of transport properties of room-temperature ionic liquids from nonpolarizable force field-based simulations has long been a challenge. The uniform charge scaling method has been widely used to improve the agreement with the experiment by incorporating the polarizability and charge transfer effects in an effective manner. While this method improves the performance of the force fields, this prescription is ad hoc in character; further, a quantitative prediction is still not guaranteed. In such cases, the nonbonded interaction parameters too need to be refined, which requires significant effort. In this work, we propose a three-step semiautomated refinement procedure based on (1) atomic site charges obtained from quantum calculations of the bulk condensed phase; (2) quenched Monte Carlo optimizer to shortlist suitable force field candidates, which are then tested using pilot simulations; and (3) manual refinement to further improve the accuracy of the force field. The strategy is designed in a sequential manner with each step improving the accuracy over the previous step, allowing the users to invest the effort commensurate with the desired accuracy of the refined force field. The refinement procedure is applied on N,N-diethyl-N-methyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide (DEME-TFSI), a front-runner as an electrolyte for electric double-layer capacitors and single-molecule-based devices. The transferability of the refined force field is tested on N,N-dimethyl-N-ethyl-N-methoxyethoxyethylammonium bis(trifluoromethanesulfonyl)imide (N_112,2O2O1-TFSI). The refined force field is found to be better at predicting both structural and transport properties compared to the uniform charge scaling procedure, which showed a discrepancy in the X-ray structure factor. The refined force field showed quantitative agreement with structural (density and X-ray structure factor) and transport properties-diffusion coefficients, ionic conductivity, and shear viscosity over a wide temperature range, building a case for the wide adoption of the procedure.

Resolving the Equal Number Density Puzzle: Molecular Picture from Simulations of LiCl(aq) and NaCl(aq)

The change in number densities of aqueous solutions of alkali chlorides should be qualitatively predictable. Typically, as cations get larger, the number density of the solution decreases. However, aqueous solutions of lithium and sodium chloride exhibit at ambient conditions practically identical number densities at equal molalities despite different ionic sizes. Here, we provide an atomistic interpretation of this experimentally observed anomalous behavior using molecular dynamics simulations. The obtained results show that the rigidity of the Li⁺ first and second solvation shells and the associated compromised hydrogen bonding result in practically equal average water densities in the local hydration regions for Li⁺ and Na⁺ despite different sizes of the cations. In addition, in more distant regions from the cations, the water densities of these two solutions also coincide. These findings thus provide an atomistic interpretation for matching number densities of LiCl and NaCl solutions. In contrast, the number density differences between NaCl and KCl solutions as well as between LiCl and KCl solutions behave in a regular fashion with lower number densities of solutions observed for larger cations.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

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.

Heterogeneous dynamics in partially disordered proteins

Importance of disordered protein regions is increasingly recognized in biology, but their characterization remains challenging due to the lack of suitable experimental and theoretical methods. NMR experiments can detect multiple timescale dynamics and structural details of disordered protein regions, but their detailed interpretation is often difficult. Here we combine protein backbone 15N spin relaxation data with molecular dynamics (MD) simulations to detect not only heterogeneous dynamics of large partially disordered proteins but also their conformational ensembles. We observed that the rotational dynamics of folded regions in partially disordered proteins is dominated by similar rigid body rotation as in globular proteins, thereby being largely independent of flexible disordered linkers. Disordered regions, on the other hand, exhibit complex rotational motions with multiple timescales below ∼30 ns which are difficult to detect from experimental data alone, but can be captured by MD simulations. Combining MD simulations and backbone 15N spin relaxation data, measured applying segmental isotopic labeling with salt-inducible split intein, we resolved the conformational ensemble and dynamics of partially disordered periplasmic domain of TonB protein from Helicobacter pylori containing 250 residues. To demonstrate the universality of our approach, it was applied also to the partially disordered region of chicken Engrailed 2. Our results pave the way in understanding how TonB transfers energy from inner membrane to the outer membrane receptors in Gram-negative bacteria, as well as the function of other proteins with disordered domains.

Structural Effects of Cation Binding to DPPC Monolayers

Ions at the two sides of the plasma membrane maintain the transmembrane potential, participate in signaling, and affect the properties of the membrane itself. The extracellular leaflet is particularly enriched in phosphatidylcholine lipids and under the influence of Na⁺, Ca²⁺, and Cl^- ions. In this work, we combined molecular dynamics simulations performed using state-of-the-art models with vibrational sum frequency generation (VSFG) spectroscopy to study the effects of these key ions on the structure of dipalmitoylphosphatidylcholine. We used lipid monolayers as a proxy for membranes, as this approach enabled a direct comparison between simulation and experiment. We find that the effects of Na⁺ are minor. Ca²⁺, on the other hand, strongly affects the lipid headgroup conformations and induces a tighter packing of lipids, thus promoting the liquid condensed phase. It does so by binding to both the phosphate and carbonyl oxygens via direct and water-mediated binding modes, the ratios of which depend on the monolayer packing. Clustering analysis performed on simulation data revealed that changes in area per lipid or CaCl₂ concentration both affect the headgroup conformations, yet their effects are anticorrelated. Cations at the monolayer surface also attract Cl^-, which at large CaCl₂ concentrations penetrates deep to the monolayer. This phenomenon coincides with a radical change in the VSFG spectra of the phosphate group, thus indicating the emergence of a new binding mode.

Interfacial and bulk properties of concentrated solutions of ammonium nitrate

We conducted molecular dynamics (MD) simulations to calculate the density and surface tension of concentrated ammonium nitrate (AN) solutions up to the solubility limit of ammonium nitrate in water, by combining the SPC/E, SPCE/F and TIP4P/2005 water models with OPLS model for ammonium and nitrate ions. This is the first time that the properties of concentrated solutions of nitrates, especially AN, have been studied by molecular dynamics. We effectively account for the polarisation effects by the electronic continuum correction (ECC), practically realised via rescaling of the ionic charges. We found that, the full-charge force field MD simulations overestimate the experimental results, as the ions experience repulsion from the interface and prefer to remain in the subsurface layer and the bulk solution. In contrast, reducing the ionic charges results in the behaviour that fits well with the experimental data. The nitrate anions display a greater propensity for the interface than the ammonium cations. We accurately predict both the density and the rise in the surface tension of concentrated solutions of AN, recommending TIP4P/2005 for water and the scaled-charge OPLS model (OPLS/ECC) for the ions in the solutions. We observe that, the adsorption of anions to the interface accompanies their depletion in the subsurface layer, which is preferentially occupied by cations, resulting in the formation of the electric double layer. We demonstrate the ion deficiency for up to 3 Å below the surface and establish the requirement to include the polarisability effects in the OPLS model for AN. While these results confirmed the findings of the previous studies for dilute solutions, they are new in the solubility limit. Concentrated solutions exhibit a strong effect of the abundance of solute on the coordination numbers of ions and on the degree of ion pairing. Surprisingly, ion pairing decreases significantly at the interface compared with the bulk. The present study identifies OPLS/ECC, along with TIP4P/2005, to yield accurate predictions of physical properties of concentrated AN, with precision required for industrial applications, such as a formulation of emulsion and fuel-oil explosives that now predominate the civilian use of AN. An application of this model will allow one to predict the surface properties of supersaturated solutions of AN which fall outside the capability of the present laboratory experiments but are important industrially.

Spontaneous outflow efficiency of confined liquid in hydrophobic nanopores

The suspension of nanoporous particles in a nonwetting liquid provides a unique solution to the crux of superfluid, sensing, and energy conversion, yet is challenged by the incomplete outflow of intruded liquid out of nanopores for the system reusability. We report that a continuous and spontaneous liquid outflow from hydrophobic nanopores with high and stable efficiency can be achieved by regulating the confinement of solid-liquid interactions with functionalized nanopores or/and liquids. Full-scale molecular-dynamics simulations reveal that the grafted silyl chains on nanopore wall surfaces will promote the hydrophobic confinement of liquid molecules and facilitate the molecular outflow; by contrast, the introduction of ions in the liquid weakens the hydrophobic confinement and congests the molecular outflow. Both one-step and multistep well-designed quasistatic compression experiments on a series of nanopores/nonwetting liquid material systems have been performed, and the results confirm the outflow mechanism in remarkable agreement with simulations. This study offers a fundamental understanding of the outflow of confined liquid from hydrophobic nanopores, potentially useful for devising emerging nanoporous-liquid functional systems with reliable and robust reusability.

Blockage of Store-Operated Ca2+ Influx by Synta66 is Mediated by Direct Inhibition of the Ca2+ Selective Orai1 Pore

The Ca2+ sensor STIM1 and the Ca2+ channel Orai1 that form the store-operated Ca2+ (SOC) channel complex are key targets for drug development. Selective SOC inhibitors are currently undergoing clinical evaluation for the treatment of auto-immune and inflammatory responses and are also deemed promising anti-neoplastic agents since SOC channels are linked with enhanced cancer cell progression. Here, we describe an investigation of the site of binding of the selective inhibitor Synta66 to the SOC channel Orai1 using docking and molecular dynamics simulations, and live cell recordings. Synta66 binding was localized to the extracellular site close to the transmembrane (TM)1 and TM3 helices and the extracellular loop segments, which, importantly, are adjacent to the Orai1-selectivity filter. Synta66-sensitivity of the Orai1 pore was, in fact, diminished by both Orai1 mutations affecting Ca2+ selectivity and permeation of Na+ in the absence of Ca2+. Synta66 also efficiently blocked SOC in three glioblastoma cell lines but failed to interfere with cell viability, division and migration. These experiments provide new structural and functional insights into selective drug inhibition of the Orai1 Ca2+ channel by a high-affinity pore blocker.

Two-Dimensional Infrared Spectroscopy of Aqueous Solutions of Metal Nitrates: Slowdown of Spectral Diffusion in the Presence of Divalent Cations

The hydrogen-bonded network of water can be affected both structurally and dynamically by the presence of ions. In the present study, we have considered three aqueous solutions of metal nitrates to investigate the effects of divalent cations (Mg²⁺ and Ca²⁺), compared to that of monovalent Na⁺ ions, on hydrogen-bond fluctuations and vibrational spectral diffusion through calculations of linear and two-dimensional infrared spectra of these solutions at room temperature. We have employed the methods of molecular dynamics simulations using effective polarizable models of ions combined with quantum mechanical calculations of transition variables and statistical mechanical calculations of spectral response functions of vibrational spectroscopy. Divalent cations are found to have much stronger and longer-ranged effects on the structure and dynamics of the hydrogen-bonded network than that induced by the monovalent sodium ions. The blue shifts in the calculated linear spectra are found to follow the Hofmeister trend for the cations. The 2D-IR spectral lineshape and intensity corresponding to three-pulse echo peak shift (3PEPS) experiments are calculated. The timescales of these nonlinear spectral responses and also frequency-time correlations show significant slowing down of spectral diffusion for solutions containing divalent Mg²⁺ and Ca²⁺ ions compared to the corresponding dynamics of the solution containing Na⁺ ions. Unlike NaNO₃ solution, the relaxation of frequency and dipole orientational fluctuations of anion-bound water in Mg(NO₃)₂ and Ca(NO₃)₂ solutions are found to be somewhat slower than bulk water, which can be attributed to the presence of divalent cations whose effects go beyond their first solvation shells. This is also seen in the dynamics of bulk water in these solutions which is found to be notably slower for the solutions containing divalent cations than that in the NaNO₃ solution. Unlike Mg²⁺ and Ca²⁺ ions, no specific cationic effect is observed for the Na⁺ ions.

Pressure in Molecular Simulations with Scaled Charges. 1. Ionic Systems

Charge scaling, rationalized as MDEC (molecular dynamics in electronic continuum) or ECC (electronic continuum correction), has become a widely used simple approach to how to avoid self-consistent induced dipoles yet approximately take into account the effects of electronic polarizability. It has been assumed that the continuum permittivity does not depend on density; in turn, pressure is calculated by standard formulas. In this work, we elaborate a complementary approximation of density-independent molecular polarizability and derive formulas for pressure corrections within the MDEC framework; real behavior lies between these two extremes. The pressure corrections for test ionic systems are huge and negative, leading to sizable densities in constant-pressure MDEC simulations. A comparison of MDEC results with equivalent polarizable systems gives a good pressure match for a crystal but very low MDEC pressures for ionic liquids. These results witness about the importance of a correct density dependence not only of continuum permittivity in MDEC simulations but also of polarizability in polarizable simulations.

Role of Ion-Phospholipid Interactions in Zwitterionic Phospholipid Bilayer Ion Permeation

Despite the central role of Na⁺ and K⁺ in physiological processes, it is still unclear whether they interact or alter the physical properties of simple zwitterionic phospholipid bilayers at physiologically relevant concentrations. Here we report a difference in membrane permeability between Na⁺ and K⁺, as measured with electrical impedance spectroscopy and tethered bilayer lipid membranes. We reveal that the differences in membrane permeability originate from distinct ion coordination by carbonyl oxygens at the phospholipid-water interface, altering the propensity for bilayer pore formation. Molecular dynamics simulations showed differences in the coordination of Na⁺ and K⁺ at the phospholipid-water interface of zwitterionic phospholipid bilayers. The ability of Na⁺ to conscript more phospholipids with a greater number of coordinating interactions causes a higher localized energy barrier for pore formation. These results provide evidence that ion-specific interactions at the phospholipid-water interface can modulate the physical properties of zwitterionic phospholipid bilayers.

Molecular simulation of liquid-vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.