On the hydrodynamics of bilayer membranes

https://doi.org/10.1016/j.ijnonlinmec.2015.02.006Get rights and content

Highlights

  • Derivation of membrane hydrodynamics is detailed for almost planar membranes.

  • Onsager׳s variational method is used including all dissipation sources.

  • Relaxation rates are discussed in the presence of membrane tension.

  • Legitimacy of neglecting some dissipation source is discussed.

Abstract

Using Onsager׳s variational principle, we derive in the Monge gauge the equations describing the dynamics of almost planar bilayer membranes. All dissipations sources are taken into account: intermonolayer friction, solvent viscosity and monolayer viscosity. We recover and extend the results of Seifert and Langer and we discuss in detail the effect of membrane tension on the relaxation rates. Above a threshold tension, the avoided crossing of the relaxation rates disappears and the long-time dynamics is controlled by intermonolayer friction at all scales. The flow within the solvent and the monolayers is calculated both for the relaxation of a sinusoidal membrane shape and a Gaussian bud. The two vortices localized on the membrane split into four vortices at very small scales due to the monolayers slip. We discuss, depending on the scale, the legitimacy to neglect some of the dissipation sources, proposing an approximation scheme that allows to retain only the intermonolayer friction at small scales.

Introduction

Bilayer membranes form the outer walls of biological cells and organelles [1]. They are fluid, two-dimensional (2D) structures consisting of two opposite lipid monolayers that may host proteins, carbohydrates and other biological material. Lipids molecules of opposite monolayers contact each other at the level of the extremity of their aliphatic chains while the polar heads of the lipids face the aqueous solvent that surrounds the membrane.

Bilayer membranes are complex systems from the point of view of hydrodynamics: they consist of a pair of contacting 2D fluids with soft out-of-plane elasticity [2] moving in a 3D solvent. The monolayers must be treated as compressible fluids. Indeed, the 2D density is unconstrained as the monolayers can adapt their thickness in order to keep the 3D density constant. Besides, any variation of the lipid density between the two monolayers causes a spontaneous curvature of the bilayer [3], [4]. In the early studies of membrane hydrodynamics the bilayer structure of the membrane was neglected [5], [6]. While this is a good approximations for tensionless membranes at length-scales much larger than microns, experiments and theoretical studies have shown that taking into account the bilayer structure is essential at lower length-scales [7], [8], [9], [10], [11]. This is chiefly due to the importance of the dissipation caused by intermonolayer friction.

The full equations of bilayer hydrodynamics, including curvature and density-difference elasticity, intermonolayer friction, monolayer 2D viscosity and solvent 3D viscosity, were first derived by Seifert and Langer for almost planar membranes [8]. The method employed was a careful balance of in-plane and out-of-plane stresses. Generalization to non-linear membrane deformations, within the approximations of neglecting some dissipation sources, was performed by DeSimone, Arroyo and Rahimi [12], [13], [14] using covariant elasticity and Onsager׳s variational principle [15], [16].

In this paper we apply Onsager׳s variational principle [16] to the detailed derivation of the dynamics of almost planar bilayer membranes, including all types of dissipations (2 Parametrization, 3 Conservation laws and constraints, 4 Dissipation, 5 Elastic energy, 6 Differentiation of the Rayleighian, 7 Dynamical equation, 8 Integration of the bulk dynamical equation, 9 Dynamical equations for the membrane in Fourier space). We recover the linear equations of Ref. [8] with very little modifications (Section 10). We solve in detail these equations for an arbitrary Fourier component and we calculate the associated hydrodynamic 3D flow around the membrane (Section 10). We discuss the relaxation rates below and above a threshold tension separating two distinct behaviors and we physically interpret each relaxation rate in all of the various regimes (Section 11). To study in detail the hydrodynamics of the membrane, we calculate at various length scales the relaxation of a sinusoidal membrane shape (Section 12) and a Gaussian bud (Section 13), each time describing the flow pattern. Finally, we discuss depending on the scale, the legitimacy to neglect some of the dissipation sources, proposing an approximation scheme that allows to retain only the intermonolayer friction at small scales (Section 14).

Section snippets

Parametrization

We describe the membrane shape by the height z=h(r,t), above a reference plane (x,y), of the surface separating the two monolayers, as shown in Fig. 1. We call r=(x,y) a generic point in this reference plane and R=(r,z) a point in 3D space. Because the membrane is assumed to undergo only small deformations, at first order in h the hydrodynamical equations can be derived by considering that both monolayers lie in the z=0 plane. The superscript + (resp. −) will be associated with the monolayer in

Conservation laws and constraints

There are several conservations laws and constraints that must be satisfied. Regarding the bulk solvent as an incompressible fluid, we requireαVα±=0,where α=/Rα. In the membrane, the lipid mass conservation equation is ṅ±+i(n±vi±)=0, where i=/ri and a dot indicates time derivative. It yields at first-order:ρ̇±+ivi±=0.In addition, we require a no-slip boundary condition at the interface between each monolayer and the solvent:Vi±|z=0=vi±,and a no-permeation condition of the solvent

Dissipation

As mentioned in the Introduction, there are three sources of energy dissipation. The associated dissipation functions, defined as one-half of the dissipated power [15], arePb±=B±d3RηDαβ±Dαβ±,Ps±=d2r(η2dij±dij±+λ22dii±djj±),Pi=d2rb2(v+v)2,where B+ (resp. B) is the volume defined by z>0 (resp. z<0), Dαβ±(r)=12(αVβ±+βVα±) is the rate-of-deformation tensor in the bulk solvent and dij±(rs)=12(ivj±+jvi±) is the rate-of-deformation tensor in the monolayer fluids. Pb corresponds to the

Elastic energy

In a bilayer membrane, the two monolayers are strongly held in contact by hydrophobic forces (Fig. 1), hence they are strongly coupled from the point of curvature elasticity. However, the level of interdigitation between the lipids being usually very small, the density fields are essentially uncoupled and the excess elastic energy of the membrane can be written asH=SdS[f+(ρ+,c)+f+(ρ,c)]σSp.Because the monolayers are fluid only the membrane curvature c and the excess densities ρ± appear in

Differentiation of the Rayleighian

In the Stokes approximation, i.e., neglecting all inertial effects, the dynamical equations describing the motion of the membrane within the bulk solvent can be obtained by extremalizing with respect to all the dynamical variables the total Rayleighian of the system [16], [12], [13]:R=Pb++Ps++Pi+Ps+Pb+Ḣ.Taking into account the constraints (5), (6), (7), (8) through the introduction of the Lagrange multiplier fields P±(R), ζ±(r), μi±(r), and γ±(r), we need to extremalize the functionalR=ϵ=±

Dynamical equation

Eliminating the membrane Lagrange fields, we obtain the dynamical equations for the membrane:κ˜4hσ2h+ke2(ρ+ρ)+P+2ηzVz+P+2ηzVz=0,η22vi±(η2+λ2)i·v±+b(vi±vi)+ki(ρ±±e2h)η(zVi±+iVz±)=0.Note that these equations are defined in z=0 and ∇ is the gradient operator in the (x,y) space. In addition, we have the Stokes equation for the bulkη2Vα±+αP±=0αVα±=0,and the remaining constraints, which we regroup here for the sake of clarityρ̇±+ivi±=0,Vi±|z=0=vi±,Vz±|z=0=ḣ.These

Integration of the bulk dynamical equation

The bulk Stokes equations can be dealt with by Fourier transforming along the (x,y) plane and solving in the z direction. We define the Fourier transform byX(r,t)=d2q(2π)2X(q,t)eiq·r,Y(r,z,t)=d2q(2π)2Y(q,z,t)eiq·r,for surface and bulk quantities, respectively. Note that we shall use the same letters for the functions and their Fourier transforms in order to avoid heavy notations. In terms of V(q,z) and p(q,z), the bulk Stokes equations (26), (27) readη(z2q2)Vk±+iqkP±=0,η(z2q2)Vz±+zP±=0,

Dynamical equations for the membrane in Fourier space

In Fourier space, the viscous bulk stresses in (24) read P+2ηzVz+P+2ηzVz=2ηq(A++A)=4ηqḣ(q), while those in (25) read η(zV±+Vz±)=η2i(B±+qA±)=2ηqv±(q) along q and η(zV±+Vz±)=η(qv±(q)+0)=ηqv±(q) in the direction perpendicular to q. Hence, Fourier transforming the membrane dynamical equations (24), (25) yields(σq2+κ˜q4)hkeq2(ρ+ρ)+4ηqḣ=02ηsq2v±+b(v±v)+ikq(ρ±eq2h)+2ηqv±=0η2q2v±+b(v±v)+ηqv±=0ρ̇±+iqv±=0whereηs=η2+λ2/2appears as a renormalized surface viscosity.

Final dynamical equations and methodology

The procedure to solve the dynamics of a bilayer membrane is the following. Consider at t=0 a membrane with shape h0(r) and lipid excess densities ρ0±(r). First compute ρ0(r)=ρ+(r)ρ(r) and ρ¯0(r)=ρ+(r)+ρ(r), then set ρ¯(t)=0 at all times, since, as shown above, the average density relaxes to zero in an infinitesimal time. Next, decompose h0(r) and ρ0(r) into Fourier modes and solve for each mode Eqs. (57), (58):4ηqḣ(q,t)=(σq2+κ˜q4)h(q,t)+keq2ρ(q,t),2(b+ηq+ηsq2)ρ̇(q,t)=2keq4h(q,t)kq2ρ(q,t),

Relaxation rates

Eqs. (62), (63) that give the time evolution of h(q,t) and ρ(q,t) can be expressed as(ḣ(q,t)ρ̇(q,t))=M(q)(h(q,t)ρ(q,t))with M(q) the dynamical matrix:M(q)=(σq+κ˜q34ηkeq4ηkeq4b+ηq+ηsq2kq22(b+ηq+ηsq2)).Its two eigenvalues γ1(q)>γ2(q) are the two relaxation rates of the bilayer. They exhibit the avoided crossing1 phenomenon (Fig. 2). Let us discuss the time evolution of a Fourier mode q with initial

Relaxation of a Fourier mode and hydrodynamic flow

Let us consider, at t=0, a single Fourier mode h0cos(qx) for the membrane deformation. We assume that this shape has been fixed for some time, so the lipid density is at equilibrium. Minimizing the elastic energy (15) with respect to ρ±(r) at fixed h(r) yields ρ±=e2h, hence the initial density difference is ρ0cos(qx) with ρ0=2eq2h0. Let us determine the time evolution of the membrane and the hydrodynamic flow it generates.

Relaxation of a Gaussian bud

Let us now consider, at t=0, a Gaussian bud of equation h(r,0)=h0e12r2/w2,where w is the width at 5% of the height. Again, we assume that this deformation has been fixed for some time, so the lipid density is at equilibrium, which implies ρ(r,0)=2e2h(r,0). The dynamics of the bud can easily be obtained in the following way. In Fourier components, the bud is described by(h(q,0)ρ(q,0))=πw212h0e(1/48)w2q2(12eq2).Calling P(q) the matrix of the eigenvectors of M(q) and defining D(t)=diag(eγ1t,e

Legitimacy of neglecting dissipation sources

Let us discuss (within the framework of almost planar membranes) whether or not it is legitimate to neglect the bulk viscosity η, the monolayer viscosity ηs or the intermonolayer friction coefficient b. Such approximations are encountered in the literature.

Summary

We have derived the hydrodynamical equations of bilayer membranes from Onsager׳s variational principle for almost planar membranes. Membrane tension was taken into account together with all the sources of dissipation, i.e., solvent viscosity, monolayer shear and dilatational viscosities, and intermonolayer friction. The dynamical equations we obtain confirm and generalize those of Seifert and Langer [8]. We showed that there are two regimes for the relaxation rates whether the tension σ is

Acknowledgments

I thank Marino Arroyo for valuable discussions, and I acknowledge financial support form the French Agence Nationale de la Recherche (Contract no. ANR-12-BS04-0023-MEMINT).

References (21)

  • E. Evans et al.

    Hidden dynamics in rapid changes of bilayer shape

    Chem. Phys. Lipids

    (1994)
  • W. Rawicz et al.

    Effect of chain length and unsaturation on elasticity of lipid bilayers

    Biophys. J.

    (2000)
  • H. Lodish

    Molecular Cell Biology

    (2012)
  • W. Helfrich

    Elastic properties of lipid bilayers: theory and possible experiments

    Z. NaturForsch. C

    (1973)
  • S. Svetina et al.

    Lipid bilayer elasticity and the bilayer couple interpretation of red cell shape transformation and lysis

    Stud. Biophys.

    (1985)
  • L. Miao et al.

    Budding transitions of fluid-bilayer vesiclesthe effect of area-difference elasticity

    Phys. Rev. E

    (1994)
  • L. Kramer

    Theory of light scattering from fluctuations of membranes and monolayers

    J. Chem. Phys.

    (1971)
  • F. Brochard et al.

    Frequency spectrum of flicker phenomenon in erythrocytes

    J. Phys. France

    (1975)
  • U. Seifert et al.

    Viscous modes of fluid bilayer membranes

    Europhys. Lett.

    (1993)
  • P. Sens

    Dynamics of nonequilibrium membrane bud formation

    Phys. Rev. Lett.

    (2004)
There are more references available in the full text version of this article.

Cited by (26)

  • A viscous active shell theory of the cell cortex

    2022, Journal of the Mechanics and Physics of Solids
    Citation Excerpt :

    Most of the recent works on dynamic cell shape changes have considered the mechanics of the cortex, modeled as an active and viscous surface equipped with tensions only (Yeung and Evans, 1989; Turlier et al., 2014; Mietke et al., 2019; Torres-Sánchez et al., 2019). Other works have focused on the lipid bilayer (or plasma membrane) mechanics (onto which the cortex is attached through specialized linker proteins), which is generally described as a 2-dimensional fluid layer equipped with an elastic bending rigidity (Schneider et al., 1984; Arroyo and DeSimone, 2009; Rangamani et al., 2013; Fournier, 2015; Sauer et al., 2017; Sahu et al., 2017). Lipid bilayers may indeed be regarded as a stack of two lipid monolayers of nanometric thickness that do not flow transversely, which justifies why viscous dissipation is inherently 2-dimensional in that case.

View all citing articles on Scopus
View full text