the online version of this article ( doi:10.1007/s00249 - 015 - 1017-x ) contains supplementary material , which is available to authorized users .
the presence of coexisting liquid - ordered and liquid - disordered lipid phases has been implicated in diverse biological processes with high sensitivity to perturbations in membrane composition , tension , curvature , and temperature ( lingwood et al .
lipid phase dynamics may inherently provide nanoscale clustering and facilitate lipid - raft - mediated processes ( veatch et al .
however , trace membrane additives can disturb the lipid phase dynamics and the associated biological processes .
some additives possess rotational asymmetry , preferentially localize to the phase boundary , reduce the phase miscibility , decrease the miscibility transition temperature ( tmis ) , and act as linactants , 1d analogs of surfactants ( trabelsi et al .
other membrane additives are rotationally symmetric , preferentially localize into one particular lipid phase , and may increase or decrease tmis . broadly , the effects of membrane additives are quantified by the change to the lipid phase tmis .
the difference between a sample  temperature and tmis determines fundamental equilibrium qualities such as the correlation length of the coexisting phases ( honerkamp - smith et al .
2008 ; palmieri and safran 2013a ) and the dynamics of the phase mobility ( honerkamp - smith et al . 2012 ; palmieri and safran 2013b ) .
numerous molecules have demonstrated potent abilities to alter membrane phase miscibility without necessarily containing a dependence on the rotation of the molecule within the bilayer , including cholesterol ( dietrich et al .
2001 ; zhao et al . 2007a ; levental et al . 2009 ; heberle et al .
1991 ; gray et al . 2013 ) , insecticides ( jorgensen et al . 1991 ) , n - propyl gallate ( zhao et al .
2007b ) , c - reactive protein ( senz et al . 2010 ) , and flavonoids ( ostroumova et al .
other molecules have been hypothesized to prefer the boundary between two lipid phases because of their molecular structure ; for example , hybrid lipids contain one liquid - order - preferring tail and one liquid - disorder - preferring tail ( brewster et al .
2009 ; hassan - zadeh et al . 2014 ; li and gorfe 2014 ) , and asymmetrically lipidated proteins , such as n - ras and h - ras , have shown preference for partitioning at the lipid phase boundary ( nicolini et al .
further , molecules such as the nonsteroidal anti - inflammatory drug indomethacin may contain both a net phase preference and significant molecular structure asymmetry , although the nanoscale partitioning has yet to be studied ( zhou et al .
previous studies have focused on theoretical predictions and experimental examination of systems with a significant fraction of the membrane consisting of phase - polarized particles ( > 10  mol% ) ( brewster et al . 2009 ; yamamoto et al . 2010 ; yamamoto and safran 2011 ; palmieri and safran 2013a , b ; hassan - zadeh et al .
2014 ) , whereas this manuscript focuses on small composition changes that can greatly affect membrane phase dynamics and provide a means for live cells to efficiently regulate lipid phase - dependent processes .
monte carlo simulations of 2d ising models with conserved order parameters have been used to demonstrate key membrane properties ( frazier et al . 2007 ; yethiraj and weisshaar 2007 ; honerkamp - smith et al .
2011 , 2012 ) . within the ising model , a particle s phase preference (  )
is quantified for liquid - disordered ( black ,   =  1 ) or liquid - ordered ( white ,   =  1 ) phases , although the particular phase associated with the color and sign of  is arbitrary in this symmetric system . the internal energy between two particles ( ji , j ) is equal to j for nearest neighbors and zero otherwise , and it can be calculated from the sum of all particle interactions in the hamiltonian ( h ) according to:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$h = -\sum\limits _ { \langle i , j \rangle } { j_{i , j } \sigma_{i } \sigma_{j } } .$$\end{document}h=-i , jji , jij .
when this model contains 50  % white and 50  % black particles on a square lattice , it exhibits critical behavior with \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{\text{mis } } = 2j/(k_{\text{b } } \ln ( 1 + \sqrt 2 ) ) $ $ \end{document}tmis=2j/(kbln(1 + 2 ) ) , where kb is boltzmann s constant .
isolated plasma membrane vesicles from mast cells happen to be of a composition that demonstrates near - critical behavior with tmis  =  22  c ( veatch et al . 2008 ) , corresponding to a value of j  =  0.26  kcal / mol .
to relate the complex molecular composition of near - critical membranes to the two - state ising model , each pixel from the model represents the mean composition of each state . for example
, the white particles in the ising model represent the cell membrane s average liquid - ordered phase composed of lipids with longer , more saturated acyl tails , greater sphingomyelin content , and slightly higher cholesterol content than the black particles .
this model permits the observation of complex phase behaviors without explicit incorporation of the molecular details or concentrations of the membrane , as long as the average phase segregation preferences are represented by the value of j. in this manuscript , we demonstrate the influence of trace additives in a two - component nearest neighbor model on a square lattice with conserved order parameters as a means of simulating how additives of differing molecular structure could influence lipid phase mixing .
instead of modifying the value of j to represent changes in membrane composition , a third particle type was incorporated in these simulations to provide greater structural detail of the additive s effect on phase separation and additive partitioning . in a ternary lipid membrane , for example
, the additive could be a fourth type of molecule or more of one molecule type that was already present in the membrane .
simulated additives were either rotationally symmetric gray particles or rotational asymmetric phase - polarized particles with gray values of g or polarizations of p , respectively ( fig .  1 ) . generally , the addition of a third type of particle to the system results in a change to the phase diagram , new phase transitions , and variations from critical behavior .
simulations performed with additives of g  =  0 were analogous to the blume - emery - griffiths model ( blume et al .
1971 ) , for which the phase transitions were modeled against an additive concentration with mean - field approximations . however , the additive concentrations explored here are consistent with the composition fluctuations in living and model membranes , including variations in content of cholesterol ( dietrich et al . 2001 ; pokorny et al .
2010 ) , fluorophores ( veatch et al . 2007a ; frazier et al . 2007 ) , and a variety of other biologically active molecules ( jorgensen et al .
2010 ; gray et al . 2013 ; ostroumova et al . 2014 ) . while real  membranes have been demonstrated to be near critical in composition ( veatch et al .
2007b , 2008 ) , they have not been shown to be any closer to critical than the additive - included simulations shown here .
further , the experimental method for observing plasma membrane extracts involved the incorporation of fluorescent lipid additives ( < 0.25  mol% ) ( veatch et al .
2008 ) , which altered the membrane composition and presumably its phase behavior but did not diminish the significance or confidence in the result that the membranes were near critical in composition.fig .  1 
 white , liquid - ordered (
  =  1 ) ; black , liquid - disordered (   =  1 ) ; gray particles (   =  g ) were rotationally symmetric .
phase - polarized additives , or janus additives , were rotationally asymmetric with a polarized phase preference ( ||  =  p ) 
 white , liquid - ordered (   =  1 ) ; black , liquid - disordered (   =  1 ) ; gray particles (   =  g ) were rotationally symmetric .
phase - polarized additives , or janus additives , were rotationally asymmetric with a polarized phase preference ( ||  =  p ) this work focuses on small perturbations to the well - understood two - state ising model and demonstrates how trace additives would affect the characteristic size of the co - existing single - phase domains (  ) and tmis . for 3  mol% additives , tmis and  both changed monotonically with increasing additive concentration , while tmis changed linearly .
rotationally symmetric particles of g  <  1 and rotationally asymmetric particles of all values of p caused a decrease in both tmis and . however , rotationally symmetric additives of g  >  1 resulted in an increase in both tmis and . the morphological detail provided by this model yields the distributions of additives within the system , which is unavailable from mean - field models .
the resulting distributions of the additives at the boundary between phases were quantified , and the applicability of this analysis for assessing experimental membrane additives or cellular perturbations is discussed .
non - local particle exchanges were allowed via a monte carlo algorithm conserving the system composition ( e.g. , 49.5  % white , 49.5  % black , and 1  % additive ) and equilibrated for 10 global equilibration sweeps on a 512    512 square grid with biperiodic boundary conditions .
an additional 10 global sweeps were performed over which an average of the system parameters was reported , e.g. , the correlation length .
the symmetry of this configuration with equal fractions of white and black particles resulted in the additives having the same effect , regardless of the sign of g or p. further details are provided in the supplemental material including a movie of 50 different time points of the equilibrated system.fig . 
2increasing the mole fraction of phase - polarized particles caused an increase in the phase miscibility and a decrease of t 
 mis .
a with p  =  1 , phase - polarized additives are displayed as red squares in these images , and zoomed - in regions show the individual additives in fig .
blue outlined images indicate where t  =  t 
 mis with j set such that t 
 mis of the additive - free system \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$(t_{mis}^{\prime})$$\end{document}(tmis ) is 22  c . increasing
the concentration of phase - polarized particles with p  =  1 resulted in a decrease of b 
  and c 
 t 
 mis .
t 
 mis  =  t 
 mis\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis. uncertainty of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis was approximately 0.6  c , and uncertainty of  was smaller than the symbol unless otherwise indicated increasing the mole fraction of phase - polarized particles caused an increase in the phase miscibility and a decrease of t 
 mis .
a with p  =  1 , phase - polarized additives are displayed as red squares in these images , and zoomed - in regions show the individual additives in fig .
blue outlined images indicate where t  =  t 
 mis with j set such that t 
 mis of the additive - free system \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$(t_{mis}^{\prime})$$\end{document}(tmis ) is 22  c . increasing
the concentration of phase - polarized particles with p  =  1 resulted in a decrease of b 
  and c 
 t 
 mis .
t 
 mis  =  t 
 mis\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis. uncertainty of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis was approximately 0.6  c , and uncertainty of  was smaller than the symbol unless otherwise indicated the correlation length of the phases was calculated from the azimuthal average of the 2d spatial correlation of the system configuration ( i ) where all white particles were valued at 1 , black particles at 0 , and additives at ( g/2  +  0.5 ) .
unless otherwise stated , if p    0 , then g  =  0 , and vice versa . the 2d correlation function ( c ) as a function of the position \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\left ( { \bar{r } } \right)$$\end{document}r can be calculated by2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c\left ( { \bar { r } } \right ) = \frac{{\left\langle { i\left ( { \overline{r } } \right)i\left ( { \overline{r } + \overline{r } } \right ) } \right\rangle } } { { \left\langle { i\left ( { \overline{r } } \right ) } \right\rangle^{2 } } } , \,$$\end{document}cr=irir+rir2,where  \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\langle \ , \rangle$$\end{document}  represents the average over all 2d space \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( { \bar { r } } ) $ $ \end{document}(r ) .
computationally , this was expedited by the use of fast fourier transforms ( ffts ) with3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c\left ( { \bar { r } } \right ) = \frac{{{\text{fft}}^ { - 1 } \left ( { \left| { { \text{fft}}\left ( i \right ) } \right|^{2 } } \right)}}{n},$$\end{document}cr=fft-1ffti2n , and normalization ( n ) ( veatch et al . 2012 ) .
the 1d correlation function ( c ) was computed from4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c\left ( r \right ) = \frac{1}{2\pi } \int\limits_{0}^{2\pi } { c\left ( { r,\theta } \right){\text{d}}\theta } , $ $ \end{document}cr=1202cr,d,and the correlation length (  ) was calculated by nonlinear least square fitting of cfit to c(r ) where5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c_{\text{fit } } = c_{0 } \frac{{{\text{e}}^ { - r/\xi } } } { { r^{\eta } } } $ $ \end{document}cfit = c0e - r/rand the critical exponent  equals 0.25 , which is exact for the ising model .
the critical exponents for the three - component systems and finite systems are not equal to that of the ising model ; however , additionally allowing  to change did not significantly change fitting results . changes to tmis
tmis was determined by measuring the shift in  as a function of temperature and by ascertaining the maximum of the specific heat versus temperature ( veatch et al .
the finite size effects that alter the perceived transition temperature were consistent for all of these simulations on the same size system , which permitted this simplistic calculation of tmis and its broad applicability to systems of other sizes .
a demonstration of the independence of these results on system size is presented in the supplemental material ( fig .
changes in phase miscibility with an increasing additive fraction are visibly obvious in images of the system configuration ( figs .  2 , s5 , s6 ) .
tmis decreased linearly with increasing fractions of phase - polarized particles , or janus particles , as predicted from mean - field models ( yethiraj and weisshaar 2007 ; brewster et al .
the miscibility temperature is reduced for any phase polarization ( p ) , while rotationally symmetric gray particles either increased or decreased tmis depending on their gray value ( g ) ( figs .  3 , s6 ) . increasing p for the phase - polarized additives resulted in a greater decrease in tmis until p  >  3 , upon which the additives segregated into a third phase , separate from the white and black particles , and reduced the white and black particle mixing .
gray particles with g  <  1 encouraged phase miscibility by decreasing the differences between the ordered and disordered phases , analogous to the incorporation of n - alcohols into membranes ( gray et al . 2013 )
. however , when g  >  1 , the opposite effect dominates ; the difference between the phases was increased when g  >  1 and tmis was increased , analogous to the addition of -lysin to membranes ( pokorny et al .
gray particles with g  >  1 were no longer of in - between phase and were more liquid - order preferring than white particles , i.e. ,  super-white.fig .  3 
 a 
  and b 
 t 
 mis and either increased or decreased upon additive addition depending on the additive properties .
t 
 mis  =  t 
 mis\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis  =  22  c .
measurement of  is limited by the system size , resulting in large uncertainties for situations with   >  100 pixels ( fig .
t 
 mis and  changes are symmetrical for negative values of g or p 
 
 a 
  and b 
 t 
 mis and either increased or decreased upon additive addition depending on the additive properties .
t 
 mis  =  t 
 mis\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis  =  22  c .
measurement of  is limited by the system size , resulting in large uncertainties for situations with   >  100 pixels ( fig .
t 
 mis and  changes are symmetrical for negative values of g or p 
  when p or g  > 
3 , the additives condensed and formed a distinct phase excluding white or black particles ( fig .
s6 ) , and the white - black miscibility temperature became less affected by the additives ( fig . 
3 ) . for p or g  >  3 , simulations were unable to capture a diverse set of configurations for the additive aggregation because of the improbable rearrangement of the condensed phases , analogous to the kinetically trapped configurations that are commonly experienced experimentally .
aggregates of phase - polarized particles at p  =  4 demonstrated rotational ordering of the additives ( fig .
a particle was declared to be at the phase boundary if it was adjacent to two white and two black particles .
the ratio of additives at the phase boundary to all additives with no adjacent additives was calculated ( fig . 
as p increased , the fraction of phase - polarized particles at the phase boundary increased ; as g increased , the fraction of gray particles at the phase boundary decreased .
lowering the temperature would have made phase - polarized particles more likely to be found at the phase boundary if the phase morphology had stayed consistent .
however , increasing temperature from 8.5 to 36  c resulted in 2.3    0.4 times more phase boundary ( fig .
s8 ) , which was dependent on both simulation size and additive properties , and generally resulted in a higher probability for the additives being located at the phase boundary for higher temperatures .
additive partitioning was normalized to the amount of phase boundary in order to calculate the fold enhancement of the additive to the phase boundary relative to the white or black particles ( fig . 
additives more concentrated at the phase boundary than the g  =  1 particles preferentially partitioned at the phase boundary , whereas additives less concentrated at the phase boundary than the g  =  1 particles preferentially partitioned away from the phase boundary , which is consistent with the sign of tmis ( fig .  3 ) .
although there is not a local energy difference for the g  =  0 additive if it is immersed in a single phase or at the phase boundary , there is an energy difference for the system that encourages g  <  1 additives to be concentrated at the phase boundary . at colder temperatures , the phase - polarized additives especially concentrated at the phase boundary because of the decreased entropic drive for mixing.fig .  4 
 a the probability that an additive was found at the phase boundary varied with the additive properties and temperature .
the fraction of additives at the boundary was calculated as the ratio of the number of additives surrounded by two white and two black particles versus the total number of additives , excluding additives adjacent to other additives .
b the additive locations were normalized by the amount of phase boundary in each condition to yield the fold enhancement of the additive at the phase boundary versus generic white particles .
uncertainty was calculated as the standard deviation of different time points and smaller than the symbol unless otherwise indicated 
 a the probability that an additive was found at the phase boundary varied with the additive properties and temperature .
the fraction of additives at the boundary was calculated as the ratio of the number of additives surrounded by two white and two black particles versus the total number of additives , excluding additives adjacent to other additives .
b the additive locations were normalized by the amount of phase boundary in each condition to yield the fold enhancement of the additive at the phase boundary versus generic white particles .
uncertainty was calculated as the standard deviation of different time points and smaller than the symbol unless otherwise indicated these results apply to additives in systems that are near critical compositions .
it is important to note that these results would vary greatly for compositions that are far from critical compositions because of the increased phase boundary line tension and the increased likelihood for boundary - preferring particles to localize at the boundary in such systems . in all the prior described simulations , only additives with
however , real membrane additives would likely possess both a net phase preference ( g    0 ) and a rotational phase asymmetry ( p    0 ) . to predict how changing both g and p for
a membrane additive would affect the miscibility transition temperature of the membrane , simulations were performed for additives with rotational asymmetry and a net phase preference at 1  mol% ( fig . 
the changes to the miscibility transition temperature were dependent on both the g and p value of the additive in a cumulative effect . for small values of g and p ( i.e. , either g or p    1 ) , tmis as a function of g and p was approximately equal to the sum of the change in miscibility temperature from 1  mol% gray particles of gray value g plus the change from 1  mol% polarized particles of polarization p ( fig .  5b).fig . 
5 
 a changes to the miscibility transition temperature were induced by the incorporation of 1  mol% additives with both net and rotationally asymmetric phase preferences , i.e. , dependent on both g and p values , as shown in the inset .
b for small values of g and p , t 
 mis as a function of g and p , t 
 mis ( g , p ) , is approximately equal to t 
 mis ( g , 0 )  +  t 
 mis ( 0 , p ) 
 a changes to the miscibility transition temperature were induced by the incorporation of 1  mol% additives with both net and rotationally asymmetric phase preferences , i.e. , dependent on both g and p values , as shown in the inset .
b for small values of g and p , t 
 mis as a function of g and p , t 
 mis ( g , p ) , is approximately equal to t 
 mis ( g , 0 )  +  t 
 mis ( 0 , p )
taken together , these idealized membrane additives represent different means by which molecules may alter the lipid phases in a membrane .
a membrane - associated molecule may be experimentally quantified with both a p value and g value to describe the effects of the additive molecular structure on lipid phase dynamics in a particular membrane .
g is analogous to the partition constant and quantifies the molecule s preference for one phase over the other .
p quantifies the molecule s phase asymmetry and preference to partition at the phase boundary .
the effects of a particular membrane additive are determined entirely by the membrane s additive - free miscibility temperature ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis ) or internal energy between two particles ( j ) of the near - critical membrane without a dependence on the particular molecular details of the individual membrane components .
the temperatures reported throughout this article could be alternately represented as reduced temperatures ( tr ) , where tr  =  ( t     \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis)/ \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis. upon rescaling , these results would be directly applicable to other membranes with different values of j to predict tmis for different \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis. further , experimental use of membranes of varying \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis would permit the exploration of different tr values within experimentally limited temperature ranges .
for example , increasing the fatty acid tail length in the saturated , liquid - order preferring lipid of a ternary mixture typically increases \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis and j , while decreasing the additive s g value and enhancing the additive s phase boundary preference .
similarly , the p and g values for a membrane additive in a complex cellular membrane could be predicted by measuring the additive 's effects on synthetic model membranes of greatly different molecular composition , but similar \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$t_{mis}^{\prime}$$\end{document}tmis and j. values of g and p that are greater than one represent the phase preferences for the additives that are stronger that then constituents of the phases themselves .
for example , if dipalmitoylphosphatidylcholine ( dppc ) was added to a dioleoylphosphatidylcholine ( dopc)/dimyristoylphosphatidylcholine ( dmpc)/cholesterol bilayer , it would likely have a g  > 
1 since dppc has a stronger preference for the liquid - ordered phase than dopc , dmpc , or cholesterol . similarly , it is feasible that a particular phase - preferring component of a hybrid molecule could have a stronger phase preference than average constituents of that phase , i.e. , p  > 
, the palmitoyl group of h - ras or palmitoyloleoylphosphatidylcholine ( popc ) could prefer the liquid - ordered phase in a dopc / dmpc / cholesterol bilayer even more than dopc , dmpc , or cholesterol .
the quantification of p and g values of cholesterol , for example , requires knowledge of the effects cholesterol addition has on lipid phase separation and the likelihood of cholesterol to be found on the phase boundary .
cholesterol preferentially partitions into liquid - ordered phases , thereby indicating a value of g  >  0 .
the addition of cholesterol to a membrane typically encourages phase miscibility , indicating that cholesterol has a value of g  <  1 . the value of g for cholesterol will vary with the surrounding membrane lipids ; however , quantification of cholesterol s effects on the tmis for near - critical membranes composed of dopc , dppc , and cholesterol displayed approximately 15  c reduction in tmis with the addition of 5  mol% cholesterol ( veatch et al .
2007b ) , which implies a g value of 0.3 for cholesterol in this system ( i.e. , fig . 
the focused quantification of changes in tmis with variations in  cholesterol concentration in near - critical ternary systems would provide greater certainty in the determination of g and p. similar analysis on the addition of the fluorescent lipid dii - c12 into membranes composed of dopc , dppc , and cholesterol at molar ratios of 35:35:30 demonstrates how a more complex mechanism may also be present ( veatch et al .
in this system , tmis increased linearly with increasing dii - c12 concentrations up to 0.1  mol% with reduced effects on tmis at higher dii - c12 concentrations .
this difference may be due to a particularly drastic effect of dii - c12 pushing the membrane toward or away from a critical composition ( veatch et al .
2007a )  or specific interactions between dii - c12 and the other lipids in the system .
the theoretical framework presented here provides a basis for the interpretation of experimental data for the lipid phase and boundary preference of membrane - associated molecules .
similarly , this framework provides a means of predicting the lipid - raft association and lipid phase - dependent clustering in complex systems through the extension of results from model membrane studies . this ability to quantify the phase preference of membrane additives may prove to be a valuable parameter for the therapeutic development and evaluation of the therapeutic mechanisms . 

supplementary material 1 ( avi 12307  kb ) supplementary material 2 ( docx 19959  kb )
 the addition of trace molecules into membranes can significantly alter the morphology of the co - existing liquid phases and lipid phase transition temperature . 
 membrane additives may affect lipid phase dynamics through preferentially partitioning to the boundary between lipid phases or preferentially mixing into one lipid phase . 
 the characteristic differences between these mechanisms are demonstrated here in a minimalistic nearest neighbor model to provide a framework for how slight changes to membrane composition may affect lipid - phase - dependent processes , such as lipid - raft formation , immunological signaling , and molecular sorting preceding endocytosis with coexisting liquid phases . within the low mole fractions 
 explored here ( 3  mol% ) , increasing the additive concentration linearly changed the phase miscibility temperature . 
 rotationally asymmetric janus particles reduced the miscibility transition temperature for all fractions and degree of phase polarization . 
 rotationally symmetric additives , however , either increased or decreased the phase miscibility temperature depending on the phase preference of the additive . while most experimental molecules may contain aspects of both of these idealized additives , this model provides a broad framework to quantify the effects of membrane additives in regard to lipid phase preference , lipid - raft association , and contribution to lipid phase - dependent molecular sorting.electronic supplementary materialthe online version of this article ( doi:10.1007/s00249 - 015 - 1017-x ) contains supplementary material , which is available to authorized users . 
Electronic supplementary material
Introduction
Methods
Results
Discussion
Electronic supplementary material