Stability of a melittin pore in a lipid bilayer: a molecular dynamics study

We have investigated the configuration and the stability of a single membrane pore bound by four melittin molecules and embedded in a fully hydrated bilayer lipid membrane. We used molecular dynamics simulations up to 5.8 ns. It is found that the initial tetrameric configuration decays with increasing time into a stable trimer and one monomer. This continuous transformation is accompanied by a lateral expansion of the aqueous pore exhibiting a final size comparable to experimental findings. The expansion-induced formation of an interface between the pore-lining acyl chains of the lipids and the pore water ("hydrophobic pore") is transformed into an energetically more favorable toroidal pore structure where some lipid heads are translocated from the rim to the central part of the interface ("hydrophilic pore"). The expansion of the pore is supported by the electrostatic repulsion among the alpha-helices. It is hypothesized that pore growth, and hence cell lysis, is induced by a melittin-mediated line tension of the pore.