Lecture 8. Bio-Membranes
Total Page:16
File Type:pdf, Size:1020Kb
Lecture 8. Bio-membranes Zhanchun Tu ( 涂展春 ) Department of Physics, BNU Email: [email protected] Homepage: www.tuzc.org Main contents ● Introduction ● Mathematical and physical preliminary ● Lipid membrane ● Cell membrane ● Summary and perspectives §8.1 Introduction Size and morphology of cells ● Size: several to tens of μm ● Various shapes (a) 5 cells of E. coli bacteria (b) 2 cells of yeast (c) Human red blood cell (d) Human white blood cell (e) Human sperm cell (f) Human epidermal (skin) cell (g) Human striated muscle cell (myofibril) (h) Human nerve cell Why various shapes? Animal cell The shapes of most of animal cells are determined by cytoskeleton Red blood cell No inner cellular organelles. Shape is determined by membrane. Human (normal): diameter 8μm, height 2 μm; biconcave discoid (why?) Cell membrane ●Timeline: cell-membrane bilayers awaiting a new model [Edidin (2003) Nature Reviews Molecular Cell Biology] ● Lipid structures micelle 胶囊 bilayer vesicle lipid molecule Liquid crystal phase. Cannot endure shear strain! hexagonal phase ● Cell membrane is usually in liquid crystal phase Liquid crystal phase is a necessary condition for cell as an open system Solid shell ===> cell is dead Isotropic fluid ==> no difference between inner and outside of cells in equilibrium ==> cell cannot exist as an basic unit for life 取自《从肥皂泡 到液晶生物膜》 Cancer might be related to the transition from LC to isotropic fluid Some problems we may deal with ● How to describe shapes mathematically? ● Does there exist a universal equation to govern the shapes? ● Why is human normal RBC a biconcave discoid? ● What is the mechanical function of membrane skeleton? ● To what extend membrane proteins will influence the shapes of membranes? §8.2 Mathematical & physical preliminary Curvature and torsion of a curve ● Curvature and torsion Each 3 points determine a curvature circle 1 2= r planeO123∧planeO' 234 2= 1,2,3,4 close to each other s23 ● Tangent, normal, and binormal vector t : tangent vector n: normal vector , point to O b : binormal vector , b⊥t ,b⊥ n t,n,b right-handed ● Two examples z b + t R b n = 2 2 R h/2 2R n t h h/2 =− R2h/22 1 = , =0 2R R y ● Frenet formula x s: arc length parameter Each s corresponds a point in the curve t˙ 0 0 t Different t, n, b at different points n˙ = − 0 n Different κ and τ at different points [b˙ ] [ 0 − 0][b] A curve in a surface ● Geodesic curvature & normal curvatures κ : Curvature of C' at P (roughly) g g=n⋅n ' (exactly) κ : Curvature of C'' at P (roughly) n n= n⋅N (exactly) 2 2 2 Obviously, gn= t, n: tangent and normal vectors of C N: normal vector of surface ● Geodesic torsion n': normal vector of C', such that {t,n',N} right-handed g=−N˙ ⋅n ' ● Two examples R h/2 1 = , =− = , =0 R2h/22 R2h/22 Rsin R cot 1 g=0 , n=− 2 2 = , =− R h/2 g R n R h/2 g= g=0 R2 h/22 Curvatures of a surface ● Principal curvatures l a Rotate 2 normal plane, curvature radii of 2 m r curves varies. o n 1 1 c1=− ,c2=− min{R1 } max {R2 } ● Mean and gaussian curvatures c c H = 1 2 ,K=c c 2 1 2 ● Two examples R R 1 1 c1=c2=− c1=− ,c2=0 R R c c 1 c c 1 1 H = 1 2 =− ,K=c c =0 H = 1 2 =− ,K=c c = 2 2 R 1 2 2 R 1 2 R2 Topological invariant of closed surface ● Gauss-Bonnet formula ∬KdA=41− g g=0⇒∬ KdA=4 g=1⇒∬ KdA=0 g=2⇒∬KdA=−4 Free energy ● Min F <=> equilibrium shapes F Configuration space cylinder torus ? sphere Finding min F <=> Solving δF=0 Stable: δ2F>0; unstable:δ2F<0 The meaning of variation Surface S Free energy F[S] Small deformation Surface S' Free energy F[S'] F=F [S ' ]−F [S ] δF=0 => Euler-Lagrange equation(s) describing equilibrium shapes Euler-Lagrange equation(s) <=> force balance equation(s) Variational problems on shapes in history ● Fluid films Viewed as a surface in mathematics # Soap films ---- minimal surfaces, Plateau (1803) F =∫ dA F=0⇒ H =0 Rotation axis 取自《从肥皂泡到液晶生物膜》 # Soap bubble ---- sphere, Young (1805), Laplace (1806) F =∫ dA p∫ dV p F=0⇒ H = p/2 in p "An embedded surface with constant mean out curvature in E3 must be a spherical surface" ---Alexandrov (1950s) p= pout− pin 1 − p 1 − p Sphere = Cylinder = R 2 R 取自《从肥皂泡 到液晶生物膜》 ● Solid shells # Possion (1821) F =∫ H 2dA # Schadow (1922) F=0⇒∇2 H 2 H H 2−K =0 Laplace operator # Willmore (1982) problem of surfaces Finding surfaces satisfying the above equation. ● Lipid bilayer (almost in-plane incompressible) # Spontaneous curvature energy, Helfrich (1973) k g= c 2 H c 2−kK 2 0 spontaneous curvature Analogy # Shape equation of vesicles, Ou-Yang & Helfrich (1987) F =∫ gdA∫ dA p∫ dV 2 2 F=0⇒ p−2 H 2 kc ∇ H kc 2 H c0 2H −c0 H −2K =0 k c=0⇒ p−2H =0 Young−Laplaceequation 2 2 p=0,=0,c0=0⇒∇ H 2 H H −K =0 Willmore surfaces Puzzle from the shape of RBC ● Sandwich model (Fung & Tong, 1968) To obtain the shape of biconcave discoid, Solid shell they should assume the thickness of the Fluid film membrane is nonuniform in μm scale. Solid shell Pinder's experiment (1972): the nonuniform thickness exists only in molecular (nm) scale. The thickness is uniform in large scale of μm. ● Nonuniform charge model (Lopez, 1968) Nonuniform charge distribution results in the shape of biconcave discoid. Experiment by Greet & Baker (1970): NO ● Incompressible shell model (Canham, 1970) Given the area and volume of membrane, the biconcave discoid minimizes the curvature energy ∫ H 2 dA Helfrich & Deuling (1975): the dumbbell-like shape can have the same curvature energy as the biconcave disk. But dumbbell-like shape has never observed in the experiment. Dumbbell-like Biconcave discoid 取自《从肥皂泡到液晶生物膜》 ● Spontaneous curvature model (Helfrich, 1973) Given the area and volume of membrane, the biconcave discoid minimizes the spontaneous curvature energy 2 ∫2 H c0 dA c00⇒ biconcave discoid is energetically favorable experiment Axisymmetric numerical result Can we give a analytic result from the shape equation? 2 2 p−2 H 2k c ∇ H k c2H c02 H −c0 H −2K =0 §8.3 Lipid membranes ----Soft-incompressible fluid film can endure bending but not static shear. May be asymmetric in inner side and outer one. Lipid vesicles ● Spherical vesicles 2 2 p−2 H 2k c ∇ H k c2H c02 H −c0 H −2K =0 2 2 f R≡ pR 2kc c0R−2k c c0=0 2R 1 1 f(R) f(R) H =− ,K= 2 R R p0,c00 p0,c00 f(R) p0,c00 R R R Might be related to endocytosis ● Torus [Ou-Yang (1990) PRA] =r/ {rcoscos ,rcossin ,sin } r2 cos cos 2 H =− ;K= rcos rcos 2 k 2 k c2 2 1 2 p 2 c / c 0 − r 2 2 2 3 = / 4 k c −4 k c 8 6 p c 0 c 0 cos 2 2 2 3 5 k c c0 −8 k c c0 10 6 p 2 cos 10μm 2 2 2 2 3 2 k c c0 −4 kc c0 4 2 p 3 cos =0 3 The coefficients of {1, cosj, cos2j, cos3j} should vanish! =2 2k c c c 0 2 0 [Mutz-Bensimon (1991) PRA] p=− ,=kc c0 − 2 2 ● Axisymmetric surface and Biconcave discoid 2 2 p−2 H 2k c ∇ H k c2H c02 H −c0 H −2K =0 Axisymmetric [Hu & Ou-Yang (1993) PRE] For−ec0 B0 [Naito-Okuda-OY (1993) PRE] [Evans-Fung (1972) Microvasc. Res.] describe a biconcave outline [Naito-Okuda-OY (1993) PRE] Lipid membranes with free edges ● Experiment: Opening process of lipid vesicles by Talin [Saitoh et al. (1998) PNAS] 5μm Bar: 2μm ● Previous theories – Derived from axisymmetric variation ● [Jülicher -Lipowsky (1993) PRL] ● [Zhou (2002) PhD thesis] Confine the variational problem in a subspace. Unreasonable results exist. – General case ● [Capovilla-Guven-Santiago (2002) PRE] Governing equations of edges are not expressed in the explicit forms of curvature and torsion. ● [Tu-OuYang (2003) PRE] Overcome the above shortages. ● Main results in [Tu-OuYang (2003)PRE] Free energy per area k G= c 2 H c 2−kK 2 0 Total free energy F =∫GdA∮ ds F =0⇒shape equationboundary conditions Shape equation: force balance in the normal direction 2 2 kc 2 H c02 H −c0 H −2 K −2 H 2k c ∇ H =0 Boundary conditions (curve C satisfies...) kc 2 H c0−kk n=0 Force balance equation of points in the edge along normal direction ∂ H , 2 k k −k =0 c ∂ b n g Moment balance equation of points in the edge around t Gk g=0 Force balance equation of points in the edge along b Axisymmetric analytic solutions Center of torus Cup-like open membrane Axisymmetric numerical solutions Solid squares: experiment data [Saitoh etal. (1998) PNAS] §8.4 Cell membrane ----It is beyond the lipid bilayer. How can we model it? Composite membrane model [Sackmann (2002)] ● Cell membrane = bilayer + membrane skeleton shear bending Bilayer NO YES MSK YES NO CM YES YES ● Basic assumptions – 1. CM: “smooth” surface – 2. Polymer in MSK: almost same chain length – 3. CM: in-plane isotropic, i.e. lipid crystal phase – 4. Chain length << curvature radius of CM – 5. Small deformations – 6. Free energy per area: analytical function – 7. Invariance of free energy: Strain tensor (+) (-) ● Energy density (energy per area) Assumption 1-4 => G=G2 H , K , 2 J ,Q 2 J =Tr 11 12 ,Q=det 11 12 12 22 12 22 Assumption 5-7 => up to the second order terms k G= c 2 H c 2−kK k /22 J 2−kQ 2 0 d Contribution from LB Contribution from MSK Bending energy Compress and shear energy Shape equation & in-plane strain equation F =∫GdA p∫ dV F =0 (Remark: only consider closed cell membrane) 2 2 p2k c [2 H c02 H −c0 H −2 K 2∇ H ]−2 H 2H k −k d 2 J −k ℜ: ∇ u=0 [Tu-OuYang, J.