Molecular-dynamic simulation of phospholipid bilayers: Ion distribution at the surface of neutral and charged bilayer in the liquid crystalline state

The electric field and ion distribution at the surface of neutral and charged lipid bilayers (BeCl₂ and dipalmitoyl phosphatidylcholine/dipalmitoyl phosphatidylserine (DPPC/DPPS) + KCl) were studied with molecular dynamic (MD) methods. It is shown that the contributions of lipid molecules, water and ions to the electric potential compensate each other in the region of the diffuse double layer and decrease the potential value close to zero. It is also demonstrated that the ion distribution at the charged surface is determined not only by the electrostatic ion-medium interaction. The total energy of this interaction was compared with the potential of mean ion force. It was shown that cations and anions have a different effect on the state of water molecules at the surface. The order parameter of water in the system DPPC + BeCl₂ and the Clion distribution have the extremum at the distance of 10 α atoms of the phospholipid glycerol. This position was chosen as the “electrical” interface of the electrical double layer (EDL) for all lipid systems studied. The potential of mean force of counter ions in EDL allows us to obtain the value of potential at the lipid surface suitable for experimental test of the MD data. This surface potential and surface charge density was found from MD simulation different electrolyte concentrations and DPPS content of 20, 40 and 60% in the mixture with DPPC and was shown to be in a good agreement with the Gouy-Chapman-Stern model upon fitting parameters close to their experimental values.