Molecular Dynamics Simulations of Lipid Bilayers: Major Artifacts due to Truncating Electrostatic Interactions

Abstract

We study the influence of truncating the electrostatic interactions in a fully hydrated pure dipalmitoylphosphatidylcholine (DPPC) bilayer through 20 ns molecular dynamics simulations. The computations in which the electrostatic interactions were truncated are compared to similar simulations using the Particle-Mesh Ewald (PME) technique. All examined truncation distances (1.8 to 2.5 nm) lead to major effects on the bilayer properties, such as enhanced order of acyl chains together with decreased areas per lipid. The results obtained using PME, on the other hand, are consistent with experiments. These artifacts are interpreted in terms of radial distribution functions $g(r)$ of molecules and molecular groups in the bilayer plane. Pronounced maxima or minima in g(r) appear exactly at the cutoff distance indicating that the truncation gives rise to artificial ordering between the polar phosphatidyl and choline groups of the DPPC molecules. In systems described using PME, such artificial ordering is not present.

References in corpus (2)

• Integration Schemes for Dissipative Particle Dynamics Simulations: From Softly Interacting Systems Towards Hybrid Models
• Towards Better Integrators for Dissipative Particle Dynamics Simulations

Cited by in corpus (20)

• Lessons of Slicing Membranes: Interplay of Packing, Free Area, and Lateral Diffusion in Phospholipid/Cholesterol Bilayers
• A deep potential model with long-range electrostatic interactions
• Systematic comparison of force fields for microscopic simulations of NaCl in aqueous solutions: Diffusion, free energy of hydration and structural properties
• Cationic DMPC/DMTAP Lipid Bilayers: Molecular Dynamics Study
• Sensing of Fluctuating Nanoscale Magnetic Fields Using NV Centres in Diamond
• The good, the bad and the user in soft matter simulations
• Coarse-Grained Model for Phospholipid/Cholesterol Bilayer
• Impact of Cholesterol on Voids in Phospholipid Membranes
• The $\text{u}$-series: A separable decomposition for electrostatics computation with improved accuracy
• Molecular Dynamics Study of Charged Dendrimers in Salt-Free Solution: Effect of Counterions
• Interatomic potentials: Achievements and challenges
• Lekner summations and Ewald summations for quasi-two dimensional systems
• Supported pulmonary surfactant bilayers on silica nanoparticles: Formulation, stability and impact on lung epithelial cells
• Investigating Interactions of Biomembranes and Alcohols: A Multiscale Approach
• On the efficient and accurate short-ranged simulations of uniform polar molecular liquids
• Particle-Particle, Particle-Scaling function (P3S) algorithm for electrostatic problems in free boundary conditions
• Yukawa potentials in systems with partial periodic boundary conditions I : Ewald sums for quasi-two dimensional systems
• Effect of local stress on accurate modeling of bacterial outer membranes using all-atom molecular dynamics
• Ewald sums for Yukawa potentials in quasi-two-dimensional systems
• A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations