Module 1 · The Universal Boundary
Membrane Biophysics
Every compartment considered in this course, except the membraneless ones of Module 6, is bounded by a phospholipid bilayer. Before discussing what lies inside each compartment, we must understand the surface that encloses it. The bilayer is a two-dimensional fluid — lipids diffuse within each leaflet with coefficients of order D ~ 10−12 m²/s — but it resists bending, stretching, and rupture with sharp and measurable mechanical constants.
Featured Lecture — Tom Rapoport (2/3)
Rapoport’s second iBiology lecture, How are cellular organelles shaped?, is the definitive introduction to this module’s physics. He treats ER tubule curvature (reticulons/DP1), three-way junctions, cisternal stacks, and the protein machinery that enforces high-curvature geometries against the bending cost derived below.
Cell Membrane Structure & Function
Phospholipid bilayer composition, membrane proteins, fluid mosaic model — an undergraduate primer before the Helfrich derivations.
1. The Canham–Helfrich Hamiltonian
The cornerstone of membrane mechanics is the elastic energy functional due independently to Canham (1970) and Helfrich (1973):
\[ \mathcal{H}_{\mathrm{mem}} = \int_{\mathcal{S}} dA \left[\tfrac{1}{2}\kappa(2H - C_0)^2 + \bar\kappa K + \sigma\right] \]
where H is the mean curvature, K the Gaussian curvature, C0the spontaneous curvature set by lipid composition and leaflet asymmetry, σ the surface tension, and κ ~ 10–25 kBT the bending rigidity. The Gaussian term κ̄K, via Gauss–Bonnet, is a topological invariant and only contributes when the membrane changes genus — which is precisely what happens during vesicle budding, mitochondrial fission, and nuclear envelope reassembly.
Every membrane remodelling event in the cell is a topological event. You cannot ignore the Gaussian term if you want to count vesicles: it is what distinguishes "pinching" from "pushing."
2. Thermal Fluctuation Spectrum
Fourier-decomposing small fluctuations h(r) around a flat reference membrane and integrating out the normal modes gives the thermal fluctuation spectrum
\[ \langle |h_{\mathbf{q}}|^2 \rangle = \dfrac{k_B T}{\kappa q^4 + \sigma q^2} \]
a result first derived by Helfrich in 1973 and still the most-used equation in membrane biophysics. At high q (short wavelengths) bending dominates and ⟨|hq|²⟩ ~ q−4. At long wavelengths tension dominates and the spectrum crosses over to q−2. The crossover scale q* = √(σ/κ) carries the effective persistence of a free membrane.
At q*, the bending penalty and tension penalty per unit area are equal. Below this scale, the membrane fluctuates freely and bends; above it, tension pins it flat. Flicker-spectroscopy measurements on red blood cells exploit precisely this dependence to extract κ (typically ~25 kBT for erythrocytes) from the cell’s passive shimmer.
Simulation: Helfrich Fluctuation Spectrum
Click Run to execute the Python code
Code will be executed with Python 3 on the server
3. Three Curvature Regimes
Flat (H = 0)
Large membrane domains — the bulk of the plasma membrane, ER cisternal faces, the inner-boundary membrane of mitochondria. Minimises bending cost. Stabilised by lipid rafts and cortical cytoskeleton.
Positively curved (H > 0)
Buds, transport vesicles, ER tubules, the tips of filopodia. Require curvature-sensing/generating proteins: BAR domains (sensing and generating in 10–70 nm range), ENTH domains (epsin, amphipathic helix insertion), dynamin (constricting necks). A 30 nm ER tubule represents ~15 kBT per lipid in bending energy — significant, and paid for by active protein tethering.
Saddle (K < 0)
Mitochondrial cristae junctions, fission necks, the pore-like contortions of nuclear envelope remodelling. Saddle geometry has H ~ 0 (low mean-curvature cost) but negative Gaussian curvature K < 0, contributing to the topological term. Stabilised by cone-shaped lipids with negative spontaneous curvature: cardiolipin at mitochondrial cristae junctions, PA during fission.
4. Surface-to-Volume as a Design Constraint
Click Run to execute the Python code
Code will be executed with Python 3 on the server
5. Why the Bending Rigidity Matters
A membrane with κ = 20 kBT has a persistence length of many micrometres — the scale of the cell itself. This is why a free lipid bilayer at physiological tension looks flat on cellular length scales. Every curved structure you will meet in this course — the tubular ER, the cristae, the lysosomal limiting membrane around a 250 nm organelle — is maintained against this rigidity by active machinery.
Geometry, in the cell, is always paid for.The bending energy of a sphere of radius R is exactly 8πκ (scale-invariant), so a 100 nm transport vesicle and a 10 μm plasma membrane cost the same — but the vesicle is unstable without dynamin pinching its neck, while the plasma membrane is stabilised by an actin cortex. Scale invariance means the cell cannot win through size; it must win through active enforcement.
6. Asymmetry & Spontaneous Curvature
Lipid composition is not symmetric across a bilayer. The outer leaflet of the plasma membrane is enriched in phosphatidylcholine (PC) and sphingomyelin (SM); the inner leaflet in phosphatidylethanolamine (PE), phosphatidylserine (PS), and phosphoinositides. This asymmetry is maintained by ATP-driven flippases (P4-ATPases) and is biologically essential — exposure of PS to the outer leaflet is the “eat me” signal for phagocytosis of apoptotic cells.
Cone- and inverted-cone shaped lipids contribute to spontaneous curvature C0. PE and cardiolipin (small headgroup relative to tails) drive negative curvature; lysolipids and PIP2 (large headgroup) drive positive curvature. By partitioning curvature-generating lipids across leaflets, the cell can build a membrane that naturally wants to bend in a preferred direction — a prerequisite for tubules and cristae.