Bilayer elasticity at the nanoscale: the need for new terms

Probing membrane deformation energy by KcsA potassium channel gating under varied membrane thickness and tension

This study investigated how membrane thickness and tension modify the gating of KcsA potassium channels when simultaneously varied. The KcsA channel undergoes global conformational changes upon gating: expansion of the cross-sectional area and longitudinal shortening upon opening. Thus, membranes impose differential effects on the open and closed conformations, such as hydrophobic mismatches. Here, the single-channel open probability was recorded in the contact bubble bilayer, by which variable thickness membranes under a defined tension were applied. A fully open channel in thin membranes turned to sporadic openings in thick membranes, where the channel responded moderately to tension increase. Quantitative gating analysis prompted the hypothesis that tension augmented the membrane deformation energy when hydrophobic mismatch was enhanced in thick membranes.

Socioeconomic agents as active matter in nonequilibrium Sakoda-Schelling models

How robust are socioeconomic agent-based models with respect to the details of the agents' decision rule? We tackle this question by considering an occupation model in the spirit of the Sakoda-Schelling model, historically introduced to shed light on segregation dynamics among human groups. For a large class of utility functions and decision rules, we pinpoint the nonequilibrium nature of the agent dynamics, while recovering an equilibrium-like phase separation phenomenology. Within the mean-field approximation we show how the model can be mapped, to some extent, onto an active matter field description. Finally, we consider nonreciprocal interactions between two populations and show how they can lead to nonsteady macroscopic behavior. We believe our approach provides a unifying framework to further study geography-dependent agent-based models, notably paving the way for joint consideration of population and price dynamics within a field theoretic approach.

Dynamic structure factor of undulating vesicles: finite-size and spherical geometry effects with application to neutron spin echo experiments

We consider the dynamic structure factor (DSF) of quasi-spherical vesicles and present a generalization of an expression that was originally formulated by Zilman and Granek (ZG) for scattering from isotropically oriented quasi-flat membrane plaquettes. The expression is obtained in the form of a multi-dimensional integral over the undulating membrane surface. The new expression reduces to the original stretched exponential form in the limit of sufficiently large vesicles, i.e., in the micron range or larger. For much smaller unilamellar vesicles, deviations from the asymptotic, stretched exponential equation are noticeable even if one assumes that the Seifert-Langer leaflet density mode is completely relaxed and membrane viscosity is neglected. To avoid the need for an exhaustive numerical integration while fitting to neutron spin echo (NSE) data, we provide a useful approximation for polydisperse systems that tests well against the numerical integration of the complete expression. To validate the new expression, we performed NSE experiments on variable-size vesicles made of a POPC/POPS lipid mixture and demonstrate an advantage over the original stretched exponential form or other manipulations of the original ZG expression that have been deployed over the years to fit the NSE data. In particular, values of the membrane bending rigidity extracted from the NSE data using the new approximations were insensitive to the vesicle radii and scattering wavenumber and compared very well with expected values of the effective bending modulus ([Formula: see text]) calculated from results in the literature. Moreover, the generalized scattering theory presented here for an undulating quasi-spherical shell can be easily extended to other models for the membrane undulation dynamics beyond the Helfrich Hamiltonian and thereby provides the foundation for the study of the nanoscale dynamics in more complex and biologically relevant model membrane systems.

Bilayer tension-induced clustering of the UPR sensor IRE1

The endoplasmic reticulum acts as a protein quality control center where a range of chaperones and foldases facilitates protein folding. IRE1 is a sensory transmembrane protein that transduces signals of proteotoxic stress by forming clusters and activating a cellular program called the unfolded protein response (UPR). Recently, membrane thickness variation due to membrane compositional changes have been shown to drive IRE1 cluster formation, activating the UPR even in the absence of proteotoxic stress. Here, we demonstrate a direct relationship between bilayer tension and UPR activation based on IRE1 dimer stability. The stability of the IRE1 dimer in a (50%DOPC-50%POPC) membrane at different applied bilayer tensions was analyzed via molecular dynamics simulations. The potential of mean force for IRE1 dimerization predicts a higher concentration of IRE1 dimers for both tensed and compressed ER membranes. This study shows that IRE1 may be a mechanosensitive membrane protein and establishes a direct biophysical relationship between bilayer tension and UPR activation.

Demixing of homogeneous binary lipid membranes induced by protein inclusions

We study a model of a lipid bilayer membrane described by two order parameters: the chemical composition described using the Gaussian model and the spatial configuration described with the elastic deformation model of a membrane with a finite thickness or, equivalently, for an adherent membrane. We assume and explain on physical grounds the linear coupling between the two order parameters. Using the exact solution, we calculate the correlation functions and order parameter profiles. We also study the domains that form around inclusions on the membrane. We propose and compare six distinct ways to quantify the size of such domains. Despite its simplicity, the model has many interesting features like the Fisher-Widom line and two distinct critical regions.

Multiscale lipid membrane dynamics as revealed by neutron spectroscopy

The plasma membrane is one of the principal structural components of the cell and, therefore, one of the key components of the cellular life. Because the membrane's dynamics links the membrane's structure and function, the complexity and the broad range of the membrane's motions are essential for the enormously diverse functionality of the cell membrane. Even for the main membrane component, the lipid bilayer, considered alone, the range and complexity of the lipid motions are remarkable. Spanning the time scale from sub-picosecond to minutes and hours, the lipid motion in a bilayer is challenging to study even when a broad array of dynamic measurement techniques is employed. Neutron scattering plays a special role among such dynamic measurement techniques, particularly, because it involves the energy transfers commensurate with the typical intra- and inter- molecular dynamics and the momentum transfers commensurate with intra- and inter-molecular distances. Thus, using neutron scattering-based techniques, the spatial and temporal information on the lipid motion can be obtained and analysed simultaneously. Protium vs. deuterium sensitivity and non-destructive character of the neutron probe add to the remarkable prowess of neutron scattering for elucidating the lipid dynamics. Herein we present an overview of the neutron scattering-based studies of lipid dynamics in model membranes, with a discussion of the direct relevance and implications to the real-life cell membranes. The latter are much more complex systems than simple model membranes, consisting of heterogeneous non-stationary domains composed of lipids, proteins, and other small molecules, such as carbohydrates. Yet many fundamental aspects of the membrane behavior and membrane interactions with other molecules can be understood from neutron scattering measurements of the model membranes. For example, such studies can provide a great deal of information on the interactions of antimicrobial compounds with the lipid matrix of a pathogen membrane, or the interactions of drug molecules with the plasma membrane. Finally, we briefly discuss the recently emerging field of neutron scattering membrane studies with a reach far beyond the model membrane systems.

Field-mediated interactions of passive and conformation-active particles: multibody and retardation effects

Particles in soft matter interact through the deformation field they create, as in the "cheerios" effect or the curvature-mediated interactions of membrane proteins. Using a simple model for field-mediated interactions between passive particles, or active particles that switch conformation randomly or synchronously, we derive generic results concerning multibody interactions, activity driven patterns, and retardation effects.

Membrane-Mediated Interactions Between Protein Inclusions

Integral or peripheral membrane proteins, or protein oligomers often get close to each other on cell membranes and carry out biological tasks in a collective manner. In addition to electrostatic and van der Waals interactions, those proteins also experience membrane-mediated interactions, which may be necessary for their functionality. The membrane-mediated interactions originate from perturbation of lipid membranes by the presence of protein inclusions, and have been the subject of intensive research in membrane biophysics. Here we review both theoretical and numerical studies of such interactions for membrane proteins and for nanoparticles bound to lipid membranes.

Additional contributions to elastic energy of lipid membranes: Tilt-curvature coupling and curvature gradient

Lipid bilayer membranes under biologically relevant conditions are flexible thin laterally fluid films consisting of two unimolecular layers (monolayers) each about 2 nm thick. On spatial scales much larger than the bilayer thickness, the membrane elasticity is well determined by its shape. The classical Helfrich theory considers the membrane as an elastic two-dimensional (2D) film, which has no particular internal structure. However, various local membrane heterogeneities can result in a lipids tilt relative to the membrane surface normal. On the basis of the classical elasticity theory of 3D bodies, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)10.1007/s101890070003] derived the most general energy functional, taking into account the tilt and lipid monolayer curvature. Recently, Terzi and Deserno [J. Chem. Phys. 147, 084702 (2017)10.1063/1.4990404] showed that Hamm and Kozlov's derivation was incomplete because the tilt-curvature coupling term had been missed. However, the energy functional derived by Terzi and Deserno appeared to be unstable, thereby being invalid for applications that require minimizations of the overall energy of deformations. Here, we derive a stable elastic energy functional, showing that the squared gradient of the curvature was missed in both of these works. This change in the energy functional arises from a more accurate consideration of the transverse shear deformation terms and their influence on the membrane stability. We also consider the influence of the prestress terms on the stability of the energy functional, and we show that it should be considered small and the effective Gaussian curvature should be neglected because of the stability requirements. We further generalize the theory, including the stretching-compressing deformation modes, and we provide the geometrical interpretation of the terms that were previously missed by Hamm and Kozlov. The physical consequences of the new terms are analyzed in the case of a membrane-mediated interaction of two amphipathic peptides located in the same monolayer. We also provide the expression for director fluctuations, comparing it with that obtained by Terzi and Deserno.

Spatial organization of active particles with field-mediated interactions

We consider a system of independent pointlike particles performing a Brownian motion while interacting with a Gaussian fluctuating background. These particles are in addition endowed with a discrete two-state internal degree of freedom that is subjected to a nonequilibrium source of noise, which affects their coupling with the background field. We explore the phase diagram of the system and pinpoint the role of the nonequilibrium drive in producing a nontrivial patterned spatial organization. We are able, by means of a weakly nonlinear analysis, to account for the parameter-dependence of the boundaries of the phase and pattern diagram in the stationary state.

Directed Supramolecular Organization of N-BAR Proteins through Regulation of H0 Membrane Immersion Depth

Many membrane remodeling events rely on the ability of curvature-generating N-BAR membrane proteins to organize into distinctive supramolecular configurations. Experiments have revealed a conformational switch in N-BAR proteins resulting in vesicular or tubular membrane shapes, with shallow membrane immersion of the H0 amphipathic helices of N-BAR proteins on vesicles but deep H0 immersion on tubes. We develop here a minimal elastic model of the local thinning of the lipid bilayer resulting from H0 immersion. Our model predicts that the observed conformational switch in N-BAR proteins produces a corresponding switch in the bilayer-mediated N-BAR interactions due to the H0 helices. In agreement with experiments, we find that bilayer-mediated H0 interactions oppose N-BAR multimerization for the shallow H0 membrane immersion depths measured on vesicles, but promote self-assembly of supramolecular N-BAR chains for the increased H0 membrane immersion depths measured on tubes. Finally, we consider the possibility that bilayer-mediated H0 interactions might contribute to the concerted structural reorganization of N-BAR proteins suggested by experiments. Our results indicate that the membrane immersion depth of amphipathic protein helices may provide a general molecular control parameter for membrane organization.

Nonequilibrium dynamic structure factor of a lipid bilayer in the presence of an in-plane temperature gradient

There is rapidly increasing evidence that nanoscale temperature heterogeneities are involved in important biological processes. Combining nanoheating and nanoscale thermosensors forms the basis of emerging unique methods of cell therapy, tissue engineering, and regenerative medicine. Understanding corresponding phenomena seems to require a mesoscopic nonequilibrium hydrodynamic theory. In this paper, a Langevin-type model of dynamics of phonon modes propagating along a bilayer lipid membrane in the presence of an in-plane temperature gradient is proposed. Corresponding quantitative estimates for the Brillouin components of the nonequilibrium dynamic structure factor and the equal-time longitudinal momentum-density correlation function for a lipid bilayer are obtained. The analysis reveals that for typical values of parameters of lipid bilayer, the longitudinal temperature gradient of the order of 5qK for wave numbers q from 0.01 to 1nm^{-1} induces significant asymmetry of the Brillouin components of the dynamic structure factor and long-range spatial correlations in the plane of the bilayer. The corresponding membrane temperature gradients seem to be typical or achievable for cellular processes responsible for intracellular temperature variations and such external physical impacts as high-intensity electromagnetic pulses or heating of membrane-associated nanoparticles.

Field-Embedded Particles Driven by Active Flips

Systems of independent active particles embedded into a fluctuating environment are relevant to many areas of soft-matter science. We use a minimal model of noninteracting spin-carrying Brownian particles in a Gaussian field and show that activity-driven spin dynamics leads to patterned order. We find that the competition between mediated interactions and active noise alone can yield such diverse behaviors as phase transitions and microphase separation, from lamellar up to hexagonal ordering of clusters of opposite magnetization. These rest on complex multibody interactions. We find regimes of stationary patterns, but also dynamical regimes of relentless birth and growth of lumps of magnetization opposite of the surrounding one. Our approach combines Monte Carlo simulations with analytical methods based on dynamical density functional approaches.

Novel tilt-curvature coupling in lipid membranes

On mesoscopic scales, lipid membranes are well described by continuum theories whose main ingredients are the curvature of a membrane's reference surface and the tilt of its lipid constituents. In particular, Hamm and Kozlov [Eur. Phys. J. E 3, 323 (2000)] have shown how to systematically derive such a tilt-curvature Hamiltonian based on the elementary assumption of a thin fluid elastic sheet experiencing internal lateral pre-stress. Performing a dimensional reduction, they not only derive the basic form of the effective surface Hamiltonian but also express its emergent elastic couplings as trans-membrane moments of lower-level material parameters. In the present paper, we argue, though, that their derivation unfortunately missed a coupling term between curvature and tilt. This term arises because, as one moves along the membrane, the curvature-induced change of transverse distances contributes to the area strain-an effect that was believed to be small but nevertheless ends up contributing at the same (quadratic) order as all other terms in their Hamiltonian. We illustrate the consequences of this amendment by deriving the monolayer and bilayer Euler-Lagrange equations for the tilt, as well as the power spectra of shape, tilt, and director fluctuations. A particularly curious aspect of our new term is that its associated coupling constant is the second moment of the lipid monolayer's lateral stress profile-which within this framework is equal to the monolayer Gaussian curvature modulus, κ¯_m. On the one hand, this implies that many theoretical predictions now contain a parameter that is poorly known (because the Gauss-Bonnet theorem limits access to the integrated Gaussian curvature); on the other hand, the appearance of κ¯_m outside of its Gaussian curvature provenance opens opportunities for measuring it by more conventional means, for instance by monitoring a membrane's undulation spectrum at short scales.

The effect of gramicidin inclusions on the local order of membrane components

We study the local effect of the antimicrobial peptide Gramicidin A on bilayers composed of lipids or surfactants using nuclear magnetic resonance spectroscopy and wide-angle X-ray scattering, techniques that probe the orientational and positional order of the alkyl chains, respectively. The two types of order vary with temperature and peptide concentration in complex ways which depend on the membrane composition, highlighting the subtlety of the interaction between inclusions and the host bilayer. The amplitude of the variation is relatively low, indicating that the macroscopic constants used to describe the elasticity of the bilayer are unlikely to change with the addition of peptide.

Coupling between Inclusions and Membranes at the Nanoscale

The activity of cell membrane inclusions (such as ion channels) is influenced by the host lipid membrane, to which they are elastically coupled. This coupling concerns the hydrophobic thickness of the bilayer (imposed by the length of the channel, as per the hydrophobic matching principle) but also its slope at the boundary of the inclusion. However, this parameter has never been measured so far. We combine small-angle x-ray scattering data and a complete elastic model to measure the slope for the model gramicidin channel and show that it is surprisingly steep in two membrane systems with very different elastic properties. This conclusion is confirmed and generalized by the comparison with recent results in the simulation literature and with conductivity measurements.

Probing Elastic and Viscous Properties of Phospholipid Bilayers Using Neutron Spin Echo Spectroscopy

The elastic and viscous properties of self-assembled amphiphilic membranes dictate the intricate hierarchy of their structure and dynamics ranging from the diffusion of individual molecules to the large-scale deformation of the membrane. We previously demonstrated that neutron spin echo spectroscopy measurements of model amphiphilic membranes can access the naturally occurring submicrosecond membrane motions, such as bending and thickness fluctuations. Here we show how the experimentally measured fluctuation parameters can be used to determine the inherent membrane properties and demonstrate how membrane viscosity and compressibility modulus are influenced by lipid composition in a series of simple phosphatidylcholine bilayers with different tail lengths as a function of temperature. This approach highlights the interdependence of the bilayer elastic and viscous properties and the collective membrane dynamics and opens new avenues to investigating the mechanical properties of more complex and biologically inspired systems.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Bilayer-thickness-mediated interactions between integral membrane proteins

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for noncylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. We provide error estimates of our numerical solution schemes and systematically assess their convergence properties. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.

Hydrodynamics of bilayer membranes with diffusing transmembrane proteins

We consider the hydrodynamics of lipid bilayers containing transmembrane proteins of arbitrary shape. This biologically-motivated problem is relevant to the cell membrane, whose fluctuating dynamics play a key role in phenomena ranging from cell migration, intercellular transport, and cell communication. Using Onsager's variational principle, we derive the equations that govern the relaxation dynamics of the membrane shape, of the mass densities of the bilayer leaflets, and of the diffusing proteins' concentration. With our generic formalism, we obtain several results on membrane dynamics. We find that proteins that span the bilayer increase the intermonolayer friction coefficient. The renormalization, which can be significant, is in inverse proportion to the protein's mobility. Second, we find that asymmetric proteins couple to the membrane curvature and to the difference in monolayer densities. For practically all accessible membrane tensions (σ > 10(-8) N m(-1)) we show that the protein density is the slowest relaxing variable. Furthermore, its relaxation rate decreases at small wavelengths due to the coupling to curvature. We apply our formalism to the large-scale diffusion of a concentrated protein patch. We find that the diffusion profile is not self-similar, owing to the wavevector dependence of the effective diffusion coefficient.

Line tension of multicomponent bilayer membranes

The line tension or edge energy of bilayer membranes self-assembled from binary amphiphilic molecules is studied using self-consistent-field theory (SCFT). Specifically, solutions of the SCFT equations corresponding to an infinite membrane with a circular pore, or an open membrane, are obtained for a coarse-grained model in which the amphiphilic species and hydrophilic solvents are represented by ABandED diblock copolymers and C homopolymers, respectively. The edge energy of the membrane is extracted from the free energy of the open membranes. Results for membranes composed of mixtures of symmetric and cone- or inverse cone-shaped amphiphilic molecules with neutral and/or repulsive interactions are obtained and analyzed. It is observed that an increase in the concentration of the cone-shaped species leads to a decrease of the line tension. In contrast, adding inverse cone-shaped copolymers results in an increase of the line tension. Furthermore, the density profile of the copolymers reveals that the line tension is regulated by the distribution of the amphiphiles at the bilayer edge.

The role of membrane-mediated interactions in the assembly and architecture of chemoreceptor lattices

In vivo fluorescence microscopy and electron cryo-tomography have revealed that chemoreceptors self-assemble into extended honeycomb lattices of chemoreceptor trimers with a well-defined relative orientation of trimers. The signaling response of the observed chemoreceptor lattices is remarkable for its extreme sensitivity, which relies crucially on cooperative interactions among chemoreceptor trimers. In common with other membrane proteins, chemoreceptor trimers are expected to deform the surrounding lipid bilayer, inducing membrane-mediated anisotropic interactions between neighboring trimers. Here we introduce a biophysical model of bilayer-chemoreceptor interactions, which allows us to quantify the role of membrane-mediated interactions in the assembly and architecture of chemoreceptor lattices. We find that, even in the absence of direct protein-protein interactions, membrane-mediated interactions can yield assembly of chemoreceptor lattices at very dilute trimer concentrations. The model correctly predicts the observed honeycomb architecture of chemoreceptor lattices as well as the observed relative orientation of chemoreceptor trimers, suggests a series of “gateway” states for chemoreceptor lattice assembly, and provides a simple mechanism for the localization of large chemoreceptor lattices to the cell poles. Our model of bilayer-chemoreceptor interactions also helps to explain the observed dependence of chemotactic signaling on lipid bilayer properties. Finally, we consider the possibility that membrane-mediated interactions might contribute to cooperativity among neighboring chemoreceptor trimers.

The chemotaxis system allows bacteria to respond to minute changes in chemical concentration, and serves as a paradigm for biological signal processing and the self-assembly of large protein lattices in living cells. The sensitivity of the chemotaxis system relies crucially on cooperative interactions among chemoreceptor trimers, which are organized into intricate honeycomb lattices. Chemoreceptors are membrane proteins and, hence, are expected to deform the surrounding lipid bilayer, leading to membrane-mediated interactions between chemoreceptor trimers. Using a biophysical model of bilayer-chemoreceptor interactions we show that the membrane-mediated interactions induced by chemoreceptor trimers provide a mechanism for the observed self-assembly of chemoreceptor lattices. We find that the directionality of membrane-mediated interactions between trimers complements protein-protein interactions in the stabilization of the observed honeycomb architecture of chemoreceptor lattices. Our results suggest that the symmetry of membrane protein complexes such as chemoreceptor trimers is reflected in the anisotropy of membrane-mediated interactions, yielding a general mechanism for the self-assembly of ordered protein lattices in cell membranes.

Small scale membrane mechanics

Large scale changes to lipid bilayer shapes are well represented by the Helfrich model. However, there are membrane processes that take place at smaller length scales that this model cannot address. In this work, we present a one-dimensional continuum model that captures the mechanics of the lipid bilayer membrane at the length scale of the lipids themselves. The model is developed using the Cosserat theory of surfaces with lipid orientation, or 'tilt', as the fundamental degree of freedom. The Helfrich model can be recovered as a special case when the curvatures are small and the lipid tilt is everywhere zero. We use the tilt model to study local membrane deformations in response to a protein inclusion. Parameter estimates and boundary conditions are obtained from a coarse-grained molecular model using dissipative particle dynamics (DPD) to capture the same phenomenon. The continuum model is able to reproduce the membrane bending, stretch and lipid tilt as seen in the DPD model. The lipid tilt angle relaxes to the bulk tilt angle within 5-6 nm from the protein inclusion. Importantly, for large tilt gradients induced by the proteins, the tilt energy contribution is larger than the bending energy contribution. Thus, the continuum model of tilt accurately captures behaviors at length scales shorter than the membrane thickness.