Interplay of Interfacial Viscosity, Specific-Ion, and Impurity Adsorption Determines Zeta Potentials of Phospholipid...

Ion-specific induced changes of the ζ-potential of phospholipid vesicles are commonly used to quantify the affinity of different ions to the lipid interface. The negative ζ-potential of zwitterionic net-neutral phospholipid vesicles in neat water, which changes sign and increases in solutions of NaCl or KCl, is a phenomenon consistently observed in experiments but not fully understood theoretically. Using atomistic molecular dynamics simulations in the presence of applied electric fields which drive electroosmotic flows, in combination with an electrostatic continuum model based on the modified Poisson-Boltzmann and Helmholtz-Smoluchowski equations, we study the electrokinetic and electrostatic properties as well as the specific ion affinities to the phospholipid-water interface, in order to resolve these puzzling observations. Our modified continuum equations account for the dielectric profile at the lipid-water interface, ion-specific interactions between ions and the lipid-water interface, and the interfacial viscosity profile, which are all extracted from our atomistic simulations and rather accurately predict ion-density and electrostatic-potential distributions as well as ζ-potentials in comparison with our atomistic simulations. Our continuum model can explain experimental ζ-potentials only when we assume minute amounts of surface-active anionic impurities in the aqueous solution. In fact, the amount of impurities needed to explain the experimental data increases linearly with the salt concentration, suggesting that surface-active species, which might be already present in the lab water or lipid samples, could further be introduced through the added salt.