Lipid bilayer mechanics

Lipid bilayer mechanics is the study of the physical material properties of lipid bilayers, classifying bilayer behavior with stress and strain rather than biochemical interactions. Local point deformations such as membrane protein interactions are typically modelled with the complex theory of biological liquid crystals but the mechanical properties of a homogeneous bilayer are often characterized in terms of only three mechanical elastic modulus: the area expansion modulus K_a, a bending modulus K_b and an edge energy $\Lambda$. For fluid bilayers the shear modulus is by definition zero, as the free rearrangement of molecules within plane means that the structure will not support shear stresses. These mechanical properties affect several membrane-mediated biological processes. In particular, the values of K_a and K_b affect the ability of proteins and small molecules to insert into the bilayer. Bilayer mechanical properties have also been shown to alter the function of mechanically activated ion channels.

Area expansion Modulus

Further information: Elasticity of cell membranes

See also: Elastic Modulus

Since lipid bilayers are essentially a two dimensional structure, K_a is typically defined only within the plane. Intuitively, one might expect that this modulus would vary linearly with bilayer thickness as it would for a thin plate of isotropic material. In fact this is not the case and K_a is only weakly dependent on bilayer thickness. The reason for this is that the lipids in a fluid bilayer rearrange easily so, unlike a bulk material where the resistance to expansion comes from intermolecular bonds, the resistance to expansion in a bilayer is a result of the extra hydrophobic area exposed to water upon pulling the lipids apart. Based on this understanding, a good first approximation of K_a for a monolayer is 2γ, where gamma is the surface tension of the water-lipid interface. Typically gamma is in the range of 20-50mJ/m². To calculate K_a for a bilayer it is necessary to multiply the monolayer value by two, since a bilayer is composed of two monolayer leaflets. Based on this calculation, the estimate of K_a for a lipid bilayer should be 80-200 mN/m (note: N/m is equivalent to J/m²). It is not surprising given this understanding of the forces involved that studies have shown that K_a varies strongly with solution conditions but only weakly with tail length and unsaturation.

The compression modulus is difficult to measure experimentally because of the thin, fragile nature of bilayers and the consequently low forces involved. One method utilized has been to study how vesicles swell in response to osmotic stress. This method is, however, indirect and measurements can be perturbed by polydispersity in vesicle size. A more direct method of measuring K_a is the pipette aspiration method, in which a single giant unilamellar vesicle (GUV) is held and stretched with a micropipette. More recently, atomic force microscopy (AFM) has been used to probe the mechanical properties of suspended bilayer membranes, but this method is still under development.

One concern with all of these methods is that, since the bilayer is such a flexible structure, there exist considerable thermal fluctuations in the membrane at many length scales down to sub-microscopic. Thus, forces initially applied to an unstressed membrane are not actually changing the lipid packing but are rather “smoothing out” these undulations, resulting in erroneous values for mechanical properties. This can be a significant source of error. Without the thermal correction typical values for Ka are 100-150 mN/m and with the thermal correction this would change to 220-270 mN/m.

Bending Modulus

See also: Beam (structure)

Bending modulus is defined as the energy required to deform a membrane from its intrinsic curvature to some other curvature. For an ideal bilayer the intrinsic curvature is zero, so this expression is somewhat simplified. The bending modulus, compression modulus and bilayer thickness are related by $K_b= K_a *t$ such that if two of these parameters are known the other can be calculated. This relationship derives from the fact that to bend the inner face must be compressed and the outer face must be stretched. The thicker the membrane, the more each face must deform to accommodate a given curvature (see bending moment). Many of the values for K_a in literature have actually been calculated from experimentally measured values of K_b and t. This relation holds only for small deformations, but this is generally a good approximation as most lipid bilayers can support only a few percent strain before rupturing.

Cross section view of the structures that can be formed by phospholipids in aqueous solutions

Curvature

Only certain classes of lipids can form bilayers. Two factors primarily govern whether a lipid will form a bilayer or not: solubility and shape. For a self assembled structure such as a bilayer to form, the lipid should have a low solubility in water, which can also be described as a low critical micelle concentration (CMC). Above the CMC, molecules will aggregate and form larger structures such as bilayers, micelles or inverted micelles.

The primary factor governing which structure a given lipid forms is its shape (i.e.- its intrinsic curvature). Intrinsic curvature is defined by the ratio of the diameter of the head group to that of the tail group. For two-tailed PC lipids, this ratio is nearly one so the intrinsic curvature is nearly zero. Other headgroups such as PS and PE are smaller and the resulting diacyl (two-tailed) lipids thus have a negative intrinsic curvature. Lysolipids tend to have positive spontaneous curvature because they have one rather than two alkyl chains in the tail region. If a particular lipid has too large a deviation from zero intrinsic curvature it will not form a bilayer.

Edge Energy

Schematic showing two possible conformations of the lipids at the edge of a pore. In the top image the lipids have not rearranged, so the pore wall is hydrophobic. In the bottom image some of the lipid heads have bent over, so the pore wall is hydrophilic.

Edge energy is the energy per unit length of a free edge contacting water. This can be thought of as the work needed to create a hole in the bilayer of unit length L. The origin of this energy is the fact that creating such an interface exposes some of the lipid tails to water, which is unfavorable. $\Lambda$ is also an important parameter in biological phenomena as it regulates the self-healing properties of the bilayer following electroporation or mechanical perforation of the cell membrane. Unfortunately, this property is both difficult to measure experimentally and to calculate. One of the major difficulties in calculation is that the structural properties of this edge are not known. The simplest model would be no change in bilayer orientation, such that the full length of the tail is exposed. This is a high energy conformation and, to stabilize this edge, it is likely that some of the lipids rearrange their head groups to point out in a curved boundary. The extent to which this occurs is currently unknown and there is some evidence that both hydrophobic (tails straight) and hydrophilic (heads curved around) pores can coexist.