Eduardo R.A. Lima

Specific ion adsorption and surface forces in colloid science.

Mean-field theories that include nonelectrostatic interactions acting on ions near interfaces have been found to accommodate many experimentally observed ion specific effects. However, it is clear that this approach does not fully account for the liquid molecular structure and hydration effects. This is now improved by using parametrized ionic potentials deduced from recent nonprimitive model molecular dynamics (MD) simulations in a generalized Poisson-Boltzmann equation. We investigate how ion distributions and double layer forces depend on the choice of background salt. There is a strong ion specific double layer force set up due to unequal ion specific short-range potentials acting between ions and surfaces.

Co-Ion and Ion Competition Effects: Ion Distributions Close to a Hydrophobic Solid Surface in Mixed Electrolyte Solutions

We consider within a modified Poisson-Boltzmann theory an electrolyte, with different mixtures of NaCl and NaI, near uncharged and charged solid hydrophobic surfaces. The parametrized potentials of mean force acting on Na+, Cl-, and I- near an uncharged self-assembled monolayer were deduced from molecular simulations with polarizable force fields. We study what happens when the surface presents negative charges. At moderately charged surfaces, we observe strong co-ion adsorption and clear specific ion effects at biological concentrations. At high surface charge densities, the co-ions are pushed away from the interface. We predict that Cl- ions can also be excluded from the surface by increasing the concentration of NaI. This ion competition effect (I- versus Cl-) may be relevant for ion-specific partitioning in multiphase systems where polarizable ions accumulate in phases with large surface areas.

The influence of ion binding and ion specific potentials on the double layer pressure between charged bilayers at low salt concentrations

Measurements of surface forces between double-chained cationic bilayers adsorbed onto molecularly smooth mica surfaces across different millimolar salt solutions have revealed a large degree of ion specificity [Pashley et al., J. Phys. Chem. 90, 1637 (1986)]. This has been interpreted in terms of highly specific anion binding to the adsorbed bilayers. We show here that inclusion in the double layer theory of nonspecific ion binding and ion specific nonelectrostatic potentials acting between ions and the two surfaces can account for the phenomenon. It also gives the right Hofmeister series for the double layer pressure.

Finite volume solution of the modified Poisson–Boltzmann equation for two colloidal particles

The double layer forces between spherical colloidal particles, according to the Poisson-Boltzmann (PB) equation, have been accurately calculated in the literature. The classical PB equation takes into account only the electrostatic interactions, which play a significant role in colloid science. However, there are at, and above, biological salt concentrations other non-electrostatic ion specific forces acting that are ignored in such modelling. In this paper, the electrostatic potential profile and the concentration profile of co-ions and counterions near charged surfaces are calculated. These results are obtained by solving the classical PB equation and a modified PB equation in bispherical coordinates, taking into account the van der Waals dispersion interactions between the ions and both surfaces. Once the electrostatic potential is known we calculate the double layer force between two charged spheres. This is the first paper that solves the modified PB equation in bispherical coordinates. It is also the first time that the finite volume method is used to solve the PB equation in bispherical coordinates. This method divides the calculation domain into a certain number of sub-domains, where the physical law of conservation is valid, and can be readily implemented. The finite volume method is implemented for several geometries and when it is applied to solve PB equations presents low computational cost. The proposed method was validated by comparing the numerical results for the classical PB calculations with previous results reported in the literature. New numerical results using the modified PB equation successfully predicted the ion specificity commonly observed experimentally.

Anion-Specific Partitioning in Two-Phase Finite Volume Systems: Possible Implications for Mechanisms of Ion Pumps

In two-phase finite volume systems of electroneutral phospholipids, the electrolyte concentration is different in the two phases. The partitioning is highly anion-specific, a phenomenon not accounted for by classical electrolyte theories. It is explained if ionic dispersion forces that lead to specific ion binding are taken into account. The mechanism provides a contribution to active ion pumps not previously considered.

Development of an AMBER-compatible transferable force field for poly(ethylene glycol) ethers (glymes)

AbstractAn all-atom force field consistent with the general AMBER force field (GAFF) format for poly(ethylene glycol) dimethyl ether (diglyme or G2) was developed by fitting to experimental liquid densities and dielectric constants. Not surprisingly, the new force field gives excellent agreement with experimental liquid phase densities and dielectric constants over a wide temperature range. Other dynamic and thermodynamic properties of liquid G2 such as its self-diffusion coefficient, shear viscosity, and vaporization enthalpy were also calculated and compared to experimental data. For all of the properties studied, the performance of the proposed new force field is better than that of the standard GAFF force field. The force field parameters were transferred to model two other poly(ethylene glycol) ethers: monoglyme (G1) and tetraglyme (G4). The predictive ability of the modified force field for G1 and G4 was significantly better than that of the original GAFF force field. The proposed force field provides an alternative option for the simulation of mixtures containing glymes using GAFF-compatible force fields, particularly for electrochemical applications. The accuracy of a previously published force field based on the OPLS-AA format and the accuracies of two modified versions of that force field were also examined for G1, G2, and G4. It was found that the original OPLS-AA force field is superior to the modified versions of it, and that it has a similar accuracy to the proposed new GAFF-compatible force field. Graphical abstractTransferability of an AMBER-compatible force field parameterized for G2 to other glymes

Attractive double-layer forces between neutral hydrophobic and neutral hydrophilic surfaces.

The interaction between surface patches of proteins with different surface properties has a vital role to play driving conformational changes in proteins in different salt solutions. We demonstrate the existence of ion-specific attractive double-layer forces between neutral hydrophobic and hydrophilic surfaces in the presence of certain salt solutions. This is performed by solving a generalized Poisson-Boltzmann equation for two unequal surfaces. In the calculations, we utilize parametrized ion-surface potentials and dielectric-constant profiles deduced from recent non-primitive-model molecular dynamics simulations that partially account for molecular structure and hydration effects.

Ion-specific thermodynamical properties of aqueous proteins

Ion-specific interactions between two colloidal particles are calculated using a modified Poisson-Boltzmann (PB) equation and Monte Carlo (MC)simulations. PB equations present good results of ionic concentration profiles around a macroion, especially for salt solutions containing monovalent ions. These equations include not only electrostatic interactions, but also dispersion potentials originated from polarizabilities of ions and proteins. This enables us to predict ion-specific properties of colloidal systems. We compared results obtained from the modified PB equation with those from MC simulations and integral equations. Phase diagrams and osmotic second virial coefficients are also presented for different salt solutions at different pH and ionic strengths, in agreement with the experimental results observed Hofmeister effects. In order to include the water structure and hydration effect, we have used an effective interaction obtained from molecular dynamics of each ion and a hydrophobic surface combined with PB equation. The method has been proved to be efficient and suitable for describing phenomena where the water structure close to the interface plays an essential role. Important thermodynamic properties related to protein aggregation, essential in biotechnology and pharmaceutical industries, can be obtained from the method shown here.

Membrane Potential and Ion Partitioning in an Erythrocyte Using the Poisson–Boltzmann Equation

In virtually all mammal cells, we can observe a much higher concentration of potassium ions inside the cell and vice versa for sodium ions. Classical theories ignore the specific ion effects and the difference in the thermodynamic reference states between intracellular and extracellular environments. Usually, this differential ion partitioning across a cell membrane is attributed exclusively to the active ion transport. Our aim is to investigate how much the dispersion forces contribute to active ion pumps in an erythrocyte (red blood cell) as well as the correction of chemical potential reference states between intracellular and extracellular environments. The ionic partition and the membrane potential in an erythrocyte are analyzed by the modified Poisson-Boltzmann equation, considering nonelectrostatic interactions between ions and macromolecules. Results show that the nonelectrostatic potential calculated by Lifshitz theory has only a small influence with respect to the high concentration of K(+) in the intracellular environment in comparison with Na(+).

The electrostatic behavior of the bacterial cell wall using a smoothing function to describe the charge-regulated volume charge density profile.

The Donnan potential can be observed in many biological systems due to the presence of polyelectrolytes as proteins and nucleic acids. The aim of this work was to present a useful tool to describe the fixed and charge-regulated volume charge density profile through the use of a smoothing function and to obtain the electrostatic potential profile as well as the Donnan potential of this system by solving Poisson-Boltzmann (PB) equation. When we use the smoothing function, the Donnan potential arises automatically from the solution of only one Poisson-Boltzmann equation and it is not necessary to impose this potential for treating charged system in the presence of a membrane. The electrostatic behavior across the Bacillus brevis wall considering the dependence on the ionization of the cell wall functional groups as a function of the solution pH was analyzed. An important issue was to show that potentiometric titration data could be used together with the Poisson-Boltzmann equation to predict the electrostatic behavior (e.g., zeta potential) of the bacterial cell surface.