Molecular dynamics simulation of dipalmitoylphosphatidylserine bilayer with Na+ counterions

We performed a molecular dynamics simulation of dipalmitoylphosphatidylserine (DPPS) bilayer with Na+ counterions. We found that hydrogen bonding between the NH3+ group and the phosphate group leads to a reduction in the area per headgroup when compared to the area in dipalmitoylphosphatidylcholine bilayer. The Na+ ions bind to the oxygen in the carboxyl group of serine, thus giving rise to a dipolar bilayer similar to dipalmitoylphosphatidylethanolamine bilayer. The results of the simulation show that counterions play a crucial role in determining the structural and electrostatic properties of DPPS bilayer.