# Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains

Moritz Frewein, Benjamin Kollmitzer, Peter Heftberger, Georg Pabst

► To cite this version:

Moritz Frewein, Benjamin Kollmitzer, Peter Heftberger, Georg Pabst. Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains. Soft Matter, 2016, 12, pp.3189 - 3195. 10.1039/c6sm00042h. hal-05454298

HAL Id: hal-05454298

https://hal.science/hal-05454298v1

Submitted on 12 Jan 2026

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

CrossMark

← click for updates

Cite this: *Soft Matter*, 2016, **12**, 3189

Received 7th January 2016,
Accepted 12th February 2016

DOI: 10.1039/c6sm00042h

www.rsc.org/softmatter

# Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains

Moritz Frewein,ab Benjamin Kollmitzer,c Peter Heftberger†ab and Georg Pabst*ab

We have studied the contributions of stored elastic energies in liquid-ordered (Lo) and liquid-disordered (Ld) domains to transmembrane proteins using the lateral pressure concept. In particular we applied previously reported experimental data for the membrane thickness, intrinsic curvature and bending elasticities of coexisting Lo/Ld domains to calculate whether proteins of simple geometric shapes would preferentially diffuse into Lo or Ld domains and form oligomers of a certain size. For the studied lipid mixture we generally found that proteins with convex shapes prefer sorting to Ld phases and the formation of large clusters. Lo domains in turn would be enriched in monomers of concave shaped proteins. We further observed that proteins which are symmetric with respect to the bilayer center prefer symmetric Lo or Ld domains, while asymmetric proteins favor a location in domains with Lo/Ld asymmetry. In the latter case we additionally retrieved a strong dependence on protein directionality, thus providing a mechanism for transmembrane protein orientation.

## 1 Introduction

For several decades lipid-only membranes have served as chemically well-defined mimics of biological membranes enabling detailed physicochemical and biophysical studies of diverse structural and dynamical membrane properties.¹⁻³ One aspect that has ever attracted significant scientific attention is the coupling of membrane properties to protein function.

These interactions can be divided into specific lipid-protein interactions, where lipids interact with either given protein binding sites or grooves,⁴⁻⁷ and unspecific interactions, mediated by the membranes' elastic and structural properties.⁸⁻¹⁵ Furthermore, peripheral membrane proteins may act as scaffolds for the global membrane curvature.¹²,¹⁶ For flat bilayers, hydrophobic matching is one of the most frequently discussed unspecific lipid-protein interactions, relating to the energy needed either to stretch or compress membranes near protein inclusions to alleviate thickness differences with the protein's hydrophobic length.¹¹,¹⁷,¹⁸ Alternatively, a mechanical coupling to the lateral pressure profile¹⁰,¹⁹,²⁰ or stored intrinsic lipid curvatures²¹ has been considered. Importantly, any of the above discussed interactions may affect the protein function through changes

in its conformational equilibrium and/or its preferred partitioning into a given lipid environment.

Here we focus on membrane-mediated protein-sorting into liquid-ordered (Lo) or liquid-disordered (Ld) domains of flat, tension-free bilayers. It is well-established that cholesterol-containing mixtures of high-melting and low-melting lipids display Lo/Ld phase coexistence over a broad range of compositions and temperatures.²²⁻²⁴ These systems serve commonly as models for outer plasma membranes that can be studied by an array of biophysical techniques. For example our laboratory has recently reported detailed *in situ* values for the domains' structural and elastic properties using small-angle X-ray scattering.²⁵⁻²⁷

Differences in domain thickness have been applied to explain protein sorting based on hydrophobic matching.²⁸⁻³⁰ However, it has also been demonstrated that hydrophobic matching cannot be the unique driving force for protein partitioning into Lo or Ld domains. In particular transmembrane peptides designed to match the thickness of either Lo or Ld domains were consistently reported to be primarily localized in Ld domains.³¹⁻³⁴ Moreover, single-membrane-spanning raft proteins were reported to partition into raft-like domains in vesicles prepared from plasma membranes, but not into Lo domains of a ternary lipid mixture.³⁵

Three additional factors can be considered to resolve the disparity with protein partitioning into highly ordered phases: (i) specific interactions with lipid factors such as e.g. raft gangliosides,³⁶ (ii) protein palmitoylation,³⁶,³⁷ or (iii) distinct elastic properties or lipid packing densities of a given domain.³³,³⁴

In this report we consider the latter mechanism motivated by the availability of a theoretical framework and corresponding experimental data. Specifically, the lateral pressure mechanism¹⁹,²⁰,³⁸,³⁹

<figure><img src="image_3.png" /></figure>

<footnote>a University of Graz, Institute of Molecular Biosciences, Biophysics Division, NAWI Graz, Humboldtstr. 50/III, A-8010 Graz, Austria. E-mail: georg.pabst@uni-graz.at; Tel: +43 316 380 4989</footnote>

<footnote>b BioTechMed-Graz, A-8010 Graz, Austria</footnote>

<footnote>c Infineon Technologies Austria AG, Development Center Graz, Babenbergerstr. 10, A-8010 Graz, Austria</footnote>

<footnote>† Present address: Rue de l'Orient 21, Brussels, Belgium.</footnote>
<figure><img src="image_4.png" /></figure>

allowed us to calculate energetic contributions to protein partitioning as a function of overall protein size, shape, and oligomerization state by applying experimental values for Lo/Ld domain properties such as thickness, intrinsic curvature, bending rigidity, and Gaussian modulus of curvature,²⁵⁻²⁷ which are integral parameters of the lateral pressure profile.

We found that convex-shaped proteins generally prefer Ld domains, while concave-shaped proteins would sort into Lo domains. These dependencies are amplified upon increasing protein size. For proteins with cone-like shapes no significant energy gain was found to diffuse from Lo to Ld domains or vice versa. Furthermore, we discuss the effects of lateral pressure differences in Lo and Ld domains on protein oligomerization. Here, pressures favor the aggregation of convex-shaped proteins, while concave proteins would preferentially occur as monomers.

# 2 Methods

## 2.1 Partitioning of single proteins into Lo/Ld domains

The ratio of the occupation probabilities, or molar fractions $X_{1,2}$ of two realizable states 1, 2 in a protein's phase state, which could differ in e.g. protein conformation or lipid environment, is given by the partitioning coefficient $k_p$. In thermal equilibrium, $k_p$ solely depends on the thermal energy $k_B T$ and the difference $\Delta W$ of the states' energy levels $W_{1,2}$ and is given by¹⁹

$$k_p \equiv \frac{X_1}{X_2} = \exp\left(-\frac{W_1 - W_2}{k_B T}\right) = \exp(-\Delta W / k_B T). \quad (1)$$

In what follows, we consider states of different lipid environments, corresponding to Lo and Ld. For convenience we will discuss our results with respect to the natural logarithm of $k_p$, i.e. $-\Delta W/k_B T$. Negative values of this term therefore reflect preferred partitioning into Ld phases and vice versa for positive $\ln k_p$. Furthermore, because $X_{Lo} + X_{Ld} = 1$ we can calculate the equilibrium concentration of proteins in the Lo phase using $X_{Lo} = k_p/(1 + k_p)$.

The transfer energy $\Delta W$ can depend on various contributions, e.g. hydrophobic matching,²⁸ or lateral pressures.¹⁰ Here we focus on the latter mechanism. The lateral pressure profile $p(z)$ is known to emerge from the amphiphilic properties of membrane lipids and the free energy associated with minimizing contact of the apolar regions with the aqueous phase (Fig. 1).⁴⁰,⁴¹ The lateral pressure profile is difficult to determine experimentally.⁴² Thus, either mean-field theories with a lattice model for the hydrocarbon chains¹⁹,²⁰ or molecular dynamics (MD) simulations of diverse kinds have been performed (for review see, e.g. ref. 43). Here we take an alternative approach that allows us to use experimental data.

The energy stored for a protein in a given lateral pressure field can be written as¹⁹,²⁰

$$W_{1,2} = \int_{-d_B/2}^{d_B/2} A(z)p(z)dz, \quad (2)$$

where $d_B$ is the membrane thickness and $A(z)$ the variation of the protein's cross sectional area along the bilayer normal $z$.

<figure><img src="image_5.png" /><figcaption>Fig. 1 Schematic overview of the lateral pressure profile $p(z)$ and its coupling to a membrane protein, where $z$ is the coordinate normal to the bilayer surface. For calculations the complex shape of a membrane protein is transferred into a simple rotationally symmetric body with cross sectional area $A(z)$. The molecular view on the top has been created using the CHARMM-GUI membrane builder for a mixture of distearoyl phosphatidylcholine, dioleoyl phosphatidylcholine and cholesterol in combination with a mechanosensitive channel (PDB File: 2OAR).</figcaption></figure>

Following,³⁸ we can simplify eqn (2) by expanding the protein's cross section into a Taylor series ($A(z) = a_0 + a_1^\pm|z| + a_2^\pm|z|^2 + \dots$) to

$$W_{1,2} = \sum_j a_j^\pm p_j, \quad (3)$$

where $\pm$ refers to the upper or lower monolayer, respectively and $p_j = \int_0^{d_B/2} z^j p(z) dz$ is the $j$-th moment of the pressure profile. The zero'th moment gives the surface tension, which vanishes for flat, tension-free bilayers. The first and second integral moments have been shown to be¹⁰

$$p_1 = J_0 \kappa_C \quad (4)$$

$$p_2 = 2\kappa_C J_0 h - \kappa_G, \quad (5)$$

where $J_0$ is the intrinsic lipid curvature, $\kappa_C$ the monolayer bending rigidity, $h$ the location of the neutral plane with respect to the center of the bilayer and $\kappa_G$ the monolayer Gaussian curvature modulus. $p_1$ is a measure for the lateral torque tension.⁴⁴ All these parameters are experimentally accessible, see e.g. ref. 25–27; for $\kappa_G$ we use the suggested approximation $\kappa_G \approx -0.8\kappa_C$.¹⁰

The cross sectional area of rotationally symmetric proteins $A(z) = \pi r^2(z)$ depends only on its radius $r(z)$, and diverse shapes can be modeled using $r(z) = (r_0 + |z| \tan \varphi^\pm)$, see also Fig. 2. The protein area's Taylor coefficients are then given by

$$a_1^\pm = 2\pi r_0 \tan \varphi^\pm, \quad a_2^\pm = \pi \tan^2 \varphi^\pm. \quad (6)$$

Slightly more complex shapes with smooth contour variations, see e.g. Fig. 1, can be achieved upon free variation of the Taylor coefficients. The energy for partitioning into a given domain is consequently calculated as $\Delta W = \sum_j a_j^\pm \Delta p_j$.
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 2 Parameterization of linear-shaped proteins. Note that cone-shaped proteins are achieved for $\phi^+ = -\phi^-$.</figcaption></figure>

## 2.2 Formation of protein oligomers

A first-order approximation for the aggregation of $n > 2$ proteins can be achieved mathematically by considering dense packing of congruent circles in a circle.⁴⁵ The proteins' maximum radii $r_m$ thus determine the maximum radius $R_m$ of a densely packed aggregate, see Fig. 3. A cluster's radius is then given by $R(z) = R_m - r_m + r(z)$ and its area by $A_n(z) = \pi R^2(z)$. Consistently, substituting $R_0 \equiv R(0)$ for $r_0$ in eqn (6) gives the Taylor coefficients of $A_n(z)$.

To study the influence of lateral pressure on protein clustering, we are considering the changes in the protein area $\Delta A(z) = A_n(z) - nA_1(z)$, which equals for a given $z$ the grey shaded area shown in Fig. 3, while $p(z)$ remains constant. The difference in stored energy between an $n$-mer and $n$ monomers is thus determined by $\Delta W_n = \sum_j \Delta a_j^\pm p_j$, where $\Delta a_j^\pm$ denotes the Taylor coefficients of $\Delta A(z)$. In equilibrium the partitioning coefficient is then defined as³⁹

$$k_{P,n} = \frac{X_n}{X_1^n} = \exp\left(-\frac{\Delta W_n}{k_B T}\right). \qquad (7)$$

Calculation of the protein concentration in an $n$-mer aggregate leads to

$$X_n - k_{P,n}(1 - X_n)^n = 0, \qquad (8)$$

which can be solved numerically.

<figure><img src="image_8.png" /><figcaption>Fig. 3 Modeling protein aggregation by dense packing of congruent circles in a circle. Here $r_m$ is the maximum outer radius of a protein monomer of a given shape and $R_m$ the maximum radius of a densely packed aggregate. $R_m/r_m$ ratios were taken from ref. 45.</figcaption></figure>

## 3 Results

For the present calculations we applied structural data for coexisting Lo/Ld domains in a ternary mixture of dioleoyl phosphatidylcholine (DOPC) distearoyl phosphatidylcholine (DSPC) cholesterol (Chol) reported from X-ray scattering experiments.²⁵⁻²⁷ For completeness data are summarized in Table 1. The differences between Lo and Ld structural and elastic properties, discussed in detail in our previous reports, lead to distinct values for the first and second lateral pressure moments. Most significantly, $p_2$ changes its sign from Ld to Lo, which is mainly due to the more negative intrinsic curvature of the Lo phase and its increased thickness ($h_{Lo} > h_{Ld}$).

In the following we will first present the effects of these differences on the partitioning of protein monomers of different shapes and then discuss contributions of lateral pressures in Lo and Ld to protein aggregation.

### 3.1 Shape-dependence of protein partitioning

Using the parameterization described in Fig. 2, we first calculated the partitioning coefficients for inward and outward bent proteins, varying the opening angle $\phi^+ = \phi^-$ (Fig. 4). Results are either symmetric or anti-symmetric with respect to cylindrically-shaped proteins, which do not exhibit preferred partitioning in Lo or Ld phases, because their shape does not act against the lateral strains stored in the bilayers. In turn, concave-shaped proteins prefer partitioning into Lo domains and convex-shaped proteins into Ld domains, respectively. This can be understood qualitatively in view of the change of the first moment in going from Ld to Lo domains $\Delta p_1 < 0$, which signifies that lateral pressures are redistributed from the lipid/water interface to the bilayer interior, thus favoring a location of inward-bent proteins in the Lo phase.

Cone-shaped proteins ($\phi^+ = -\phi^-$) exhibit symmetric partitioning preferences with respect to $\phi^+ = 0$ (Fig. 4) because terms linear in $\tan \phi^\pm$ compensate due to the bilayer symmetry. Quadratic terms affect a preference of these proteins for partitioning into Lo domains. However, the involved energies are only slightly above thermal energies and consequently rather insignificant.

Concerning size, only convex/concave proteins exhibit a distinct dependence. In particular we found a linear increase for the preference of sorting into either Lo or Ld domains (Fig. 5).

<table><caption>Table 1 Structural and elastic data of DOPC/DSPC/Chol (0.42 : 0.37 : 0.21) at 20 °C and the corresponding lateral pressure profile moments</caption><thead><tr><th></th><th>Ld</th><th>Lo</th></tr></thead><tbody><tr><td>J<sub>0</sub><sup>a</sup> (nm<sup>-1</sup>)</td><td>-0.12 ± 0.01</td><td>-0.20 ± 0.04</td></tr><tr><td>h<sup>b</sup> (nm)</td><td>1.68 ± 0.03</td><td>2.12 ± 0.04</td></tr><tr><td>κ<sub>C</sub><sup>c</sup> (k<sub>B</sub>T)</td><td>5.4 ± 1.2</td><td>14.8 ± 2.5</td></tr><tr><td>κ<sub>G</sub><sup>d</sup> (k<sub>B</sub>T)</td><td>-4.4 ± 1.2</td><td>-11.9 ± 2.7</td></tr><tr><td>p<sub>1</sub> (pN)</td><td>-2.6 ± 0.7</td><td>-12 ± 33</td></tr><tr><td>p<sub>2</sub> (k<sub>B</sub>T)</td><td>2.2 ± 1.2</td><td>-0.7 ± 3.5</td></tr></tbody></table>

a Calculated from ref. 25 using (DOPC/DSPC/Chol)$_{Ld}$ = 0.79 : 0.09 : 0.12 and (DOPC/DSPC/Chol)$_{Lo}$ = 0.05 : 0.65 : 0.30. b Derived from the position of the carbon glycerol groups reported in ref. 26. c Taken from ref. 27.

d Calculated using $\kappa_G/\kappa_C = -0.80 \pm 0.05$.¹⁰
<figure><img src="image_9.png" /><figcaption>Fig. 4 Effect of the opening angle $\phi^+$ on the partitioning of proteins into Lo ($-\Delta W/k_B T > 0$) and Ld ($-\Delta W/k_B T < 0$) domains displayed by a DOPC/DSPC/Chol mixture ($r_0 = 2$ nm). The blue line describes the results for proteins changing their shape from convex to concave forms; the green line for cone-shaped proteins. Shaded areas indicate uncertainties due to experimental errors (Table 1).</figcaption></figure>

Cone-shaped proteins in turn do not change their preferred sorting to Lo domains with size. This is again due to the neutralization of linear terms.

So far we have described proteins with a plane of symmetry in center of the bilayer. For $|\phi^+| \neq |\phi^-|$ we found similar tendencies for Lo/Ld partitioning as for symmetric proteins. This can be generalized in terms of the angle $\alpha$ describing the protein's bending direction (Fig. 2). Proteins with $\alpha < \pi$ prefer sorting to Lo phases, while shapes with $\alpha > \pi$ would diffuse into Ld domains.

## 3.2 Influence of curvature

Using $r(z) = \sqrt{(a_0 + a_1^+|z| + a_2^+|z|^2)/\pi}$ allows the generation of smooth protein contours. In the following we will restrict – due to symmetry – the presentation of our results to the upper half of the protein's shape. Interestingly, proteins with significantly different shapes have equal partitioning energies for constant

<figure><img src="image_10.png" /><figcaption>Fig. 5 Effect of the protein size ($r_0$) on partitioning of proteins into Lo and Ld domains of a DOPC/DSPC/Chol mixture for proteins with concave ($\phi^- = \phi^+ = 0.1$ rad; blue) and convex ($\phi^- = \phi^+ = -0.05$ rad; green) shapes.</figcaption></figure>

$r'(0) = \left( \frac{dr(z)}{dz} \right)_{z=0}$. This can be quantified upon comparison to eqn (6) through the angle

$$ \theta = \arctan \left( \frac{a_1}{2\pi r(0)} \right) = \arctan \left( \frac{a_1}{2\sqrt{a_0}} \right). \quad (9) $$

Since $a_0$ does not contribute to $W$ in tension free bilayers, varying $a_0$, while keeping $a_1$ and $a_2$ constant, yields equal $k_B$ values. Furthermore, even minor variations of $\theta$ induced by small changes of $a_2$, result in comparable partitioning probabilities for Lo domains at significant different protein shapes (Fig. 6).

If $r'(0)$ is allowed to vary significantly, we find strikingly different partitioning coefficients for proteins of similar dimensions. To this end, let us consider proteins under the constraint $a_1 + a_2 \cdot h = \text{const}$. That is, the proteins are tied to the same cross-sections at $z = 0$ and at $z = h$. We observed that concave-shaped proteins, see Fig. 7A, have a decreased (increased) propensity to partition into Lo domains if their contour is inward (outward) bent. Likewise, convex-shaped proteins increase (decrease) their preference for Ld phases for inward (outward) bent contours, see Fig. 7B.

## 3.3 Asymmetric domains

So far we have restricted our analysis to the simple picture of phase-separated, but symmetric membranes. Natural plasma membranes exhibit, however, a considerable degree of lipid asymmetry.46 Most recently protocols have become available, which enable a characterization of asymmetric model membranes with a number of biophysical techniques.47–50 It can be anticipated therefore that experimental data such as those reported in Table 1 will become available for phase separated asymmetric bilayers in due time. In the meantime it is instructive to estimate lateral pressure effects on protein partitioning assuming that Lo and Ld monolayer domains have the same properties as in symmetric bilayers.

Besides considering transbilayer correlation or anti-correlation of Lo and Ld phases, our calculations also included a variation of protein symmetry with respect to $z = 0$. Fig. 8 shows the

<figure><img src="image_11.png" /><figcaption>Fig. 6 Inward-bent proteins with different contours, but similar partitioning probabilities into Lo domains ($a_0 = 4\pi \text{ nm}^2$, $a_1^+ = 1$ nm). The inset shows the results for different $a_2^+$ values.</figcaption></figure>

<figure><img src="image_12.png" /></figure>
<figure><img src="image_13.png" /><figcaption>Fig. 7 Effect of inward or outward bending of the protein's contour for concave (panel A) and convex (panel B) proteins ($r_0$ = 2 nm).</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 8 Stored lateral strain energy for proteins of different shapes in symmetric or asymmetric Lo and Ld domains. Data have been calculated using $r_0$ = 2 nm and ($\varphi^+$; $\varphi^-$) values of shape 1: (0.1; 0.1), shape 2: (0.1; 0.05), shape 3: (0.1; -0.1), shape 4: (0.02; -0.05) and shape 5: (-0.05; -0.05). The partitioning of a protein of given shape in a specific lipid environment is calculated through $-\Delta W/k_B T$ (see (eqn 1)). The overall preferred lipid environment for a given protein shape is given by the lowest $W/k_B T$-value.</figcaption></figure>

stored elastic energies W in the different lipid environments for selected protein shapes. The preferred lipid environment exhibits the lowest stored elastic energy value for a given protein.

Intriguingly, we found that symmetric proteins (shapes 1 and 5 of Fig. 8) would sort either to symmetric Lo (concave-shaped) or Ld (convex-shaped) domains, while asymmetric proteins (shapes 2-4) favor anti-correlated domains. The sorting to anti-correlated domains depends strongly on protein orientation, however. Asymmetric proteins with larger diameters on upper membrane boundary ($r(h) > r(-h)$) prefer sorting to $Ld^{inner}/Lo^{outer}$ domains, while $Lo^{inner}/Ld^{outer}$ is even less energetically favorable than symmetric Lo or Ld domains. Hence, transmembrane proteins may flip horizontally within asymmetric membranes due to the lateral pressure field in order to lower their free energy. Thus, besides commonly considered contributions, such as the overall charge distribution of polar amino acid residues and ionic membrane lipids,⁵¹ also lateral membrane pressures provide a means to orient membrane proteins.

### 3.4 Protein oligomerization

Already early concepts for complex membrane organization considered membrane rafts as platforms for protein assembly.⁵² Recent super-resolution microscopy experiments further indicated

<figure><img src="image_15.png" /><figcaption>Fig. 9 Cluster formation of convex-shaped proteins ($\varphi^+ = -0.05$) in Ld domains (blue) and concave-shaped proteins ($\varphi^- = 0.1$) in Lo domains (red), using $r_0$ = 2 nm. Panel A shows the energetic gain for cluster formation and panel B the relative concentration of aggregates of a given size $n$.</figcaption></figure>

the formation of large membrane protein clusters.⁵³⁻⁵⁵ Here, we consider in a highly simplified way the contributions of membrane lateral pressure to this effect using the methodology described in Section 2.2.

Our previous calculations showed the preference of concave-shaped proteins for Lo and convex-shaped proteins for Ld domains. Hence, we focus on the question whether the proteins would tend to form clusters in their preferred lipid environment or not. In particular we considered the formation of trimers to heptamers (Fig. 3).

Our results demonstrate that the lateral pressure distribution in Lo domains would drive concave proteins towards monomeric forms, while convex proteins in Ld domains would tend to aggregate into large clusters (Fig. 9A). Furthermore, the fraction of clustered proteins is significantly higher for convex proteins at all aggregate sizes (Fig. 9B).

This behavior can be understood in terms of the packing differences between concave and convex proteins. Because the opening angle $\varphi(=\varphi^+ = \varphi^-)$ is assumed to be the same for the monomer and aggregate, $\Delta A(z) = \Delta a_1 \cdot |z| = 2\pi \tan \varphi(R_0 - nr_0) \cdot |z|$ (eqn (6)). Aggregate stability requires $\Delta W_n > 0$. Since $p_1 < 0$ for both considered lipid phases (Table 1), aggregates are stable if $\Delta a_1$ is negative as well, which is equivalent to $R_0/r_0 < n$. Clusters composed of convex protein monomers have the highest packing density at the center of the bilayer, i.e. they are stable if $R_m/r_m < n$. This is always achieved for $n > 2$ (see ratios in Fig. 3). Oligomers of concave proteins in turn have the highest packing density at the lipid/water interface ($\pm d_B/2$) and are consequently more loosely packed at the bilayer midplane. This increases their $R_0/r_0$ ratio, with respect to convex aggregates making the clusters formed of concave proteins less stable. Specifically for the $r_0$ and $\varphi$ values used in Fig. 9 these oligomers are unstable.

## 4 Discussion

We have studied, based on the availability of experimental data, the influence of lateral pressures on the sorting and cluster formation of transmembrane proteins. Several assumptions were made to perform these calculations.

Firstly, all considered proteins were of simple geometric shape with smooth surfaces. For more complex shapes higher

<figure><img src="image_16.png" /></figure>
lateral pressure moments would need to be defined and measured. Alternatively, MD simulations could be applied in combination with crystallographic data for membrane proteins, e.g. ref. 43. However, uncertainties due to limitations of simulation box size and inaccuracies in MD force fields⁵⁶ as well as unknowns of exact protein conformation in a given lipid environment⁵⁷ would be still rather significant. Yet, our calculations using somewhat more complex protein shapes with smoothly curved contours show that trends of preferential partitioning are conserved (Fig. 6 and 7), i.e. outward-bent proteins prefer Ld and inward-bent proteins Lo phases. Thus, the overall tendencies for protein sorting due to lateral pressures are captured by our simplifications. The strength of these tendencies of course depends on the exact shape and would require exact knowledge of protein conformation and p(z).

Secondly, our calculations are based on a single lipid mixture of DOPC/DSPC/Chol. It is highly conceivable that changing lipid composition may influence the here observed tendencies considerably due to changes in *h*, *J*₀, *κ*<sub>C</sub> and *κ*<sub>G</sub>. For example, we found previously that increasing temperature leads to a redistribution of cholesterol from Lo to Ld domains,²⁶ which due to its large negative intrinsic curvature²⁵ will significantly affect the domain's elastic properties. Importantly, present calculations do not consider protein-induced modifications of structural and elastic properties of Lo and Ld phases. Previous studies on single lipid membranes demonstrated that proteins may shift these properties significantly and correlate strongly with protein concentration (see e.g. ref. 13). Distinct effects on Lo and Ld domains are presently unknown, but would warrant further research.

Thirdly, we neglected contributions from specific lipid/protein correlations⁴⁻⁷,³⁶ as well as unspecific interactions such as hydrophobic matching²⁸⁻³⁰ or protein diffusion barriers at the domain boundaries which may act as local sinks for the proteins. The energies involved in these interactions are not trivial to determine. However, for single membrane-spanning peptides hydrophobic matching was not found to contribute to protein partitioning.³⁴ The same authors report, however, that hydrophobic matching affects protein aggregation.

Despite all approximations and limitations discussed in the above paragraphs, the lateral pressure fields stored in Lo and Ld domains provide a fundamental contribution to protein sorting into a given lipid environment, which needs to be considered in a comprehensive picture of lipid/protein interactions in complex membranes. Our results demonstrate the preference of outward-bent proteins to the more loosely packed Ld domains, where they would tend to form clusters of large size. Our calculations did not result in an optimal aggregation size, however (Fig. 9). We speculate that this will be strongly determined by specific protein/protein or lipid/protein interactions. Inward-bent proteins, in turn preferentially locate in the more dense Lo domains in monomeric form if only contributions from lateral pressures are considered.

Qualitatively, a complementary view can be taken by considering the proteins to be rigid bodies, whose shapes affect lipid packing in their vicinity. Dan and Safran²¹ considered the

free energy contributions of interfacial packing mismatches and found them to be dominated by differences in intrinsic lipid curvatures. Assuming that lipids do not demix in the vicinity of the protein inclusion this would mean that Lo phases, having a more negative *J*₀ than Ld phases (Table 1), favor concave proteins and *vice versa* for convex proteins, which agrees with our findings.

It is interesting to compare our results to experimental findings. For example, the multitransmembrane strand protein perfringolysin O (PFO) was found to prefer sorting to Lo domains.³⁰ Although authors have attributed this to hydrophobic matching, we note that PFO forms a multimeric barrel with an overall concave shape, which according to our results favors lateral pressures in Lo domains. The pentameric nicotinic acetylcholine receptor (nAChR) in turn was found to lack preference for Lo domains.⁵⁸ Structural studies suggest a cone-like structure of nAChR's transmembrane domain,⁵⁹ for which shapes our calculations do not yield significant contributions for Lo partitioning. However, lateral pressures in asymmetric Ld<sup>inner</sup>/Lo<sup>outer</sup> domains strongly favor the partitioning of such proteins (Fig. 8). This observation matches with the recent findings by Perillo *et al*.⁶⁰ who suggested the specific interactions of nAChR with outer leaflet sphingomyelin. Thus, partitioning of both proteins could also be rationalized in terms of lateral pressures, although – in view of the many assumptions involved in our calculations – we explicitly refrain from stating that this is the only contribution driving this behavior. Lateral pressures, however, do not contribute to the sorting of single-membrane-spanning proteins due to their nearly cylindrical shape (Fig. 4) and small size (Fig. 5) within the membrane's interior.

For asymmetric proteins in correlated or anti-correlated Lo/Ld domains we found an additional strong coupling to the protein's preferred orientation. Thus, besides influencing protein sorting, lateral pressures also contribute to the direction of transmembrane proteins. Our results consequently allude to the importance of directional membrane-mediated protein sorting and encourage further research along these directions.

## Acknowledgements

We thank Saša Svetina for helpful discussions. This work was supported by the Austrian Science Fund FWF, Project no. P24459-B20.

## References

1. *Handbook of Biological Physics*, ed. R. Lipowsky and E. Sackmann, Elsevier, Amsterdam, 1995.

2. O. G. Mouritsen, *Life – as a matter of fat: The emerging science of lipidomics*, Springer, Berlin, 2005.

3. *Liposomes, Lipid Bilayers and Model Membranes: From Basic Research to Application*, ed. G. Pabst, N. Kučerka, M.-P. Nieh and J. Katsaras, CRC Press, Boca Raton, FL, 2014.

4. A. G. Lee, *Biochim. Biophys. Acta*, 2004, **1666**, 62–87.

5. D. Marsh, *Biochim. Biophys. Acta*, 2008, **1778**, 1545–1575.
<figure><img src="image_17.png" /></figure>

6 Y. D. Paila, S. Tiwari and A. Chattopadhyay, *Biochim. Biophys. Acta*, 2009, **1788**, 295-302.

7 A. W. Smith, *Biochim. Biophys. Acta*, 2012, **1818**, 172-177.

8 S. May, *Curr. Opin. Colloid Interface Sci.*, 2000, **5**, 244-249.

9 J. A. Killian, *FEBS Lett.*, 2003, **555**, 134-138.

10 D. Marsh, *Biophys. J.*, 2007, **93**, 3884-3899.

11 D. Marsh, *Biophys. J.*, 2008, **94**, 3996-4013.

12 T. Baumgart, B. R. Capraro, C. Zhu and S. L. Das, *Annu. Rev. Phys. Chem.*, 2011, **62**, 483-506.

13 G. Pabst, *Advances in Planar Lipid Bilayers and Liposomes*, Academic Press, Burlington, MA, 2013, vol. 18, pp. 81-109.

14 S. Aimon, A. Callan-Jones, A. Berthaud, M. Pinot, G. E. S. Toombes and P. Bassereau, *Dev. Cell*, 2014, **28**, 212-218.

15 B. Božič, S. L. Das and S. Svetina, *Soft Matter*, 2015, **11**, 2479-2487.

16 J. Zimmerberg and M. M. Kozlov, *Nat. Rev. Mol. Cell Biol.*, 2006, **7**, 9-19.

17 N. Dan, P. Pincus and S. A. Safran, *Langmuir*, 1993, **9**, 2768-2771.

18 J. A. Lundbæk, S. A. Collingwood, H. I. Ingólfsson, R. Kapoor and O. S. Andersen, *J. R. Soc., Interface*, 2010, **7**, 373-395.

19 R. S. Cantor, *J. Phys. Chem. B*, 1997, **101**, 1723-1725.

20 R. S. Cantor, *Biophys. J.*, 1999, **76**, 2625-2639.

21 N. Dan and S. A. Safran, *Biophys. J.*, 1998, **75**, 1410-1414.

22 B. Ramstedt and J. P. Slotte, *Biochim. Biophys. Acta*, 2006, **1758**, 1945-1956.

23 D. Marsh, *Biochim. Biophys. Acta*, 2009, **1788**, 2114-2123.

24 D. G. Ackerman and G. W. Feigenson, *Essays Biochem.*, 2015, **57**, 33-42.

25 B. Kollmitzer, P. Heftberger, M. Rappolt and G. Pabst, *Soft Matter*, 2013, **9**, 10877-10884.

26 P. Heftberger, B. Kollmitzer, A. A. Rieder, H. Amenitsch and G. Pabst, *Biophys. J.*, 2015, **108**, 854-862.

27 B. Kollmitzer, P. Heftberger, R. Podgornik, J. F. Nagle and G. Pabst, *Biophys. J.*, 2015, **108**, 2833-2842.

28 J. A. Lundbæk, O. S. Andersen, T. Werge and C. Nielsen, *Biophys. J.*, 2003, **84**, 2080-2089.

29 H.-J. Kaiser, A. Orlowski, T. Róg, T. K. M. Nyholm, W. Chai, T. Feizi, D. Lingwood, I. Vattulainen and K. Simons, *Proc. Natl. Acad. Sci. U. S. A.*, 2011, **108**, 16628-16633.

30 Q. Lin and E. London, *J. Biol. Chem.*, 2013, **288**, 1340-1352.

31 B. Y. van Duyl, D. T. S. Rijkers, B. d. Kruijff and J. A. Killian, *FEBS Lett.*, 2002, **523**, 79-84.

32 M. E. Fastenberg, H. Shogomori, X. Xu, D. A. Brown and E. London, *Biochemistry*, 2003, **42**, 12376-12390.

33 A. Vidal and T. J. McIntosh, *Biophys. J.*, 2005, **89**, 1102-1108.

34 L. V. Schäfer, D. H. d. Jong, A. Holt, A. J. Rzepiela, A. H. d. Vries, B. Poolman, J. A. Killian and S. J. Marrink, *Proc. Natl. Acad. Sci. U. S. A.*, 2011, **108**, 1343-1348.

35 H. J. Kaiser, D. Lingwood, I. Levental, J. L. Sampaio, L. Kalvodova, L. Rajendran and K. Simons, *Proc. Natl. Acad. Sci. U. S. A.*, 2009, **106**, 16645-16650.

36 D. H. d. Jong, C. A. Lopez and S. J. Marrink, *Faraday Discuss.*, 2013, **161**, 347-363.

37 J. Greaves and L. H. Chamberlain, *J. Cell Biol.*, 2007, **176**, 249-254.

38 R. S. Cantor, *Chem. Phys. Lipids*, 1999, **101**, 45-56.

39 R. S. Cantor, *Biophys. J.*, 2002, **82**, 2520-2525.

40 J. N. Israelachvili, *Intermolecular and surface forces*, Academic Press, Burlington, MA, 3rd edn, 2011.

41 A. Ben Shaul, in *Handbook of Biological Physics*, ed. R. Lipowsky and E. Sackmann, Elsevier, Amsterdam, 1995, pp. 359-401.

42 R. H. Templer, S. J. Castle, A. R. Curran, G. Rumbles and D. R. Klug, *Faraday Discuss.*, 1998, **111**, 41-53.

43 O. H. Ollila and I. Vattulainen, in *Molecular Simulations and Biomembranes*, ed. M. S. P. Sansom and P. C. Biggin, RSC Publishing, London, UK, 2010, pp. 26-55.

44 J. M. Seddon and R. H. Templer, in *Structure and dynamics of membranes*, ed. R. Lipowsky and E. Sackmann, North-Holland, Amsterdam, 1995, pp. 97-160.

45 R. L. Graham, B. D. Lubachevsky, K. J. Nurmela and P. R. J. Östergård, *Discrete Math.*, 1998, **181**, 139-154.

46 M. S. Bretscher, *Science*, 1973, **181**, 622-629.

47 S. Pautot, B. J. Frisken and D. A. Weitz, *Proc. Natl. Acad. Sci. U. S. A.*, 2003, **100**, 10718-10721.

48 V. Kiessling, J. M. Crane and L. K. Tamm, *Biophys. J.*, 2006, **91**, 3313-3326.

49 H. T. Cheng and E. London, *Biophys. J.*, 2011, **100**, 2671-2678.

50 D. Marquardt, B. Geier and G. Pabst, *Membranes*, 2015, **5**, 180-196.

51 M. Bogdanov, W. Dowhan and H. Vitrac, *Biochim. Biophys. Acta*, 2014, **1843**, 1475-1488.

52 D. Lingwood and K. Simons, *Science*, 2010, **327**, 46-50.

53 J. J. Sieber, K. I. Willig, C. Kutzner, C. Gerding-Reimers, B. Harke, G. Donnert, B. Rammner, C. Eggeling, S. W. Hell, H. Grubmüller and T. Lang, *Science*, 2007, **317**, 1072-1076.

54 A. Anderluh, E. Klotzsch, A. W. A. F. Reismann, M. Brameshuber, O. Kudlacek, A. H. Newman, H. H. Sitte and G. J. Schütz, *J. Biol. Chem.*, 2014, **289**, 4387-4394.

55 S. K. Saka, A. Honigmann, C. Eggeling, S. W. Hell, T. Lang and S. O. Rizzoli, *Nat. Commun.*, 2014, **5**, 4509.

56 A. Botan, F. Favela-Rosales, P. F. J. Fuchs, M. Javanainen, M. Kanduč, W. Kulig, A. Lamberg, C. Loison, A. Lyubartsev, M. S. Miettinen, L. Monticelli, J. Määttä, O. H. S. Ollila, M. Retegan, T. Róg, H. Santuz and J. Tynkkynen, *J. Phys. Chem. B*, 2015, 15075-15088.

57 P. J. Stansfeld, J. E. Goose, M. Caffrey, E. P. Carpenter, J. L. Parker, S. Newstead and M. S. P. Sansom, *Structure*, 2015, **23**, 1350-1361.

58 V. Bermúdez, S. S. Antollini, G. A. Fernández Nievas, M. I. Aveldaño and F. J. Barrantes, *J. Lipid Res.*, 2010, **51**, 2629-2641.

59 R. J. Hilf and R. Dutzler, *Nature*, 2008, **452**, 375-379.

60 V. L. Perillo, D. A. Peñalva, A. J. Vitale, F. J. Barrantes and S. S. Antollini, *Arch. Biochem. Biophys.*, 2015, **591**, 76-86.