Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations

Lipid adhesion forces can be measured using several experimental techniques, but none of these techniques provide insight on the atomic level. Therefore, we performed extensive nonequilibrium molecular dynamics simulations of a phospholipid membrane in the liquid-crystalline phase out of which individual lipid molecules were pulled. In our method, as an idealization of the experimental setups, we have simply attached a harmonic spring to one of the lipid headgroup atoms. Upon retraction of the spring, the force needed to drag the lipid out of the membrane is recorded. By simulating different retraction rates, we were able to investigate the high pull rate part of the dynamical spectrum of lipid adhesion forces. We find that the adhesion force increases along the unbinding path, until the point of rupture is reached. The maximum value of the adhesion force, the rupture force, decreases as the pull rate becomes slower, and eventually enters a friction-dominated regime. The computed bond lengths depend on the rate of rupture, and show some scatter due to the nonequilibrium nature of the experiment. On average, the bond length increases from approximately 1.7 nm to 2.3 nm as the rates go down. Conformational analyses elucidate the detailed mechanism of lipid-membrane bond rupture. We present results of over 15 ns of membrane simulations. Implications for the interpretation and understanding of experimental rupture data are discussed.