Atomistic Monte Carlo simulation of lipid membranes

Daniel Wüstner, Heinz Sklenar

Research output: Contribution to journalJournal articleResearchpeer-review

Abstract

Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.
Original languageEnglish
Article numberijms15021767
JournalInternational Journal of Molecular Sciences
Volume15
Issue number2
Pages (from-to)1767-803
Number of pages37
ISSN1661-6596
DOIs
Publication statusPublished - 24. Jan 2014

Fingerprint

lipids
membranes
Monte Carlo method
simulation
closures
molecules
electron density profiles
cholesterol
assemblies
torsion
monomers
sampling
entropy
molecular dynamics
fluids
approximation
energy

Keywords

  • Monte Carlo; phospholipid bilayer; cholesterol; diffusion; coordinate transformation; entropy; sampling

Cite this

@article{85c5424b383f4c55bcd1a26d6c03f3d2,
title = "Atomistic Monte Carlo simulation of lipid membranes",
abstract = "Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.",
keywords = "Monte Carlo; phospholipid bilayer; cholesterol; diffusion; coordinate transformation; entropy; sampling",
author = "Daniel W{\"u}stner and Heinz Sklenar",
note = "Invited manuscript belonging to Special Issue on 'Computational Modelling of Biological Membranes'",
year = "2014",
month = "1",
day = "24",
doi = "10.3390/ijms15021767",
language = "English",
volume = "15",
pages = "1767--803",
journal = "International Journal of Molecular Sciences",
issn = "1661-6596",
publisher = "MDPI",
number = "2",

}

Atomistic Monte Carlo simulation of lipid membranes. / Wüstner, Daniel; Sklenar, Heinz.

In: International Journal of Molecular Sciences, Vol. 15, No. 2, ijms15021767, 24.01.2014, p. 1767-803.

Research output: Contribution to journalJournal articleResearchpeer-review

TY - JOUR

T1 - Atomistic Monte Carlo simulation of lipid membranes

AU - Wüstner, Daniel

AU - Sklenar, Heinz

N1 - Invited manuscript belonging to Special Issue on 'Computational Modelling of Biological Membranes'

PY - 2014/1/24

Y1 - 2014/1/24

N2 - Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.

AB - Biological membranes are complex assemblies of many different molecules of which analysis demands a variety of experimental and computational approaches. In this article, we explain challenges and advantages of atomistic Monte Carlo (MC) simulation of lipid membranes. We provide an introduction into the various move sets that are implemented in current MC methods for efficient conformational sampling of lipids and other molecules. In the second part, we demonstrate for a concrete example, how an atomistic local-move set can be implemented for MC simulations of phospholipid monomers and bilayer patches. We use our recently devised chain breakage/closure (CBC) local move set in the bond-/torsion angle space with the constant-bond-length approximation (CBLA) for the phospholipid dipalmitoylphosphatidylcholine (DPPC). We demonstrate rapid conformational equilibration for a single DPPC molecule, as assessed by calculation of molecular energies and entropies. We also show transition from a crystalline-like to a fluid DPPC bilayer by the CBC local-move MC method, as indicated by the electron density profile, head group orientation, area per lipid, and whole-lipid displacements. We discuss the potential of local-move MC methods in combination with molecular dynamics simulations, for example, for studying multi-component lipid membranes containing cholesterol.

KW - Monte Carlo; phospholipid bilayer; cholesterol; diffusion; coordinate transformation; entropy; sampling

U2 - 10.3390/ijms15021767

DO - 10.3390/ijms15021767

M3 - Journal article

VL - 15

SP - 1767

EP - 1803

JO - International Journal of Molecular Sciences

JF - International Journal of Molecular Sciences

SN - 1661-6596

IS - 2

M1 - ijms15021767

ER -