Molecular motion in cell membranes: analytic study of fence-hindered random walks

Particle-environment interactions in arbitrary dimensions: A unifying analytic framework to model diffusion with inert spatial heterogeneities

Inert interactions between randomly moving entities and spatial disorder play a crucial role in quantifying the diffusive properties of a system, with examples ranging from molecules advancing along dendritic spines to antipredator displacements of animals due to sparse vegetation. Despite the ubiquity of such phenomena, a general framework to model the movement explicitly in the presence of spatial heterogeneities is missing. Here, we tackle this challenge and develop an analytic theory to model inert particle-environment interactions in domains of arbitrary shape and dimensions. We use a discrete space formulation, which allows us to model the interactions between an agent and the environment as perturbed dynamics between lattice sites. Interactions from spatial disorder, such as impenetrable and permeable obstacles or regions of increased or decreased diffusivity, as well as many others, can be modelled using our framework. We provide exact expressions for the generating function of the occupation probability of the diffusing particle and related transport quantities such as first-passage, return, and exit probabilities and their respective means. We uncover a surprising property, the disorder indifference phenomenon of the mean first-passage time in the presence of a permeable barrier in quasi-1D systems. We demonstrate the widespread applicability of our formalism by considering three examples that span across scales and disciplines. (1) We explore an enhancement strategy of transdermal drug delivery. (2) We represent the movement decisions of an animal undergoing thigomotaxis, the tendency to remain at the peripheries of its enclosure, using a spatially disordered environment. (3) We illustrate the use of spatial heterogeneities to model inert interactions between particles by modeling the search for a promoter region on the DNA by transcription factors during gene transcription.

Renewal equations for single-particle diffusion through a semipermeable interface

Diffusion through semipermeable interfaces has a wide range of applications, ranging from molecular transport through biological membranes to reverse osmosis for water purification using artificial membranes. At the single-particle level, one-dimensional diffusion through a barrier with constant permeability κ_{0} can be modeled in terms of so-called snapping out Brownian motion (BM). The latter sews together successive rounds of partially reflected BMs that are restricted to either the left or right of the barrier. Each round is killed (absorbed) at the barrier when its Brownian local time exceeds an exponential random variable parameterized by κ_{0}. A new round is then immediately started in either direction with equal probability. It has recently been shown that the probability density for snapping out BM satisfies a renewal equation that relates the full density to the probability densities of partially reflected BM on either side of the barrier. Moreover, generalized versions of the renewal equation can be constructed that incorporate non-Markovian, encounter-based models of absorption. In this paper we extend the renewal theory of snapping out BM to single-particle diffusion in bounded domains and higher spatial dimensions. In each case we show how the solution of the renewal equation satisfies the classical diffusion equation with a permeable boundary condition at the interface. That is, the probability flux across the interface is continuous and proportional to the difference in densities on either side of the interface. We also consider an example of an asymmetric interface in which the directional switching after each absorption event is biased. Finally, we show how to incorporate an encounter-based model of absorption for single-particle diffusion through a spherically symmetric interface. We find that, even when the same non-Markovian model of absorption applies on either side of the interface, the resulting permeability is an asymmetric time-dependent function with memory. Moreover, the permeability functions tend to be heavy tailed.

A probabilistic model of diffusion through a semi-permeable barrier

Diffusion through semi-permeable structures arises in a wide range of processes in the physical and life sciences. Examples at the microscopic level range from artificial membranes for reverse osmosis to lipid bilayers regulating molecular transport in biological cells to chemical and electrical gap junctions. There are also macroscopic analogues such as animal migration in heterogeneous landscapes. It has recently been shown that one-dimensional diffusion through a barrier with constant permeability κ 0 is equivalent to snapping out Brownian motion (BM). The latter sews together successive rounds of partially reflecting BMs that are restricted to either the left or the right of the barrier. Each round is killed when its Brownian local time exceeds an exponential random variable parameterized by κ 0 . A new round is then immediately started in either direction with equal probability. In this article, we use a combination of renewal theory, Laplace transforms and Green’s function methods to show how an extended version of snapping out BM provides a general probabilistic framework for modelling diffusion through a semi-permeable barrier. This includes modifications of the diffusion process away from the barrier (e.g. stochastic resetting) and non-Markovian models of membrane absorption that kill each round of partially reflected BM. The latter leads to time-dependent permeabilities.

Numerical calculation on a two-step subdiffusion behavior of lateral protein movement in plasma membranes

A two-step subdiffusion behavior of lateral movement of transmembrane proteins in plasma membranes has been observed by using single-molecule experiments. A nested double-compartment model where large compartments are divided into several smaller ones has been proposed in order to explain this observation. These compartments are considered to be delimited by membrane-skeleton "fences" and membrane-protein "pickets" bound to the fences. We perform numerical simulations of a master equation using a simple two-dimensional lattice model to investigate the heterogeneous diffusion dynamics behavior of transmembrane proteins within plasma membranes. We show that the experimentally observed two-step subdiffusion process can be described using fence and picket models combined with decreased local diffusivity of transmembrane proteins in the vicinity of the pickets. This allows us to explain the two-step subdiffusion behavior without explicitly introducing nested double compartments.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

We review recent advances in multiscale modeling of the biomechanical characteristics of red blood cells (RBCs) in hematological diseases, and their relevance to the structure and dynamics of defective RBCs. We highlight examples of successful simulations of blood disorders including malaria and other hereditary disorders, such as sickle-cell anemia, spherocytosis, and elliptocytosis.

Modeling of band-3 protein diffusion in the normal and defective red blood cell membrane

We employ a two-component red blood cell (RBC) membrane model to simulate lateral diffusion of band-3 proteins in the normal RBC and in the RBC with defective membrane proteins. The defects reduce the connectivity between the lipid bilayer and the membrane skeleton (vertical connectivity), or the connectivity of the membrane skeleton itself (horizontal connectivity), and are associated with the blood disorders of hereditary spherocytosis (HS) and hereditary elliptocytosis (HE) respectively. Initially, we demonstrate that the cytoskeleton limits band-3 lateral mobility by measuring the band-3 macroscopic diffusion coefficients in the normal RBC membrane and in a lipid bilayer without the cytoskeleton. Then, we study band-3 diffusion in the defective RBC membrane and quantify the relation between band-3 diffusion coefficients and percentage of protein defects in HE RBCs. In addition, we illustrate that at low spectrin network connectivity (horizontal connectivity) band-3 subdiffusion can be approximated as anomalous diffusion, while at high horizontal connectivity band-3 diffusion is characterized as confined diffusion. Our simulations show that the band-3 anomalous diffusion exponent depends on the percentage of protein defects in the membrane cytoskeleton. We also confirm that the introduction of attraction between the lipid bilayer and the spectrin network reduces band-3 diffusion, but we show that this reduction is lower than predicted by the percolation theory. Furthermore, we predict that the attractive force between the spectrin filament and the lipid bilayer is at least 20 times smaller than the binding forces at band-3 and glycophorin C, the two major membrane binding sites. Finally, we explore diffusion of band-3 particles in the RBC membrane with defects related to vertical connectivity. We demonstrate that in this case band-3 diffusion can be approximated as confined diffusion for all attraction levels between the spectrin network and the lipid bilayer. By comparing the diffusion coefficients measured in horizontal vs. vertical defects, we conclude that band-3 mobility is primarily controlled by the horizontal connectivity.

Consequences of animal interactions on their dynamics: emergence of home ranges and territoriality

Animal spacing has important implications for population abundance, species demography and the environment. Mechanisms underlying spatial segregation have their roots in the characteristics of the animals, their mutual interaction and their response, collective as well as individual, to environmental variables. This review describes how the combination of these factors shapes the patterns we observe and presents a practical, usable framework for the analysis of movement data in confined spaces. The basis of the framework is the theory of interacting random walks and the mathematical description of out-of-equilibrium systems. Although our focus is on modelling and interpreting animal home ranges and territories in vertebrates, we believe further studies on invertebrates may also help to answer questions and resolve unanswered puzzles that are still inaccessible to experimental investigation in vertebrate species.

Dynamics of the serine chemoreceptor in the Escherichia coli inner membrane: a high-speed single-molecule tracking study

We investigated the mobility of the polar localized serine chemoreceptor, Tsr, labeled by the fluorescent protein Venus in the inner membrane of live Escherichia coli cells at observation rates up to 1000 Hz. A fraction (7%) of all Tsr molecules shows free diffusion over the entire cell surface with an average diffusion coefficient of 0.40 ± 0.01 μm(2) s(-1). The remaining molecules were found to be ultimately confined in compartments of size 290 ± 15 nm and showed restricted diffusion at an inner barrier found at 170 ± 10 nm. At the shortest length-scales (<170 nm), all Tsr molecules diffuse equally. Disruption of the cytoskeleton and rounding of the cells resulted in an increase in the mobile fraction of Tsr molecules and a fragmenting of the previously polar cluster of Tsr consistent with a curvature-based mechanism of Tsr cluster maintenance.

Lateral diffusion in a discrete fluid membrane with immobile particles

Due to the coupling between the plasma membrane and the actin cytoskeleton, membrane molecules such as receptor proteins can become immobilized by binding to cytoskeletal structures. We investigate the effect of immobile membrane molecules on the diffusion of mobile ones by modeling the membrane as a two-dimensional (2D) fluid composed of hard particles and performing event-driven molecular dynamics simulations at a particle density where the system is in an isotropic liquid state. We show that the diffusion coefficient sharply decreases with increasing immobile fraction, dropping by a factor of ∼ 3 as the fraction of immobile particles increases from 0 to 0.1, in a system-size dependent manner. By combining our results with earlier calculations, we estimate that a factor-of-∼ 20 reduction in diffusion coefficients in live cell membranes, a puzzling finding in cell biology, can be accounted for when less than ∼ 22% of the particles in our model system is immobilized. Furthermore, we investigate the effects of confinement induced by a correlated distribution of immobile particles by calculating the distribution of the time it takes for particles to escape from a corral. In the regime where the particles can always escape from the corral, it is found that the escape times follow an exponential distribution, and the mean escape time grows exponentially with the density of obstacles at the corral boundary, increasing by a factor of 3-5 when immobile particles cover 50% of the boundary, and is approximately proportional to the area of the corral. We believe that our findings will be useful in interpreting (1) single molecule observations of membrane molecules and (2) results of particle based simulations that explore the effect of fluid dynamics on molecular transport in a 2D fluid.

Analytic solutions for some reaction-diffusion scenarios

Motivated currently by the problem of coalescence of receptor clusters in mast cells in the general subject of immune reactions, and formerly by the investigation of exciton trapping and sensitized luminescence in molecular systems and aggregates, we present analytic expressions for survival probabilities of moving entities undergoing diffusion and reaction on encounter. Results we provide cover several novel situations in simple 1-d systems as well as higher-dimensional counterparts along with a useful compendium of such expressions in chemical physics and allied fields. We also emphasize the importance of the relationship of discrete sink term analysis to continuum boundary condition studies.

Insulin regulates Glut4 confinement in plasma membrane clusters in adipose cells

Insulin-stimulated delivery of glucose transporter-4 (GLUT4) to the plasma membrane (PM) is the hallmark of glucose metabolism. In this study we examined insulin’s effects on GLUT4 organization in PM of adipose cells by direct microscopic observation of single monomers tagged with photoswitchable fluorescent protein. In the basal state, after exocytotic delivery only a fraction of GLUT4 is dispersed into the PM as monomers, while most of the GLUT4 stays at the site of fusion and forms elongated clusters (60–240 nm). GLUT4 monomers outside clusters diffuse freely and do not aggregate with other monomers. In contrast, GLUT4 molecule collision with an existing cluster can lead to immediate confinement and association with that cluster. Insulin has three effects: it shifts the fraction of dispersed GLUT4 upon delivery, it augments the dissociation of GLUT4 monomers from clusters ∼3-fold and it decreases the rate of endocytic uptake. All together these three effects of insulin shift most of the PM GLUT4 from clustered to dispersed states. GLUT4 confinement in clusters represents a novel kinetic mechanism for insulin regulation of glucose homeostasis.

Dynamic organizing principles of the plasma membrane that regulate signal transduction: commemorating the fortieth anniversary of Singer and Nicolson's fluid-mosaic model

The recent rapid accumulation of knowledge on the dynamics and structure of the plasma membrane has prompted major modifications of the textbook fluid-mosaic model. However, because the new data have been obtained in a variety of research contexts using various biological paradigms, the impact of the critical conceptual modifications on biomedical research and development has been limited. In this review, we try to synthesize our current biological, chemical, and physical knowledge about the plasma membrane to provide new fundamental organizing principles of this structure that underlie every molecular mechanism that realizes its functions. Special attention is paid to signal transduction function and the dynamic aspect of the organizing principles. We propose that the cooperative action of the hierarchical three-tiered mesoscale (2-300 nm) domains--actin-membrane-skeleton induced compartments (40-300 nm), raft domains (2-20 nm), and dynamic protein complex domains (3-10 nm)--is critical for membrane function and distinguishes the plasma membrane from a classical Singer-Nicolson-type model.

Bacterial secretion and the role of diffusive and subdiffusive first passage processes

By funneling protein effectors through needle complexes located on the cellular membrane, bacteria are able to infect host cells during type III secretion events. The spatio-temporal mechanisms through which these events occur are however not fully understood, due in part to the inherent challenges in tracking single molecules moving within an intracellular medium. As a result, theoretical predictions of secretion times are still lacking. Here we provide a model that quantifies, depending on the transport characteristics within bacterial cytoplasm, the amount of time for a protein effector to reach either of the available needle complexes. Using parameters from Shigella flexneri we are able to test the role that translocators might have to activate the needle complexes and offer semi-quantitative explanations of recent experimental observations.

Insights into cell membrane microdomain organization from live cell single particle tracking of the IgE high affinity receptor FcϵRI of mast cells

Current models propose that the plasma membrane of animal cells is composed of heterogeneous and dynamic microdomains known variously as cytoskeletal corrals, lipid rafts and protein islands. Much of the experimental evidence for these membrane compartments is indirect. Recently, live cell single particle tracking studies using quantum dot-labeled IgE bound to its high affinity receptor FcϵRI, provided direct evidence for the confinement of receptors within micrometer-scale cytoskeletal corrals. In this study, we show that an innovative time-series analysis of single particle tracking data for the high affinity IgE receptor, FcϵRI, on mast cells provides substantial quantitative information about the submicrometer organization of the membrane. The analysis focuses on the probability distribution function of the lengths of the jumps in the positions of the quantum dots labeling individual IgE FcϵRI complexes between frames in movies of their motion. Our results demonstrate the presence, within the micrometer-scale cytoskeletal corrals, of smaller subdomains that provide an additional level of receptor confinement. There is no characteristic size for these subdomains; their size varies smoothly from a few tens of nanometers to a over a hundred nanometers. In QD-IGE labeled unstimulated cells, jumps of less than 70 nm predominate over longer jumps. Addition of multivalent antigen to crosslink the QD-IgE-FcϵRI complexes causes a rapid slowing of receptor motion followed by a long tail of mostly jumps less than 70 nm. The reduced receptor mobility likely reflects both the membrane heterogeneity revealed by the confined motion of the monomeric receptor complexes and the antigen-induced cross linking of these complexes into dimers and higher oligomers. In both cases, the probability distribution of the jump lengths is well fit, from 10 nm to over 100 nm, by a novel power law. The fit for short jumps suggests that the motion of the quantum dots can be modeled as diffusion in a fractal space of dimension less than two.

Confining domains lead to reaction bursts: reaction kinetics in the plasma membrane

Confinement of molecules in specific small volumes and areas within a cell is likely to be a general strategy that is developed during evolution for regulating the interactions and functions of biomolecules. The cellular plasma membrane, which is the outermost membrane that surrounds the entire cell, was considered to be a continuous two-dimensional liquid, but it is becoming clear that it consists of numerous nano-meso-scale domains with various lifetimes, such as raft domains and cytoskeleton-induced compartments, and membrane molecules are dynamically trapped in these domains. In this article, we give a theoretical account on the effects of molecular confinement on reversible bimolecular reactions in a partitioned surface such as the plasma membrane. By performing simulations based on a lattice-based model of diffusion and reaction, we found that in the presence of membrane partitioning, bimolecular reactions that occur in each compartment proceed in bursts during which the reaction rate is sharply and briefly increased even though the asymptotic reaction rate remains the same. We characterized the time between reaction bursts and the burst amplitude as a function of the model parameters, and discussed the biological significance of the reaction bursts in the presence of strong inhibitor activity.

Membrane mechanisms for signal transduction: the coupling of the meso-scale raft domains to membrane-skeleton-induced compartments and dynamic protein complexes

Virtually all biological membranes on earth share the basic structure of a two-dimensional liquid. Such universality and peculiarity are comparable to those of the double helical structure of DNA, strongly suggesting the possibility that the fundamental mechanisms for the various functions of the plasma membrane could essentially be understood by a set of simple organizing principles, developed during the course of evolution. As an initial effort toward the development of such understanding, in this review, we present the concept of the cooperative action of the hierarchical three-tiered meso-scale (2-300 nm) domains in the plasma membrane: (1) actin membrane-skeleton-induced compartments (40-300 nm), (2) raft domains (2-20 nm), and (3) dynamic protein complex domains (3-10nm). Special attention is paid to the concept of meso-scale domains, where both thermal fluctuations and weak cooperativity play critical roles, and the coupling of the raft domains to the membrane-skeleton-induced compartments as well as dynamic protein complexes. The three-tiered meso-domain architecture of the plasma membrane provides an excellent perspective for understanding the membrane mechanisms of signal transduction.

Hierarchical organization of the plasma membrane: investigations by single-molecule tracking vs. fluorescence correlation spectroscopy

Single-molecule tracking and fluorescence correlation spectroscopy (FCS) applied to the plasma membrane in living cells have allowed a number of unprecedented observations, thus fostering a new basic understanding of molecular diffusion, interaction, and signal transduction in the plasma membrane. It is becoming clear that the plasma membrane is a heterogeneous entity, containing diverse structures on nano-meso-scales (2-200 nm) with a variety of lifetimes, where certain membrane molecules stay together for limited durations. Molecular interactions occur in the time-dependent inhomogeneous two-dimensional liquid of the plasma membrane, which might be a key for plasma membrane functions.

Effect of heterogeneity on the characterization of cell membrane compartments: I. Uniform size and permeability

Observations of the motion of individual molecules in the membrane of a number of different cell types have led to the suggestion that the outer membrane of many eukaryotic cells may be effectively partitioned into microdomains. A major cause of this suggested partitioning is believed to be due to the direct/indirect association of the cytosolic face of the cell membrane with the cortical cytoskeleton. Such intimate association is thought to introduce effective hydrodynamic barriers into the membrane that are capable of frustrating molecular Brownian motion over distance scales greater than the average size of the compartment. To date, the standard analytical method for deducing compartment characteristics has relied on observing the random walk behavior of a labeled lipid or protein at various temporal frequencies and different total lengths of time. Simple theoretical arguments suggest that the presence of restrictive barriers imparts a characteristic turnover to a plot of mean squared displacement versus sampling period that can be interpreted to yield the average dimensions of the compartment expressed as the respective side lengths of a rectangle. In the following series of articles, we used computer simulation methods to investigate how well the conventional analytical strategy coped with heterogeneity in size, shape, and barrier permeability of the cell membrane compartments. We also explored questions relating to the necessary extent of sampling required (with regard to both the recorded time of a single trajectory and the number of trajectories included in the measurement bin) for faithful representation of the actual distribution of compartment sizes found using the SPT technique. In the current investigation, we turned our attention to the analytical characterization of diffusion through cell membrane compartments having both a uniform size and permeability. For this ideal case, we found that (i) an optimum sampling time interval existed for the analysis and (ii) the total length of time for which a trajectory was recorded was a key factor.

Computer simulations of protein diffusion in compartmentalized cell membranes

The diffusion of proteins in the cell membrane is investigated using computer simulations of a two-dimensional model. The membrane is assumed to be divided into compartments, with adjacent compartments separated by a barrier of stationary obstacles. Each compartment contains traps represented by stationary attractive disks. Depending on their size, these traps are intended to model either smaller compartments or binding sites. The simulations are intended to model the double-compartment model, which has been used to interpret single molecule experiments in normal rat kidney cells, where five regimes of transport are observed. The simulations show, however, that five regimes are observed only when there is a large separation between the sizes of the traps and large compartments, casting doubt on the double compartment model for the membrane. The diffusive behavior is sensitive to the concentration and size of traps and the strength of the barrier between compartments suggesting that the diffusion of proteins can be effectively used to characterize the structure of the membrane.

Time series analysis of particle tracking data for molecular motion on the cell membrane

Biophysicists use single particle tracking (SPT) methods to probe the dynamic behavior of individual proteins and lipids in cell membranes. The mean squared displacement (MSD) has proven to be a powerful tool for analyzing the data and drawing conclusions about membrane organization, including features like lipid rafts, protein islands, and confinement zones defined by cytoskeletal barriers. Here, we implement time series analysis as a new analytic tool to analyze further the motion of membrane proteins. The experimental data track the motion of 40 nm gold particles bound to Class I major histocompatibility complex (MHCI) molecules on the membranes of mouse hepatoma cells. Our first novel result is that the tracks are significantly autocorrelated. Because of this, we developed linear autoregressive models to elucidate the autocorrelations. Estimates of the signal to noise ratio for the models show that the autocorrelated part of the motion is significant. Next, we fit the probability distributions of jump sizes with four different models. The first model is a general Weibull distribution that shows that the motion is characterized by an excess of short jumps as compared to a normal random walk. We also fit the data with a chi distribution which provides a natural estimate of the dimension d of the space in which a random walk is occurring. For the biological data, the estimates satisfy 1 < d < 2, implying that particle motion is not confined to a line, but also does not occur freely in the plane. The dimension gives a quantitative estimate of the amount of nanometer scale obstruction met by a diffusing molecule. We introduce a new distribution and use the generalized extreme value distribution to show that the biological data also have an excess of long jumps as compared to normal diffusion. These fits provide novel estimates of the microscopic diffusion constant. Previous MSD analyses of SPT data have provided evidence for nanometer-scale confinement zones that restrict lateral diffusion, supporting the notion that plasma membrane organization is highly structured. Our demonstration that membrane protein motion is autocorrelated and is characterized by an excess of both short and long jumps reinforces the concept that the membrane environment is heterogeneous and dynamic. Autocorrelation analysis and modeling of the jump distributions are powerful new techniques for the analysis of SPT data and the development of more refined models of membrane organization. The time series analysis also provides several methods of estimating the diffusion constant in addition to the constant provided by the mean squared displacement. The mean squared displacement for most of the biological data shows a power law behavior rather the linear behavior of Brownian motion. In this case, we introduce the notion of an instantaneous diffusion constant. All of the diffusion constants show a strong consistency for most of the biological data.

Diffusion in a fluid membrane with a flexible cortical cytoskeleton

We calculate the influence of a flexible network of long-chain proteins, which is anchored to a fluid membrane, on protein diffusion in this membrane. This is a model for the cortical cytoskeleton and the lipid bilayer of the red blood cell, which we apply to predict the influence of the cytoskeleton on the diffusion coefficient of a mobile band 3 protein. Using the pressure field that the cytoskeleton exerts on the membrane, from the steric repulsion between the diffusing protein and the cytoskeletal filaments, we define a potential landscape for the diffusion within the bilayer. We study the changes to the diffusion coefficient on removal of one type of anchor proteins, e.g., in several hemolytic anemias, as well as for isotropic and anisotropic stretching of the cytoskeleton. We predict an overall increase of the diffusion for a smaller number of anchor proteins and increased diffusion for anisotropic stretching in the direction of the stretch, because of the decrease in the spatial frequency as well as in the height of the potential barriers.