Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers
2009 (English)In: Journal of Chemical Theory and Computation, ISSN 1549-9618, E-ISSN 1549-9626, Vol. 5, no 8, 2125-2134 p.Article in journal (Refereed) Published
The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers; is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na+ ions than K+ ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na+ and K+ ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to sigma(MO) (M represents Na+/K+), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control. procedure is provided.
Place, publisher, year, edition, pages
2009. Vol. 5, no 8, 2125-2134 p.
IdentifiersURN: urn:nbn:se:kth:diva-18665DOI: 10.1021/ct9000763ISI: 000268720100017ScopusID: 2-s2.0-73549108322OAI: oai:DiVA.org:kth-18665DiVA: diva2:336712
QC 201005252010-08-052010-08-052010-12-17Bibliographically approved