Emma Falck

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

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-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.

Coarse-Grained Model for Phospholipid / Cholesterol Bilayer

We construct a coarse-grained (CG) model for dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayers and apply it to large-scale simulation studies of lipid membranes. Our CG model is a two-dimensional representation of the membrane, where the individual lipid and sterol molecules are described by pointlike particles. The effective intermolecular interactions used in the model are systematically derived from detailed atomic-scale molecular dynamics simulations using the Inverse Monte Carlo technique, which guarantees that the radial distribution properties of the CG model are consistent with those given by the corresponding atomistic system. We find that the coarse-grained model for the DPPC/cholesterol bilayer is substantially more efficient than atomistic models, providing a speedup of approximately eight orders of magnitude. The results are in favor of formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the experimental phase diagram of the system. We also explore the limits of the coarse-grained model, and discuss the general validity and applicability of the present approach.

Impact of cholesterol on voids in phospholipid membranes

Free volume pockets or voids are important to many biological processes in cell membranes. Free volume fluctuations are a prerequisite for diffusion of lipids and other macromolecules in lipid bilayers. Permeation of small solutes across a membrane, as well as diffusion of solutes in the membrane interior are further examples of phenomena where voids and their properties play a central role. Cholesterol has been suggested to change the structure and function of membranes by altering their free volume properties. We study the effect of cholesterol on the properties of voids in dipalmitoylphosphatidylcholine (DPPC) bilayers by means of atomistic molecular dynamics simulations. We find that an increasing cholesterol concentration reduces the total amount of free volume in a bilayer. The effect of cholesterol on individual voids is most prominent in the region where the steroid ring structures of cholesterol molecules are located. Here a growing cholesterol content reduces the number of voids, completely removing voids of the size of a cholesterol molecule. The voids also become more elongated. The broad orientational distribution of voids observed in pure DPPC is, with a 30% molar concentration of cholesterol, replaced by a distribution where orientation along the bilayer normal is favored. Our results suggest that instead of being uniformly distributed to the whole bilayer, these effects are localized to the close vicinity of cholesterol molecules.

Coarse-grained model for phospholipid/cholesterol bilayer employing inverse Monte Carlo with thermodynamic constraints

The authors introduce a coarse-grained (CG) model for a lipid membrane comprised of phospholipids and cholesterol at different molar concentrations, which allows them to study systems that are approximately 100 nm in linear size. The systems are studied in the fluid phase above the main transition temperature. The effective interactions for the CG model are extracted from atomic-scale molecular dynamics simulations using the inverse Monte Carlo (IMC) technique, an approach similar to the one the authors used earlier to construct another CG bilayer model [T. Murtola et al., J. Chem. Phys. 121, 9156 (2004)]. Here, the authors improve their original CG model by employing a more accurate description of the molecular structure for the phospholipid molecules. Further, they include a thermodynamic constraint in the IMC procedure to yield area compressibilities in line with experimental data. The more realistic description of the molecular structure of phospholipids and a more accurate representation of the interaction between cholesterols and phospholipid tails are shown to improve the behavior of the model significantly. In particular, the new model predicts the formation of denser transient regions in a pure phospholipid system, a finding that the authors have verified through large scale atomistic simulations. They also find that the model predicts the formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the original model and the experimental phase diagram. However, the domains observed here are much more distinct compared to the previous model. Finally, the authors also explore the limitations of the model, discussing its advantages and disadvantages.

Influence of hydrodynamics on many-particle diffusion in 2D colloidal suspensions

Abstract.We study many-particle diffusion in 2D colloidal suspensions with full hydrodynamic interactions through a novel mesoscopic simulation technique. We focus on the behaviour of the effective scaled tracer and collective-diffusion coefficients

Dynamics and scaling of two-dimensional polymers in a dilute solution

D_{\rm T}(\rho) / D_0

Probing biomembranes with positrons.

and

Conformational analysis of lipid molecules by self-organizing maps

D_{\rm C}(\rho) / D_0

Dynamics and Scaling of Polymers in a Dilute Solution: Analytical Treatment in Two and Higher Dimensions

, respectively, where D0 is the single-particle diffusion coefficient, as a function of the density of the colloids

Insights into activation and RNA binding of trp RNA-binding attenuation protein (TRAP) through all-atom simulations.

\rho