Phosphatidylserine (PS) lipids are important signaling molecules and the most common negatively charged lipids in eukaryotic membranes. The signaling can be often regulated by calcium, but its interactions with PS headgroups are not fully understood. Classical molecular dynamics (MD) simulations can potentially give detailed description of lipid-ion interactions, but the results strongly depend on the used force field. Here, we apply the electronic continuum correction (ECC) to the Amber Lipid17 parameters of 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipid to improve its interactions with K+, Na+, and Ca2+ ions. The partial charges of the headgroup, glycerol backbone, and carbonyls of POPS, bearing a unit negative charge, were scaled with a factor of 0.75, derived for monovalent ions, and the Lennard-Jones σ parameters of the same segments were scaled with a factor of 0.89. The resulting ECC-POPS model gives more realistic interactions with Na+ and Ca2+ cations than the original Amber Lipid17 parameters when validated using headgroup order parameters and the "electrometer concept". In ECC-lipids simulations, populations of complexes of Ca2+ cations with more than two PS lipids are negligible, and interactions of Ca2+ cations with only carboxylate groups are twice more likely than with only phosphate groups, while interactions with carbonyls almost entirely involve other groups as well. Our results pave the way for more realistic MD simulations of biomolecular systems with anionic membranes, allowing signaling processes involving PS and Ca2+ to be elucidated.

Groningen Biomolecular Sciences and Biotechnology Institute and The Zernike Institute for Advanced Materials University of Groningen 9747 AG Groningen The Netherlands

Institute of Biotechnology University of Helsinki Helsinki FI 00014 Finland

Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo nám 542 2 CZ 16610 Prague 6 Czech Republic

NMR Group Institut for Physics Martin Luther University Halle Wittenberg 06120 Halle Germany

Citace poskytuje Crossref.org

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters via Ab Initio Molecular Dynamics Simulations To Accurately Model Glycosaminoglycan Electrostatic Interactions

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

Overlay databank unlocks data-driven analyses of biomolecules for all

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

Modulation of Anionic Lipid Bilayers by Specific Interplay of Protons and Calcium Ions

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