Particle-based simulations of bilayer membranes: self-assembly, structural analysis, and shock-wave damage

详细信息    查看全文

• 作者：Martin O. Steinhauser ; Tanja Schindler
• 关键词：Molecular dynamics simulation ; Biological membrane ; Particle simulation ; Shock ; wave physics ; Cell research ; Multiscale modeling
• 刊名：Computational Particle Mechanics
• 出版年：2017
• 出版时间：January 2017
• 年：2017
• 卷：4
• 期：1
• 页码：69-86
• 全文大小：
• 刊物类别：Theoretical and Applied Mechanics; Computational Science and Engineering; Classical Continuum Physic
• 刊物主题：Theoretical and Applied Mechanics; Computational Science and Engineering; Classical and Continuum Physics;
• 出版者：Springer International Publishing
• ISSN：2196-4386
• 卷排序：4

文摘

We report on the results of particle-based, coarse-grained molecular dynamics simulations of amphiphilic lipid molecules in aqueous environment where the membrane structures at equilibrium are subsequently exposed to strong shock waves, and their damage is analyzed. The lipid molecules self-assemble from unbiased random initial configurations to form stable bilayer membranes, including closed vesicles. During self-assembly of lipid molecules, we observe several stages of clustering, starting with many small clusters of lipids, gradually merging together to finally form one single bilayer membrane. We find that the clustering of lipids sensitively depends on the hydrophobic interaction \(h_\mathrm{c}\) of the lipid tails in our model and on temperature T of the system. The self-assembled bilayer membranes are quantitatively analyzed at equilibrium with respect to their degree of order and their local structure. We also show that—by analyzing the membrane fluctuations and using a linearized theory— we obtain area compression moduli \(K_\mathrm{A}\) and bending stiffnesses \(\kappa _\mathrm{B}\) for our bilayer membranes which are within the experimental range of in vivo and in vitro measurements of biological membranes. We also discuss the density profile and the pair correlation function of our model membranes at equilibrium which has not been done in previous studies of particle-based membrane models. Furthermore, we present a detailed phase diagram of our lipid model that exhibits a sol–gel transition between quasi-solid and fluid domains, and domains where no self-assembly of lipids occurs. In addition, we present in the phase diagram the conditions for temperature T and hydrophobicity \(h_\mathrm{c}\) of the lipid tails of our model to form closed vesicles. The stable bilayer membranes obtained at equilibrium are then subjected to strong shock waves in a shock tube setup, and we investigate the damage in the membranes due to their interaction with shock waves. Here, we find a transition from self-repairing membranes (reducing their damage after impact) and permanent (irreversible) damage, depending on the shock front speed. The here presented idea of using coarse-grained (CG) particle models for soft matter systems in combination with the investigation of shock-wave effects in these systems is a quite new approach.