Heterogeneous nanoscopic lipid diffusion in the live cell membrane and its dependency on cholesterol

Cholesterol plays a unique role in the regulation of membrane organization and dynamics by modulating the membrane phase transition at the nanoscale. Unfortunately, due to their small sizes and dynamic nature, the effects of cholesterol-mediated membrane nanodomains on membrane dynamics remain elusive. Here, using ultrahigh-speed single-molecule tracking with advanced optical microscope techniques, we investigate the diffusive motion of single phospholipids in the live cell plasma membrane at the nanoscale and its dependency on the cholesterol concentration. We find that both saturated and unsaturated phospholipids undergo anomalous subdiffusion on the length scale of 10–100 nm. The diffusion characteristics exhibit considerable variations in space and in time, indicating that the nanoscopic lipid diffusion is highly heterogeneous. Importantly, through the statistical analysis, apparent dual-mobility subdiffusion is observed from the mixed diffusion behaviors. The measured subdiffusion agrees well with the hop diffusion model that represents a diffuser moving in a compartmentalized membrane created by the cytoskeleton meshwork. Cholesterol depletion diminishes the lipid mobility with an apparently smaller compartment size and a stronger confinement strength. Similar results are measured with temperature reduction, suggesting that the more heterogeneous and restricted diffusion is connected to the nanoscopic membrane phase transition. Our conclusion supports the model that cholesterol depletion induces the formation of gel-phase, solid-like membrane nanodomains. These nanodomains undergo restricted diffusion and act as diffusion obstacles to the membrane molecules that are excluded from the nanodomains. This work provides the experimental evidence that the nanoscopic lipid diffusion in the cell plasma membrane is heterogeneous and sensitive to the cholesterol concentration and temperature, shedding new light on the regulation mechanisms of nanoscopic membrane dynamics.     

PEG molecular weight and lateral diffusion of PEG-ylated lipids in magnetically aligned bicelles

Lateral diffusion coefficients of PEG-ylated lipids with three different molecular weight PEG groups (1000, 2000 and 5000) were measured in magnetically-aligned bicelles using the stimulated echo (STE) pulsed field gradient (PEG) (1)H nuclear magnetic resonance (NMR) method. At concentrations below the PEG "mushroom-to-brush" transition, all three PEG-ylated lipids exhibited unrestricted lateral diffusion, with lateral diffusion coefficients comparable to those of corresponding non-PEG-ylated lipids and independent of PEG molecular weight. At concentrations above this transition, lateral diffusion slowed progressively with increasing concentration of PEG-ylated lipid as a result of surface crowding. As well, the lateral diffusion coefficients exhibited a pronounced decrease with increasing PEG group molecular weight and a diffusion-time dependence indicative of obstructed diffusion. We conclude that, while lateral diffusion of PEG-ylated lipids within lipid bilayers is determined primarily by the hydrophobic anchoring group, when crowding at the lipid bilayer surface becomes significant, the size of the extra-membranous domain, in this case the PEG group, can influence lateral diffusion, leading to decreased diffusivity with increasing size and producing obstructed diffusion at high crowding. These findings imply that similar considerations will pertain to lateral diffusion of membrane proteins with large extra-membranous domains.     

Lateral diffusion in an archipelago. The effect of mobile obstacles.

Lateral diffusion of mobile proteins and lipids (tracers) in a membrane is hindered by the presence of proteins (obstacles) in the membrane. If the obstacles are immobile, their effect may be described by percolation theory, which states that the long-range diffusion constant of the tracers goes to zero when the area fraction of obstacles is greater than the percolation threshold. If the obstacles are themselves mobile, the diffusion constant of the tracers depends on the area fraction of obstacles and the relative jump rate of tracers and obstacles. This paper presents Monte Carlo calculations of diffusion constants on square and triangular lattices as a function of the concentration of obstacles and the relative jump rate. The diffusion constant for particles of various sizes is also obtained. Calculated values of the concentration-dependent diffusion constant are compared with observed values for gramicidin and bacteriorhodopsin. The effect of the proteins as inert obstacles is significant, but other factors, such as protein-protein interactions and perturbation of lipid viscosity by proteins, are of comparable importance. Potential applications include the diffusion of proteins at high concentrations (such as rhodopsin in rod outer segments), the modulation of diffusion by release of membrane proteins from cytoskeletal attachment, and the diffusion of mobile redox carriers in mitochondria, chloroplasts, and endoplasmic reticulum.     

Obstructed diffusion in phase-separated supported lipid bilayers: a combined atomic force microscopy and fluorescence recovery after photobleaching approach

Proteins and other macromolecules are believed to hinder molecular lateral diffusion in cellular membranes. We have constructed a well-characterized model system to better understand how obstacles in lipid bilayers obstruct diffusion. Fluorescence recovery after photobleaching was used to measure the lateral diffusion coefficient in single supported bilayers composed of mixtures of 1,2-dilauroylphosphotidylcholine (DLPC) and 1,2-distearoylphosphotidylcholine (DSPC). Because these lipids are immiscible and phase separate at room temperature, a novel quenching technique allowed us to construct fluid DLPC bilayers containing small disk-shaped gel-phase DSPC domains that acted as obstacles to lateral diffusion. Our experimental setup enabled us to analyze the same samples with atomic force microscopy and exactly characterize the size, shape, and number of gel-phase domains before measuring the obstacle-dependent diffusion coefficient. Lateral obstructed diffusion was found to be dependent on obstacle area fraction, size, and geometry. Analysis of our results using a free area diffusion model shows the possibility of unexpected long-range ordering of fluid-phase lipids around the gel-phase obstacles. This lipid ordering has implications for lipid-mediated protein interactions in cellular membranes.     

PEG molecular weight and lateral diffusion of PEG-ylated lipids in magnetically aligned bicelles

Lateral diffusion coefficients of PEG-ylated lipids with three different molecular weight PEG groups (1000, 2000 and 5000) were measured in magnetically-aligned bicelles using the stimulated echo (STE) pulsed field gradient (PEG) ¹H nuclear magnetic resonance (NMR) method. At concentrations below the PEG “mushroom-to-brush” transition, all three PEG-ylated lipids exhibited unrestricted lateral diffusion, with lateral diffusion coefficients comparable to those of corresponding non-PEG-ylated lipids and independent of PEG molecular weight. At concentrations above this transition, lateral diffusion slowed progressively with increasing concentration of PEG-ylated lipid as a result of surface crowding. As well, the lateral diffusion coefficients exhibited a pronounced decrease with increasing PEG group molecular weight and a diffusion-time dependence indicative of obstructed diffusion. We conclude that, while lateral diffusion of PEG-ylated lipids within lipid bilayers is determined primarily by the hydrophobic anchoring group, when crowding at the lipid bilayer surface becomes significant, the size of the extra-membranous domain, in this case the PEG group, can influence lateral diffusion, leading to decreased diffusivity with increasing size and producing obstructed diffusion at high crowding. These findings imply that similar considerations will pertain to lateral diffusion of membrane proteins with large extra-membranous domains.     

A milling crowd model for local and long-range obstructed lateral diffusion. Mobility of excimeric probes in the membrane of intact erythrocytes.

A new model for lateral diffusion, the milling crowd model (MC), is proposed and is used to derive the dependence of the monomeric and excimeric fluorescence yields of excimeric membrane probes on their concentration. According to the MC model, probes migrate by performing spatial exchanges with a randomly chosen nearest neighbor (lipid or probe). Only nearest neighbor probes, one of which is in the excited state, may form an excimer. The exchange frequency, and hence the local lateral diffusion coefficient, may then be determined from experiment with the aid of computer simulation of the excimer formation kinetics. The same model is also used to study the long-range lateral diffusion coefficient of probes in the presence of obstacles (e.g., membrane proteins). The dependence of the monomeric and excimeric fluorescence yields of 1-pyrene-dodecanoic acid probes on their concentration in the membranes of intact erythrocytes was measured and compared with the prediction of the MC model. The analysis yields an excimer formation rate for nearest neighbor molecules of approximately 1 X 10(7) s-1 and an exchange frequency of approximately greater than 2 X 10(7) s-1, corresponding to a local diffusion coefficient of greater than 3 X 10(-8) cm2 s-1. This value is several times larger than the long-range diffusion coefficient for a similar system measured in fluorescence photobleaching recovery experiments. The difference is explained by the fact that long-range diffusion is obstructed by dispersed membrane proteins and is therefore greatly reduced when compared to free diffusion. The dependence of the diffusion coefficient on the fractional area covered by obstacles and on their size is derived from MC simulations and is compared to those of other theories lateral diffusibility.     

Anomalous diffusion due to binding: a Monte Carlo study.

In classical diffusion, the mean-square displacement increases linearly with time. But in the presence of obstacles or binding sites, anomalous diffusion may occur, in which the mean-square displacement is proportional to a nonintegral power of time for some or all times. Anomalous diffusion is discussed for various models of binding, including an obstruction/binding model in which immobile membrane proteins are represented by obstacles that bind diffusing particles in nearest-neighbor sites. The classification of binding models is considered, including the distinction between valley and mountain models and the distinction between singular and nonsingular distributions of binding energies. Anomalous diffusion is sensitive to the initial conditions of the measurement. In valley models, diffusion is anomalous if the diffusing particles start at random positions but normal if the particles start at thermal equilibrium positions. Thermal equilibration leads to normal diffusion, or to diffusion as normal as the obstacles allow.     

Accelerating membrane insertion of peripheral proteins with a novel membrane mimetic model

Characterizing atomic details of membrane binding of peripheral membrane proteins by molecular dynamics (MD) has been significantly hindered by the slow dynamics of membrane reorganization associated with the phenomena. To expedite lateral diffusion of lipid molecules without sacrificing the atomic details of such interactions, we have developed a novel membrane representation, to our knowledge, termed the highly mobile membrane-mimetic (HMMM) model to study binding and insertion of various molecular species into the membrane. The HMMM model takes advantage of an organic solvent layer to represent the hydrophobic core of the membrane and short-tailed phospholipids for the headgroup region. We demonstrate that using these components, bilayer structures are formed spontaneously and rapidly, regardless of the initial position and orientation of the lipids. In the HMMM membrane, lipid molecules exhibit one to two orders of magnitude enhancement in lateral diffusion. At the same time, the membrane atomic density profile of the headgroup region produced by the HMMM model is essentially identical to those obtained for full-membrane models, indicating the faithful representation of the membrane surface by the model. We demonstrate the efficiency of the model in capturing spontaneous binding and insertion of peripheral proteins by using the membrane anchor (γ-carboxyglutamic-acid-rich domain; GLA domain) of human coagulation factor VII as a test model. Achieving full insertion of the GLA domain consistently in 10 independent unbiased simulations within short simulation times clearly indicates the robustness of the HMMM model in capturing membrane association of peripheral proteins very efficiently and reproducibly. The HMMM model will provide significant improvements to the current all-atom models by accelerating lipid dynamics to examine protein-membrane interactions more efficiently.     

Phospholipids undergo hop diffusion in compartmentalized cell membrane

The diffusion rate of lipids in the cell membrane is reduced by a factor of 5-100 from that in artificial bilayers. This slowing mechanism has puzzled cell biologists for the last 25 yr. Here we address this issue by studying the movement of unsaturated phospholipids in rat kidney fibroblasts at the single molecule level at the temporal resolution of 25 micros. The cell membrane was found to be compartmentalized: phospholipids are confined within 230-nm-diameter (phi) compartments for 11 ms on average before hopping to adjacent compartments. These 230-nm compartments exist within greater 750-nm-phi compartments where these phospholipids are confined for 0.33 s on average. The diffusion rate within 230-nm compartments is 5.4 microm2/s, which is nearly as fast as that in large unilamellar vesicles, indicating that the diffusion in the cell membrane is reduced not because diffusion per se is slow, but because the cell membrane is compartmentalized with regard to lateral diffusion of phospholipids. Such compartmentalization depends on the actin-based membrane skeleton, but not on the extracellular matrix, extracellular domains of membrane proteins, or cholesterol-enriched rafts. We propose that various transmembrane proteins anchored to the actin-based membrane skeleton meshwork act as rows of pickets that temporarily confine phospholipids.     

Lateral diffusion in an archipelago. Dependence on tracer size.

In a pure fluid-phase lipid, the dependence of the lateral diffusion coefficient on the size of the diffusing particle may be obtained from the Saffman-Delbrück equation or the free-volume model. When diffusion is obstructed by immobile proteins or domains of gel-phase lipids, the obstacles yield an additional contribution to the size dependence. Here this contribution is examined using Monte Carlo calculations. For random point and hexagonal obstacles, the diffusion coefficient depends strongly on the size of the diffusing particle, but for fractal obstacles--cluster-cluster aggregates and multicenter diffusion-limited aggregates--the diffusion coefficient is independent of the size of the diffusing particle. The reason is that fractals have no characteristic length scale, so a tracer sees on average the same obstructions, regardless of its size. The fractal geometry of the excluded area for tracers of various sizes is examined. Percolation thresholds are evaluated for a variety of obstacles to determine how the threshold depends on tracer size and to compare the thresholds for compact and extended obstacles.     

Sources of anomalous diffusion on cell membranes: a Monte Carlo study

A stochastic random walk model of protein molecule diffusion on a cell membrane was used to investigate the fundamental causes of anomalous diffusion in two-dimensional biological media. Three different interactions were considered: collisions with fixed obstacles, picket fence posts, and capture by, or exclusion from, lipid rafts. If motion is impeded by randomly placed, fixed obstacles, we find that diffusion can be highly anomalous, in agreement with previous studies. In contrast, collision with picket fence posts has a negligible effect on the anomalous exponent at realistic picket fence parameters. The effects of lipid rafts are more complex. If proteins partition into lipid rafts there is a small to moderate effect on the anomalous exponent, whereas if proteins are excluded from rafts there is a large effect on the anomalous exponent. In combination, these mechanisms can explain the level of anomaly in experimentally observed membrane diffusion, suggesting that anomalous diffusion is caused by multiple mechanisms whose effects are approximately additive. Finally, we show that the long-range diffusion rate, D(macro), estimated from fluorescence recovery after photobleaching studies, can be much smaller than D(micro), the small-scale diffusion rate, and is highly sensitive to obstacle densities and other impeding structures.     

Two-Dimensional Continuum Percolation Threshold for Diffusing Particles of Nonzero Radius

Lateral diffusion in the plasma membrane is obstructed by proteins bound to the cytoskeleton. The most important parameter describing obstructed diffusion is the percolation threshold. The thresholds are well known for point tracers, but for tracers of nonzero radius, the threshold depends on the excluded area, not just the obstacle concentration. Here thresholds are obtained for circular obstacles on the continuum. Random obstacle configurations are generated by Brownian dynamics or Monte Carlo methods, the obstacles are immobilized, and the percolation threshold is obtained by solving a bond percolation problem on the Voronoi diagram of the obstacles. The percolation threshold is expressed as the diameter of the largest tracer that can cross a set of immobile obstacles at a prescribed number density. For random overlapping obstacles, the results agree with the known analytical solution quantitatively. When the obstacles are soft disks with a 1/r¹² repulsion, the percolating diameter is ∼20% lower than for overlapping obstacles. A percolation model predicts that the threshold is highly sensitive to the tracer radius. To our knowledge, such a strong dependence has so far not been reported for the plasma membrane, suggesting that percolation is not the factor controlling lateral diffusion. A definitive experiment is proposed.     

Supported double membranes

Planar model membranes, like supported lipid bilayers and surface-tethered vesicles, have been proven to be useful tools for the investigation of complex biological functions in a significantly less complex membrane environment. In this study, we introduce a supported double membrane system that should be useful for studies that target biological processes in the proximity of two lipid bilayers such as the periplasm of bacteria and mitochondria or the small cleft between pre- and postsynaptic neuronal membranes. Large unilamellar vesicles (LUV) were tethered to a preformed supported bilayer by a biotin–streptavidin tether. We show from single particle tracking (SPT) experiments that these vesicle are mobile above the plane of the supported membrane. At higher concentrations, the tethered vesicles fuse to form a second continuous bilayer on top of the supported bilayer. The distance between the two bilayers was determined by fluorescence interference contrast (FLIC) microscopy to be between 16 and 24 nm. The lateral diffusion of labeled lipids in the second bilayer was very similar to that in supported membranes. SPT experiments with reconstituted syntaxin-1A show that the mobility of transmembrane proteins was not improved when compared with solid supported membranes.     

The rate of formation of surface membrane ether lipids in Ehrlich ascites tumor cells: kinetic considerations

A previous investigation has shown that O-alkyl phospholipids are present in the surface membrane of Ehrlich ascites tumor cells. In the present investigation it was shown that 90% or more of [1-3H]hexadecanol injected intraperitoneally into mice bearing Ehrlich ascites tumors is taken up by the neoplastic cells in less than 15 min. Near maximum formation of surface membrane O-alkyl phospholipids requires approximately 8 h. The rate of accumulation of O-alkyl phospholipids is very similar both for the whole cell and for the surface membrane. Further examination of the data revealed that the conversion of hexadecanol into O-alkyl glycerophospholipids can be described by a simple model in which O-alkyl lipids appear at a single rate constant of 0.25 to 0.35 per hour and disappear at a rate of 0.02 per hour or less. These rate constants were obtained initially by stochastic analysis and validated both by deterministic methods and by compartmental analysis using the SAAM computer program. The method of kinetic analysis described may find broader application in providing comparative rate constants for the in vivo turnover of O-alkyl lipids in both normal and neoplastic tissues. The advantage of a stochastic approach is that kinetic data may be obtained with fewer assumptions relating to pool structure or specific models.     

Transbilayer (flip-flop) lipid motion and lipid scrambling in membranes

This paper reviews the current knowledge on the various mechanisms for transbilayer, or flip-flop, lipid motion in model and cell membranes, enzyme-assisted lipid transfer by flippases, floppases and scramblases is briefly discussed, while non-catalyzed lipid flip-flop is reviewed in more detail. Transbilayer lipid motion may occur as a result of the insertion of foreign molecules (detergents, lipids, or even proteins) in one of the membrane leaflets. It may also be the result of the enzymatic generation of lipids, e.g. diacylglycerol or ceramide, at one side of the membrane. Transbilayer motion rates decrease in the order diacylglycerol ? ceramide ? phospholipids. Ceramide, but not diacylglycerol, can induce transbilayer motion of other lipids, and bilayer scrambling. Transbilayer lipid diffusion and bilayer scrambling are defined as two conceptually and mechanistically different processes. The mechanism of scrambling appears to be related to local instabilities caused by the non-lamellar ceramide molecule, or by other molecules that exhibit a relatively slow flip-flop rate, when asymmetrically inserted or generated in one of the monolayers in a cell or model membrane.     

STED nanoscopy reveals molecular details of cholesterol- and cytoskeleton-modulated lipid interactions in living cells

Details about molecular membrane dynamics in living cells, such as lipid-protein interactions, are often hidden from the observer because of the limited spatial resolution of conventional far-field optical microscopy. The superior spatial resolution of stimulated emission depletion (STED) nanoscopy can provide new insights into this process. The application of fluorescence correlation spectroscopy (FCS) in focal spots continuously tuned down to 30 nm in diameter distinguishes between free and anomalous molecular diffusion due to, for example, transient binding of lipids to other membrane constituents, such as lipids and proteins. We compared STED-FCS data recorded on various fluorescent lipid analogs in the plasma membrane of living mammalian cells. Our results demonstrate details about the observed transient formation of molecular complexes. The diffusion characteristics of phosphoglycerolipids without hydroxyl-containing headgroups revealed weak interactions. The strongest interactions were observed with sphingolipid analogs, which showed cholesterol-assisted and cytoskeleton-dependent binding. The hydroxyl-containing headgroup of gangliosides, galactosylceramide, and phosphoinositol assisted binding, but in a much less cholesterol- and cytoskeleton-dependent manner. The observed anomalous diffusion indicates lipid-specific transient hydrogen bonding to other membrane molecules, such as proteins, and points to a distinct connectivity of the various lipids to other membrane constituents. This strong interaction is different from that responsible for forming cholesterol-dependent, liquid-ordered domains in model membranes.     

Formation of Raft-Like Assemblies within Clusters of Influenza Hemagglutinin Observed by MD Simulations

The association of hemagglutinin (HA) with lipid rafts in the plasma membrane is an important feature of the assembly process of influenza virus A. Lipid rafts are thought to be small, fluctuating patches of membrane enriched in saturated phospholipids, sphingolipids, cholesterol and certain types of protein. However, raft-associating transmembrane (TM) proteins generally partition into Ld domains in model membranes, which are enriched in unsaturated lipids and depleted in saturated lipids and cholesterol. The reason for this apparent disparity in behavior is unclear, but model membranes differ from the plasma membrane in a number of ways. In particular, the higher protein concentration in the plasma membrane may influence the partitioning of membrane proteins for rafts. To investigate the effect of high local protein concentration, we have conducted coarse-grained molecular dynamics (CG MD) simulations of HA clusters in domain-forming bilayers. During the simulations, we observed a continuous increase in the proportion of raft-type lipids (saturated phospholipids and cholesterol) within the area of membrane spanned by the protein cluster. Lateral diffusion of unsaturated lipids was significantly attenuated within the cluster, while saturated lipids were relatively unaffected. On this basis, we suggest a possible explanation for the change in lipid distribution, namely that steric crowding by the slow-diffusing proteins increases the chemical potential for unsaturated lipids within the cluster region. We therefore suggest that a local aggregation of HA can be sufficient to drive association of the protein with raft-type lipids. This may also represent a general mechanism for the targeting of TM proteins to rafts in the plasma membrane, which is of functional importance in a wide range of cellular processes.     

The dynamical Matryoshka model: 2. Modeling of local lipid dynamics at the sub-nanosecond timescale in phospholipid membranes

Biological membranes are generally formed by lipids and proteins. Often, the membrane properties are studied through model membranes formed by phospholipids only. They are molecules composed by a hydrophilic head group and hydrophobic tails, which can present a panoply of various motions, including small localized movements of a few atoms up to the diffusion of the whole lipid or collective motions of many of them. In the past, efforts were made to measure these motions experimentally by incoherent neutron scattering and to quantify them, but with upcoming modern neutron sources and instruments, such models can now be improved. In the present work, we expose a quantitative and exhaustive study of lipid dynamics on DMPC and DMPG membranes, using the Matryoshka model recently developed by our group. The model is confronted here to experimental data collected on two different membrane samples, at three temperatures and two instruments. Despite such complexity, the model describes reliably the data and permits to extract a series of parameters. The results compare also very well to other values found in the literature.     

Methods to measure the lateral diffusion of membrane lipids and proteins

In this chapter, we discuss methods to measure lateral mobility of membrane lipids and proteins using techniques based on the light microscope. These methods typically sample lateral mobility in very small, micron-sized regions of the membrane so that they can be used to measure diffusion in regions of single cells. The methods are based on fluorescence from the molecules of interest or from light scattered from particles attached to single or small groups of membrane lipids or proteins. Fluorescence recovery after photobleaching (FRAP), fluorescence correlation spectroscopy (FCS) and Single particle tracking (SPT) are presented in that order. FRAP and FCS methodologies are described for a dedicated wide field microscope although many confocal microscopes now have software permitting these measurement to be made; nevertheless, the principles of the measurement are the same for a wide field or confocal microscope. SPT can be applied to trace the movements of single fluorescent molecules in membranes but this aspect will not be treated in detail.     

Diffusion measurements of water, ubiquinone and lipid bilayer inside a cylindrical nanoporous support: A stimulated echo pulsed-field gradient MAS-NMR investigation

Stimulated echo pulsed-field gradient ¹H magic angle spinning NMR has been used to investigate the mobility of water, ubiquinone and tethered phospholipids, components of a biomimetic model membrane. The diffusion constant of water corresponds to an isotropic motion in a cylinder. When the lipid bilayer is obtained after the fusion of small unilamellar vesicles, the extracted value of lipid diffusion indicates unrestricted motion. The cylindrical arrangement of the lipids permits a simplification of data analysis since the normal bilayer is perpendicular to the gradient axis. This feature leads to a linear relation between the logarithm of the attenuation of the signal intensity and a factor depending on the gradient strength, for lipids covering the inner wall of aluminium oxide nanopores as well as for lipids adsorbed on a polymer sheet rolled into a cylinder. The effect of the bilayer formation on water diffusion has also been observed. The lateral diffusion coefficient of ubiquinone is in the same order of magnitude as the lipid lateral diffusion coefficient, in agreement with its localization within the bilayer.