Taiki Shigematsu

Effects of Stretching Speed on Mechanical Rupture of Phospholipid/Cholesterol Bilayers: Molecular Dynamics Simulation

Rupture of biological cell membrane under mechanical stresses is critical for cell viability. It is triggered by local rearrangements of membrane molecules. We investigated the effects of stretching speed on mechanical rupture of phospholipid/cholesterol bilayers using unsteady molecular dynamics simulations. We focused on pore formation, the trigger of rupture, in a 40 mol% cholesterol-including bilayer. The unsteady stretching was modeled by proportional and temporal scaling of atom positions at stretching speeds from 0.025 to 30 m/s. The effects of the stretching speed on the critical areal strain, where the pore forms, is composed of two regimes. At low speeds (<1.0 m/s), the critical areal strain is insensitive to speed, whereas it significantly increases at higher speeds. Also, the strain is larger than that of a pure bilayer, regardless of the stretching speeds, which qualitatively agrees with available experimental data. Transient recovery of the cholesterol and phospholipid molecular orientations was evident at lower speeds, suggesting the formation of a stretch-induced interdigitated gel-like phase. However, this recovery was not confirmed at higher speeds or for the pure bilayer. The different responses of the molecular orientations may help explain the two regimes for the effect of stretching speed on pore formation.

Molecular dynamics simulations of pore formation in stretched phospholipid/cholesterol bilayers.

Molecular dynamics (MD) simulations of pore formation in stretched dipalmitoylphosphatidylcholine (DPPC) bilayers containing different concentrations of cholesterol (0, 20, 40, and 60 mol%) are presented. The stretched bilayers were simulated by constant NPZA||T MD simulations with various constant areas. The effects of the cholesterol concentration on pore formation are examined in terms of the critical areal strain where the pore is formed, the processes of pore formation, and the change in molecular orientation of the DPPC molecules by analyzing the order parameters and radial distribution functions of the DPPC molecules. With increasing cholesterol concentration, the critical areal strain initially increases, peaks at 40 mol%, and then decreases, which agrees well with the available experimental data. For the bilayers containing cholesterol, DPPC molecules become disordered at low areal strains, whereas the order slightly increases when the areal strain exceeds a certain value depending on the cholesterol concentration. For 40 mol% cholesterol, the two monolayers in the bilayer interpenetrate under high areal strains, inducing an increase of the order parameters and the peak positions of the radial distribution function compared with their states at low areal strains, indicating the formation of an interdigitated gel-phase-like structure. The transient increasing of the order of the molecular orientations may inhibit water penetration into the bilayer, resulting in increased critical areal strain in the phospholipid/cholesterol bilayers.

Stretch-Induced Interdigitation of a Phospholipid/Cholesterol Bilayer

The interdigitated gel (LβI) phase is one of the membrane phases of phospholipid molecules in which the hydrophobic tails of the phospholipid molecules penetrate the opposite leaflet of the bilayer. Recent molecular dynamics (MD) simulations have shown that interdigitation can take place as a phase transition from the liquid-ordered (Lo) phase to the LβI phase under stretching. However, there is still no conclusive experimental evidence for this process, so its existence remains controversial. In this study, to explain the transition from energy balance, we propose a free-energy model. The model consists of three energy components: the elastic deformation energy, surface energy at the bilayer-water interface, and interphase boundary energy. To determine the parameters of the model, we perform MD simulations of a stretched 1,2-dipalmitoyl- sn-glycero-3-phosphocholine/cholesterol bilayer. The phase diagrams from our model are in good agreement with those obtained from MD simulations. The energy balance among the components in the stretched bilayer quantitatively explains the stretch-induced transition. In the model, increasing the system size to that used in experiments shows that interdigitation is favorable for rigid bilayers under stretching or in alcohol solutions. These results suggest that the stretch-induced interdigitation might be observed in microscopic experiments.