Abstract
Gramicidin A (gA) is a short β-helical peptide known to form conducting channels in lipid membranes because of transbilayer dimerization. The gA conducting dimer, being shorter than the lipid bilayer thickness, deforms the membrane in its vicinity, and the bilayer elastic energy contributes to the gA dimer formation energy. Likewise, membrane incorporation of a gA monomer, which is shorter than the lipid monolayer thickness, creates a void, thereby forcing surrounding lipid molecules to tilt to fill it. The energy of membrane deformation was calculated in the framework of the continuum elasticity theory, taking into account splay, tilt, lateral stretching/compression, Gaussian splay deformations, and external membrane tension. We obtained the interaction energy profiles for two gA monomers located either in the same or in the opposite monolayers. The profiles demonstrated the long-range attraction and short-range repulsion behavior of the monomers resulting from the membrane deformation. Analysis of the profile features revealed conditions under which clusters of gA monomers would not dissipate because of diffusion. The calculated dependence of the dimer formation and decay energy barriers on the membrane elastic properties was in good agreement with the available experimental data and suggested an explanation for a hitherto contentious phenomenon.