Bending Stiffness of Lipid Bilayer Calculator

The bending modulus kappa is the energy cost of curving a membrane, defined by the Helfrich energy density (kappa/2)(c1 + c2 - c0)^2. Fluid phosphatidylcholine bilayers sit around 0.8 to 1.0 times 10^-19 J, which is about 20 to 25 kBT at room temperature, and this calculator returns it in both zeptojoules and kBT.

The polymer brush model of Rawicz and colleagues gives kappa = K_A h^2 / 24, where K_A is the area compressibility modulus and h is the hydrocarbon core thickness of the bilayer. The factor 24 reflects two monolayers that are mechanically coupled but can still slide against each other at the midplane.

The prefactor beta encodes how the two monolayers share bending strain. Fully coupled monolayers acting like a single solid plate give beta = 12, the polymer brush case gives beta = 24, and freely sliding uncoupled monolayers give beta = 48. Switching from polymer brush to uncoupled therefore halves the computed kappa.

The area compressibility modulus K_A measures how hard it is to stretch a membrane's area, with tension = K_A times the fractional area strain. Fluid phosphatidylcholine bilayers measure roughly 230 to 270 mN/m, and cholesterol-rich mixtures can exceed 600 mN/m.

Divide kappa in joules by kBT, where kB is 1.380649 times 10^-23 J/K and T is absolute temperature. At 25 degrees Celsius kBT is 4.1164 times 10^-21 J, so a bilayer with kappa = 8.66 times 10^-20 J has a bending modulus of about 21 kBT.

Bending a flat membrane patch into a closed sphere costs 8 pi kappa of Helfrich bending energy, independent of the vesicle radius, because the curvature energy density falls exactly as fast as the surface area grows. For a typical fluid bilayer this is around 500 kBT, far above thermal energy, which is why vesicles do not form spontaneously without extra driving forces.

Enter the hydrocarbon core thickness, which is roughly 2.5 to 3.5 nm for common fluid phospholipids and is typically 1 nm smaller than the head-to-head distance. Using the full head-to-head thickness inflates kappa by roughly 50% because of the squared dependence.

Cholesterol thickens and condenses saturated and mono-unsaturated bilayers, raising both K_A and h, so kappa can rise by a factor of two to four. Its effect on highly poly-unsaturated bilayers such as those with docosahexaenoic acid chains is much weaker because the disordered chains resist condensation.

The main methods are micropipette aspiration of giant unilamellar vesicles, thermal shape-fluctuation (flicker) spectroscopy, X-ray diffuse scattering on stacked bilayers, neutron spin echo, and pulling membrane tethers with optical tweezers. Values from different methods can disagree by a factor of two, which is one reason the coupling model matters.

Bending a plate stretches the outer leaflet and compresses the inner one, and the strain at a given curvature grows linearly with distance from the neutral midplane. Integrating that strain energy across the thickness produces a term proportional to h^2 times the stretching modulus, exactly the K_A h^2 form used here.

The bending modulus of the red blood cell membrane, including its spectrin skeleton, is usually quoted at roughly 2 to 7 times 10^-19 J, which is several times stiffer than a bare fluid phospholipid bilayer because the cytoskeleton contributes extra rigidity that this lipid-only formula does not capture.

The bending modulus in joules changes only weakly with temperature over a narrow fluid-phase range, but its value in kBT units falls as temperature rises simply because kBT grows. Near a gel-to-fluid phase transition kappa itself drops sharply, and the simple K_A h^2 form should not be applied across that transition.