Internal configuration and electric potential in planar negatively charged lipid head group region in contact with ionic solution

Abstract The lipid bilayer composed of negatively charged lipid 1-palmitoyl-3-oleoyl- sn -glycero-3-phosphatidylserine (POPS) in contact with an aqueous solution of monovalent salt ions was studied theoretically by using the mean-field modified Langevin–Poisson–Boltzmann (MLPB) model. The MLPB results were tested by using molecular dynamic (MD) simulations. In the MLPB model the charge distribution of POPS head groups is theoretically described by the negatively charged surface which accounts for negatively charged phosphate groups, while the positively charged amino groups and negatively charged carboxylate groups are assumed to be fixed on the rod-like structures with rotational degree of freedom. The spatial variation of relative permittivity, which is not considered in the well-known Gouy–Chapman (GC) model or in MD simulations, is thoroughly derived within a strict statistical mechanical approach. Therefore, the spatial dependence and magnitude of electric potential within the lipid head group region and its close vicinity are considerably different in the MLPB model from the GC model. The influence of the bulk salt concentration and temperature on the number density profiles of counter-ions and co-ions in the lipid head group region and aqueous solution along with the probability density function for the lipid head group orientation angle was compared and found to be in qualitative agreement in the MLPB and MD models.