Membrane-mediated interactions

TL;DR

This paper reviews how inclusions like proteins influence membrane behavior through shape and fluctuation constraints, highlighting generic, non-chemical interactions in biological membranes.

Contribution

It provides a comprehensive overview of membrane-mediated interactions arising from shape and thickness constraints, emphasizing generic physical mechanisms.

Findings

Long-range interactions from curvature and fluctuation constraints

Short-range interactions from hydrophobic mismatch

Use of coarse-grained models for analysis

Abstract

Interactions mediated by the cell membrane between inclusions, such as membrane proteins or antimicrobial peptides, play important roles in their biological activity. They also constitute a fascinating challenge for physicists, since they test the boundaries of our understanding of self-assembled lipid membranes, which are remarkable examples of two-dimensional complex fluids. Inclusions can couple to various degrees of freedom of the membrane, resulting in different types of interactions. In this chapter, we review the membrane-mediated interactions that arise from direct constraints imposed by inclusions on the shape of the membrane. These effects are generic and do not depend on specific chemical interactions. Hence, they can be studied using coarse-grained soft matter descriptions. We deal with long-range membrane-mediated interactions due to the constraints imposed by inclusions on membrane curvature and on its fluctuations. We also discuss the shorter-range interactions that arise from the constraints on membrane thickness imposed by inclusions presenting a hydrophobic mismatch with the membrane.

Figures

Tables

Table 1. Table 1: Summary of the power laws obtained for the leading-order terms of the two types of membrane-mediated interactions, as a function of the separation d 𝑑 d between the inclusions, in the regime of small deformations of a flat membrane and distant inclusions. Different inclusion geometries are considered. In the case labeled “disks+anisotropy”, the anisotropy can be either in the inclusion shape (e.g. ellipsoidal  Park96 ) or in the constraint it imposes (e.g. an anisotropic local curvature  Dommersnes99 ).

Dominant term in		Fluctuation-induced	Interaction due to the
the Hamiltonian	Geometry	interaction	ground-state deformation –
in Eq. 7			Vanishes if up-down symmetric
Bending rigidity $\kappa$	Disks	$1/d^4$ Goulian93etc ; Park96 ; Dommersnes99	$1/d^4$ Goulian93etc ; Park96 ; Dommersnes99
	Disks	$1/d^4$ Park96 ; Dommersnes99	$1/d^2$ Park96 ; Dommersnes99
	+anisotropy
	Distant rods	$1/d^4$ Golestanian96a ; Golestanian96b
Tension $\sigma$	Disks	$1/d^8$ Noruzifar09 ; Lin11	Exponentially suppressed
	Disks	$1/d^8$ Noruzifar09 ; Lin11	$1/d^4$ Stamou00 ; Yolcu12b
	+anisotropy
	Distant rods	$1/d^4$ Golestanian96a ; Golestanian96b

Equations

$$H=\int dA\left[\frac{\kappa }{2}{\left(c_1+c_2\right)}^2+\bar{\kappa }{c}_1{c}_2\right],$$

$$H[h]=\int d\mathbf{r}\left\{\frac{\kappa }{2}{\left[{\nabla }^{2}h(\mathbf{r})\right]}^2+\bar{\kappa }\det [{\partial }_i{\partial }_{j}h(\mathbf{r})]\right\},$$

$$U_1(d)=4\pi \kappa ({\alpha }_1^2+{\alpha }_2^2)\frac{a^4}{d^4},$$

$$U_2(d)=-6k_{\mathrm{B}}T\frac{a^4}{d^4}.$$

$$|\alpha |=\sqrt{\frac{3}{4\pi }\frac{k_{\mathrm{B}}T}{\kappa }}.$$

$$U_3(d)=-8\pi \kappa \frac{a^4}{d^2}\left\{2J^2\cos (2({\theta }_1+{\theta }_2))+JK\left[\cos (2{\theta }_1)+\cos (2{\theta }_2)\right]\right\},$$

$$H[h]=\int d\mathbf{r}\left\{\frac{\kappa }{2}{\left[{\nabla }^{2}h(\mathbf{r})\right]}^2+\frac{\sigma }{2}{\left[\nabla h(\mathbf{r})\right]}^2\right\},$$

$$H[h]=\frac{\sigma }{2}\int d\mathbf{r}{\left[\nabla h(\mathbf{r})\right]}^2.$$

$$U_4(d)=-9k_{\mathrm{B}}T\frac{a^8}{d^8}.$$

$$2b\frac{\partial ({\rho }^{+}-{\rho }^{-})}{\partial t}=-k{q}^2({\rho }^{+}-{\rho }^{-})+2ke{q}^{4}h,$$

$$4\eta q\frac{\partial h}{\partial t}=-(\sigma q^2+\tilde{\kappa }{q}^4)h+ke{q}^2({\rho }^{+}-{\rho }^{-})+F(\mathbf{q},t),$$

$$H=\int dxdy\left[\frac{K_a}{2d_0^2}{u}^2+\frac{\gamma }{4}(\nabla u{)}^2+\frac{\kappa }{8}{\left({\nabla }^{2}u\right)}^2\right].$$

$$H=\int dxdy\left[\frac{K_a}{2d_0^2}{u}^2+\frac{\kappa c_0}{2}{\nabla }^{2}u+\frac{\kappa }{2d_0}\left(c_0-c_0^{\prime }{\Sigma }_0\right)u{\nabla }^{2}u+\frac{\kappa }{8}{\left({\nabla }^{2}u\right)}^2\right],$$

$$H=\int dxdy\left[\frac{\sigma }{d_0}u+\frac{K_a}{2d_0^2}{u}^2+\frac{K_a^{\prime }}{2}(\nabla u{)}^2+\frac{K_a^{\prime \prime }}{2}({\nabla }^{2}u{)}^2\right],$$

$$K_a^{\prime }$$

$$K_a^{\prime \prime }$$

$$H=\int dxdy\left\{\frac{K_a}{2d_0^2}{u}^2+\gamma \left[\frac{u}{d_0}+\frac{(\nabla u{)}^2}{8}\right]+\frac{\kappa }{8}{\left({\nabla }^{2}u\right)}^2\right\},$$

Full text

1 Introduction

Although membrane proteins were traditionally described as free to diffuse in the cell membrane Singer72 , it was soon acknowledged that the lipid bilayer can influence their organization and thus have an impact on many aspects of their activity Sackmann84 . Hence, interactions between proteins and the host membrane, as well as the resulting protein-protein interactions, have become fundamental topics in biophysics.

Membrane inclusions such as proteins can couple to various degrees of freedom of the membrane (curvature, thickness, composition, tilt, etc.), thus giving rise to several types of membrane-mediated interactions. It is noteworthy that these interactions are often non-specific, i.e. they do not involve the formation of chemical bonds between the various components. Thus, understanding these interactions calls for a description of the membrane as a self-assembled system whose properties are collectively determined, and not merely given by the chemical properties of the molecules involved Jensen04 . Over the last few decades, it has become clear that the concepts developed in soft matter physics to describe self-organized systems are extremely useful in this context, and that coarse-grained effective models such as the Helfrich model of membrane elasticity Helfrich73 can yield valuable insight.

In this chapter, we review the membrane-mediated interactions between inclusions such as membrane proteins that arise from direct constraints imposed by these inclusions on the shape of the membrane. Our point of view is mostly theoretical, in agreement with the history of this research field, but we also discuss the numerical and experimental results that are available. For clarity, we treat separately the effects that result from the coupling of the inclusions with membrane curvature and those that arise from their coupling with membrane thickness. Note however that a given inclusion can couple to both of these degrees of freedom. The first case, presented in Section 2, leads to interactions with a much larger range than the characteristic size of the inclusions, which will be referred to as “long-range interactions”. Such effects can be described starting from the coarse-grained Helfrich model Helfrich73 . The second case, discussed in Section 3, yields a much shorter-range interaction, and requires more detailed effective models of the membrane.

Other types of membrane-mediated interactions, arising from other underlying membrane degrees of freedom such as lipid composition and tilt, will not be discussed in detail. Besides, important applications such as the crystallization of membrane proteins and the interaction between constituents of such crystals, are outside of the scope of this chapter.

2 Long-range membrane-mediated interactions

Inclusions such as proteins are generally more rigid than the membrane. Therefore, they effectively impose constraints on the shape of the membrane, especially on its curvature, which plays a crucial part in membrane elasticity. These constraints in turn yield long-range membrane-mediated interactions between inclusions.

We will review the first theoretical predictions of these interactions, before moving on to further results in the analytically tractable regime of distant inclusions embedded in almost-flat membranes, including anisotropy, multi-body effects, and dynamics. Extensions to other geometries will then be discussed, including the compelling but tricky regime of large deformations, where numerical simulations provide useful insight. Finally, we will examine the available experimental results.

2.1 First predictions

2.1.1 Seminal paper

The existence of long-range membrane-mediated forces between inclusions in lipid membranes was first predicted in Ref. Goulian93etc . The curvature elasticity of the membrane was described by the tensionless Helfrich Hamiltonian Helfrich73 . For an up-down symmetric membrane, it reads:

[TABLE]

where $\kappa$ is the bending rigidity of the membrane and $\bar{\kappa}$ is its Gaussian bending rigidity, while $c_{1}$ and $c_{2}$ denote the local principal curvatures of the membrane, and $A$ its area. This elastic energy penalizes curvature. For small deformations of the membrane around a planar shape, Eq. 1 can be approximated by

[TABLE]

where $h(\bm{r})$ is the height of the membrane at position $\bm{r}=(x,y)\in\mathbb{R}^{2}$ with respect to a reference plane, and $(i,\,j)\in\{x,y\}^{2}$. The Hamiltonian in Eq. 2 is massless and features a translation symmetry ($h\rightarrow h+C$ where $C$ is independent of position) that is broken in a ground-state configuration, yielding Goldstone modes. The associated long-range correlations give rise to long-range membrane-mediated interactions. Neglecting the effect of the membrane tension $\sigma$, as in Eqs. 1 and 2, is legitimate below the length scale $\sqrt{\kappa/\sigma}$. Note that the simplified Hamiltonian in Eq. 2 is quadratic in the field $h$, i.e. the field theory is Gaussian.

In Ref. Goulian93etc , inclusions are characterized by bending rigidities different from those of the membrane bulk. A zone with slightly different rigidities can represent a phase-separated lipid domain, while a very rigid zone can represent a protein. Both regimes (perturbative and strong-coupling) are discussed, in the geometry of two identical circular domains of radius $a$ at large separation $d\gg a$ (see Fig. 1). An interaction potential proportional to $1/d^{4}$ is obtained in both regimes.

In the perturbative regime, the interaction depends on the perturbations of $\kappa$ and $\bar{\kappa}$ in the inclusions and on the value of $\kappa$ in the membrane, as well as on $k_{\mathrm{B}}T$.

Besides, a low-temperature interaction is obtained for rigid inclusions that impose a contact angle with the membrane, e.g. cone-shaped inclusions Goulian93etc ; Fournier97comment :

[TABLE]

where $\alpha_{1}$ and $\alpha_{2}$ are the contact angles imposed by inclusion 1 and inclusion 2 (see Fig. 1). This interaction is obtained by calculating the membrane shape that minimizes the membrane curvature energy in Eq. 2 in the presence of the inclusions. It arises from the ground-state membrane deformation due to the inclusions, and vanishes for up-down symmetric inclusions. It is repulsive. Note that this interaction does not depend on the Gaussian bending rigidity of the membrane Fournier97comment , in contrast with the perturbative case Goulian93etc . Indeed, the Gaussian curvature energy term only depends on the topology of the membrane and on boundary conditions Fournier97comment . Hence, in most subsequent studies of the membrane-mediated forces between rigid membrane inclusions, the Gaussian curvature term in Eq. 2 is discarded.

Another interaction, which is attractive and originates from the thermal fluctuations of the membrane shape, was predicted as well between rigid inclusions Goulian93etc ; Golestanian96b :

[TABLE]

Importantly, this fluctuation-induced interaction is independent of elastic constants and of contact angles. It exists even for up-down symmetric inclusions (imposing $\alpha_{1}=0$ and $\alpha_{2}=0$) that do not deform the ground-state membrane shape.

Multipole expansions valid for $a\ll d$ were used to calculate these interactions for rigid inclusions. Details on these expansions are presented in Refs. Golestanian96b ; Fournier15 . Only the leading-order terms in $a/d$ were obtained in Ref. Goulian93etc . This method was recently pushed further, yielding higher-order terms in $a/d$ Fournier15 .

2.1.2 Point-like approach

Ref. Park96 extended the study of of Ref. Goulian93etc . Membrane elasticity was described by Eq. 2 as in Ref. Goulian93etc , but different membrane-inclusion couplings were considered. Rigid inclusions were treated through a coupling Hamiltonian favoring a relative orientation of their main axis and of the normal of the membrane. The membrane-mediated interaction was calculated in the limit of very small inclusions, where the ultraviolet cutoff of the theory $\Lambda$ appears. The radius $a$ of the inclusions was related to $\Lambda$ through $\Lambda=2/a$ Park96 , yielding agreement with the results of Goulian93etc : the total interaction energy obtained is the sum of $U_{1}$ and $U_{2}$ (Eqs. 3 and 4).

This opened the way to direct point-like descriptions of membrane inclusions. In Ref. Netz97 , a perturbative approach was taken, where the coupling with the membrane and the inclusions was assumed to be linear or quadratic in the local mean curvature at the point location of the inclusion. In Ref. Kim98 , the insertion energy of a protein in the membrane was approximated by a term proportional to the Gaussian curvature of the membrane at the insertion point. Then, in Refs. Dommersnes99 ; Dommersnes99b , inclusions were modeled as more general local constraints on the membrane curvature tensor. Considering inclusions as point-like is justified in the case of membrane proteins, since their typical radius is comparable to membrane thickness, which is neglected when the membrane is considered as a surface, as in Eq. 2. This description simplifies the calculation of membrane-mediated interactions, by eliminating the need for a multipole expansion. In practice, one writes the partition function of the membrane described by the elastic energy in Eq. 2 (discarding Gaussian curvature), modeling inclusions as point curvature constraints Dommersnes99 ; Dommersnes99b . For one inclusion imposing a local isotropic curvature $c$ in $\bm{r}_{0}$, these constraints read $\partial^{2}_{x}h(\bm{r}_{0})=\partial^{2}_{y}h(\bm{r}_{0})=c$ and $\partial_{x}\partial_{y}h(\bm{r}_{0})=0$. Then, the part of the free energy that depends on the distance $d$ between the inclusions is the sum of $U_{1}$ and $U_{2}$ (Eqs. 3 and 4), where the effective radius $a$ of the point-like inclusions appears through the cutoff $\Lambda=2/a$, and the effective contact angle is $\alpha=ac$.

Refs. Yolcu11 ; Yolcu12 formalized the connection between the original description of inclusions as rigid objects Goulian93etc , and the more convenient point-like description. The effective field theory formalism developed in Refs. Yolcu11 ; Yolcu12 for membranes (see also Ref. Yolcu12b for fluid interfaces, and Ref. Yolcu14 for a review) considers inclusions as point-like particles, and captures their structure and the boundary conditions they impose via localized coupling terms. In practice, a series of generic scalar localized terms consistent with the symmetries of the system is added to the curvature energy describing the bare membrane. Each term in the series is polynomial in the derivatives of the membrane height $h$, taken at the point position of the inclusion. The coefficients of each term of the series are then obtained by matching observables, such as the ground-state membrane shape responding to an imposed background, between the full model with extended inclusions and the effective field theory Yolcu12 . These Wilson coefficients are analogous to charges, polarizabilities etc. of the inclusions and describe the interplay between the membrane and the inclusions, by encoding the long-range effects of short-range coupling Yolcu14 . Membrane-mediated interactions can be obtained from this effective field theory. It gives back the leading terms in $a/d$ obtained previously, with a generalization to inclusions with different radii, and yields higher-order corrections Yolcu11 ; Yolcu12 . This general and powerful method could be extended to complex inclusions with specific Wilson coefficients, and also enables general derivation of scaling laws through power counting. However, one should bear in mind that its existing application to rigid disk-shaped inclusions a priori yields results specific to this particular model of the inclusions. In particular, the discrepancy obtained with previous point-like approaches on certain higher-order terms Yolcu12 should be regarded as a different result obtained for a different model, since previous point-like approaches did not aim to fully mimic rigid disk-shaped inclusions. Note that higher-order terms were recently calculated in the framework of extended disks Fournier15 , showing agreement with Yolcu12 and pushing the expansion further.

2.1.3 Two types of interactions

The long-range membrane-mediated interaction between rigid inclusions comprises two leading-order terms that both depend on the fourth power of $a/d$ (Eqs. 3 and 4) Goulian93etc . Subsequent works Park96 ; Dommersnes99 ; Dommersnes99b ; Yolcu11 ; Yolcu12 demonstrated that the total interaction is the sum of these two terms, one coming from the ground-state deformation of the membrane by the inclusions (Eq. 3) and the second one arising from entropic effects (Eq. 4). However, it should be noted that the separation of these two terms is mostly of formal interest, since the ground-state shape, which is obtained by minimizing the Hamiltonian of the system, may not be of much practical relevance. In practice, one may be able to measure experimentally the average shape of a membrane, but in general it would not coincide with the ground-state one, except in the regime of small deformations. In this regime, which has been the focus of most theoretical work, the membrane Hamiltonian is quadratic (Eq. 2): then, the separation of the two terms makes sense. Let us now discuss each of these two terms.

The first term, $U_{1}$ (Eq. 3), arises from the interplay of the ground-state deformations of the membrane due to the presence of each of the inclusions, and it was first obtained in Ref. Goulian93etc by taking the (fictitious) zero-temperature limit. It also corresponds to the membrane-mediated interaction within a mean-field approximation.

The second term, $U_{2}$ (Eq. 4), is a fluctuation-induced or entropic effect, which exists even if both inclusions impose vanishing contact angles. Remarkably, in the case of rigid inclusions, the only energy scale involved is $k_{\mathrm{B}}T$: this interaction is universal. It arises from the constraints imposed by the inclusions on the thermal fluctuations of the shape of the membrane, which is a field with long-range correlations. It is analogous to the Casimir force in quantum electrodynamics (see e.g. Park96 ; Golestanian96c ; Golestanian96b ; Dommersnes99 ; Dommersnes99b ), which arises from the constraints imposed by non-charged objects (e.g. metal plates) on the quantum fluctuations of the electromagnetic field Casimir48 ; Milonni . This fluctuation-induced interaction is thus often termed “Casimir” or “Casimir-like”. In Ref. Helfrich01 , the fluctuation-induced force between membrane inclusions was recovered from the entropy loss associated to the suppression of fluctuation modes, thus reinforcing the formal analogy with the Casimir force. Fluctuation-induced forces analogous to the Casimir force exist in several other soft matter systems, where thermal fluctuations play an important part Kardar99 ; Gambassi09 . They were first discussed by Fisher and de Gennes in the context of critical binary mixtures Fisher78 . This “critical Casimir” force has been measured experimentally between a colloid and a surface immersed in a critical binary mixture Hertlein08 . Interestingly, such critical Casimir forces have been predicted to exist in membranes close to a critical point in lipid composition, and that they are very long-range, with power laws up to $(a/d)^{1/4}$ Machta12 . Their sign depends on the boundary conditions imposed by the inclusions Machta12 , as in the three-dimensional critical case Gambassi09 .

Let us now compare the magnitude of these two types of interactions. For two identical inclusions imposing the same contact angle $\alpha$, the interactions in Eq. 3 and 4 have the same modulus if

[TABLE]

Using the typical value $\kappa\approx 25\,k_{\mathrm{B}}T$ gives $|\alpha|\approx 6^{\circ}$: for larger contact angles, the mean-field repulsion dominates over the fluctuation-induced attraction.

2.2 Further developments on distant inclusions embedded in almost-flat membranes

2.2.1 In-plane anisotropy

Until now, we discussed the simple case of two inclusions with isotropic (i.e. disk-shaped) in-plane cross-section, which was the first case investigated Goulian93etc . However, real membrane inclusions, such as proteins, have various shapes. Fig. 2 shows a schematic of the different cases at stake: those in panels a and b were discussed above, and those in panels c and d will be discussed here.

In Ref. Park96 , the case of anisotropic cross-sections was treated through a coupling between membrane curvature and symmetric traceless tensor order parameters constructed from the main direction of the inclusion cross-section, integrated over the surface of the inclusion cross-section. The interaction energies obtained are anisotropic, and depend on $d$ as $1/d^{4}$ for up-down symmetric inclusions that interact only through the fluctuation-induced interaction (see Fig. 2c), just as in the case of isotropic cross-sections. However, inclusions that break the up-down symmetry of the membrane feature an anisotropic interaction with a stronger $1/d^{2}$ power law. Its angle dependence is $\cos(2(\theta_{1}+\theta_{2}))$, where $\theta_{i}$ is the angle between the main in-plane axis of inclusion $i$ and the line joining the two inclusion centers (Fig. 2c). This orientation dependence is that of a quadrupole-quadrupole interaction Golestanian96a ; Stamou00 , and the interaction energy is minimized whenever $\theta_{1}+\theta_{2}=0$ (or equivalently $\theta_{1}+\theta_{2}=\pi$). This interaction is attractive for a wide range of relative orientations, while the analogous interaction between inclusions with an isotropic cross-section is always repulsive (see Eq. 3).

The in-plane anisotropic case of rigid up-down symmetric rods imposing vanishing contact angles to the membrane on their edges was treated in Refs. Golestanian96a ; Golestanian96b . Only the fluctuation-induced interaction is then at play (as in Fig. 2c). In this study, thin rods were considered in the limit of vanishing width, and in the “distant” regime where their length $L$ is much smaller than their separation $d$. The opposite case $d\ll L$ will be discussed in Sec. 2.4.2. The power law obtained is in $1/d^{4}$, as in the case of isotropic cross sections (Eq. 4), and the only energetic scale involved in this fluctuation-induced force is $k_{\mathrm{B}}T$. The angular dependence of the interaction is $\cos^{2}[2(\theta_{1}+\theta_{2})]$, yielding energy minima for $\theta_{1}+\theta_{2}=0$ and $\pi/2$.

Anisotropic cross-sections were revisited within the point-like approach in Refs. Dommersnes99 ; Dommersnes02 . In this model, inclusions couple to the membrane by locally imposing a generic curvature tensor, with eigenvalues (principal curvatures) denoted by $K+J$ and $K-J$. The interaction between two such identical inclusions then reads, to leading order in $a/d$ Dommersnes99 ; Dommersnes02 :

[TABLE]

where $\theta_{i}$ are angles between the line joining the inclusion centers and their axis of smallest principal curvature (see Fig. 2d). This term vanishes for isotropic inclusions ($J=0$), consistently with Refs. Goulian93etc ; Park96 . Furthermore, in the fully anisotropic case $K=0$, corresponding to a saddle, the power law and the angular dependence both agree with the up-down symmetry breaking and anisotropic cross-section case of Ref. Park96 . Eq. 6 shows that in the generic case where $J$ and $K$ are nonzero, the angular degeneracy of the lowest-energy state is lifted, and (assuming without loss of generality that $K$ and $J$ have the same sign) the inclusions tend to align their axis of smallest principal curvature along the line joining their centers. Their interaction is then attractive Dommersnes99 . This interaction (Eq. 6) was recovered in Ref. Yolcu11 (with different angle notations), and generalized to inclusions with different radii.

Subleading terms in $1/d^{4}$ were also calculated in Refs. Dommersnes99 and Yolcu11 , featuring different results (as for the subleading terms in the isotropic case). One should keep in mind that the models at stake are different, since Ref. Dommersnes99 considers fully point-like inclusions while Ref. Yolcu11 models disk-shaped ones with finite radius through the effective field theory. While the agreement of these models on the leading-order term is a nice sign of robustness, there is no reason to expect an exact agreement at all orders.

Ref. Dommersnes99 also investigated the fluctuation-induced interaction, but its leading-order term was found not to be modified with respect to the isotropic case (Eq. 4). This is at variance with the anisotropy obtained in Refs. Golestanian96a ; Golestanian96b for the flat rods, but one should keep in mind that the point-like saddles do not correspond to the limit of the distant flat rods.

2.2.2 Multi-body effects and aggregation

A crucial and biologically relevant question is how long-range membrane-mediated interactions drive the collective behavior of inclusions, in particular aggregation. One would be tempted to start by summing the pairwise potentials discussed above, but these long-range membrane-mediated interactions are not pairwise additive. Non-pairwise additivity is a general feature of fluctuation-induced interactions. For instance, the existence of a three-body effect in the van der Waals–London interaction was demonstrated in Ref. Axilrod . The interaction due to the ground-state membrane deformation is not additive either. Indeed, if one considers inclusions that impose boundary conditions to the membrane on their edges, a shape minimizing the energy in the presence of one inclusions will generically not satisfy the boundary conditions imposed by the other one, yielding non-additivity Yolcu14 .

Three-body and four-body long-range membrane-mediated interactions were first calculated within a perturbative height-displacement model, breaking up-down symmetry but retaining in-plane anisotropy, in Ref. Park96 . The distance dependence of the three-body term involves terms in $1/(d_{12}^{2}d_{23}^{2})$ where $d_{ij}$ is the distance between particles $i$ and $j$. These interactions were also investigated in Ref. Netz97 , in a different perturbative approach, considering in particular inclusions that favor a given average curvature, and then in Ref. Kim98 in a point-like framework, but this particular calculation was recently shown to miss some contributions Yolcu11 .

In Ref. Dommersnes99 , the multi-body interactions and the aggregation of point-like inclusions locally imposing a curvature tensor were investigated. This generic model can include both up-down symmetry breaking and in-plane anisotropy depending on the curvature tensor imposed. The leading three-body interaction was found to involve terms in $1/(d_{12}^{2}d_{23}^{2})$, as in Ref. Park96 , and to vanish for inclusions imposing a zero curvature tensor Dommersnes99 . Monte-Carlo simulations including the full multi-body interactions were performed, allowing to study the phase diagram of the system (see Fig. 3). Polymer-like linear aggregates were obtained for sufficient values of $K$ and $J$, as predicted from the leading pairwise term (Eq. 6). A gas phase was found for small $J$, consistent with the fact that for isotropic inclusions ($J=0$) that break the up-down symmetry ($K\neq 0$), the leading pairwise interaction is repulsive (Eq. 3). Finally, for small $K$ and large $J$, aggregates were obtained, some of which had an “egg-carton” structure. This is made possible by the angular degeneracy of the lowest-energy state for $K=0$ in the leading pairwise term (Eq. 6). Multi-body interactions were shown to be quantitatively important, but the effect of the fluctuation-induced interaction (Eq. 4) was found to be negligible Dommersnes99 . The analytical calculation of multi-body effects was performed in this framework in Ref. Dommersnes02 , where the “egg-carton” aggregates were also further studied and related to experimentally-observed structures.

Coarse-grained molecular-dynamics simulations of the highly anisotropic curvature-inducing N-BAR domain proteins adhering on membranes have demonstrated linear aggregation of these proteins on the membrane. This is a first self-assembly step, which then yields the formation of meshes enabling budding Simunovic13 . This is qualitatively in good agreement with the predictions of Ref. Dommersnes99 .

The influence of the long-range elastic repulsion between isotropic inclusions that break the up-down symmetry of the membrane on their aggregation was also discussed in Ref. Destainville08 , but within a less specific framework including other types of interactions. In this work, this repulsive interaction (Eq. 3) plays the role of an energetic barrier to aggregation.

In Ref. Weikl01 , the collective behavior of inclusions locally penalizing local curvature (either only mean curvature or also Gaussian curvature) was studied using a mean-field theory for the inclusion concentration and Monte-Carlo simulations. Since the inclusions considered retain both up-down symmetry and in-plane isotropy, the only membrane-mediated interaction at play is an attractive fluctuation-induced one similar to that in Eq. 4. Direct interactions were also included. Aggregation was found to occur even for vanishing direct interactions, provided that the rigidity of the inclusions was sufficient Weikl01 . Hence, fluctuation-induced interactions may be relevant for aggregation, at least in the absence of other, stronger, interactions. Note that Eq. (4) shows that the amplitude of fluctuation-induced interactions is quite small. For instance, $d=4a$ yields $U_{2}\approx 0.02\,k_{\mathrm{B}}T$ (all the results discussed so far are strictly relevant only for $d\gg a$).

In Ref. Yolcu11 , the general effective field theory framework was used in the case of in-plane isotropic inclusions. The leading-order and next-order three-body interaction terms due to the ground-state membrane deformation between up-down symmetry-breaking inclusions were obtained, as well as the leading three-body and four-body fluctuation-induced interactions.

2.2.3 Membrane tension

Until now, we have focused on the regime where bending rigidity dominates over membrane tension. This is appropriate for length scales below $\sqrt{\kappa/\sigma}$. As $\sigma$ is in the range $10^{-6}-10^{-8}$ N/m for floppy membranes, while $\kappa\simeq 10^{-19}$ J, this length scale is then of order 1 $\mu$m. However, membrane tensions can span several orders of magnitude Rawicz00 depending on external conditions (e.g. osmotic pressure), so it is relevant to go beyond $\sqrt{\kappa/\sigma}$. For small deformations around a planar shape, the quadratic Hamiltonian of a membrane including tension reads

[TABLE]

where notations are the same as in Eq. 2, and where the Gaussian curvature term has been discarded. Note that, in a self-assembled membrane not submitted to external actions, each lipid adopts an equilibrium area. Hence, a membrane has no intrinsic surface tension (contrary to a liquid-gas interface), and stretching the membrane has an energy cost quadratic in the area variation. However, one usually considers a patch of membrane in contact with a reservoir made up by the rest of the membrane, so the tension term in Eq. 7 can be interpreted as arising from the chemical potential of this reservoir.

For length scales much larger than $\sqrt{\kappa/\sigma}$, tension dominates and Eq. 7 can be simplified into

[TABLE]

This case applies to a tense membrane at large scales, but also to a liquid interface (neglecting gravity). From a formal point of view, techniques similar to those employed in the bending-dominated case can be used, since the Hamiltonian is also quadratic with a single term.

Let us first focus on inclusions that do not break the up-down symmetry of the membrane. In Refs. Golestanian96a ; Golestanian96b , the fluctuation-induced interaction between two distant up-down symmetric rigid thin rods embedded in such a surface was calculated. It was found to be similar to the analogous bending-dominated case (see above), with the same $1/d^{4}$ power law, but with a different angular dependence.

Refs. Lehle07 ; Noruzifar09 considered the tension-dominated case of ellipsoidal colloids trapped at a fluid interface. In the case where the colloid height fluctuations are included but their contact line with the fluid is pinned, long-range fluctuation-induced interactions were obtained. This case is analogous to that of rigid in-plane anisotropic membrane inclusions preserving the up-down symmetry. Interestingly, the power law obtained was found to depend on whether or not in-plane orientational fluctuations of the colloids were allowed. If they are not allowed, the result of Refs. Golestanian96a ; Golestanian96b with the $1/d^{4}$ power law is recovered in the limit of full anisotropy. If they are allowed, a weaker anisotropic interaction with $1/d^{8}$ power law is obtained Noruzifar09 . This strong dependence of the power law of fluctuation-induced forces on boundary conditions was confirmed in Ref. Yolcu12b through the effective field theory method, in the specific case of in-plane isotropic (disk-shaped) rigid inclusions Yolcu11 ; Yolcu14 . In the case of membranes, the physical case should allow orientational fluctuations of the inclusions, and hence the $1/d^{8}$ power law should be considered. It is attractive and reads:

[TABLE]

Hence, we expect a crossover between a $1/d^{4}$ power law (Eq. 4) and a $1/d^{8}$ power law (Eq. 9) as the tension becomes more important.

In Ref. Lin11 , a scattering-matrix approach analogous to the one developed for the study of Casimir forces Genet03 ; Emig08 ; Rahi09 was developed, and applied to the full Hamiltonian in Eq. 7 including both tension and bending. The focus was on disk-shaped elastic inclusions preserving the up-down symmetry, and on their fluctuation-induced interaction. The results obtained in the case of rigid inclusions were consistent with Eq. 4 in the bending-rigidity–dominated regime, and with Eq. 9 in the tension-dominated regime. Moreover, the crossover between these two regimes was studied numerically. The method developed in Ref. Lin11 can potentially deal with more general cases, involving multiple complex inclusions. It appears to be complementary to the effective field theory method of Refs. Yolcu11 ; Yolcu14 , and was more straightforward in the transition regime where both tension and bending are relevant Lin11 .

Let us now focus on the interaction due to the ground-state deformation of the membrane. Ref. Weikl98 studied the case of conical inclusions breaking up-down symmetry but retaining in-plane isotropy, and considered the full Hamiltonian in Eq. 7. They showed that for non-vanishing tension, this interaction has a sign that depends on the relative orientation of the cones with respect to the membrane plane (i.e. on the signs of the angles they impose), contrary to the vanishing-tension case (see Eq. 3). Furthermore, at long distances between inclusions, the interaction is exponentially cut off with a decay length $\sqrt{\kappa/\sigma}$ (it involves Bessel K functions). This property was confirmed in Ref. Evans03 . Hence, at long distances, the fluctuation-induced force in Eq. 9 should dominate over the force due to the ground-state deformation. Conversely, in the case of colloids or inclusions with anisotropic cross-sections, Refs. Stamou00 and Yolcu12b demonstrated the existence of a long-range interaction due to the ground-state deformation of the membrane. The leading term of this interaction is anisotropic and decays as $1/d^{4}$.

In Ref. Simunovic15 , the effect of tension on the aggregation of the highly anisotropic curvature-inducing N-BAR domain proteins adhering on membranes was investigated through coarse-grained molecular-dynamics simulations. Increasing tension was shown to weaken the tendency of these proteins to linear aggregation, in agreement with the predicted weakening of the ground-state membrane-mediated interaction.

2.2.4 Summary of the interaction laws

Table 1 presents a summary of the power laws of the leading-order term of the membrane-mediated interactions in the various situations discussed until now.

2.2.5 External forces and torques

Until now, we have only discussed cases where inclusions couple to the membrane shape through its curvature, either explicitly or implicitly (e.g. through rigidity). This is the relevant case in the absence of external forces or torques. External forces can yield local constraints directly on the height of the membrane, e.g. quadratic ones in the case of local trapping or linear ones in the case of local pulling Netz97 . More specifically, inclusions may experience direct mechanical constraints if they are attached to the cytoskeleton, and torques in the presence of electrical fields because of their dipole moments Dommersnes99b . In these cases, one expects membrane-mediated interactions to be enhanced, because the ground-state deformations will generically be stronger than in the case where inclusions can freely reorient to minimize them, and because the constraints imposed on fluctuations will be stronger too.

The case of inclusions subject to external torques was studied in Ref. Dommersnes99b , for point-like inclusions setting a curvature tensor, in in-plane isotropic case. Both external fields strong enough to effectively pin the orientations of the inclusions, and finite external fields that set a preferred orientation, were considered. In both cases, membrane-mediated forces are strongly enhanced, even more in the strong-field case. A logarithmic fluctuation-induced interaction was obtained, as well as an interaction due to the ground-state deformation which either scales as $1/d^{2}$ if the preferred orientations are the same for both inclusions, or logarithmically if they are different. Interestingly, these interactions depend on the relative orientation of the preferred curvatures set by the inclusions, while in the torque-free case, the interaction only depends on their absolute values (see Eq. 3) Dommersnes99b .

In Ref. Lehle07 , colloids at a fluid interface were considered, with different types of boundary conditions. In the case were the position of the colloids is considered to be frozen (both in height and in orientation), strong logarithmic fluctuation-induced interactions are obtained.

2.2.6 Fluctuations of the interactions

Until now, we have discussed the average values at thermal equilibrium of membrane-mediated forces. Thermal fluctuations already play an important part since they are the physical origin of fluctuation-induced forces. But membrane-mediated forces themselves fluctuate as the shape of the membrane fluctuates. The fluctuations of these forces have been studied in Ref. Bitbol10 , using the stress tensor of the membrane Capovilla02 ; Fournier07 . This approach is inspired from those used previously for the fluctuations of Casimir forces Barton91 , and of Casimir-like forces between parallel plates imposing Dirichlet boundary conditions on a thermally fluctuating scalar field Bartolo02 .

The case of two point-like membrane inclusions that locally impose a curvature tensor was studied in Ref. Bitbol10 , for in-plane isotropic inclusions but including the up-down symmetry breaking case. Integrating the stress tensor on a contour surrounding one of the two inclusions allowed to calculate the force exerted on an inclusion by the rest of the system, in any shape of the membrane Bitbol11_force . The average of the force obtained gives back the known results Eqs. 3 and 4 that were obtained from the free energy in previous works. The variance of the force was also calculated, showing that the membrane-mediated force is dominated by its fluctuations. The distance dependence of the fluctuations, present in the sub-leading term of the variance, was also discussed. Interestingly, it shares a common physical origin with the fluctuation-induced (Casimir-like) force Bitbol10 .

2.3 Dynamics

Fundamental interactions, e.g., electrostatic ones, are usually considered as instantaneous, in the sense that they propagate at a velocity much higher than that of the particles experiencing them. This is not the case for membrane-mediated interactions, as the spreading of membrane deformations involves slow dissipative phenomena. The dynamics of membrane-mediated interactions is a promising subject for future research. Studying out-of-equilibrium membrane-mediated interactions intrinsically requires taking into account the dynamics of the membrane. Taking care both of the motion of the membrane and of that of the inclusions is very difficult. Hence, the first theoretical study in this direction to our knowledge, Ref. Fournier14 , considered two immobile inclusions that simultaneously change conformation, i.e., that simultaneously create a source of deformation, and therefore trigger a time-dependent interaction as the membrane deformation spreads dissipatively.

In Ref. Fournier14 , inclusions were modeled as simple point-like sources of mean curvature that are triggered simultaneously at $t=0$. One could imagine cylindrical integral proteins such as ion channels transforming into conical ones upon receiving a chemical signal. The time-dependent Hamiltonian of these inclusions is $\mathcal{H}_{\mathrm{inc}}(t)=\theta(t)\sum_{i}B_{i}\nabla^{2}h(\bm{r}_{i})$, with $\theta(t)$ the Heaviside step function, $B_{i}$ the curving strength, and $\bm{r}_{i}$ the position of inclusion $i$. The dynamical reaction of the membrane to such a perturbation was studied in Ref. Fournier14 .

As shown in the pioneering works of Refs. Seifert93 ; Evans94 , the dominant dissipation mechanism at short length scales is the friction between the two monolayers of the membrane. The corresponding dissipated power per unit area is $b(\bm{v}^{+}-\bm{v}^{-})^{2}$, where $\bm{v}^{\pm}$ are the velocities of the two lipid monolayers (the monolayers are denoted by $+$ and $-$) and $b\approx 10^{9}~{}\mathrm{J\,s}/\mathrm{m}^{4}$ is the intermonolayer friction coefficient. In addition, the membrane is subjected to viscous forces from the bulk solvent, of viscosity $\eta\approx 10^{-3}~{}\mathrm{J\,s}/\mathrm{m}^{3}$, and each monolayer behaves as a compressible fluid with elastic energy density $\frac{1}{2}k(\rho^{\pm}\pm e\nabla^{2}h)^{2}$. In this expression, $\rho^{\pm}$ are the monolayer relative excess densities (normalized by their equilibrium density), measured on the membrane mid-surface, $e\approx 1~{}\mathrm{nm}$ is the distance between this surface and the neutral surface of the monolayers (where density and curvature effects are decoupled) and $k\approx 0.1~{}\mathrm{J}/\mathrm{m}^{2}$. For most practical purposes, the two-dimensional viscosities of the monolayers can be neglected Fournier15b .

Taking into account all these effects, Ref. Fournier14 showed that the relaxation dynamics of a Fourier mode $\{h(q,t),\rho^{\pm}(q,t)\}$ in the membrane with two identical triggered inclusions is given, to linear order, by a set of two first-order dynamical equations:

[TABLE]

where $F(\bm{q},t)$ is the Fourier transform of $-\delta\mathcal{H}_{\mathrm{inc}}/\delta h(\bm{r},t)$, $\sigma$ is the membrane tension, and $\tilde{\kappa}=\kappa+2ke^{2}$ the bending rigidity at frozen lipid density Seifert93 . Solving these linear differential equations for time evolution and integrating over the Fourier modes $\bm{q}$ yields the time-dependent membrane deformation produced by one or more inclusions. Then the force $f(t)$ exerted by one inclusion on the other can be obtained by integrating the membrane stress tensor Capovilla02 ; Fournier07 ; Bitbol11_stress around one inclusion.

Two striking behaviors were observed in Ref. Fournier14 (see Fig. 4): (i) the force $f(t)$ reaches a maximum $f_{\mathrm{m}}$ and then decreases to the equilibrium force $f_{\mathrm{eq}}$. (ii) While $f_{\mathrm{eq}}$ decreases exponentially with the separation $d$ between the inclusions, the maximum force $f_{\mathrm{m}}$ decreases as a power-law $\sim d^{-3}$ until it reaches $f_{\mathrm{eq}}$. Hence $f_{\mathrm{m}}$ is long-ranged. Although these results were obtained with a simplified Hamiltonian for the inclusions, it is likely that the general trends observed will also apply to more realistic cases. It should be straightforward to extend the model of Ref. Fournier14 to inclusions that trigger at different times, but considering the movement of the inclusions at the same time as the movement of the membrane would be more challenging.

2.4 Other geometries

Until now, we focused on the case of inclusions with separation $d$ larger than their characteristic size, embedded in a membrane with small deformations around the flat shape. This is the case that has attracted the most attention in the literature, because of its relevance for proteins embedded in the membrane, and because of its technical tractability. We now move on to other geometries.

2.4.1 Spherical vesicle

Ref. Dommersnes98 focused on the membrane-mediated interaction arising from the ground-state deformation between two disk-shaped inclusions embedded in the membrane of a spherical vesicle, and imposing contact angles. The case of the spherical vesicle is practically relevant both in biology and in in-vitro experiments. The energy of the membrane was considered to be dominated by bending rigidity, which requires the length scales at play (in particular the vesicle radius) to be small with respect to $\sqrt{\kappa/\sigma}$. The covariant Helfrich Hamiltonian (Eq. 1 with no Gaussian curvature term) was adapted to small deformations with respect to a sphere.

The interaction was evaluated thanks to an expansion of the energy-minimizing profile of the membrane, and it was found to be strongly enhanced with respect to the flat-geometry interaction (Eq. 3) at length scales where the spherical shape of the vesicle is relevant. At sufficient angular separation, the effective power law of the interaction is $\sim 1/d^{1/3}$ Dommersnes98 . This sheds light on the strong impact of the underlying geometry of the membrane on membrane-mediated forces. Qualitatively, in a flat membrane, the interaction is weaker because the curvature energy in Eq. 1 can be minimized quite well between the inclusions (with an almost perfect saddle that has very little curvature energy), which is not possible in the spherical geometry. Similarly, in the case of external torques (Sec. 2.2.5), the imposed orientations did not allow for this low-energy saddle, thus enhancing the interaction.

2.4.2 Close parallel rods

We already discussed the case of rigid rods of length $L$, at a distance $d\gg L$ Golestanian96a ; Golestanian96b , which is close to the point-like case. The opposite regime $d\ll L$ is also relevant biologically, since it can model semi-flexible polymers adsorbed on the membrane. In Ref. Podgornik95 , the effect of the reduction of the membrane fluctuations by the presence of a semiflexible (wormlike) polymer was discussed. An effective nematic interaction was found between different segments of the polymer, and it was shown that this interaction can yield an orientational ordering transition.

Let us first consider rods that do not break the up-down symmetry of the membrane. The case of such stiff parallel rods in the limit $d\ll L$ (see Fig. 5a) embedded in a membrane with energy dominated either by bending rigidity (Eq. 2) or by tension (Eq. 8) was studied in Ref. Golestanian96c . A constant scale-free Casimir-like interaction per unit length is then expected Li92 , and indeed the Casimir-like interaction potential is then proportional to $-k_{\mathrm{B}}T\,L/d$ Golestanian96c . This interaction is much stronger than the one between point-like objects (Eq. 4), because the constraints imposed on fluctuation modes are much stronger in the geometry of parallel close rods. Ref. Golestanian96c further showed that such rods tend to bend toward one another below a certain critical distance, and that their interaction is screened by out-of-plane fluctuations if the rigidity of the polymer is finite.

This situation was further studied in Ref. Bitbol11_rods . Rods were modeled as constraints imposed on the membrane curvature along a straight line, allowing to define four types of rods, according to whether the membrane can twist along the rod and/or curve across it (see Fig. 5b-c for two examples of these rod types). The numerical prefactors of the potential in $L/d$ were obtained for interactions between the different types of rods, and they were all found to be attractive, provided that the rods are rigid, i.e. that they impose $\partial_{y}\partial_{y}h=0$ along them, with the notations of Fig. 5. However, repulsion was obtained between objects imposing completely antagonistic conditions (i.e. a rigid rod only imposing $\partial_{y}\partial_{y}h=0$ along it, see Fig. 5b, and a non-rigid “ribbon” only imposing $\partial_{x}\partial_{x}h=0$ along it), which is reminiscent of results obtained in critical binary mixtures Gambassi09 . In addition, the interaction energy was studied numerically versus $d/L$, thanks to a discretization scheme SinRonia12 , showing the transition between the asymptotic behaviors at large $d/L$ Golestanian96a and at small $d/L$ Bitbol11_rods were recovered. Finally, the bending and coming into contact of the rods due to the fluctuation-induced interaction was discussed: it was predicted to occur below a certain value of $d$ Bitbol11_rods .

The $L\gg d$ geometry gives insight into what happens between two generic inclusions that are very close to one another, through the proximity force approximation Derjaguin34 . This approximation was used in the case of disk-shaped inclusions in Refs. Lin11 ; Bitbol11_rods , showing that the fluctuation-induced interaction potential then scales as $1/d^{1/2}$.

In Ref. Weikl03 , the interaction due to the ground-state deformation between parallel rigid cylinders adsorbed on a membrane and interacting with it through an adhesion energy was studied. The membrane was assumed to be in the regime of small deformations, but both tension and bending were accounted for (see Eq. 7), and the geometry where $d\ll L$ was considered. The interaction due to the ground-state deformation was calculated explicitly in this effectively one-dimensional case. It was found to be repulsive for a pair of cylinders adhering to the same side of the membrane, and attractive for cylinders adhering to opposite sides (and hence imposing an opposite curvature). This is at variance with the point-like case, where the interaction only depends on the modulus of the curvatures imposed (see Eq. 3). The dependence in $d$ is in $\tanh(d/\sqrt{\kappa/\sigma})$ in the first case, and in $\coth(d/\sqrt{\kappa/\sigma})$ in the second one Weikl03 .

2.4.3 Large deformation regime

All cases discussed until now focused on small deformations. Then, the Hamiltonian of the membrane is quadratic, and the field theory is Gaussian. This provides tractability, both to solve the Euler-Lagrange equations that give the ground-state shape, which are then linear, and to compute thermodynamical quantities such as the free energy. Here, we will discuss the biologically relevant but much trickier regime of large deformations.

In Ref. Muller05 ; Muller05b , the covariant membrane stress and torque tensors associated to the full Helfrich Hamiltonian Capovilla02 were used to determine formal expressions of the forces between objects adsorbed on fluid membranes that are due to the ground-state deformation of the membrane. These expressions are valid without assuming small deformations, but the ground-state shape needs to be determined in order to obtain a more explicit expression. This is not an easy task in the large-deformation regime. Equilibrium shapes in the large deformation regime were further investigated in Ref. Muller07 , allowing to plot the force between cylinders, in the case of a fixed adhesion area between them and the membrane. The direction of the force and its asymptotic exponential decay at large $d/\sqrt{\kappa/\sigma}$ were found to remain the same as in the small-deformation regime Weikl03 . This situation was also investigated numerically in Ref. Mkrtchyan10 in the case of cylinders interacting with the membrane through an adhesion energy, yielding phase diagrams of the system.

In Ref. Gosselin11 , the entropic contribution to the membrane-mediated interaction between two long cylinders adsorbed on the same side of a membrane was studied in the regime of large deformations, in the case of a fixed adhesion area between the cylinders and the membrane. The free energy of the system was calculated by assuming Gaussian fluctuations around the ground-state shape. Interestingly, this entropic contribution enhances the ground-state repulsion between the two cylinders Gosselin11 , while the fluctuation-induced interaction between identical rods in the small-deformation regime is attractive Golestanian96c ; Bitbol11_rods . This is presumably a non-trivial effect coming from the non-linearities at play in the large deformations. It would be interesting to go beyond the approximation of Gaussian fluctuations around the ground-state shape.

Solving the shape/Euler-Lagrange equation for membranes beyond the domain of small deformations is technically very hard for most geometries, and incorporating fluctuations too, but numerical simulations can provide further insight. The coarse-grained molecular-dynamics membrane simulations without explicit solvent description of Ref. Reynwar07 showed that the elastic interaction between two isotropic curvature-inducing membrane inclusions (quasi-spherical caps) can become attractive at short separations, provided that the inclusions induce a strong enough curvature. Recall that the interaction due to the ground-state deformation, which is dominant with respect to the fluctuation-induced one for large enough curvatures imposed by inclusions, is always repulsive in the regime of small deformations (see Eq. 3). This hints at highly non-trivial effects of the large-deformation regime. The attractive membrane-mediated interaction was found to be able to yield aggregation of the caps and vesiculation of the membrane Reynwar07 (see Fig. 6). The case of curved phase-separated lipid domains was explored in Ref. Yuan11 through coarse-grained molecular-dynamics simulations. The interaction between domains was found to be attractive, but the angles imposed by the domains were smaller than those yielding attraction in Ref. Reynwar07 .

A numerical minimization via Surface Evolver of the Helfrich Hamiltonian Eq. 1 for a membrane with two in-plane isotropic curvature-inducing inclusions was presented in Ref. Reynwar11 , and forces were calculated by studying infinitesimal displacements. A change of sign of the membrane-mediated interaction due to the ground-state deformation of the membrane was obtained, consistently with Ref. Reynwar07 . The repulsive interaction, agreeing quantitatively with Eq. 3 at large $d/a$ and for small deformations, turned attractive for $d/a$ of order one, provided that the curvature imposed by the inclusions (and hence the membrane deformation) was large enough. The separation $d$ is defined as the center-to-center distance projected on a reference plane, while $a$ is the real radius of the inclusions, so that in the large deformation regime where inclusions are very tilted, it is possible to have $d<2a$. Attraction occurs in this regime, which is inaccessible to the small-deformation approach. Recently, Ref. Schweitzer15 studied anisotropic protein scaffolds, modeling e.g. BAR proteins, in the large-deformation regime, through similar numerical minimization methods: strongly anisotropic attractive interactions were obtained.

Ref. Bahrami12 presented a Monte-Carlo simulation of spherical nanoparticles adsorbed on a spherical vesicle modeled as a triangulated surface. Aggregation of the nanoparticles and inward tubulation of the vesicle were observed, implying strong attractive interactions. Note however that adhesion might have a strong impact on these structures Saric13 . A similar coarse-grained description of a membrane vesicle was used in Ref. Ramakrishnan13 to investigate the collective effects of anisotropic curvature-inducing inclusions, modeling e.g. BAR proteins. Vesicles were strongly deformed by the numerous inclusions, with sheet-like shapes or tubulation depending on inclusion concentration, and aggregation and nematic ordering of these inclusions were observed.

2.5 Experimental studies

While membrane-mediated interactions have been the object of significant theoretical and numerical attention, quantitative experimental tests of the theoretical predictions remain scarce to this day. A very active research area in biophysics deals with the morphological changes of the cell (invagination Roemer07 , vesiculation Reynawar07 , etc.) under the action of various proteins (see McMahon05 for a recent review). However, many other ingredients than membranes and inclusions are at play in these biological systems, for instance the cytoskeleton, out-of-equilibrium events, etc., which makes it hard to isolate membrane-mediated interactions. Biomimetic lipid membranes such as giant unilamellar vesicles Angelova86 are a good model system to study such effects. In principle, the inclusions could be real proteins, but these molecules have complex shapes, which makes it difficult to test predictions of models developed for simple geometries. Many studies have focused on the simpler and more easily controlled system of colloids adhering to membranes (see Ref. Saric13 for a review), and some have investigated interactions between phase-separated membrane domains Semrau09 . However, even in these simpler cases, membrane-mediated interactions may involve other effects, such as adhesion of the colloids, variability of contact angles imposed by domains, etc.

An experimental study of the aggregation of spherical colloidal particles adhering to biomimetic lipid membranes was presented in Ref. Koltover99 . The observed aggregation of two particles was deemed consistent with a short-range (e.g. exponential) attractive force, and no signature of a longer-range force was obtained. Note that theoretical studies predict a mostly repulsive membrane-mediated force in this geometry, except at very high deformations and small distances. Surprisingly, triplets were observed to form chains, and a linear ring-like aggregate was observed around the waist of a vesicle. Linear chain-like arrangements were also obtained in simulations of a very similar situation in Refs. Yue12 ; Saric12 , for certain sizes of particles and adhesion regimes. Ref. Saric12 used a scaling argument to show that this was not due to membrane-mediated interactions, but to the adhesion of the particles to the membrane, as a linear aggregate yields a higher adhesion area than a compact one. Such a phenomenon would thus not arise in the case of inclusions Saric13 .

Apart from proteins and colloids, another source of membrane deformation is the presence of phase separated (liquid-ordered/liquid disordered) domains, which can be partially budded. Contrast between the domains is obtained in fluorescence microscopy by adding a dye which partitions into one phase Rozovsky05 or by selectively labeling one lipid species Ursell09 . Selective deuteration can also be used to induce contrast in small-angle nuclear scattering Heberle13 . In Ref. Semrau09 , the stability of partially budded domains was interpreted as a signature of repulsive interactions, since flat ones rapidly fused. The strength of this interaction was evaluated by measuring the distribution of inter-domain distance, and then by evaluating the effective spring constant of the confining potential. It was found to be consistent with the membrane-mediated interaction arising from the ground-state deformation of a tension-free membrane in the small-deformation and large-separation regime (Eq. 3). In Ref. Ursell09 , a good agreement was obtained between the observed in-plane distribution of the domains and the predictions of the elastic theory in the presence of tension Weikl98 (see Figure 7).

3 Short-range membrane-mediated interactions

In Part 1, we dealt with long-range membrane-mediated interactions between inclusions, which arise from the curvature constraints imposed by rigid inclusions. There exist several other ways in which inclusions can couple to the surrounding membrane and thus interact with other inclusions through the membrane, but these effects are generally short-ranged. The study of these interactions was in fact initiated before that of their long-range counterparts Goulian93etc . Membrane proteins were shown experimentally to tend to immobilize neighboring lipids Jost73 . A membrane-mediated attraction between proteins was predicted to arise due to this local ordering Marcelja76 , and to decay exponentially above the correlation length of the membrane order parameter Schroder77 . Proteins can locally perturb the thickness of the membrane due to this local ordering, but they may also couple preferentially to one component of a lipid mixture Owicki78 .

Here, we are going to focus on the coupling of proteins to membrane thickness. Intrinsic membrane proteins can have a hydrophobic mismatch with the membrane: their hydrophobic thickness is slightly different from that of the unperturbed membrane. Hydrophobic mismatch is ubiquitous, and has important biological consequences, since the activity of many membrane proteins has been shown to depend on membrane thickness Killian98 . As proteins are more rigid than membranes, the membrane generically deforms in the vicinity of the protein, in order to match its thickness and avoid exposing part of the hydrophobic chains of lipids to water. This local deformation of the membrane thickness yields a membrane-mediated interaction between two such proteins.

Membrane thickness deformations are not included in the traditional Helfrich description of the membrane Helfrich73 . Describing them is tricky since they occur on the nanometer scale, which corresponds to the limit of validity of usual continuum theories where only long-scale terms are kept. Let us focus on these models before moving on to the actual interactions.

3.1 Models for local membrane thickness deformations

3.1.1 Early models

The idea that the membrane hydrophobic thickness must locally match that of an intrinsic protein was first used in theoretical descriptions of lipid-protein interactions that focused on the thermodynamic phase behavior of the lipid-protein system and on protein aggregation. In Ref. Mouritsen84 , a thermodynamic model called the “mattress model” was proposed in order to describe the phase diagrams of lipid bilayers containing proteins with a hydrophobic mismatch.

More detailed theoretical investigations of local membrane thickness deformations and of resulting membrane-mediated interactions were motivated by experimental results on the antimicrobial peptide gramicidin. In lipid membranes, two gramicidin monomers, one on each side of the bilayer, can associate to form a dimer, which acts as an ion channel. While isolated monomers do not deform the membrane, the dimeric channel generically possesses a hydrophobic mismatch with the membrane Kelkar07 . Conductivity measurements yield the formation rate and lifetime of the channel, which are directly influenced by membrane properties Kolb77 ; Elliott83 ; Goulian98 . Hence, gramicidin constitutes a very convenient experimental system to probe the effects of local membrane thickness deformations.

The first attempt to explain the dependence of gramicidin channel lifetime on the membrane thickness was provided by Ref. Elliott83 . It is based on the idea that the relevant membrane energy variation upon dimer breaking is mostly due to membrane tension, which pulls apart the monomers in a membrane with hydrophobic thickness larger than that of the dimer. The resulting estimate of the gap between the two monomers in the transition state is $\delta\simeq 1.8$ nm Elliott83 . However, this is far larger than the separation required for the breaking of the hydrogen bonds that stabilize the dimer Kelkar07 , which is of order 1 Å. Hence, this first model was not complete.

3.1.2 Huang’s model

The first full continuum model describing membrane thickness deformations was proposed in Ref. Huang86 . The Hamiltonian per unit area of the membrane was written by analogy with a smectic A liquid crystal, in which the elongated molecules organize in layers with the molecules oriented along the layers’ normal. These two systems present the same symmetries. The most important energetic terms in smectic A liquid crystals correspond to compression of the layers, and to splay distortion, i.e. curvature orthogonal to the layers dG_LC . In addition, the contribution of the “surface tension” of the membrane was included Huang86 . Restricting to symmetric deformations of the two monolayers, the effective Hamiltonian $H$ of the membrane reads Huang86

[TABLE]

In this expression, $u$ denotes the thickness excess of the membrane relative to its equilibrium thickness $d_{0}$ (see Fig. 8), $K_{a}$ is the stretching modulus of the membrane, $d_{0}$ its equilibrium thickness, $\gamma$ its “surface tension”, and $\kappa$ an elastic constant associated to splay. Finally, $x$ and $y$ denote Cartesian coordinates on the mid-plane of the membrane.

Ref. Huang86 assimilated $\gamma$ to the tension of a Plateau border and $\kappa$ to the Helfrich bending modulus, which may be questioned (see below). The corresponding typical values allowed to neglect the contribution of the “tension” term. By minimizing the resulting membrane Hamiltonian, analytical expressions were obtained for the membrane deformation profiles close to a mismatched protein such as the gramicidin channel, obtaining a decay length of a few nanometers. This model yields a satisfactory agreement with the experimental results of Ref. Elliott83 .

3.1.3 Models based on the work of Dan, Pincus and Safran

Refs. Dan93 ; ArandaEspinoza96 proposed another construction of the membrane Hamiltonian associated to thickness deformations. The energy per lipid molecule in each monolayer of the membrane was written for small deformations as a generic second-order expansion in the variation of area per lipid and in the local “curvature” of the monolayer thickness (different from the curvature of the shape of the membrane involved in the Helfrich model, which disregards thickness). Incompressibility of the lipids was used to relate the monolayer thickness and and the area per lipid. Using the same notations as in Eq. 12, and restricting again to up-down symmetric deformations of the membrane, the membrane Hamiltonian of Ref. ArandaEspinoza96 reads:

[TABLE]

where $c_{0}$ is the spontaneous curvature of a monolayer, while $c^{\prime}_{0}$ denotes its derivative with respect to the area per molecule, and $\Sigma_{0}$ the equilibrium area per lipid.

The main difference between this model and that of Ref. Huang86 is that the effect of monolayer spontaneous curvature is included in Eq. 13. It was shown in Refs. Dan93 ; ArandaEspinoza96 that this ingredient can yield oscillations in the membrane deformation profile, and in the resulting interaction potential between two mismatched proteins. Note that no “tension” term is included in Eq. 13, but the “tension” term in Eq. 12 was neglected in all the calculations of Ref. Huang86 too.

The model of Refs. Dan93 ; ArandaEspinoza96 was generalized in Refs. Brannigan06 ; Brannigan07 , where results of coarse-grained molecular-dynamics simulations for mismatched proteins in lipid membranes were also presented. The deformations of the average shape of the membrane (i.e., those usually described by the Helfrich model), and the small-scale protrusions were accounted for, as well as the symmetric thickness deformations Brannigan06 ; Brannigan07 . The effect of Gaussian curvature was also included in Ref. Brannigan07 , and Ref. Watson11 added the effect of tilt.

The model of Refs. Dan93 ; ArandaEspinoza96 ; Brannigan06 ; Brannigan07 was further generalized in Ref. Bitbol12_grami , where an additional term, proportional to the squared gradient of thickness, was included in the initial expression of the energy per lipid molecule in each monolayer of the membrane. Physically, this term should involve a microscopic interfacial tension contribution, associated to variations of the area per lipid. Note that this is different from the Plateau border tension discussed and discarded in Ref. Huang86 , as in a Plateau border, molecules can move along the surface and exchange with the bulk, yielding a smaller tension. Macroscopic membrane tension was also incorporated explicitly in Ref. Bitbol12_grami , through a chemical potential $\mu$ set by the rest of the membrane on the patch considered: $\sigma=-2\mu/\Sigma_{0}$ then plays the part of an externally applied tension. The Helfrich Hamiltonian with tension Eq. 7 was recovered from this model for average height deformations. In the case where the average shape of the membrane is flat, and integrating out anti-symmetric thickness deformations to focus on symmetric ones, the Hamiltonian reads:

[TABLE]

plus omitted boundary terms (see Ref. Bitbol12_grami ), with

[TABLE]

and the same notations as in Eqs. 12 and 13, and where the new contribution $k^{\prime}_{a}$ with respect to Eq. 13 arises from the term proportional to the squared gradient of the thickness $u$. (The definition of $u$ in Ref. Bitbol12_grami is slightly different from that of Refs. Dan93 ; ArandaEspinoza96 ; Brannigan06 ; Brannigan07 , but it does not affect the present discussion.)

The predictions of the model of Ref. Bitbol12_grami were compared with numerical profiles of membrane thickness close to a mismatched protein Brannigan06 ; Brannigan07 ; West09 , and with experimental data regarding gramicidin lifetime Elliott83 and formation rate Goulian98 . This analysis yielded consistent results for the term stemming from the gradient of the area per molecule, and its order of magnitude was found to be of order of the contribution of the interfacial tension between water and the hydrophobic part of the membrane. In addition, the presence of this new term allowed to explain for the first time a systematic dependence in previous numerical data.

3.1.4 Isotropic cross-section

The first models of short-range interactions between transmembrane proteins assumed that the proteins are coupled to a local order parameter describing the internal state of the membrane, either the conformational/chain-packing properties of the lipids, or the bilayer thickness $u$ Marcelja76 ; Schroder77 . Both are equivalent for a fully incompressible membrane hydrophobic core. In Refs. Owicki78 ; Owicki79 , a generic Landau–Ginzburg expansion of the free energy density in terms of $u$ and its first gradient was used to investigate the energy of a hexagonal lattice of embedded proteins imposing a value $u_{0}$ of the order parameter, i.e., a fixed hydrophobic mismatch, on their edge. Approximating the Wigner-Seitz cell of the lattice by a circle, which yields cylindrical symmetry, the authors derived a monotonically attractive short-range interaction caused by the overlap of the membrane regions deformed by the inclusions.

As discussed in Sec. 3.1, several models based on the thickness order parameter $u$ have been developed. They have been used to study membrane-mediated interactions. These models essentially introduced terms involving the second-order derivative of $u$, based on the (recently questioned Bitbol12_grami ) expectation that the term proportional to $(\bm{\nabla}u)^{2}$ was negligible. In particular, Ref. Huang86 introduced a term proportional to $(\nabla^{2}u)^{2}$, by analogy with the splay term for smectic liquid crystals. Later, Ref. Dan93 introduced additional terms, linear in $\nabla^{2}u$ and in $u\nabla^{2}u$, which arise from the spontaneous curvature of the monolayers and its dependence on the area per lipid. This initiated a series of works Dan93 ; Dan94 ; ArandaEspinoza96 aiming to estimate the elastic energy of a hexagonal lattice of proteins with hydrophobic mismatch. These works showed that the interaction potential can be non-monotonous, with short-distance repulsion and a minimum energy at finite separation. These effects can arise from the spontaneous curvature term, but also from a fixed contact angle between the membrane hydrophobic-hydrophilic interface and the inclusion, thereafter referred to as “slope”. The associated multi-body effects were investigated in Ref. Harroun99_theo through a Monte-Carlo simulation of inclusions fixing both the membrane thickness and its slope, in a membrane described by the elastic energy in Eq. 12. This study also demonstrated the interest of the structure factor to test the models. Another term involving second-order derivatives of the thickness profile $u$, proportional to its Gaussian curvature, was included in Ref. Brannigan07 , improving the agreement with coarse-grained molecular-dynamics numerical simulations. Note that oscillations in the interaction potential were observed in the coarse-grained molecular-dynamics simulations of Ref. Neder11 .

The term proportional to $(\bm{\nabla}u)^{2}$ in the elastic energy density was originally discarded on the grounds that it originates from a negligible microscopic surface tension assimilated to that of a Plateau border Huang86 . However, it was recently shown by us to also originate form gradients of lipid density, and therefore to contribute significantly to the elastic Hamiltonian Bitbol12_grami . Note in addition that the term in $u\nabla^{2}u$ introduced in Refs. Dan93 ; ArandaEspinoza96 contributes to the $(\bm{\nabla}u)^{2}$ term once integrated by parts.

In the end, these models converge towards the most general quadratic expansion in terms of $u$ and its first and second-order derivatives Bitbol12_grami ; Bories_tbp . In standard statistical field theory, it is justified to neglect higher-order gradients, because the focus is on large-scale physics and the coarse-graining length is much larger than the range of the microscopic interactions ChaikinLubensky . However, here, such arguments do not hold since the distortions around proteins relax on a length comparable with the bilayer thickness. Therefore, in practice, one should rather rely on comparison with experiments and simulations to determine how many terms to include in the expansion. Our current understanding is that all linear and quadratic terms involving derivatives of $u$ up to second order should contribute, and that the best strategy is to try to fit the parameters of the elastic Hamiltonian and of the protein–membrane coupling using experimental or numerical data Bories_tbp .

The focus of this chapter is on membrane-mediated interactions arising from direct constraints on the membrane shape (mean shape and thickness). Hence, we will not discuss in detail the role of the underlying lipid tilt degree of freedom Watson11 in membrane-mediated interactions Fournier98 ; May99 ; May00 ; May02 ; Bohinc03 ; Kozlovsky04 . However, tilt certainly plays a part in these interactions. For instance, proteins with no hydrophobic mismatch but with a hour-glass shape Fournier98 ; May99 may induce a membrane deformation due to the boundary conditions they impose on lipid tilt. A legitimate question, though, is how necessary it is to include this degree of freedom. Statistical physics allows to integrate out virtually any degree of freedom ChaikinLubensky . The resulting effective elasticity for the remaining degrees of freedom takes into account the underlying distortion energy of the removed ones. For instance, integrating out the tilt degree of freedom in the presence of an hour-glass shaped inclusion would produce an effective boundary energy depending on the inclusion thickness and on its angle with the membrane. What is not clear is how many orders in the derivatives of $u$, both in the bulk and in the boundary energy, one would have to introduce in order to properly account for the removed degrees of freedom. Future works in this direction could be interesting.

3.1.5 Anisotropic cross-section

While most theoretical studies of short-range membrane-mediated interactions have considered cylinder-shaped inclusions, actual membrane proteins have various shapes. As in the case of long-range interactions, in-plane anisotropy may result in directional membrane-mediated interactions, which may impact the formation of multi-protein complexes.

In Ref. Haselwandter13 , an analytical method was developed to study membrane-mediated interactions between in-plane anisotropic mismatched inclusions. The effective Hamiltonian $H$ associated to membrane thickness deformations was expressed as

[TABLE]

where we have used the notations defined in Eq. 12. This model is based on that of Ref. Huang86 (see Eq. 12), but includes an additional “tension” term in $u/d_{0}$. Such a term is also included in Ref. Bitbol12_grami (see Eq. 14), but without the assumption that its prefactor is related to that of the squared thickness gradient term. This assumption should be viewed as a simplifying hypothesis, given the contribution of monolayer curvature to the squared thickness gradient term Dan93 ; ArandaEspinoza96 ; Brannigan06 ; Brannigan07 and the difference between externally applied tension and interfacial tension Bitbol12_grami (see Sec. 3.1.3).

In Ref. Haselwandter13 , the solution of the Euler-Lagrange equation associated with Eq. 17 in the case of a single cylinder-shaped inclusion was expressed using Fourier-Bessel series. Then, using an ansatz introduced in Ref. Weikl98 in the context of long-range membrane-mediated interactions, the ground-state shape of the membrane in the presence of two inclusions was written as a sum of two such series. The coefficients of the successive terms of these series can be chosen in order to match the boundary conditions imposed by both inclusions, using expansions in $a/d<1$.

This method was extended to weakly anisotropic inclusions, modeling mechanosensitive channels of large conductance (MscL) in Ref. Haselwandter13 . The in-plane cross-section of these pentameric proteins was described as a circle perturbed by a small-amplitude sinusoidal, with fifth-oder symmetry. Boundary conditions along the edge of these proteins were expressed perturbatively in the amplitude of the sinusoidal, allowing to use the method described above. The resulting anisotropic membrane-mediated interaction features an energy barrier to dimerization, and demonstrates that the tip-on orientation is more favorable than the face-on one, except at very short distances. Gating of the MscL channel was also studied in Ref. Haselwandter13 , by modeling open and closed channels as having different diameters and hydrophobic thicknesses Ursell07 . The impact of having different oligomeric states of MscL on these interactions and on gating by tension (see Ref. Phillips09 ) was studied in Ref. Kahraman14 .

The method developed in Ref. Haselwandter13 was used in Ref. Haselwandter14 to study the effect of membrane-mediated interactions on the self-assembly and architecture of bacterial chemoreceptor lattices. Chemotaxis enables bacteria to perform directed motion in gradients of chemicals. The chemoreceptors that bind to these chemicals are transmembrane proteins that organize into large honeycomb lattices of trimers of dimers at the poles of bacteria Briegel12 . In Ref. Haselwandter14 , it was shown that membrane-mediated interactions between chemoreceptor trimers of dimers, modeled as inclusions with three-fold symmetry, correctly predict the structure of the arrays observed in experiments. Indeed, at short distances, the face-on relative orientation of the trimers is favored by these anisotropic interactions. In addition, the collective structure of the honeycomb lattice, studied approximately through the pairwise nearest-neighbor interactions, was shown to be more favorable than other types of aggregates at realistic densities of proteins. Gateway states to this lattice were also predicted, and it was shown that membrane-mediated interactions may contribute to the cooperativity of chemotactic signaling.

3.2 Numerical studies at the microscopic scale

Continuum models account for the microscopic degrees of freedom (i.e. the positions and conformations of all molecules involved) in a coarse-grained way, via effective terms in the elastic energy and the associated prefactors. However, even in the absence of a mesoscopic deformation due to hydrophobic mismatch, the presence of an inclusion constrains the configurations accessible to the lipid chains that surround it Marcelja76 ; Sintes97 ; Lague00 ; May00 . Further insight can thus be gained by treating such microscopic degrees of freedom explicitly, in particular those describing the conformation of the lipid chains. Recent advances in numerical simulations have made such approaches possible. Here, we give a brief overview of such studies. Note that numerical studies focusing on larger-scale features were mentioned earlier.

Refs. Lague00 ; Lague01 used the lateral density-density response function of the alkyl chains, obtained by molecular dynamics simulations of lipid bilayers, to determine the interaction between “smooth” (no anchoring) hard cylinders inserted into the bilayer. Three values were considered for the cylinder radius. For the largest one (9 Å, comparable to that of the gramicidin pore, for instance), the long-range interaction is repulsive for all the lipids studied (DMPC, DPPC, POPC and DOPC), with an additional short-range attraction for DMPC. This study does not discuss how the interaction might vary with the concentration of inclusions. Other studies followed suit Janosi10 ; Kik10 ; Yoo13b ; Dunton14 .

A complete description should in principle combine the effects of hydrophobic mismatch and of these changes in chain order Marcelja99 ; Bohinc03 . Such a complete model is currently lacking, due to the theoretical difficulties but also due to the dearth of experimental data that could be used to test and validate it. As in the case of lipid tilt (see Section 3.1.4), one can wonder how integrating out these underlying degrees of freedom would affect an effective model written in terms of $u$, what effective boundary conditions non-mismatched inclusions would then impose, and whether such a model would be sufficient.

3.3 Experimental studies

It has proven very difficult to directly measure the interactions between membrane inclusions.

3.3.1 Electron microscopy

First among such attempts were freeze-fracture electron microscopy (FFEM) studies Lewis83 ; Chen73 ; James73 ; Abney87 that analyzed the spatial distribution of inclusions to determine their radial distribution function $g(r)$. The data was then described using liquid state theories Pearson83 ; Pearson84 ; Braun87 in terms of a hard-core model with an additional interaction, either repulsive or attractive depending on the system.

These pioneering results were not followed by more systematic investigations, probably due to the intrinsic difficulty of the technique. It is also very difficult to check whether the distribution function observed in the sample after freezing still corresponds to that at thermal equilibrium.

3.3.2 Atomic force microscopy

It has been known for a long time that atomic force microscopy (AFM) can resolve lateral structures down to the nanometer scale Oesterhelt00 , but data acquisition used to be relatively slow. This changed with the introduction of high-speed AFM Ando01 , which allows taking “snapshots” of the system and determining the radial distribution function. The latter gives access to [OK?] the interaction potential between inclusions, as illustrated by Ref. Casuso10 for ATP-synthase c-rings in purple membranes (see Fig. 9).

3.3.3 Small-angle scattering

A promising way of studying membrane-mediated interactions is through small-angle radiation (X-ray or neutrons) scattering from oriented samples, as demonstrated by Refs. He95 ; He96 ; Yang99 . This non-invasive technique is very well adapted to measurements of membrane-mediated interactions since the wavelength used is of the same order of magnitude as the typical length scales over which one must probe the system (nanometers). One can thus measure the structure factor of the two-dimensional system formed by the inclusions in the membrane and obtain the interaction potential between them.

This strategy was recently used to study alamethicin pores in DMPC membranes Constantin07 , inorganic particles contained in bilayers of a synthetic surfactant Constantin08 ; Constantin10b and gramicidin pores in several types of membranes Constantin09 .

4 Conclusion

Membrane-mediated interactions between inclusions constitute a very rich topic. Their study gives insight into the behavior of complex two-dimensional biological membranes. In particular, these interactions may have important impacts on membrane protein aggregation, and on the formation of specific biologically functional assemblies. Interestingly, inclusions can also serve as membrane probes, since membrane-mediated interactions are in part determined by the properties of the host membrane.

The field of long-range membrane-mediated interactions has been dominated by theory, yielding interesting theoretical developments such as the fluctuation-induced interaction, the general effective field theory and scattering approaches, and the questions currently raised by the dynamics of these interactions. Some experimental and numerical studies have enriched this field, and we hope for further progress allowing for more quantitative comparison with theory.

The study of short-range membrane thickness deformations was motivated by quantitative experiments on gramicidin. Work on these deformations and on the associated membrane-mediated interactions has led to several developments of the theoretical description of membrane elasticity at the nanoscale. Importantly, the small-scale deformations involved are at the limit of the domain of validity of standard coarse-grained continuum theories, making comparison to precise experimental and numerical data even more crucially important.

An interesting fundamental feature of membrane-mediated interactions is the existence of many-body effects, arising from the interplay of the deformations caused by each of the inclusions. It would thus be particularly interesting to vary the concentration of inclusions in experiments.

Acknowledgments

A.-F. B. acknowledges support from the Human Frontier Science Program. D. C. and J.-B. F. acknowledge support by ANR Grant MEMINT (2012-BS04-0023).