This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.107.116335v1 94/9/3565    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Google Scholar Articles by Lee, S.-J. Articles by Baker, N. A. PubMed PubMed Citation Articles by Lee, S.-J. Articles by Baker, N. A.

Molecular Dynamics Simulations of Asymmetric NaCl and KCl Solutions Separated by Phosphatidylcholine Bilayers: Potential Drops and Structural Changes Induced by Strong Na⁺-Lipid Interactions and Finite Size Effects

Sun-Joo Lee ^*, Yuhua Song  and Nathan A. Baker ^*

^* Department of Biochemistry and Molecular Biophysics, Center for Computational Biology, Washington University in St. Louis, St. Louis, Missouri; and Department of Biomedical Engineering, University of Alabama at Birmingham, Birmingham, Alabama

Correspondence: Address reprint requests to Nathan A. Baker, Tel.: 314-362-2040 E-mail: baker{at}ccb.wustl.edu.

Differences of ionic concentrations across lipid bilayers are some of the primary energetic driving forces for cellular electrophysiology. While macroscopic models of asymmetric ionic solutions are well-developed, their connection to ion, water, and lipid interactions at the atomic scale are much more poorly understood. In this study, we used molecular dynamics to examine a system of two chambers of equal ionic strength, but differing amounts of NaCl and KCl, separated by a lipid bilayer. Our expectation was that the net electrostatic potential difference between the two chambers should be small or zero. Contrary to our expectation, a large potential difference (–70 mV) slowly evolved across the two water chambers over the course of our 172-ns simulation. This potential primarily originated from strong Na⁺ binding to the carbonyls of the phosphatidylcholine lipids. This ion adsorption also led to significant structural and mechanical changes in the lipid bilayer. We discuss this surprising result in the context of indirect experimental evidence for Na⁺ interaction with bilayers as well as potential caveats in current biomembrane simulation methodology, including force-field parameters and finite size effects.

This article has been cited by other articles:

				E. J. Denning and T. B. Woolf Double Bilayers and Transmembrane Gradients: A Molecular Dynamics Study of a Highly Charged Peptide Biophys. J., October 1, 2008; 95(7): 3161 - 3173. [Abstract] [Full Text] [PDF]