Effect of Na+ and Ca2+ Ions on a Lipid Langmuir Monolayer: An Atomistic Description by Molecular Dynamics Simulations

Studying the effect of alkali and alkaline‐earth metal cations on Langmuir monolayers is relevant from biophysical and nanotechnological points of view. In this work, the effect of Na⁺ and Ca²⁺ on a model of an anionic Langmuir lipid monolayer of dimyristoylphosphatidate (DMPA^?) is studied by molecular dynamics simulations. The influence of the type of cation on lipid structure, lipid–lipid interactions, and lipid ordering is analyzed in terms of electrostatic interactions. It is found that for a lipid monolayer in its solid phase, the effect of the cations on the properties of the lipid monolayer can be neglected. The influence of the cations is enhanced for the lipid monolayer in its gas phase, where sodium ions show a high degree of dehydration compared with calcium ions. This loss of hydration shell is partly compensated by the formation of lipid–ion–lipid bridges. This difference is ascribed to the higher charge‐to‐radius ratio q/r for Ca²⁺, which makes ion dehydration less favorable compared to Na⁺. Owing to the different dehydration behavior of sodium and calcium ions, diminished lipid–lipid coordination, lipid–ion coordination, and lipid ordering are observed for Ca²⁺ compared to Na⁺. Furthermore, for both gas and solid phases of the lipid Langmuir monolayers, lipid conformation and ion dehydration across the lipid/water interface are studied.     

金属离子(Na~+、K~+、Ca~(2+)、Mg~(2+)、Zn~(2+))与鸟嘌呤异构体配合物的稳定性

在B3LYP/6-311++G**水平上用极化连续介质模型(PCM)系统研究了金属离子(M+/2+=Na+,K+,Ca2+,Mg2+,Zn2+)和十三种鸟嘌呤异构体形成的配合物GnxM+/2+(n为鸟嘌呤异构体的编号,x表示M+/2+与鸟嘌呤异构体的结合位点)在气(g)液(a)两相中的稳定性顺序.着重探讨了液相中配合物的稳定性差异,并且从溶质-溶剂效应、结合能、形变能及异构体的相对能量等几个方面分析了造成稳定顺序发生变化的原因.报道了溶液中这五种金属离子与鸟嘌呤异构体结合形成的六种基态配合物:aG1N2,N3Na+,aG1N2,N3K+,aG1O6,N7Ca2+,aG1N2,N3Mg2+(aG1O6,N7Mg2+),aG2N3,N9Zn2+.可以看出,除了在Zn2+配合物中鸟嘌呤异构体为G2外,构成其余四种金属离子配合物的鸟嘌呤异构体都是G1,但结合位点不同.同时对气相中各类配合物稳定性也进行了系统的排序,并报道了几种较稳定的配合物,如:gG3N1,O6K+,gG5N1,O6K+,gG3N1,O6Ca2+/Mg2+,gG4O6,N7Ca2+/Mg2+.     

Comparison of structural,thermodynamic, kinetic and mass transport properties of Mg2+ ion models commonly used in biomolecular simulations

The prevalence of Mg²⁺ ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg²⁺ ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion‐centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg²⁺ model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion–water and ion–solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single‐site Mg²⁺ ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12‐6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12‐6‐4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg²⁺ ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next‐generation models that target specific binding sites in biomolecules. © 2015 Wiley Periodicals, Inc.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure and properties of L-arginine and zwitterionic L-arginine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. The effect of metal ions and water on the structure of L-arginine is examined. The effects of metal ions (Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+)) and water on structures of Arg x M(H2O)m , m = 0, 1 complexes have been determined theoretically by employing the density functional theories (DFT) and using extended basis sets. Of the three stable complexes investigated, the relative stability of the gas-phase complexes computed with DFT methods (with the exception of K(+) systems) suggests metallic complexes of the neutral L-arginine to be the most stable species. The calculations of monohydrated systems show that even one water molecule has a profound effect on the relative stability of individual complexes. Proton dissociation enthalpies and Gibbs energies of arginine in the presence of the metal cations Li(+), Na(+), K(+), Mg(2+), Ca(2+), Ni(2+), Cu(2+), and Zn(2+) were also computed. Its gas-phase acidity considerably increases upon chelation. Of the Lewis acids investigated, the strongest affinity to arginine is exhibited by the Cu(2+) cation. The computed Gibbs energies DeltaG(o) are negative, span a rather broad energy interval (from -150 to -1500 kJ/mol), and are appreciably lowered upon hydration.     

Specific Ion and Concentration Effects in Acetate Solutions with Na+, K+ and Cs+

How salt ions affect solutes and the water beyond the solvation shell is not well understood. Molecular dynamics simulations of alkali-acetate solutions were analysed here in order to examine if, and how, different cations and solute concentrations affect the water structure and the interactions between water and acetates. The results revealed that water structure is perturbed to more than 1 nm away from the acetates and that this effect is more pronounced in physiological than in molar electrolyte concentrations. Analysis of simulations of a soluble protein revealed that the water orientation is perturbed to at least 1.5 nm from the protein structure. Furthermore, modifications to the orientation of water around carboxylate side chains were shown to depend on the local environment on the protein surface, and could extend to well over 1 nm, which may have an effect on protein dynamics during MD simulations in small water boxes.     

Effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water coordination on the structure of glycine and zwitterionic glycine

Interactions between metal ions and amino acids are common both in solution and in the gas phase. Here, the effect of metal ions and water on the structure of glycine is examined. The effect of metal ions (Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) and water on structures of Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (m = 0, 2, 5) complexes have been determined theoretically by employing the hybrid B3LYP exchange-correlation functional and using extended basis sets. Selected calculations were carried out also by means of CBS-QB3 model chemistry. The interaction enthalpies, entropies, and Gibbs energies of eight complexes Gly.Mn+ (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) were determined at the B3LYP density functional level of theory. The computed Gibbs energies DeltaG degrees are negative and span a rather broad energy interval (from -90 to -1100 kJ mol(-1)), meaning that the ions studied form strong complexes. The largest interaction Gibbs energy (-1076 kJ mol(-1)) was computed for the NiGly2+ complex. Calculations of the molecular structure and relative stability of the Gly.Mn+(H2O)m and GlyZwitt.Mn+(H2O)m (Mn+ = Li+, Na+, K+, Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+; m = 0, 2, and 5) systems indicate that in the complexes with monovalent metal cations the most stable species are the NO coordinated metal cations in non-zwitterionic glycine. Divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ prefer coordination via the OO bifurcated bonds of the zwitterionic glycine. Stepwise addition of two and five water molecules leads to considerable changes in the relative stability of the hydrated species. Addition of two water molecules at the metal ion in both Gly.Mn+ and GlyZwitt.Mn+ complexes reduces the relative stability of metallic complexes of glycine. For Mn+ = Li+ or Na+, the addition of five water molecules does not change the relative order of stability. In the Gly.K+ complex, the solvation shell of water molecules around K+ ion has, because of the larger size of the potassium cation, a different structure with a reduced number of hydrogen-bonded contacts. This results in a net preference (by 10.3 kJ mol(-1)) of the GlyZwitt.K+H2O5 system. Addition of five water molecules to the glycine complexes containing divalent cations Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+ results in a net preference for non-zwitterionic glycine species. The computed relative Gibbs energies are quite high (-10 to -38 kJ mol(-1)), and the NO coordination is preferred in the Gly.Mn+(H2O)5 (Mn+ = Mg2+, Ca2+, Ni2+, Cu2+, and Zn2+) complexes over the OO coordination.     

Molecular dynamics simulations on the protonated 222. H+ and 222.2H+ cryptands in water:Endo versusexo conformations

MD simulations on the 222 cryptand, monoprotonated 222.H⁺ and diprotonated 222.2H⁺ in theendo-endo andexo-exo conformations have been performedin vacuo and in a bath of water molecules. It is found that intrinsicallyendo protonation is favoured overexo protonation due to internal N-H⁺ ... O hydrogen bonding which makes the cage more rigid. On the other hand,endo protonated forms display hydrophobic hydration compared toexo forms. For the monoprotonated 222. H⁺ endo conformer, one water molecule is hydrogen bonded inside the cage thereby forming a water cryptate. From the hydration pattern found previously for the neutral 222 cryptand and for its cation complexes, we suggest mechanisms (not involvingexo toendo conversions) for the protonation of 222 and for acid catalysed decomplexation of cryptates in theendo-endo form.     

Correlated ab initio quantum chemical study of the interaction of the Na+, Mg2+, Ca2+, and Zn2+ ions with the tautomers of cytosine

The interaction of the metal ions Na(+), Mg(2+), Ca(2+), and Zn(2+) with cytosine have been reinvestigated at the density functional, M?ller-Plesset, and coupled cluster levels of theory, including hitherto unstudied tautomeric forms. It has been found that the interaction of the metal ion has a varying and often significant effect on the stabilities of the various tautomers, in some cases making most stable rare tautomeric forms. The results have been analyzed with respect to method and role of ion in binding, and confirm that, as has been found for the base cytosine tautomers, B3LYP does not give energetics consistent with highly accurate post-SCF methods for their interaction with these metal ions.     

Structural properties of hydroxyl-substituted alkyl benzenesulfonates at the water/vapor and water/decane interfaces

Molecular dynamics simulations have been performed to investigate the structural properties of hydroxyl-substituted alkyl benzenesulfonate monolayers formed at the water/vapor and water/decane interfaces.We report a detailed study of the interfacial properties—liquid density profile,hydrogen bond structure,surfactant aggregate structure and order parameter—of the novel anionic surfactant,sodium 2-hydroxy-3-decyl-5-octylbenzenesulfonate(C10C8OHphSO3Na).Simulation results show that:with increasing number of s...     

端基改性碳纳米管膜分离Li+/Mg2+的分子动力学模拟

通过向"扶手椅"型(10,10)碳纳米管一端添加不同数量的COO-和NH+3修饰基团建立连续的碳纳米管膜模型,利用分子动力学模拟的方法研究了80 MPa水压力下Li+和Mg2+在膜中的通量和密度分布并计算了两种离子进入修饰碳纳米管的平均力势。结果表明,恰当的修饰基团添加使(10,10)碳纳米管能够有效分离Li+和Mg2+。带电基团与离子间静电作用力所产生的束缚和排斥作用使离子在纳米管内通量下降,Mg2+在修饰纳米管中的通量均为0,添加8个COO-以及4个NH+3基团均能完全阻挡两离子通过,在添加1个COO-和1个NH+3基团的情况下,Li+通量达到最大,具有最佳分离效果。因此,添加特定带电修饰基团可有效改善较大直径碳纳米管膜对Li+和Mg2+的分离性能,修饰基团电荷性质和数量对分离效果影响很大。     

Studies on Synthetic Inorganic Ion Exchangers-III Quantitative Separations of Mg2+ -Pb2+, Zn2+ - Pb2+ Cu2+ - Pb2+, Al3+ - Pb2+, Zn2+ - Cd2+, and Mg2+ - Cd2+ on Lead Antimonate Columns

Abstract

A new inorganic ion exchanger, lead antimonate has been synthesized having an Pb:Sb ratio of 1:5 and cation exchange capacity of 1.46 mequiv./g. It is fairly stable in water and dilute solutions of acids, bases and salts. Ion distribution studies on twenty metal ions have been determined on this gel at pH 1,2,3 and 5. The following mixtures have been separated: Mg²⁺ - Pb²⁺, Zn²⁺ - Pb²⁺, Zn²⁺ - Pb²⁺, Cu²⁺ - Pb²⁺, Al³⁺ - Pb²⁺, Zn²⁺ - Cd²⁺ and Mg²⁺ - Cd²⁺. Mg²⁺ and Al³⁺ were removed with 0.4 M ammonium nitrate, Cu²⁺ and Zn²⁺ with 0.4 M ammonium nitrate + 0.1M nitric acid (1:1), Pb²⁺ with 0.5M nitric acid and Cd²⁺ with 0.25M nitric acid. A tentative structure of this material is proposed on the basis of chemical analysis, pH titrations, thermogravimetry and IR spectrophotometry.     

Comprehensive study on the solvation of mono- and divalent metal cations: Li+, Na+, K+, Be2+, Mg2+ and Ca2+

Hydration of mono- and divalent metal ions (Li(+), Na(+), K(+), Be(2+), Mg(2+) and Ca(2+)) has been studied using the DFT (B3LYP), second-order M?ller-Plesset (MP2) and CCSD(T) perturbation theory as well as the G3 quantum chemical methods. Double-zeta and triple-zeta basis sets containing both (multiple) polarization and diffuse functions were applied. Total and sequential binding energies are evaluated for all metal-water clusters containing 1-6 water molecules. Total binding energies predicted at lower levels of theory are compared with those from the high level G3 calculations, whereas the sequential binding energies are compared with available experimental values. An increase in the quality of the basis set from double-zeta to triple-zeta has a significant effect on the sequential binding energies, irrespective of the geometries used. Within the same group (I or II), the sequential binding energy predictions at the MP2 and B3LYP vary appreciably. We noticed that, for each addition of a water molecule, the change of the M-O distance in metal-water clusters is higher at the B3LYP than at the MP2 level. The charge of the metal ion decreases monotonically as the number of water molecules increase in the complex.     

Different binding thermodynamics of Ni2+, Cu2+, and Zn2+ to bacitracin A1 determined by capillary electrophoresis

Thermodynamics of the binding of Ni(2+), Cu(2+) and Zn(2+) to bacitracin A(1) was studied by capillary electrophoresis measuring the peptide effective mobility at different pH in the presence of increasing concentration of the three ligands. The affinity follows the order Ni(2+) > Cu(2+) > Zn(2+), with association constant values of (2.3 +/- 0.1)x10(4), (4.9 +/- 0.2)x10(3), and (1.5 +/- 0.1)x10(3) M(-1), respectively. The only model able to rationalize mobility data implies that metal ion binds to the P(0) peptide form. Moreover, mobility values indicated a change of bacitracin A(1) acidic properties on Ni(2+) and Cu(2+) binding, with a shift of the pK(a) of N-terminal Ile-1 from 7.6 to about 5 and of the pK(a) of the delta-amino group of D-Orn-7 from 9.7 to about 7. Even though on Zn(2+) binding a shift of the N-terminal Ile-1 pK(a) was observed, restrictions in the pH range suitable for investigation, due to precipitation phenomena, did not allow establish if the shift of D-Orn-7 lateral chain pK(a) also occurred. Nonetheless, if present, the shift should be limited to the 7.8-9.7 range. Mobility data indicated that the Stokes radius of the complexes is ca. 3 A lower than that of the free peptide. The present results indicate that metal-ion binding to bacitracin A(1) is more complex than previously assumed.     

Modeling electrified metal/water interfaces from ab initio molecular dynamics: Structure and Helmholtz capacitance

The structure and potential distribution of electric double layers (EDLs) are of close relevance to the performance of electrode materials. In the past years, despite tremendous efforts devoted to this topic, an atomistic picture of the EDL is still lacking, let alone understanding on how the EDL structure is related to the dielectric property of interface water. In this article, we briefly review the recent progress in modeling electrified metal/water interfaces using ab initio molecular dynamics (AIMD). The ab initio methods for EDL modeling is firstly summarized, and then we discuss the structures of interface water on metal electrodes at different potential conditions. Moreover, we illustrate the potential-dependent behavior of chemisorbed water on Pt(111) surface and its relationship with the peak of the differential Helmholtz capacitance observed by experiment. At last, we give some perspective for future development in ab initio modeling of electrochemical interfaces.     

Ionophoric properties of atropine: complexation and transport of Na+, K+, Mg2+ and Ca2+ ions across a liquid membrane

The activity of atropine on the complexation and transport of Na(+), K(+), Mg(2+) and Ca(2+) ions across a liquid membrane was investigated using a spectrophotometric method. Atropine is a natural drug that blocks muscarinic receptors. It is a competitive antagonist of the action of acetylcholine and other muscarinic agonists. Atropine is shown to extract Na(+), K(+), Mg(2+) and Ca(2+) ions from an aqueous phase into an organic one with a preference for Ca(2+) ions. According to a kinetic study, divalent cations (Mg(2+) and Ca(2+)) are more rapidly transported than monovalent ones (Na(+) and K(+)). In both complexation and transport, the flux of the ions increases with the increase of atropine concentration. Atropine might act on the membrane permeability; its complexation and ionophoric properties shed new lights on its therapeutic properties.     

Segmental dynamics of PEO/LiClO4 complex crystals and their influence on the Li+-ion transportation in crystal lattices: A 13C solid-state NMR approach

The structure and chain dynamics of polyethylene oxide (PEO)/LiClO(4) complex crystals were studied by employing state-of-the-art solid-state (13)C NMR techniques. Remarkable helical jump motions of the PEO segments in the complex crystals at ambient temperatures were clearly demonstrated. The jump motions are believed to induce the transportation of the coordinated Li(+) ions along the crystallographic c axis, providing a novel mechanism of ionic conductivity of the complex crystals. In addition, this work also shows that solid-state high-resolution (13)C NMR spectroscopy can be a powerful and general tool for elucidating phase structures, dynamics, and subsequently the conduction mechanism of crystalline polymer electrolytes.     

Systematic comparison of force fields for microscopic simulations of NaCl in aqueous solutions: diffusion, free energy of hydration, and structural properties

In this article we compare different force fields that are widely used (Gromacs, Charmm-22/x-Plor, Charmm-27, Amber-1999, OPLS-AA) in biophysical simulations containing aqueous NaCl. We show that the uncertainties of the microscopic parameters of, in particular, sodium, and, to a lesser extent, chloride, translate into large differences in the computed radial-distribution functions. This uncertainty reflect the incomplete experimental knowledge of the structural properties of ionic aqueous solutions at finite molarity. We discuss possible implications on the computation of potential of mean force and effective potentials.     

高锰酸钾-鲁米诺反应体系中碱土金属离子Mg2+,Ca2+,Sr2+,Ba2+化学发光行为研究

将被认为没有化学发光活性的第二主族(碱土金属)离子Mg^2 ,Ca^2 ,Sr^2 ,Ba^2 溶液注入到已充分反应的高锰酸钾与鲁米诺混合液中时,又发生了新的化学发光反应,并检测到强的化学发光信号．在对有关反应的动力学性质、化学发光光谱、紫外可见光谱及其它一系列实验研究的基础上,提出了可能的化学发光反应机理．同时,优化了反应条件,评价了这一反应用于Mg^2 ,Ca^2 ,Sr^2 ,Ba^2 分析的可行性．