Biology: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 moduli: the area expansion modulus K_a, a bending modulus K_b and an edge energy [math]\displaystyle{ \Lambda }[/math]. 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: Biology:Elasticity of cell membranes

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

Bending modulus is defined as the energy required to deform a membrane from its natural 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 [math]\displaystyle{ K_b= K_a t^2 }[/math] 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. [math]\displaystyle{ \Lambda }[/math] 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.

Finite element modeling

FE Modeling is a powerful tool for testing the mechanical deformation and equilibrium configuration of lipid membranes. In this context membranes are treated under the thin-shell theory where the bending behavior of the membrane is described by the Helfrich bending model which considers the bilayer as being a very thin object and interprets it as a two-dimensional surface. This consideration imply that Kirchhoff-Love plate theory can be applied to lipid bilayers to determine their stress-deformation behavior. Furthermore, in the FE approach a bilayer surface is subdivided into discrete elements, each described by the above 2D mechanics.

Under these considerations the weak form virtual work for the entire system is described as the sum of the contribution of all the discrete elements work components.

[math]\displaystyle{ G =\sum_{e=1}^{n_{\mathrm{el}}}\left(G_{\mathrm{int}}^{e}-G_{\mathrm{ext}}^{e}\right) }[/math]

For each discrete element the virtual work is determined by the force vector [math]\displaystyle{ \mathbf{f} }[/math]and the displacement vector [math]\displaystyle{ \mathbf{x} }[/math] each for an applied stress [math]\displaystyle{ \sigma }[/math] and a bending momentum [math]\displaystyle{ M }[/math]

[math]\displaystyle{ G_{\mathrm{int}}^{e}=\delta \mathbf{x}_{e}^{\mathrm{T}}\left(\mathbf{f}_{\mathrm{int} \sigma}^{e}+\mathbf{f}_{\mathrm{int} M}^{e}\right) }[/math]

The FE force vectors due to the applied bilayer stress [math]\displaystyle{ \sigma^{\alpha \beta} }[/math] are given as

[math]\displaystyle{ \mathbf{f}_{\text {int } \sigma}^{e}=\int\sigma^{\alpha \beta} \mathbf{N}_{, \alpha}^{\mathrm{T}} \boldsymbol{a}_{\beta} \mathrm{d} a }[/math]

Here [math]\displaystyle{ \mathbf{N}_{ \alpha}^{\mathrm{T}} }[/math] is the displacement state function at point [math]\displaystyle{ \alpha }[/math], and [math]\displaystyle{ \boldsymbol{a}_{\beta} }[/math] the tangent vector to the bilayer surface at point [math]\displaystyle{ \beta }[/math]

The above individual element vectors for the internal force [math]\displaystyle{ \mathbf{f}_{\sigma}^{e} }[/math] and the internal work [math]\displaystyle{ \mathbf{G}^{e} }[/math]can be expressed in a global assembly to obtain a discretized weak form as follows:

[math]\displaystyle{ \delta \mathbf{x}^{\mathrm{T}} \mathbf{f}(\mathbf{x}, \mathbf{q})+\delta \mathbf{q}^{\mathrm{T}} \mathbf{g}(\mathbf{x})=0 }[/math]

In the above equation [math]\displaystyle{ \mathbf{x} }[/math] is the deformation in each discrete element while [math]\displaystyle{ \mathbf{q} }[/math] is the Lagrange multiplier associated with area-incompressibility. The discretized weak form is satisfied when [math]\displaystyle{ \mathbf{f}=0 }[/math] and [math]\displaystyle{ \mathbf{g}=0 }[/math]. The resulting nonlinear equations are solved using Newton ‘s method. This allows the predictions of equilibrium shapes that lipid membranes adopt under different stimuli.

Most analysis are done for lipid membranes with uniform properties (isotropic), while this is partly true for simple membranes containing a single or few lipid species, this description can only approximate the mechanical response of more complex lipid bilayers which can contain several domains of segregated lipids having distinct material properties or intermembrane proteins as in the case of cellular membranes. Other complex cases requiring surface flow analysis, pH, and temperature dependency would require a more discrete model of a bilayer such as MD simulations. FE methods can predict the equilibrium conformations of a lipid bilayer in response to external forces as shown in the following cases.

Tethering

In this scenario a point in the bilayer surface is pulled with a force normal to the surface plane, this leads to the elongation of a thin tether of bilayer material. The rest of the bilayer surface is subject to a tension S reflecting the pulling force of the continuous bilayer. In such this case a finer mesh is applied near the pulling force area to have a more accurate prediction of the bilayer deformation. Tethering is an important mechanical event for cellular lipid bilayers, by this action membranes are able to mediate docking into substrates or components of the cytoskeleton.

Budding

Lipid bilayer budding is a commonplace phenomenon in living cells and relates to the transport of metabolites in the form of vesicles. During this process, a lipid bilayer is subject to internal hydrostatic stresses, in combination with strain restrictions along a bilayer surface, this can lead to elongation of areas of the lipid bilayer by elastic shear or viscous shear. This eventually leads to a deformation of a typical spherical bilayer into different budding shapes. Such shapes are not restricted to be symmetrical along their axes, but they can have different degrees of asymmetry. In FE analysis this results in budding equilibrium shapes like, elongate plates, tubular buds and symmetric budding.