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 Ka, a bending modulus Kb and an edge energy
Λ
See also: Elastic Modulus. 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.[4] 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.[5] 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[6] but only weakly with tail length and unsaturation.[7]
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.[6] A more direct method of measuring Ka is the pipette aspiration method, in which a single giant unilamellar vesicle (GUV) is held and stretched with a micropipette.[8] More recently, atomic force microscopy (AFM) has been used to probe the mechanical properties of suspended bilayer membranes,[9] 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.[7] 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.
See also: Beam (structure). 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
Kb=Kat2
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).[5] 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. [11]
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.
Λ
FE Modeling is a powerful tool for testing the mechanical deformation and equilibrium configuration of lipid membranes.[13] 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.[14] [15]
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.
G
| nel | |
| =\sum | |
| e=1 |
| e | |
| \left(G | |
| int |
| e | |
| -G | |
| ext |
\right)
For each discrete element the virtual work is determined by the force vector
f
x
\sigma
M
| e | |
| G | |
| int |
=\delta
| T | |
| x | |
| e |
| e | |
| \left(f | |
| int\sigma |
| e | |
| +f | |
| intM |
\right)
The FE force vectors due to the applied bilayer stress
\sigma\alpha
| e | |
| f | |
| int\sigma |
=\int\sigma\alpha
| T | |
| N | |
| ,\alpha |
\boldsymbol{a}\betada
Here
| T | |
| N | |
| \alpha |
\alpha
\boldsymbol{a}\beta
\beta
The above individual element vectors for the internal force
| e | |
| f | |
| \sigma |
Ge
\deltaxTf(x,q)+\deltaqTg(x)=0
In the above equation
x
q
f=0
g=0
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.
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.
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.