pith. sign in

arxiv: cond-mat/0211650 · v1 · submitted 2002-11-28 · ❄️ cond-mat.soft · physics.bio-ph

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

classification ❄️ cond-mat.soft physics.bio-ph
keywords bilayerelectrostaticinteractionsmolecularsimulationsartifactsartificialdppc
0
0 comments X
read the original 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.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.