Regulation of protein mobility in cell membranes: a dynamic corral model

Multi-physics interactions drive VEGFR2 relocation on endothelial cells

Vascular Endothelial Growth Factor Receptor-2 (VEGFR2) is a pro-angiogenic receptor, expressed on endothelial cells (ECs). Although biochemical pathways that follow the VEGFR2 activation are well established, knowledge about the dynamics of receptors on the plasma membrane remains limited. Ligand stimulation induces the polarization of ECs and the relocation of VEGFR2, either in cell protrusions or in the basal aspect in cells plated on ligand-enriched extracellular matrix (ECM). We develop a mathematical model in order to simulate the relocation of VEGFR2 on the cell membrane during the mechanical adhesion of cells onto a ligand-enriched substrate. Co-designing the in vitro experiments with the simulations allows identifying three phases of the receptor dynamics, which are controlled respectively by the high chemical reaction rate, by the mechanical deformation rate, and by the diffusion of free receptors on the membrane. The identification of the laws that regulate receptor polarization opens new perspectives toward developing innovative anti-angiogenic strategies through the modulation of EC activation.

Membrane Compartmentalization Reducing the Mobility of Lipids and Proteins within a Model Plasma Membrane

The cytoskeleton underlying cell membranes may influence the dynamic organization of proteins and lipids within the bilayer by immobilizing certain transmembrane (TM) proteins and forming corrals within the membrane. Here, we present coarse-grained resolution simulations of a biologically realistic membrane model of asymmetrically organized lipids and TM proteins. We determine the effects of a model of cytoskeletal immobilization of selected membrane proteins using long time scale coarse-grained molecular dynamics simulations. By introducing compartments with varying degrees of restraints within the membrane models, we are able to reveal how compartmentalization caused by cytoskeletal immobilization leads to reduced and anomalous diffusional mobility of both proteins and lipids. This in turn results in a reduced rate of protein dimerization within the membrane and of hopping of membrane proteins between compartments. These simulations provide a molecular realization of hierarchical models often invoked to explain single-molecule imaging studies of membrane proteins.

Elucidating the general principles of cell adhesion with a coarse-grained simulation model

Coarse-grained simulation of interplay between cell adhesion and cell signaling.

Across the great divide: the plant cell surface continuum

The plant cell wall, plasma membrane and cytoskeleton exist as a cell surface continuum. This interconnection of organelles forms the interface between the plant cell and the external environment and is important for detecting the presence of a diverse range of stimuli. A plethora of plasma membrane microdomains with putative roles in membrane localized enzymatic or signalling processes have been described. While regulation of cell wall composition is defined by proteins within the plasma membrane, the cell wall has been shown to have an anchoring role on plasma membrane proteins which affects their lateral mobility. This interplay between plasma membrane and cell wall components is necessary for plasma membrane microdomain function. Actin and microtubule cytoskeletons are also involved in maintenance and function of the cell surface continuum. Investigation of the interactions between organellar components of this mechanism are important if we are to understand how cells respond to external signals.

Molecular force transduction by ion channels: diversity and unifying principles

Cells perceive force through a variety of molecular sensors, of which the mechanosensitive ion channels are the most efficient and act the fastest. These channels apparently evolved to prevent osmotic lysis of the cell as a result of metabolite accumulation and/or external changes in osmolarity. From this simple beginning, nature developed specific mechanosensitive enzymes that allow us to hear, maintain balance, feel touch and regulate many systemic variables, such as blood pressure. For a channel to be mechanosensitive it needs to respond to mechanical stresses by changing its shape between the closed and open states. In that way, forces within the lipid bilayer or within a protein link can do work on the channel and stabilize its state. Ion channels have the highest turnover rates of all enzymes, and they can act as both sensors and effectors, providing the necessary fluxes to relieve osmotic pressure, shift the membrane potential or initiate chemical signaling. In this Commentary, we focus on the common mechanisms by which mechanical forces and the local environment can regulate membrane protein structure, and more specifically, mechanosensitive ion channels.

Characterization of differential protein tethering at the plasma membrane in response to epidermal growth factor signaling

Physical tethering of membrane proteins to the cortical actin cytoskeleton provides functional organization to the plasma membrane and contributes to diverse cellular processes including cell signaling, vesicular trafficking, endocytosis, and migration. For these processes to occur, membrane protein tethering must be dynamically regulated in response to environmental cues. In this study, we describe a novel biochemical scheme for isolating the complement of plasma membrane proteins that are physically tethered to the actin cytoskeleton. We utilized this method in combination with tandem liquid chromatography/mass spectrometry (LC-MS/MS) to demonstrate that cytoskeletal tethering of membrane proteins is acutely regulated by epidermal growth factor (EGF) in normal human kidney (HK2) cells. Our results indicate that several proteins known to be involved in EGF signaling, as well as other proteins not traditionally associated with this pathway, are tethered to the cytoskeleton in dynamic fashion. Further analysis of one hit from our proteomic survey, the receptor phosphotyrosine phosphatase PTPRS, revealed a correlation between cytoskeletal tethering and endosomal trafficking in response to EGF. This finding parallels previous indications that PTPRS is involved in the desensitization of EGFR and provides a potential mechanism to coordinate localization of these two membrane proteins in the same compartment upon EGFR activation.

Neural membrane signaling platforms

Throughout much of the history of biology, the cell membrane was functionally defined as a semi-permeable barrier separating aqueous compartments, and an anchoring site for proteins. Little attention was devoted to its possible regulatory role in intracellular molecular processes and neuron electrical signaling. This article reviews the history of membrane studies and the current state of the art. Emphasis is placed on natural and artificial membrane studies of electric field effects on molecular organization, especially as these may relate to impulse propagation in neurons. Implications of these studies for new designs in artificial intelligence are briefly examined.

Adaptive coarse-grained Monte Carlo simulation of reaction and diffusion dynamics in heterogeneous plasma membranes

Background

An adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.

Results

Coarse-graining of dimerization, rather than of diffusion, may lead to computational error. It is shown that the ACGMC method is an effective technique to minimize error in diffusion-reaction processes and is superior to the microscopic kinetic Monte Carlo simulation in terms of computational cost while retaining accuracy. The low computational cost enables sensitivity analysis calculations. Sensitivity analysis indicates that it may be possible to retain clusters of receptors over the time scale of minutes under suitable conditions and the cluster lifetime may depend on both receptor density and cluster size.

Conclusions

The ACGMC method is an ideal platform to resolve large length and time scales in heterogeneous biological systems well beyond the plasma membrane and the EGFR system studied here. Our results demonstrate that cluster size must be considered in conjunction with receptor density, as they synergistically affect EGFR cluster lifetime. Further, the cluster lifetime being of the order of several seconds suggests that any mechanisms responsible for EGFR aggregation must operate on shorter timescales (at most a fraction of a second), to overcome dissipation and produce stable clusters observed experimentally.

A dynamic corral model of receptor trafficking at a synapse

The postsynaptic density (PSD) is a cytoskeletal specialization within the postsynaptic membrane of a neuron that helps to concentrate and organize neurotransmitter receptors at a chemical synapse. The total number of receptors within the PSD, which is a major factor in determining the physiological strength or weight of a synapse, fluctuates due to the surface diffusion of receptors into and out of the PSD, and the interactions of receptors with scaffolding proteins and cytoskeletal elements within the PSD. In this article, we present a stochastic model of protein receptor trafficking at the PSD that takes into account these various processes. The PSD is treated as a stochastically gated corral, which contributes a source of extrinsic or environmental noise that supplements the intrinsic noise arising from small receptor numbers. Using a combination of stochastic analysis and Monte Carlo simulations, we determine the time-dependent variation in the mean and variance of synaptic receptor numbers for a variety of initial conditions that simulate fluorescence recovery after photobleaching experiments, and indicate how such data might be used to infer certain properties of the PSD.

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.

Investigating membrane protein dynamics in living cells

Live cell imaging is a powerful tool for understanding the function and regulation of membrane proteins. In this review, we briefly discuss 4 fluorescence-microscopy-based techniques for studying the transport dynamics of membrane proteins: fluorescence-correlation spectroscopy, image-correlation spectroscopy, fluorescence recovery after photobleaching, and single-particle and (or) molecule tracking. The advantages and limitations of each approach are illustrated using recent studies of an ion channel and cell adhesion molecules.

Molecular transport in a crowded volume created from vertically aligned carbon nanofibres: a fluorescence recovery after photobleaching study

Rapid and selective molecular exchange across a barrier is essential for emulating the properties of biological membranes. Vertically-aligned carbon nanofibre (VACNF) forests have shown great promise as membrane mimics, owing to their mechanical stability, their ease of integration with microfabrication technologies and the ability to tailor their morphology and surface properties. However, quantifying transport through synthetic membranes having micro- and nanoscale features is challenging. Here, fluorescence recovery after photobleaching (FRAP) is coupled with finite difference and Monte Carlo simulations to quantify diffusive transport in microfluidic structures containing VACNF forests. Anomalous subdiffusion was observed for FITC (hydrodynamic radius of 0.54 nm) diffusion through both VACNFs and SiO(2)-coated VACNFS (oxVACNFs). Anomalous subdiffusion can be attributed to multiple FITC-nanofibre interactions for the case of diffusion through the VACNF forest. Volume crowding was identified as the cause of anomalous subdiffusion in the oxVACNF forest. In both cases the diffusion mode changes to a time-independent, Fickian mode of transport that can be defined by a crossover length (R(CR)). By identifying the space-and time-dependent transport characteristics of the VACNF forest, the dimensional features of membranes can be tailored to achieve predictable molecular exchange.

Neural membrane microdomains as computational systems: Toward molecular modeling in the study of neural disease

Several studies indicate that the lipid biological membrane contains discrete regions known as rafts or microdomains. These structures range in size from approximately 50 to 70nm to nearly a mum and play important roles in cell signaling. In the neuron, computational models suggest that transiently polarized microdomain ethenes may regulate ion-channel dynamics and control impulse propagation. Thus the microdomain is nominated as the fundamental unit of nervous system signaling. Based on this model, the article proposes a first-approximation design for a supported-membrane device which would mimic microdomain properties. The basic architecture would consist of an electrically addressable biotemplated nanowire crossing an artificial membrane corralled in a vertical carbon nanofiber barrier. Advantages and disadvantages of model components are discussed at length. It is proposed that artificial devices of this type would be medically useful in simulating membrane states correlated with neural disease. This possibility is examined with reference to the A-current potassium channel, implicated in epilepsy.

Implicit solvent simulation models for biomembranes

Fully atomic simulation strategies are infeasible for the study of many processes of interest to membrane biology, biophysics and biochemistry. We review various coarse-grained simulation methodologies with special emphasis on methods and models that do not require the explicit simulation of water. Examples from our own research demonstrate that such models have potential for simulating a variety of biologically relevant phenomena at the membrane surface.

Dynamic simulations of membranes with cytoskeletal interactions

We describe a simulation algorithm for the dynamics of elastic membrane sheets over long length and time scales. Our model includes implicit hydrodynamic coupling between membrane and surrounding solvent and allows for arbitrary external forces acting on the membrane surface. In particular, the methodology is well suited to studying membranes in interaction with cytoskeletal filaments. We present results for the thermal undulations of a lipid bilayer attached to a regular network of spectrin filaments as a model for the red blood cell membrane. The dynamic fluctuations of the bilayer over the spectrin network are quantified and used to predict the macroscopic diffusion constant of band 3 on the surface of the red blood cell. We find that thermal undulations likely play a role in the mobility of band 3 in the plane of the erythrocyte membrane.

Modes of remodeling in the cortical cytoskeleton of vascular endothelial cells

The cortical cytoskeleton of vascular endothelial cells plays an important role in responding to mechanical stimuli and controlling the distribution of cell surface proteins. Here, we have used atomic force microscopy to visualize the dynamics of cortical cytoskeleton in living bovine pulmonary artery endothelial cells. We demonstrate that the cortical cytoskeleton, organized as a complex polygonal mesh, is highly dynamic and shows two modes of remodeling: intact-boundary-mode where mesh element boundaries remain intact but move at approximately 0.08 microm/min allowing the mesh element to change shape, and altered-boundary-mode where new mesh boundaries form and existing ones disappear.

Dynamics of pinned membranes with application to protein diffusion on the surface of red blood cells

We present a theoretical treatment and simulation algorithm for the dynamics of Helfrich elastic membrane surfaces in the presence of general harmonic perturbations and hydrodynamic coupling to the surrounding solvent. In the limit of localized and strong interactions, this harmonic model can be used to pin the membrane to intracellular/intercellular structures. We consider the case of pinning to the cytoskeleton and use such a model to estimate the macroscopic diffusion constant for band 3 protein on the surface of human erythrocytes. Comparison to experimental results suggests that thermal undulations of the membrane surface should play a significant role in protein mobility on the red blood cell.

One-dimensional transport with dynamic disorder

We study the mean quenching time distribution and its moments in a one-dimensional N-site donor-bridge-acceptor system where all sites are coupled to a two-state jump bath for arbitrary disorder and an arbitrary ratio kappa identical with <k>/R of the bath jump rate R and the average hopping rate <k>. When kappaN approximately 1, the quenching time distribution has long power-law tails even when the waiting times are exponentially distributed. These disappear for kappaN<<1 where the hopping rate self-averages on the bath relaxation time scale. In the absence of disorder or for small kappa, the mean quenching time scales linearly with N. Otherwise, we observe a power law, approximately N1+gamma, with a crossover to linear scaling (gamma=0) for large N. Distributions of particle position, its second moment, velocity and diffusion coefficient are computed in the infinite N limit. For times longer than R-1, the dynamic disorder self-averages and the average position, velocity, and diffusion coefficient scale linearly in time.

Regulation of protein mobility via thermal membrane undulations

The in-plane diffusivelike motion of membrane bound proteins on the surface of cells is considered. We suggest, on the basis of theoretical arguments and simulation, that thermally excited undulations of the lipid bilayer may serve as a mechanism for proteins to hop between adjacent regions on the cell surface separated by barriers composed of internal cellular structure (e.g., the cytoskeleton). We specifically investigate the mobility of band 3 dimer on the surface of red blood cells where the spectrin cytoskeletal meshwork defines a series of "corrals" on the cell surface known to hinder protein motion. Previous models of this system have postulated that the cytoskeleton must deform to allow passage of membrane bound proteins out of these corral regions and have ignored fluctuations of the bilayer. Our model provides a complementary mechanism and we posit that the mobility of real proteins in real cells is likely the result of several mechanisms acting in parallel.

Confined diffusion without fences of a g-protein-coupled receptor as revealed by single particle tracking

Single particle tracking is a powerful tool for probing the organization and dynamics of the plasma membrane constituents. We used this technique to study the micro -opioid receptor belonging to the large family of the G-protein-coupled receptors involved with other partners in a signal transduction pathway. The specific labeling of the receptor coupled to a T7-tag at its N-terminus, stably expressed in fibroblastic cells, was achieved by colloidal gold coupled to a monoclonal anti T7-tag antibody. The lateral movements of the particles were followed by nanovideomicroscopy at 40 ms time resolution during 2 min with a spatial precision of 15 nm. The receptors were found to have either a slow or directed diffusion mode (10%) or a walking confined diffusion mode (90%) composed of a long-term random diffusion and a short-term confined diffusion, and corresponding to a diffusion confined within a domain that itself diffuses. The results indicate that the confinement is due to an effective harmonic potential generated by long-range attraction between the membrane proteins. A simple model for interacting membrane proteins diffusion is proposed that explains the variations with the domain size of the short-term and long-term diffusion coefficients.

Lateral diffusion of membrane proteins in the presence of static and dynamic corrals: suggestions for appropriate observables

We consider the possibility of inferring the nature of cytoskeletal interaction with transmembrane proteins via optical experiments such as single-particle tracking (SPT) and near-field scanning optical microscopy (NSOM). In particular, we demonstrate that it may be possible to differentiate between static and dynamic barriers to diffusion by examining the time-dependent variance and higher moments of protein population inside cytoskeletal "corrals." Simulations modeling Band 3 diffusion on the surface of erythrocytes provide a concrete demonstration that these statistical tools might prove useful in the study of biological systems.