Binding affinities and stoichiometries of Na+ and Ca2+ ions to phospholipid bilayers are of paramount significance in the properties and functionality of cellular membranes. Current estimates of binding affinities and stoichiometries of cations are, however, inconsistent due to limitations in the available experimental and computational methods. In this work, we improve the description of the binding details of Na+ and Ca2+ ions to a 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer by implicitly including electronic polarization as a mean field correction, known as the electronic continuum correction (ECC). This is applied by scaling the partial charges of a selected state-of-the-art POPC lipid model for molecular dynamics simulations. Our improved ECC-POPC model reproduces not only the experimentally measured structural parameters for the ion-free membrane, but also the response of lipid headgroup to a strongly bound cationic amphiphile, as well as the binding affinities of Na+ and Ca2+ ions. With our new model, we observe on the one side negligible binding of Na+ ions to POPC bilayer, while on the other side stronger interactions of Ca2+ primarily with phosphate oxygens, which is in agreement with the previous interpretations of the experimental spectroscopic data. The present model results in Ca2+ ions forming complexes with one to three POPC molecules with almost equal probabilities, suggesting more complex binding stoichiometries than those from simple models used to interpret the NMR data previously. The results of this work pave the way to quantitative molecular simulations with realistic electrostatic interactions of complex biochemical systems at cellular membranes.

Department of Physical Chemistry Institute of Chemical Technology 166 28 Prague 6 Czech Republic

Department of Physics Tampere University of Technology P O Box 692 FI 33101 Tampere Finland

Institute of Biotechnology University of Helsinki 00100 Helsinki Finland

Institute of Organic Chemistry and Biochemistry Academy of Sciences of the Czech Republic 117 20 Prague 6 Czech Republic

References provided by Crossref.org

NMR-Derived Salt Bridges in Insulin Analogue: Resolving Artifactual Overbinding in Molecular Dynamics via Charge Scaling

Molecular dynamics simulations unveil the aggregation patterns and salting out of polyarginines at zwitterionic POPC bilayers in solutions of various ionic strengths

Effective Inclusion of Electronic Polarization Improves the Description of Electrostatic Interactions: The prosECCo75 Biomolecular Force Field

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Stealthy Player in Lipid Experiments? EDTA Binding to Phosphatidylcholine Membranes Probed by Simulations and Monolayer Experiments

Curvature Matters: Modeling Calcium Binding to Neutral and Anionic Phospholipid Bilayers

Ionic Strength and Solution Composition Dictate the Adsorption of Cell-Penetrating Peptides onto Phosphatidylcholine Membranes

Accurate Simulations of Lipid Monolayers Require a Water Model with Correct Surface Tension

Changes in the Local Conformational States Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Seven Force Fields with and without NBFIX and ECC Corrections

Functionalization of the Parylene C Surface Enhances the Nucleation of Calcium Phosphate: Combined Experimental and Molecular Dynamics Simulations Approach