The membrane that surrounds cells provides a selective barrier that allows some molecules through, but blocks the path of others. A cell’s membrane is made up of two layers of molecules with oily tails, and is therefore known as a bilayer. Many proteins are dotted within and on the inner and outer surfaces of the bilayer: some act as channels that control what goes in and out of the cell, while others protrude outside the cell so that they can sense changes in the environment.

Membrane proteins can move and interact within the bilayer, and various models have emerged to try to explain this dynamic system. These models are based on the membrane having some fluidity but also having regions where there is more structure, and typically describe the proteins as ordered clusters floating in an otherwise disordered fluid membrane. However, many researchers now think some proteins that pass through both layers of the bilayer (i.e., transmembrane proteins) make membranes more ordered, with a possibly gel-like state. However, it is not clear how transmembrane proteins can move and assemble together within such a relatively rigid membrane.

To investigate this, Katira, Mandadapu, Vaikuntanathan et al. carried out computer simulations using a model of a simple bilayer membrane. This membrane can exist in an ordered state, where the oily tails are neatly aligned, or a disordered state, where they are irregularly packed. Virtual ‘heating’ of the membrane caused it to shift from an ordered to a disordered state. When a simple transmembrane protein favoring the disordered state was inserted into the ordered state of the modeled membrane, disordered regions formed locally around the protein and the protein was able to diffuse within the membrane.

Modeling what would happen if two transmembrane proteins approached each other revealed that a consequence of the order–disorder transition is a strong attractive force that assembles the proteins together. Katira, Mandadapu, Vaikuntanathan et al. named this new phenomenon the 'orderphobic effect'. The forces arising from this effect were much greater than those currently believed to contribute to the assembly of membrane protein complexes, such as those generated by the elasticity of the membrane. This means that the orderphobic effect may be responsible for generating the protein clusters commonly seen in cell membranes.

Future work should next explore the opposite effect, where proteins favoring the ordered state are inserted into the disordered state of a membrane. This is expected to cause clustering of such proteins and thus large ordered regions in an otherwise disordered membrane.

This paper presents implications of first-order order–disorder phase transitions in lipid bilayers. The fluid mosaic model (Singer and Nicolson, 1972) and the lipid raft hypothesis (Simons and Ikonen, 1997; Munro, 2003) have guided intuition on how proteins diffuse and assemble in biological membranes—ordered clusters floating in an otherwise disordered fluid membrane (Simons and Toomre, 2000; Lingwood and Simons, 2010). However, recent advances show that a significant proportion of the membrane is liquid-ordered (Swamy et al., 2006; Owen et al., 2012; Polozov et al., 2008), with coexistence between the liquid-ordered and disordered phases. This coexistence suggests that effects of an order–disorder transition might be at play in the assembly of proteins. This possibility is studied here by examining the effects mediated by the simplest related order–disorder transition, that between solid-ordered and liquid-disordered phases.

Specifically, with molecular simulation, we study a coarse-grained model of a hydrated one-component bilayer and proteins that are added to the membrane. The model membrane exhibits two distinct phases—a solid-ordered phase and a liquid-disordered phase—and a first-order transition between them. We find that a transmembrane protein in the ordered bilayer can induce effects that resemble pre-melting (Lipowsky, 1982; 1984; Limmer and Chandler, 2014). In particular, within the otherwise ordered membrane phase, mesoscopic disordered domains surround proteins that favor disordered states. Importantly, the boundary of the domains resembles a stable, fluctuating order–disorder interface. The dynamic equilibrium established at the boundary allows the protein and its surrounding domain to diffuse. Moreover, because the interface has a finite stiffness, neighboring proteins can experience a membrane-induced force of adhesion, an attractive force that is distinctly stronger and can act over significantly larger lengths than those that can arise from simple elastic deformations of the membrane (Dan et al., 1993; Goulian et al., 1993; Phillips et al., 2009; Kim et al., 1998; Haselwandter and Phillips, 2013).

This force between transmembrane proteins is analogous to forces of interaction between hydrated hydrophobic objects. In particular, extended hydrophobic surfaces in water can nucleate vapor–liquid-like interfaces. In the presence of such interfaces, hydrophobic objects cluster to reduce the net interfacial free energy. This microscopic pre-transition effect manifesting the liquid–vapor phase transition can occur at ambient conditions (Chandler, 2005; Lum et al., 1999; Willard and Chandler, 2008; Stillinger, 1973; ten Wolde and Chandler, 2002; Mittal and Hummer, 2008; Patel et al., 2011; 2012). In the transmembrane case, we show here that a protein favoring the disordered phase creates a similar pre-transition effect. In this case it manifests the order–disorder transition of a lipid bilayer. Like the raft hypothesis, therefore, clusters do indeed form, but the mechanism for their assembly and mobility emerge as consequences of order–disorder interfaces in an otherwise ordered phase. We refer to this phenomenon as the 'orderphobic effect'.

While considering the effect with one specific order–disorder transition, one should bear in mind its generic nature. The orderphobic effect should be a general consequence of a first-order transition, whether the transition is between solid-ordered and liquid-disordered phases as considered explicitly herein, or between liquid-ordered and liquid-disordered phases as in multicomponent membrane systems. More is said on this point in the Implications section of this paper.

The order–disorder transition is a first-order phase transition

We choose the MARTINI model of hydrated dipalmitoyl phosphatidylcholine (DPPC) lipid bilayers (Marrink et al., 2007) to illustrate the orderphobic effect. See Materials and methods. This membrane model exhibits an ordered phase and a disordered phase. Figure 1A contrasts configurations from the two phases, and it shows our estimated phase boundary between the two phases. The ordered phase has regular tail packing compared to the disorganized tail arrangement of the disordered phase. A consequence of the regular tail packing is that hydrophobic thickness of the ordered phase, ${\displaystyle {\mathcal{𝒟}}_{\mathrm{o}}}$ is larger than that of the disordered phase, ${\displaystyle {\mathcal{𝒟}}_{\mathrm{d}}}$. Correspondingly, the area per lipid in the ordered phase is smaller than that in the disordered phase.

Figure 1

Download asset Open asset

Rendering the end particles of all the lipid chains in one of the two monolayers provides a convenient visual representation that distinguishes the two phases. These tail-end particles appear hexagonally-packed in the ordered phase and randomly arranged in the disordered phase. Regions that appear empty in this rendering are in fact typically filled by non tail-end particles or by tail-end particles from the other lipid monolayer.

To quantify the distinctions between the two phases, we consider a local rotational-invariant (Nelson, 2002; Halperin and Nelson, 1978; Frenkel et al., 1980), ${\displaystyle {\varphi }_l=|(1/6){\sum }_{j\in \mathrm{n}\mathrm{n}(l)}\exp (6i{\theta }_{lj}){|}^2}$, where ${\displaystyle {\theta }_{lj}}$ is the angle between an arbitrary axis and a vector connecting tail-end particle ${\displaystyle l}$ to tail-end particle ${\displaystyle j}$, and the summation is over the six nearest neighbors of particle ${\displaystyle l}$. The equilibrium average, ${\displaystyle ⟨{\varphi }_l⟩}$, is 1 for a perfect hexagonal packing, and it is 1/6 or smaller in the absence of bond-orientation correlations. Small periodically replicated samples of the hydrated DPPC membrane exhibit hysteretic changes in area per lipid and in ${\displaystyle ⟨{\varphi }_l⟩}$ during heating and cooling. See Appendix, and Marrink et al. (2005) and Rodgers et al. (2012). To establish whether the first-order-like behavior persists to large scales and thus actually manifests a phase transition, we consider larger systems and the behavior of the interface that separates the ordered and disordered phases.

Figure 1B shows coexistence for a system size of N = 3900 lipids with an interface between the two phases. To analyze interfacial fluctuations, we first identify the location of the interface at each instant. This location is found with a two-dimensional version of the three-dimensional constructions described in Limmer and Chandler (2014) and Willard and Chandler (2010). Specifically, and as discussed in Materials and methods, the interface is the line in the plane of the bilayer with an intermediate coarse-grained value of the orientational-order density,

(1) ${\displaystyle \varphi (\mathbf{r})=\sum _l{\varphi }_l\delta (\mathbf{r}-{\mathbf{r}}_l).}$

where ${\displaystyle {\mathbf{r}}_l}$ is the position of the ${\displaystyle l}$th tail-end particle projected onto a plane parallel to that of the bilayer, ${\displaystyle \mathbf{r}}$ is a two-dimensional vector specifying a position in that plane, and ${\displaystyle \delta (\mathbf{r})}$ is Dirac’s delta function. We focus on this field rather than the tail-end number density, ${\displaystyle \rho (\mathbf{r})={\sum }_l\delta (\mathbf{r}-{\mathbf{r}}_l)}$, because the difference between the two phases is larger for typical orientational-order than for typical tail-end density.

A director density field, ${\displaystyle u(\mathbf{r})={\sum }_l{u}_l\delta (\mathbf{r}-{\mathbf{r}}_l)}$, could also be used to distinguish disordered regions from ordered regions. ${\displaystyle u_l}$ would specify the degree to which the hydrophobic chain of lipid ${\displaystyle l}$ is perpendicular to the average plane of the membrane. A field of this form would be useful for systems where liquid-ordered behavior occurs in the absence of solid-ordered behavior. Multicomponent membranes, for example, can exist in solid-ordered, liquid-ordered, and liquid-disordered states. For constructing the order–disorder interface of the simple one-component membrane considered here, however, ${\displaystyle u(\mathbf{r})}$ offers little more information than ${\displaystyle \varphi (\mathbf{r})}$.

Figure 1C shows the Fourier spectrum of the height fluctuations of this interface, ${\displaystyle ⟨|\delta h_k{|}^2⟩}$. Two different system sizes are studied, with the larger system having approximately double the interface length of the smaller system. The Fourier component ${\displaystyle \delta h_k}$ is related to the height fluctuation ${\displaystyle \delta h(x)}$ as ${\displaystyle \delta h(x)={\sum }_k\delta h_k\exp (ikx)}$ where ${\displaystyle x}$ is a point along the horizontal in Figure 1B. Here, ${\displaystyle 0⩽x⩽L}$, and ${\displaystyle L}$ is the box length. With periodic boundary conditions, ${\displaystyle k=2\pi m/L}$, ${\displaystyle m=0,\pm 1,\pm 2,\cdots }$. According to capillarity theory for crystal–liquid interfaces (Nozières, 1992; Fisher et al., 1982), ${\displaystyle ⟨|\delta h_k{|}^2⟩\sim k_{\mathrm{B}}T/L\gamma k^2}$ for small ${\displaystyle k}$, with ${\displaystyle k_{\mathrm{B}}}$ being Boltzmann’s constant.

Given the proportionality with ${\displaystyle 1/k^2}$ at small ${\displaystyle k}$ (i.e., wavelengths larger than 10 nm), comparison of the proportionality constants from simulation and capillarity theory determines the interfacial stiffness (Camley et al., 2010), yielding ${\displaystyle \gamma =11.5\pm 0.46}$ pN. This value is significantly larger than the prior estimate of interfacial stiffness for this model, ${\displaystyle 3\pm 2}$ pN (Marrink et al., 2005). That prior estimate was obtained from simulations of coarsening of the ordered phase.

Because the ordered phase has a hexagonal packing, the interfacial stiffness depends on the angle between the interface and the lattice of the ordered phase. For a hexagonal lattice, there are three symmetric orientations for which the interfacial stiffnesses are equal. We will see that for the model we have simulated there appears to be only little angle dependence. Irrespective of that angle dependence, the stability of the interface and the quantitative consistency with capillary scaling provide our evidence for the order–disorder transition being a first-order transition in the model we have simulated.

The system sizes we have considered contain up to 10⁷ particles, allowing for membranes with N ≈ 10⁴ lipids, and requiring 10 μs to equilibrate. As such, our straightforward simulations are unable to determine whether the ordered phase is hexatic or crystal because correlation functions that would distinguish one from the other (Nelson et al., 1982) require equilibrating systems at least 10 times larger (Bernard and Krauth, 2011). Similarly, we are unable to determine the range of conditions for which the membranes organize with ripples and with tilted lipids (Sirota et al., 1988; Smith et al., 1990). Presumably, the ordered domain of the phase diagram in Figure 1A partitions into several subdomains coinciding with one or more of these possibilities. With advanced sampling techniques (Frenkel and Smit, 2001), free energy functions of characteristic order parameters can be computed to estimate the positions of boundaries between these various ordered behaviors. Here, we do not pursue this additional level of detail in the phase diagram because the additional boundaries refer to continuous transitions (Sirota et al., 1988). It is only the first-order transition, with its discontinuous change between ordered and disordered phases, that supports coexistence with a finite interfacial stiffness, and it is this stiffness that results in the orderphobic effect, which we turn to now.

Transmembrane proteins can disfavor the ordered membrane

A disordering (i.e., orderphobic) transmembrane protein is one that solvates more favorably in the disordered phase than in the ordered phase. The disordering effect of the protein could be produced by specific side chain structures. See Appendix. Here, in the main text, we consider a simpler mechanism. In particular, we have chosen to focus on the size of the protein’s hydrophobic thickness and the extent to which that thickness matches the thickness of the membrane’s hydrophobic layer (Killian, 1998; Sharpe et al., 2010). See Figure 2.

Figure 2

Download asset Open asset

The membrane’s hydrophobic layer is thicker in the ordered state than in the disordered state. For instance, at zero lateral pressure and 294 K in the model DPPC membrane, we find that the average thicknesses of the hydrophobic layers in the ordered and disordered states are ${\displaystyle {\mathcal{𝒟}}_{\mathrm{o}}=3.1}$ nm and ${\displaystyle {\mathcal{𝒟}}_{\mathrm{d}}=2.6}$ nm, respectively. A transmembrane protein with hydrophobic thickness of size ${\displaystyle \ell \approx 2.6}$ nm will therefore favor the structure of the disordered phase. If the protein is large enough, it can melt the ordered phase near the protein and result in the formation of an order–disorder interface.

An orderphobic protein nucleates a fluctuating order–disorder interface

Figure 3A shows a configuration of the instantaneous interface that forms around the orderphobic protein shown in Figure 2B. The interface is identified as described above. A video of its dynamics is provided as Video 1. As is common in crystal–liquid interfaces, the interface nucleated by an orderphobic protein may exhibit hexagonal faceting (Nozières, 1992), remnants of which can be observed in Figure 3A.

Figure 3

Download asset Open asset

Video 1

Download asset

The mean interface is a circle of radius ${\displaystyle {\mathcal{ℛ}}_0}$. Fourier analysis of fluctuations about that circle yields a spectrum of components. To the extent that these fluctuations obey statistics of capillary wave theory for a circular interface, the mean-square fluctuation for the ${\displaystyle k}$th component is ${\displaystyle ⟨|\delta {\mathcal{ℛ}}_k{|}^2⟩=k_{\mathrm{B}}T/2\pi \gamma k^2{\mathcal{ℛ}}_0}$, where ${\displaystyle k=m/{\mathcal{ℛ}}_0}$ and ${\displaystyle m=\pm 1}$, ${\displaystyle \pm 2,\cdots }$, and ${\displaystyle \gamma }$ is the order–disorder interfacial stiffness, neglecting the dependence on the angle between the interface and the lattice. The discrete values of ${\displaystyle k}$ reflect periodic boundary conditions going full circle around the model protein.

In Figure 3B, we use the interfacial stiffness from the free interface (${\displaystyle \gamma =11.5}$ pN) separating coexisting ordered and disordered phases with the capillary theory expression, and its corresponding spectrum, to compare with the spectrum of the protein-induced interface. The agreement between the theory, the free interface and the protein-induced interface is good, and it improves as the radius of the orderphobic protein increases and the wave vector ${\displaystyle k}$ decreases. This agreement indicates that the orderphobic protein does indeed nucleate an interface manifesting the order–disorder transition. The deviations of the fluctuations of the free interface from capillary wave theory occur for ${\displaystyle k\gtrsim 0.8}$ nm⁻¹, corresponding to wavelengths ${\displaystyle 2\pi /k\lesssim 7}$ nm, and a mean interface radius ${\displaystyle {\mathcal{ℛ}}_0\lesssim 1}$ nm. Indeed, Figure 2 suggests that even a small protein of radius 0.5 nm, which supports an interface of radius ${\displaystyle {\mathcal{ℛ}}_0\approx 1.2}$ nm, is sufficient to induce an order–disorder interface with fluctuations consistent with capillary theory.

The orderphobic effect generates forces of assembly and facilitates protein mobility

Figure 4 shows three snapshots from a typical trajectory initiated with two orderphobic proteins of radius 1.5 nm separated by a distance of 14 nm. Each induces a disordered region in its vicinity, with soft interfaces separating the ordered and disordered regions. The free energy of the separated state is approximately ${\displaystyle \gamma (P_1+P_2)}$, where ${\displaystyle P_i}$ is the perimeter of the order–disorder interface around protein ${\displaystyle i}$. On average, ${\displaystyle ⟨P_i⟩=2\pi {\mathcal{ℛ}}_0}$. After a few hundred nanoseconds, a fluctuation occurs where the two interfaces combine. While the single large interface remains intact, the finite tension of the interface pulls the two proteins together. Eventually, the tension pulls the two proteins together with a final perimeter, ${\displaystyle P_{\mathrm{f}}}$, that is typically much smaller than ${\displaystyle P_1+P_2}$. A video of its dynamics is provided as Video 2.

Figure 4

Download asset Open asset

Video 2

Download asset

After the separated interfaces join, the assembly process occurs on the time scale of microseconds. This time is required for the proteins to push away lipids that lie in the path of the assembling proteins. Given this time scale, a reversible work calculation of the binding free energy would best control both the distance and the number of lipids between the proteins. Moreover, the evident role of interfacial fluctuations indicates that the transition state ensemble for assembly must involve an interplay between inter-protein separations and lipid ordering as well as lipid concentration.

While we leave the study of reversible work surfaces and transition state ensembles to future work, it seems already clear that the net driving force for assembly is large compared to thermal energies. For example, with a model orderphobic protein radius of 1.5 nm, we find ${\displaystyle \gamma (⟨P_{\mathrm{f}}⟩-2⟨P_1⟩)\approx -30k_{\mathrm{B}}T}$. The range over which the force acts is given by the average radius of the two interfaces, ${\displaystyle 2{\mathcal{ℛ}}_0}$. This range is further amplified by the width of the interface, which is of ${\displaystyle \mathcal{𝒪}({\mathcal{ℛ}}_0)}$ for one-dimensional interfaces in two-dimensional systems (Kardar, 2007). The typical range is ≈10 to 30 nm. In comparison, given the elastic moduli of the membranes we consider, elastic responses will generate attractive forces between transmembrane proteins that are much smaller in strength and range, typically ${\displaystyle -5k_{\mathrm{B}}T}$ and 1 nm, respectively (Haselwandter and Phillips, 2013; de Meyer et al., 2008). Moreover, similarly weak and short ranged forces are found from solvation theory that accounts for linear response in microscopic detail while not accounting for the possibility of an underlying phase transition (Lagüe et al., 1998).

As in the hydrophobic effect (Chandler, 2005), the strength and range of the orderphobic force leverages the power of a phase transition, depending in this case on the ability of the orderphobic protein to induce a disordered layer in its vicinity. This ability depends upon the proximity to the membrane’s phase transition, and, for the simple protein models considered in this paper, it depends upon the protein’s radius and hydrophobic mismatch with the membrane. The spatial extent of the disordered region increases with proximity to phase coexistence as shown in Figure 5A.

Figure 5

Download asset Open asset

Furthermore, Figure 5B shows that the strength of the effect is maximal for a hydrophobic thickness equal to that of the disordered phase, and it decreases as the hydrophobic thickness approaches that of the ordered phase. In the case of zero mismatch (i.e., ${\displaystyle \ell ={\mathcal{𝒟}}_{\mathrm{o}}}$) the value of the order parameter in the vicinity of the protein is consistent with that of a pure bilayer in the ordered state. Therefore, the model proteins with zero mismatch do not induce a disordered region, and the orderphobic effect vanishes. See Figure 5B and D.

Figure 4 also shows that the orderphobic effect produces excess mobility, by proteins melting order in a surrounding microscopic layer and by facilitating the motions of neighboring proteins. This finding explains how protein mobility and reorganization can be relatively facile in the so-called 'gel' phases of membranes. Further information on this phenomenon is provided in Appendix. Our prediction of enhanced lipid mobility surrounding orderphobic proteins may be amenable to experimental tests by single molecule tracking techniques (Eggeling et al., 2009).

In this Appendix, we examine hysteresis in the bilayer system, show results for mean-square displacements as functions of time, and provide evidence that a small transmembrane ${\displaystyle \alpha }$-helix protein can be orderphobic if the side chains of the protein disrupt the packing of lipid tails.

The lipid bilayer system exhibits hysteresis

Appendix 1—figure 1 shows the change in area per lipid with temperature while heating and cooling a bilayer. There are finite jumps in area per lipid as the system transitions between the two phases, suggesting a first-order phase transition. Hysteresis occurs because ordering from the metastable disordered phase is much slower than disordering from the metastable ordered phase. Due to the difference in time scales, when contrasting melting and freezing from heating and cooling runs, the melting points from heating runs as shown in Appendix 1—figure 1A provide the more accurate estimates of the actual phase boundaries. Systematic errors due to small system size and heating rate have not been estimated.

Appendix 1—figure 1

Download asset Open asset

The pre-melting layer has a higher mobility than the ordered phase

In Appendix 1—figure 2 we show the mean squared displacements of lipids in the bulk ordered phase, bulk disordered phase, and the pre-melting layer induced by an orderphobic protein. As discussed in the main text, an orderphobic protein increases the mobility of the lipids in its vicinity. Note that the center of mass of the membrane fluctuates in time. These fluctuations affect the absolute positions of lipid molecules, but they are irrelevant to the issue of lipid mobility. Therefore, the mean-square displacements considered in Appendix 1—figure 2 are for tail-end particle positions relative to the instantaneous position of the membrane’s center of mass. That is to say, for ${\displaystyle ⟨|{\overline{\mathbf{r}}}_l(t)-{\overline{\mathbf{r}}}_l(0){|}^2⟩}$, where ${\displaystyle {\overline{\mathbf{r}}}_l(t)}$ is the position at time ${\displaystyle t}$ of the ${\displaystyle l}$th tail-end particle less that of the membrane’s center of mass.

Appendix 1—figure 2

Download asset Open asset

For large enough times, ${\displaystyle t}$, the mean-square displacements are asymptotic to ${\displaystyle 4Dt}$, where ${\displaystyle D}$ is the self-diffusion constant. For the disordered liquid phase, we see from Appendix 1—figure 2 that the asymptotic region is reached within 10² ns, and that ${\displaystyle D\approx 4\times {10}^{-7}}$ cm²/s. In contrast, for the ordered phase, the diffusive asymptotic limit is not reached until 10⁴ ns, and ${\displaystyle D}$${\displaystyle \approx 2\times {10}^{-9}}$ cm²/s. The mobility of lipids within the disordered layer surrounding the orderphobic protein is an order of magnitude larger than that of lipids beyond that region and in the ordered phase.

Experimental results for lipid diffusion constants in disordered and ordered bilayers are ${\displaystyle D\approx 3\times {10}^{-8}}$ cm²/s and ${\displaystyle D\approx 2\times {10}^{-10}}$ cm²/s, respectively (Korlach et al., 1999). The simulation is in harmony with experiment for the two-order of magnitude difference between the ordered- and disordered-phase values of ${\displaystyle D}$. Of course, absolute values of ${\displaystyle D}$ are beyond the scope of what can be predicted from our simulations because coarse graining omits degrees of freedom that would increase friction and decrease ${\displaystyle D}$.

A model ${\displaystyle \alpha }$-helix is orderphobic

Here, we consider the MARTINI model for KALP23—a polypeptide chain with 23 residues, consisting of alanines and leucines flanked by lysines (Schäfer et al., 2011). This molecule has a hydrophobic thickness of ${\displaystyle \ell \approx 3.0}$ nm, which means that it has essentially no hydrophobic-length mismatch with the ordered bilayer. Nevertheless, it is orderphobic because its side chains perturb the ordered lipid phase to an extent that a pre-melting layer is formed around the protein. This behavior is demonstrated with the aid of Appendix 1—figure 3. The panels render configurations from a simulation in which we have placed this model protein in the MARTINI model for the ordered DPPC bilayer system considered in the main text. The pre-melting layer that forms around the protein causes the protein to tilt so as to keep its full hydrophobic length in contact with hydrophobic tails of the lipids. The interface separating its disordered domain from the surrounding ordered phase remains stable throughout a molecular dynamics trajectory running for more than 1 µs.

Appendix 1—figure 3

Download asset Open asset

An orderphilic protein in the disordered phase nucleates an ordered domain

In a disordered bilayer, we embed a model protein with a hydrophobic thickness approximately equal to that of the ordered phase, ${\displaystyle \ell =3.3}$ nm. This embedding illustrates the orderphilic case. As a control, we also embed a model protein with a hydrophobic thickness approximately equal to that of the disordered phase, ${\displaystyle \ell =2.3}$ nm. For both cases, we then calculate the average of the bond orientational order density ${\displaystyle ⟨\varphi (r)⟩}$, and the number density ${\displaystyle ⟨\rho (r)⟩}$. We also compute the director density, ${\displaystyle ⟨u(r)⟩}$ as defined earlier, with ${\displaystyle u_l=(3/2){\cos }^2({\theta }_l)-(1/2)}$. Here, ${\displaystyle {\theta }_l}$ is the angle between the membrane normal and the ${\displaystyle l}$th chain’s orientation. The director density for ${\displaystyle \ell =3.3}$ nm as seen in Appendix 1—figure 4 changes approximately sigmoidally, connecting its values of 0.92 and 0.62 in the ordered and disordered phases, respectively. Further, we have computed the instantaneous interfaces from the coarse graining of each of these three fields ${\displaystyle \rho (\mathbf{r})}$, ${\displaystyle \varphi (\mathbf{r})}$ and ${\displaystyle u(\mathbf{r})}$. In the orderphilic case, interfaces are apparent for all three fields. The director field provides the clearest pictures, which are illustrated in Appendix 1—figures 4,5. In the control case, interfaces do not appear.

Appendix 1—figure 4

Download asset Open asset

Appendix 1—figure 5

Download asset Open asset

Distinctions between interfaces for each of the fields are microscopic effects worthy of future study. For the bulk interface between the ordered and disordered phases, the distinction disappears because there is only one order–disorder transition for the membrane model we have considered. In more complicated membrane models, those with mixtures of components exhibiting both liquid-ordered and solid-ordered phases, the different fields offer different information for both large and small length scales.

1. 1. Bernard EP
   2. Krauth W

   (2011) Two-step melting in two dimensions: first-order liquid-hexatic transition

   Physical Review Letters 107:155704.

   https://doi.org/10.1103/PhysRevLett.107.155704

   • PubMed
   • Google Scholar
2. 1. Camley BA
   2. Esposito C
   3. Baumgart T
   4. Brown FL

   (2010) Lipid bilayer domain fluctuations as a probe of membrane viscosity

   Biophysical Journal 99:L44–L46.

   https://doi.org/10.1016/j.bpj.2010.07.007

   • PubMed
   • Google Scholar
3. Book

   1. Chandler D

   (1987)

   Introduction to Modern Statistical Mechanics

   Oxford University Press.

   • Google Scholar
4. 1. Chandler D

   (2005) Interfaces and the driving force of hydrophobic assembly

   Nature 437:640–647.

   https://doi.org/10.1038/nature04162

   • PubMed
   • Google Scholar
5. 1. Dan N
   2. Pincus P
   3. Safran SA

   (1993) Membrane-induced interactions between inclusions

   Langmuir 9:2768–2771.

   https://doi.org/10.1021/la00035a005

   • Google Scholar
6. 1. de Meyer FJ
   2. Venturoli M
   3. Smit B

   (2008) Molecular simulations of lipid-mediated protein-protein interactions

   Biophysical Journal 95:1851–1865.

   https://doi.org/10.1529/biophysj.107.124164

   • PubMed
   • Google Scholar
7. 1. Domański J
   2. Marrink SJ
   3. Schäfer LV

   (2012) Transmembrane helices can induce domain formation in crowded model membranes

   Biochimica Et Biophysica Acta (BBA) - Biomembranes 1818:984–994.

   https://doi.org/10.1016/j.bbamem.2011.08.021

   • Google Scholar
8. 1. Eggeling C
   2. Ringemann C
   3. Medda R
   4. Schwarzmann G
   5. Sandhoff K
   6. Polyakova S
   7. Belov VN
   8. Hein B
   9. von Middendorff C
   10. Schönle A
   11. Hell SW

   (2009) Direct observation of the nanoscale dynamics of membrane lipids in a living cell

   Nature 457:1159–1162.

   https://doi.org/10.1038/nature07596

   • PubMed
   • Google Scholar
9. 1. Fasshauer D
   2. Sutton RB
   3. Brunger AT
   4. Jahn R

   (1998) Conserved structural features of the synaptic fusion complex: SNARE proteins reclassified as Q- and R-SNAREs

   Proceedings of the National Academy of Sciences of the United States of America 95:15781–15786.

   https://doi.org/10.1073/pnas.95.26.15781

   • PubMed
   • Google Scholar
10. 1. Fisher MPA
    2. Fisher DS
    3. Weeks JD

    (1982) Agreement of capillary-wave theory with exact results for the interface profile of the two-dimensional Ising model

    Physical Review Letters 48:368.

    https://doi.org/10.1103/PhysRevLett.48.368

    • Google Scholar
11. 1. Fratti RA
    2. Jun Y
    3. Merz AJ
    4. Margolis N
    5. Wickner W

    (2004) Interdependent assembly of specific regulatory lipids and membrane fusion proteins into the vertex ring domain of docked vacuoles

    The Journal of Cell Biology 167:1087–1098.

    https://doi.org/10.1083/jcb.200409068

    • Google Scholar
12. 1. Frenkel D
    2. Hanson FE
    3. McTague JP

    (1980)

    Ordering in Strongly Fluctuating Condensed Matter Systems

    285–292, Phase Transitions and Orientational Order in a Two Dimensional Lennard-Jones System, Ordering in Strongly Fluctuating Condensed Matter Systems, Boston, MA, Springer US, 10.1007/978-1-4684-3626-6_20.

    • Google Scholar
13. Book

    1. Frenkel D
    2. Smit B

    (2001)

    Understanding Molecular Simulation: From Algorithms to Applications

    Academic Press.

    • Google Scholar
14. 1. Goulian M
    2. Bruinsma R
    3. Pincus P

    (1993) Long-range forces in heterogeneous fluid membranes

    Europhysics Letters 22:145–150.

    https://doi.org/10.1209/0295-5075/22/2/012

    • Google Scholar
15. 1. Halperin BI
    2. Nelson DR

    (1978) Theory of two-dimensional melting

    Physical Review Letters 41:121–124.

    https://doi.org/10.1103/PhysRevLett.41.121

    • Google Scholar
16. 1. Haselwandter CA
    2. Phillips R

    (2013) Directional interactions and cooperativity between mechanosensitive membrane proteins

    EPL 101:68002.

    https://doi.org/10.1209/0295-5075/101/68002

    • Google Scholar
17. 1. James JR
    2. Vale RD

    (2012) Biophysical mechanism of T-cell receptor triggering in a reconstituted system

    Nature 487:64–69.

    https://doi.org/10.1038/nature11220

    • PubMed
    • Google Scholar
18. Book

    1. Kardar M

    (2007) Statistical Physics of Fields

    Cambridge: Cambridge University Press.

    https://doi.org/10.1017/CBO9780511815881

    • Google Scholar
19. 1. Killian JA

    (1998) Hydrophobic mismatch between proteins and lipids in membranes

    Biochimica Et Biophysica Acta (BBA) - Reviews on Biomembranes 1376:401–416.

    https://doi.org/10.1016/S0304-4157(98)00017-3

    • Google Scholar
20. 1. Kim KS
    2. Neu J
    3. Oster G

    (1998) Curvature-mediated interactions between membrane proteins

    Biophysical Journal 75:2274–2291.

    https://doi.org/10.1016/S0006-3495(98)77672-6

    • PubMed
    • Google Scholar
21. 1. Korlach J
    2. Schwille P
    3. Webb WW
    4. Feigenson GW

    (1999) Characterization of lipid bilayer phases by confocal microscopy and fluorescence correlation spectroscopy

    Proceedings of the National Academy of Sciences of the United States of America 96:8461–8466.

    https://doi.org/10.1073/pnas.96.15.8461

    • Google Scholar
22. 1. Lagüe P
    2. Zuckermann MJ
    3. Roux B

    (1999) Protein inclusion in lipid membranes: A theory based on the hypernetted chain integral equation

    Faraday Discussions 111:165–172.

    https://doi.org/10.1039/a807109h

    • Google Scholar
23. 1. Limmer DT
    2. Chandler D

    (2014) Premelting, fluctuations, and coarse-graining of water-ice interfaces

    The Journal of Chemical Physics 141:18C505.

    https://doi.org/10.1063/1.4895399

    • PubMed
    • Google Scholar
24. 1. Lingwood D
    2. Simons K

    (2010) Lipid rafts as a membrane-organizing principle

    Science 327:46–50.

    https://doi.org/10.1126/science.1174621

    • PubMed
    • Google Scholar
25. 1. Lipowsky R

    (1982) Critical surface phenomena at first-order bulk transitions

    Physical Review Letters 49:1575–1578.

    https://doi.org/10.1103/PhysRevLett.49.1575

    • Google Scholar
26. 1. Lipowsky R

    (1984) Surface-induced order and disorder: Critical phenomena at first-order phase transitions (invited)

    Journal of Applied Physics 55:2485–2490.

    https://doi.org/10.1063/1.333703

    • Google Scholar
27. 1. Lum K
    2. Chandler D
    3. Weeks JD

    (1999) Hydrophobicity at small and large length scales

    The Journal of Physical Chemistry B 103:4570–4577.

    https://doi.org/10.1021/jp984327m

    • Google Scholar
28. 1. Machta BB
    2. Veatch SL
    3. Sethna JP

    (2012) Critical Casimir forces in cellular membranes

    Physical Review Letters 109:138101.

    https://doi.org/10.1103/PhysRevLett.109.138101

    • PubMed
    • Google Scholar
29. 1. Marrink SJ
    2. Risselada J
    3. Mark AE

    (2005) Simulation of gel phase formation and melting in lipid bilayers using a coarse grained model

    Chemistry and Physics of Lipids 135:223–244.

    https://doi.org/10.1016/j.chemphyslip.2005.03.001

    • Google Scholar
30. 1. Marrink SJ
    2. Risselada HJ
    3. Yefimov S
    4. Tieleman DP
    5. de Vries AH

    (2007) The MARTINI force field: coarse grained model for biomolecular simulations

    The Journal of Physical Chemistry. B 111:7812–7824.

    https://doi.org/10.1021/jp071097f

    • PubMed
    • Google Scholar
31. 1. Milovanovic D
    2. Honigmann A
    3. Koike S
    4. Göttfert F
    5. Pähler G
    6. Junius M
    7. Müllar S
    8. Diederichsen U
    9. Janshoff A
    10. Grubmüller H
    11. Risselada HJ
    12. Eggeling C
    13. Hell SW
    14. van den Bogaart G
    15. Jahn R

    (2015) Hydrophobic mismatch sorts SNARE proteins into distinct membrane domains

    Nature Communications 6:5984.

    https://doi.org/10.1038/ncomms6984

    • PubMed
    • Google Scholar
32. 1. Mittal J
    2. Hummer G

    (2008) Static and dynamic correlations in water at hydrophobic interfaces

    Proceedings of the National Academy of Sciences of the United States of America 105:20130–20135.

    https://doi.org/10.1073/pnas.0809029105

    • PubMed
    • Google Scholar
33. 1. Munro S

    (2003) Lipid Rafts

    Cell 115:377–388.

    https://doi.org/10.1016/S0092-8674(03)00882-1

    • Google Scholar
34. 1. Nag K
    2. Perez-Gil J
    3. Ruano ML
    4. Worthman LA
    5. Stewart J
    6. Casals C
    7. Keough KMWKM

    (1998) Phase transitions in films of lung surfactant at the air-water interface

    Biophysical Journal 74:2983–2995.

    https://doi.org/10.1016/S0006-3495(98)78005-1

    • PubMed
    • Google Scholar
35. 1. Nelson DR
    2. Rubinstein M
    3. Spaepen F

    (1982) Order in two-dimensional binary random arrays

    Philosophical Magazine A 46:105–126.

    https://doi.org/10.1080/01418618208236211

    • Google Scholar
36. Book

    1. Nelson DR

    (2002)

    Defects and Geometry in Condensed Matter Physics

    Cambridge University Press.

    • Google Scholar
37. 1. Nozières P

    (1992)

    Solids Far From Equilibrium

    1–154, Shape and growth of crystals, Solids Far From Equilibrium, Cambridge, Cambridge University Press.

    • Google Scholar
38. 1. Owen DM
    2. Williamson DJ
    3. Magenau A
    4. Gaus K

    (2012) Sub-resolution lipid domains exist in the plasma membrane and regulate protein diffusion and distribution

    Nature Communications 3:1256–11257.

    https://doi.org/10.1038/ncomms2273

    • PubMed
    • Google Scholar
39. 1. Patel AJ
    2. Varilly P
    3. Jamadagni SN
    4. Acharya H
    5. Garde S
    6. Chandler D

    (2011) Extended surfaces modulate hydrophobic interactions of neighboring solutes

    Proceedings of the National Academy of Sciences of the United States of America 108:17678–17683.

    https://doi.org/10.1073/pnas.1110703108

    • Google Scholar
40. 1. Patel AJ
    2. Varilly P
    3. Jamadagni SN
    4. Hagan MF
    5. Chandler D
    6. Garde S

    (2012) Sitting at the edge: how biomolecules use hydrophobicity to tune their interactions and function

    The Journal of Physical Chemistry. B 116:2498–2503.

    https://doi.org/10.1021/jp2107523

    • PubMed
    • Google Scholar
41. 1. Phillips R
    2. Ursell T
    3. Wiggins P
    4. Sens P

    (2009) Emerging roles for lipids in shaping membrane-protein function

    Nature 459:379–385.

    https://doi.org/10.1038/nature08147

    • PubMed
    • Google Scholar
42. 1. Polozov IV
    2. Bezrukov L
    3. Gawrisch K
    4. Zimmerberg J

    (2008) Progressive ordering with decreasing temperature of the phospholipids of influenza virus

    Nature Chemical Biology 4:248–255.

    https://doi.org/10.1038/nchembio.77

    • PubMed
    • Google Scholar
43. 1. Pronk S
    2. Páll S
    3. Schulz R
    4. Larsson P
    5. Bjelkmar P
    6. Apostolov R
    7. Shirts MR
    8. Smith JC
    9. Kasson PM
    10. van der Spoel D
    11. Hess B
    12. Lindahl E

    (2013) GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

    Bioinformatics 29:845–854.

    https://doi.org/10.1093/bioinformatics/btt055

    • PubMed
    • Google Scholar
44. 1. Qi SY
    2. Groves JT
    3. Chakraborty AK

    (2001) Synaptic pattern formation during cellular recognition

    Proceedings of the National Academy of Sciences of the United States of America 98:6548–6553.

    https://doi.org/10.1073/pnas.111536798

    • PubMed
    • Google Scholar
45. 1. Risselada HJ
    2. Marrink SJ

    (2008) The molecular face of lipid rafts in model membranes

    Proceedings of the National Academy of Sciences of the United States of America 105:17367–17372.

    https://doi.org/10.1073/pnas.0807527105

    • Google Scholar
46. 1. Rodgers JM
    2. Sørensen J
    3. de Meyer FJ
    4. Schiøtt B
    5. Smit B

    (2012) Understanding the phase behavior of coarse-grained model lipid bilayers through computational calorimetry

    The Journal of Physical Chemistry. B 116:1551–1569.

    https://doi.org/10.1021/jp207837v

    • PubMed
    • Google Scholar
47. Book

    1. Rowlinson JS
    2. Widom B

    (1982)

    Molecular Theory of Capillarity

    Oxford: Clarendon Press.

    • Google Scholar
48. 1. Różycki B
    2. Boura E
    3. Hurley JH
    4. Hummer G

    (2012) Membrane-elasticity model of Coatless vesicle budding induced by ESCRT complexes

    PLoS Computational Biology 8:e1002736.

    https://doi.org/10.1371/journal.pcbi.1002736

    • PubMed
    • Google Scholar
49. 1. Schäfer LV
    2. de Jong DH
    3. Holt A
    4. Rzepiela AJ
    5. de Vries AH
    6. Poolman B
    7. Killian JA
    8. Marrink SJ

    (2011) Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes

    Proceedings of the National Academy of Sciences of the United States of America 108:1343–1348.

    https://doi.org/10.1073/pnas.1009362108

    • PubMed
    • Google Scholar
50. 1. Sharpe HJ
    2. Stevens TJ
    3. Munro S

    (2010) A comprehensive comparison of transmembrane domains reveals organelle-specific properties

    Cell 142:158–169.

    https://doi.org/10.1016/j.cell.2010.05.037

    • PubMed
    • Google Scholar
51. 1. Simons K
    2. Ikonen E

    (1997) Functional rafts in cell membranes

    Nature 387:569–572.

    https://doi.org/10.1038/42408

    • PubMed
    • Google Scholar
52. 1. Simons K
    2. Toomre D

    (2000) Lipid rafts and signal transduction

    Nature Reviews. Molecular Cell Biology 1:31–39.

    https://doi.org/10.1038/35036052

    • PubMed
    • Google Scholar
53. 1. Singer SJ
    2. Nicolson GL

    (1972) The fluid mosaic model of the structure of cell membranes

    Science 175:720–731.

    https://doi.org/10.1126/science.175.4023.720

    • PubMed
    • Google Scholar
54. 1. Sirota EB
    2. Smith GS
    3. Safinya CR
    4. Plano RJ
    5. Clark NA

    (1988) X-ray scattering studies of aligned, stacked surfactant membranes

    Science 242:1406–1409.

    https://doi.org/10.1126/science.242.4884.1406

    • PubMed
    • Google Scholar
55. 1. Smith GS
    2. Sirota EB
    3. Safinya CR
    4. Plano RJ
    5. Clark NA

    (1990) X-ray structural studies of freely suspended ordered hydrated DMPC multimembrane films

    The Journal of Chemical Physics 92:4519–4529.

    https://doi.org/10.1063/1.457764

    • Google Scholar
56. 1. Stein A
    2. Weber G
    3. Wahl MC
    4. Jahn R

    (2015) Helical extension of the neuronal SNARE complex into the membrane

    Nature 460:525–528.

    https://doi.org/10.1038/nature08156

    • Google Scholar
57. 1. Stillinger FH

    (1973) Structure in aqueous solutions of nonpolar solutes from the standpoint of scaled-particle theory

    Journal of Solution Chemistry 2:141–158.

    https://doi.org/10.1007/BF00651970

    • Google Scholar
58. 1. Swamy MJ
    2. Ciani L
    3. Ge M
    4. Smith AK
    5. Holowka D
    6. Baird B
    7. Freed JH

    (2006) Coexisting domains in the plasma membranes of live cells characterized by spin-label ESR spectroscopy

    Biophysical Journal 90:4452–4465.

    https://doi.org/10.1529/biophysj.105.070839

    • PubMed
    • Google Scholar
59. 1. ten Wolde PR
    2. Chandler D

    (2002) Drying-induced hydrophobic polymer collapse

    Proceedings of the National Academy of Sciences of the United States of America 99:6539–6543.

    https://doi.org/10.1073/pnas.052153299

    • PubMed
    • Google Scholar
60. 1. Veatch SL
    2. Keller SL

    (2005) Miscibility phase diagrams of giant vesicles containing sphingomyelin

    Physical Review Letters 94:148101.

    https://doi.org/10.1103/PhysRevLett.94.148101

    • PubMed
    • Google Scholar
61. 1. Wickner W
    2. Schekman R

    (2008)

    Membrane fusion

    Nature Structural & Molecular Biology 15:658–664.

    • PubMed
    • Google Scholar
62. 1. Willard AP
    2. Chandler D

    (2008) The role of solvent fluctuations in hydrophobic assembly

    The Journal of Physical Chemistry. B 112:6187–6192.

    https://doi.org/10.1021/jp077186+

    • PubMed
    • Google Scholar
63. 1. Willard AP
    2. Chandler D

    (2010) Instantaneous liquid interfaces

    The Journal of Physical Chemistry. B 114:1954–1958.

    https://doi.org/10.1021/jp909219k

    • PubMed
    • Google Scholar
64. 1. Zick M
    2. Stroupe C
    3. Orr A
    4. Douville D
    5. Wickner WT

    (2014) Membranes linked by trans-SNARE complexes require lipids prone to non-bilayer structure for progression to fusion

    eLife 3:e01879.

    https://doi.org/10.7554/eLife.01879

    • PubMed
    • Google Scholar

Author details

1. Shachi Katira

   Department of Chemistry, University of California, Berkeley, Berkeley, United States

   Contribution

   SK, SK, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

   Contributed equally with

   Kranthi K Mandadapu and Suriyanarayanan Vaikuntanathan

   Competing interests

   The author declares that no competing interests exist.

2. Kranthi K Mandadapu

   1. Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, United States
   2. Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, United States

   Contribution

   KKM, KM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

   Contributed equally with

   Shachi Katira and Suriyanarayanan Vaikuntanathan

   Competing interests

   The author declares that no competing interests exist.

3. Suriyanarayanan Vaikuntanathan

   Department of Chemistry, University of Chicago, Chicago, United States

   Contribution

   SV, SV, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

   Contributed equally with

   Shachi Katira and Kranthi K Mandadapu

   Competing interests

   The author declares that no competing interests exist.

4. Berend Smit

   1. Department of Chemistry, University of California, Berkeley, Berkeley, United States
   2. Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, United States
   3. Laboratory of Molecular Simulation, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne, Sion, Switzerland

   Contribution

   BS, BS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

   Competing interests

   The author declares that no competing interests exist.

5. David Chandler

   Department of Chemistry, University of California, Berkeley, Berkeley, United States

   Contribution

   DC, DC, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

   For correspondence

   chandler@berkeley.edu

   Competing interests

   The author declares that no competing interests exist.

• 3,225

  views

• 764

  downloads

• 59

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.