Lipid bilayer mechanics

Since lipid bilayers are essentially a two dimensional structure, Ka 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 Ka 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 Ka 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/m2. To calculate Ka 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 Ka for a lipid bilayer should be 80-200 mN/m (note: N/m is equivalent to J/m2). It is not surprising given this understanding of the forces involved that studies have shown that Ka 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. 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 K_b= K_a t^2 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. == 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). == Edge energy ==

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. == 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. 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. == References ==