Diffusing proteins on a fluctuating membrane: analytical theory and simulations

Diffusion coefficients and MSD measurements on curved membranes and porous media

We study some geometric aspects that influence the transport properties of particles that diffuse on curved surfaces. We compare different approaches to surface diffusion based on the Laplace-Beltrami operator adapted to predict concentration along entire membranes, confined subdomains along surfaces, or within porous media. Our goal is to summarize, firstly, how diffusion in these systems results in different types of diffusion coefficients and mean square displacement measurements, and secondly, how these two factors are affected by the concavity of the surface, the shape of the possible barriers or obstacles that form the available domains, the sinuosity, tortuosity, and constrictions of the trajectories and even how the observation plane affects the measurements of the diffusion. In addition to presenting a critical and organized comparison between different notions of MSD, in this review, we test the correspondence between theoretical predictions and numerical simulations by performing finite element simulations and illustrate some situations where diffusion theory can be applied. We briefly reviewed computational schemes for understanding surface diffusion and finally, discussed how this work contributes to understanding the role of surface diffusion transport properties in porous media and their relationship to other transport processes.

Kinetic relaxation of giant vesicles validates diffusional softening in a binary lipid mixture

The stiffness of biological membranes determines the work required by cellular machinery to form and dismantle vesicles and other lipidic shapes. Model membrane stiffness can be determined from the equilibrium distribution of giant unilamellar vesicle surface undulations observable by phase contrast microscopy. With two or more components, lateral fluctuations of composition will couple to surface undulations depending on the curvature sensitivity of the constituent lipids. The result is a broader distribution of undulations whose complete relaxation is partially determined by lipid diffusion. In this work, kinetic analysis of the undulations of giant unilamellar vesicles made of phosphatidylcholine-phosphatidylethanolamine mixtures validates the molecular mechanism by which the membrane is made 25% softer than a single-component one. The mechanism is relevant to biological membranes, which have diverse and curvature-sensitive lipids.

Capsule self-oscillating gels showing cell-like nonthermal membrane/shape fluctuations

A primary interest in cell membrane and shape fluctuations is establishing experimental models reflecting only nonthermal active contributions. Here we report a millimeter-scaled capsule self-oscillating gel model mirroring the active contribution effect on cell fluctuations. In the capsule self-oscillating gels, the propagating chemical signals during a Belousov-Zhabotinsky (BZ) reaction induce simultaneous local deformations in the various regions, showing cell-like shape fluctuations. The capsule self-oscillating gels do not fluctuate without the BZ reaction, implying that only the active chemical parameter induces the gel fluctuations. The period and amplitude depend on the gel layer thickness and the concentration of the chemical substrate for the BZ reaction. Our results allow for a solid experimental platform showing actively driven cell-like fluctuations, which can potentially contribute to investigating the active parameter effect on cell fluctuations.

Brownian particles in periodic potentials: Coarse-graining versus fine structure

We study the motion of an overdamped particle connected to a thermal heat bath in the presence of an external periodic potential in one dimension. When we coarse-grain, i.e., bin the particle positions using bin sizes that are larger than the periodicity of the potential, the packet of spreading particles, all starting from a common origin, converges to a normal distribution centered at the origin with a mean-squared displacement that grows as 2D^{*}t, with an effective diffusion constant that is smaller than that of a freely diffusing particle. We examine the interplay between this coarse-grained description and the fine structure of the density, which is given by the Boltzmann-Gibbs (BG) factor e^{-V(x)/k_{B}T}, the latter being nonnormalizable. We explain this result and construct a theory of observables using the Fokker-Planck equation. These observables are classified as those that are related to the BG fine structure, like the energy or occupation times, while others, like the positional moments, for long times, converge to those of the large-scale description. Entropy falls into a special category as it has a coarse-grained and a fine structure description. The basic thermodynamic formula F=TS-E is extended to this far-from-equilibrium system. The ergodic properties are also studied using tools from infinite ergodic theory.

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

Modeling Receptor Motility along Advecting Lipid Membranes

This work aims to overview multiphysics mechanobiological computational models for receptor dynamics along advecting cell membranes. Continuum and statistical models of receptor motility are the two main modeling methodologies identified in reviewing the state of the art. Within the former modeling class, a further subdivision based on different biological purposes and processes of proteins’ motion is recognized; cell adhesion, cell contractility, endocytosis, and receptor relocations on advecting membranes are the most relevant biological processes identified in which receptor motility is pivotal. Numerical and/or experimental methods and approaches are highlighted in the exposure of the reviewed works provided by the literature, pertinent to the topic of the present manuscript. With a main focus on the continuum models of receptor motility, we discuss appropriate multiphyisics laws to model the mass flux of receptor proteins in the reproduction of receptor relocation and recruitment along cell membranes to describe receptor–ligand chemical interactions, and the cell’s structural response. The mass flux of receptor modeling is further supported by a discussion on the methodology utilized to evaluate the protein diffusion coefficient developed over the years.

Binding of a Brownian nanoparticle to a thermally fluctuating membrane surface

We investigate the Brownian dynamics of a nanoparticle bound to a thermally undulating elastic membrane. The ligand-functionalized nanoparticle is assumed to interact monovalently with the receptor expressed on the membrane. In order to resolve the nanoparticle transient motion subject to the instantaneous membrane configuration in a consistent manner, we employ a set of coupled Langevin equations that simultaneously incorporate the hydrodynamic effects, ligand-receptor binding interaction, intramembrane elastic forces, and thermal fluctuations. We show that the presence of a deformable, elastic fluid membrane not only affects the dynamics of a bound nanoparticle but also alters the effective binding potential felt by the nanoparticle. In contrast to a nanoparticle bound to a flat surface, the oscillatory characteristics of the nanoparticle velocity autocorrelation function are suppressed and transition to an anticorrelated long-time tail. Moreover, the nanoparticle position fluctuation becomes more coherent with that of the membrane binding site, and the width of the distribution of the nanoparticle distance from the membrane decreases with increasing membrane bending rigidity. By introducing a locally harmonic, bistable potential as an effective potential for the ligand-receptor pair, the rate of nanoparticle transitioning between two bound states is facilitated by membrane undulations as a result of stronger positional variations associated with the nanoparticle.

CryoEM reveals how the complement membrane attack complex ruptures lipid bilayers

The membrane attack complex (MAC) is one of the immune system's first responders. Complement proteins assemble on target membranes to form pores that lyse pathogens and impact tissue homeostasis of self-cells. How MAC disrupts the membrane barrier remains unclear. Here we use electron cryo-microscopy and flicker spectroscopy to show that MAC interacts with lipid bilayers in two distinct ways. Whereas C6 and C7 associate with the outer leaflet and reduce the energy for membrane bending, C8 and C9 traverse the bilayer increasing membrane rigidity. CryoEM reconstructions reveal plasticity of the MAC pore and demonstrate how C5b6 acts as a platform, directing assembly of a giant β-barrel whose structure is supported by a glycan scaffold. Our work provides a structural basis for understanding how β-pore forming proteins breach the membrane and reveals a mechanism for how MAC kills pathogens and regulates cell functions.

Diffusion and Directionality of Charged Nanoparticles on Lipid Bilayer Membrane

Diffusion dynamics of charged nanoparticles on the lipid membrane is of essential importance to cellular functioning. Yet a fundamental insight into electrostatics-mediated diffusion dynamics of charged nanoparticles on the membrane is lacking and remains to be an urgent issue. Here we present the computational investigation to uncover the pivotal role of electrostatics in the diffusion dynamics of charged nanoparticles on the lipid membrane. Our results demonstrate diffusive behaviors and directional transport of a charged nanoparticle, significantly depending on the sign and spatial distribution of charges on its surface. In contrast to the Fickian diffusion of neutral nanoparticles, randomly charged nanoparticles undergo superdiffusive transport with directionality. However, the dynamics of uniformly charged nanoparticles favors Fickian diffusion that is significantly enhanced. Such observations can be explained in term of electrostatics-induced surface reconstruction and fluctuation of lipid membrane. We finally present an analytical model connecting surface reconstruction and local deformation of the membrane. Our findings bear wide implications for the understanding and control of the transport of charged nanoparticles on the cell membrane.

Suppressing membrane height fluctuations leads to a membrane-mediated interaction among proteins

Membrane-induced interactions can play a significant role in the spatial distribution of membrane-bound proteins. We develop a model that combines a continuum description of lipid bilayers with a discrete particle model of proteins to probe the emerging structure of the combined membrane-protein system. Our model takes into account the membrane's elastic behavior, the steric repulsion between proteins, and the quenching of membrane shape fluctuations due to the presence of the proteins. We employ coupled Langevin equations to describe the dynamics of the system. We show that coupling to the membrane induces an attractive interaction among proteins, which may contribute to the clustering of proteins in biological membranes. We investigate the lateral protein diffusion and find that it is reduced due to transient fluctuations in membrane shape.

Curvature correction to the mobility of fluid membrane inclusions

Using rigorous low-Reynolds-number hydrodynamic theory on curved surfaces, we provide, via a Stokeslet-type approach, a general and concise expression for the leading-order curvature correction to the canonical, planar, Saffman-Delbrück value of the diffusion constant for a small inclusion embedded in an arbitrarily (albeit weakly) curved fluid membrane. In order to demonstrate the efficacy and utility of this general result, we apply our theory to the specific case of calculating the diffusion coefficient of a locally curvature inducing membrane inclusion. By including both the effects of inclusion and membrane elasticity, as well as their respective thermal shape fluctuations, excellent agreement is found with recently published experimental data on the surface tension dependent mobility of membrane bound inclusions.

Membrane stiffness is modified by integral membrane proteins

The ease with which a cell membrane can bend and deform is important for a wide range of biological functions. Peripheral proteins that induce curvature in membranes (e.g. BAR domains) have been studied for a number of years. Little is known, however, about the effect of integral membrane proteins on the stiffness of a membrane (characterised by the bending rigidity, Kc). We demonstrate by computer simulation that adding integral membrane proteins at physiological densities alters the stiffness of the membrane. First we establish that the coarse-grained MARTINI forcefield is able to accurately reproduce the bending rigidity of a small patch of 1500 phosphatidyl choline lipids by comparing the calculated value to both experiment and an atomistic simulation of the same system. This enables us to simulate the dynamics of large (ca. 50 000 lipids) patches of membrane using the MARTINI coarse-grained description. We find that altering the lipid composition changes the bending rigidity. Adding integral membrane proteins to lipid bilayers also changes the bending rigidity, whilst adding a simple peripheral membrane protein has no effect. Our results suggest that integral membrane proteins can have different effects, and in the case of the bacterial outer membrane protein, BtuB, the greater the density of protein, the larger the reduction in stiffness.

Anomalous Dynamical Behavior of Freestanding Graphene Membranes

We report subnanometer, high-bandwidth measurements of the out-of-plane (vertical) motion of atoms in freestanding graphene using scanning tunneling microscopy. By tracking the vertical position over a long time period, a 1000-fold increase in the ability to measure space-time dynamics of atomically thin membranes is achieved over the current state-of-the-art imaging technologies. We observe that the vertical motion of a graphene membrane exhibits rare long-scale excursions characterized by both anomalous mean-squared displacements and Cauchy-Lorentz power law jump distributions.

Hydrodynamic coupling of particle inclusions embedded in curved lipid bilayer membranes

We develop theory and computational methods to investigate particle inclusions embedded within curved lipid bilayer membranes. We consider the case of spherical lipid vesicles where inclusion particles are coupled through (i) intramembrane hydrodynamics, (ii) traction stresses with the external and trapped solvent fluid, and (iii) intermonolayer slip between the two leaflets of the bilayer. We investigate relative to flat membranes how the membrane curvature and topology augment hydrodynamic responses. We show how both the translational and rotational mobility of protein inclusions are effected by the membrane curvature, ratio of intramembrane viscosity to solvent viscosity, and intermonolayer slip. For general investigations of many-particle dynamics, we also discuss how our approaches can be used to treat the collective diffusion and hydrodynamic coupling within spherical bilayers.

Hydrodynamics of bilayer membranes with diffusing transmembrane proteins

We consider the hydrodynamics of lipid bilayers containing transmembrane proteins of arbitrary shape. This biologically-motivated problem is relevant to the cell membrane, whose fluctuating dynamics play a key role in phenomena ranging from cell migration, intercellular transport, and cell communication. Using Onsager's variational principle, we derive the equations that govern the relaxation dynamics of the membrane shape, of the mass densities of the bilayer leaflets, and of the diffusing proteins' concentration. With our generic formalism, we obtain several results on membrane dynamics. We find that proteins that span the bilayer increase the intermonolayer friction coefficient. The renormalization, which can be significant, is in inverse proportion to the protein's mobility. Second, we find that asymmetric proteins couple to the membrane curvature and to the difference in monolayer densities. For practically all accessible membrane tensions (σ > 10(-8) N m(-1)) we show that the protein density is the slowest relaxing variable. Furthermore, its relaxation rate decreases at small wavelengths due to the coupling to curvature. We apply our formalism to the large-scale diffusion of a concentrated protein patch. We find that the diffusion profile is not self-similar, owing to the wavevector dependence of the effective diffusion coefficient.

Mobility Measurements Probe Conformational Changes in Membrane Proteins due to Tension

The function of membrane-embedded proteins such as ion channels depends crucially on their conformation. We demonstrate how conformational changes in asymmetric membrane proteins may be inferred from measurements of their diffusion. Such proteins cause local deformations in the membrane, which induce an extra hydrodynamic drag on the protein. Using membrane tension to control the magnitude of the deformations, and hence the drag, measurements of diffusivity can be used to infer-via an elastic model of the protein-how conformation is changed by tension. Motivated by recent experimental results [Quemeneur et al., Proc. Natl. Acad. Sci. U.S.A. 111, 5083 (2014)], we focus on KvAP, a voltage-gated potassium channel from Aeropyrum pernix. The conformation of KvAP is found to change considerably due to tension, with its "walls," where the protein meets the membrane, undergoing significant angular strains. The torsional stiffness is determined to be 26.8k(B)T per radian at room temperature. This has implications for both the structure and the function of such proteins in the environment of a tension-bearing membrane.

Seeing the Forest in Lieu of the Trees: Continuum Simulations of Cell Membranes at Large Length Scales

Biological membranes exhibit long-range spatial structure in both chemical composition and geometric shape, which gives rise to remarkable physical phenomena and important biological functions. Continuum models that describe these effects play an important role in our understanding of membrane biophysics at large length scales. We review the mathematical framework used to describe both composition and shape degrees of freedom, and present best practices to implement such models in a computer simulation. We discuss in detail two applications of continuum models of cell membranes: the formation of microemulsion and modulated phases, and the effect of membrane-mediated interactions on the assembly of membrane proteins.

Shape matters in protein mobility within membranes

The lateral mobility of proteins within cell membranes is usually thought to be dependent on their size and modulated by local heterogeneities of the membrane. Experiments using single-particle tracking on reconstituted membranes demonstrate that protein diffusion is significantly influenced by the interplay of membrane curvature, membrane tension, and protein shape. We find that the curvature-coupled voltage-gated potassium channel (KvAP) undergoes a significant increase in protein mobility under tension, whereas the mobility of the curvature-neutral water channel aquaporin 0 (AQP0) is insensitive to it. Such observations are well explained in terms of an effective friction coefficient of the protein induced by the local membrane deformation.

Diffusion of a particle quadratically coupled to a thermally fluctuating field

We study the diffusion of a Brownian particle quadratically coupled to a thermally fluctuating field. In the weak-coupling limit, a path-integral formulation allows us to compute the effective diffusion coefficient in the cases of an active particle, which tends to suppress field fluctuations, and of a passive particle, which only undergoes field fluctuations. We show that the behavior is similar to what was previously found for a linear coupling: an active particle is always slowed down, whereas a passive particle is slowed down in a slow field and accelerated in a fast field. Numerical simulations show a good agreement with the analytical calculations. The examples of a membrane protein coupled to the curvature or composition of the membrane are discussed, with a focus on the room for anomalous diffusion.

Instabilities and diffusion in a hydrodynamic model of a fluid membrane coupled to a thin active fluid layer

We construct a coarse-grained effective two-dimensional (2d hydrodynamic theory as a theoretical model for a coupled system of a fluid membrane and a thin layer of a polar active fluid in its ordered state that is anchored to the membrane. We show that such a system is prone to generic instabilities through the interplay of nonequilibrium drive, polar order and membrane fluctuation. We use our model equations to calculate diffusion coefficients of an inclusion in the membrane and show that their values depend strongly on the system size, in contrast to their equilibrium values. Our work extends the work of S. Sankararaman and S. Ramaswamy (Phys. Rev. Lett., 102, 118107 (2009)) to a coupled system of a fluid membrane and an ordered active fluid layer. Our model is broadly inspired by and should be useful as a starting point for theoretical descriptions of the coupled dynamics of a cell membrane and a cortical actin layer anchored to it.

Effect of hydroperoxides on red blood cell membrane mechanical properties

We investigate the effect of oxidative stress on red blood cell membrane mechanical properties in vitro using detailed analysis of the membrane thermal fluctuation spectrum. Two different oxidants, the cytosol-soluble hydrogen peroxide and the membrane-soluble cumene hydroperoxide, are used, and their effects on the membrane bending elastic modulus, surface tension, strength of confinement due to the membrane skeleton, and 2D shear elastic modulus are measured. We find that both oxidants alter significantly the membrane elastic properties, but their effects differ qualitatively and quantitatively. While hydrogen peroxide mainly affects the elasticity of the membrane protein skeleton (increasing the membrane shear modulus), cumene hydroperoxide has an impact on both membrane skeleton and lipid bilayer mechanical properties, as can be seen from the increased values of the shear and bending elastic moduli. The biologically important implication of these results is that the effects of oxidative stress on the biophysical properties, and hence the physiological functions, of the cell membrane depend on the nature of the oxidative agent. Thermal fluctuation spectroscopy provides a means of characterizing these different effects, potentially in a clinical milieu.

In-plane dynamics of membranes with immobile inclusions

Cell membranes are anchored to the cytoskeleton via immobile inclusions. We investigate the effect of such anchors on the in-plane dynamics of a fluid membrane and mobile inclusions (proteins) embedded in it. The immobile particles lead to a decreased diffusion coefficient of mobile ones and suppress the correlated diffusion of particle pairs. Because of the long-range, quasi-two-dimensional nature of membrane flows, these effects become significant at a low area fraction (below 1%) of immobile inclusions.

Continuum simulations of biomembrane dynamics and the importance of hydrodynamic effects

Traditional particle-based simulation strategies are impractical for the study of lipid bilayers and biological membranes over the longest length and time scales (microns, seconds and longer) relevant to cellular biology. Continuum-based models developed within the frameworks of elasticity theory, fluid dynamics and statistical mechanics provide a framework for studying membrane biophysics over a range of mesoscopic to macroscopic length and time regimes, but the application of such ideas to simulation studies has occurred only relatively recently. We review some of our efforts in this direction with emphasis on the dynamics in model membrane systems. Several examples are presented that highlight the prominent role of hydrodynamics in membrane dynamics and we argue that careful consideration of fluid dynamics is key to understanding membrane biophysics at the cellular scale.

Perturbative path-integral study of active- and passive-tracer diffusion in fluctuating fields

We study the effective diffusion constant of a Brownian particle linearly coupled to a thermally fluctuating scalar field. We use a path-integral method to compute the effective diffusion coefficient perturbatively to lowest order in the coupling constant. This method can be applied to cases where the field is affected by the particle (an active tracer) and cases where the tracer is passive. Our results are applicable to a wide range of physical problems, from a protein diffusing in a membrane to the dispersion of a passive tracer in a random potential. In the case of passive diffusion in a scalar field, we show that the coupling to the field can, in some cases, speed up the diffusion corresponding to a form of stochastic resonance. Our results on passive diffusion are also confirmed via a perturbative calculation of the probability density function of the particle in a Fokker-Planck formulation of the problem. Numerical simulations on simplified systems corroborate our results.

Microphase separation in nonequilibrium biomembranes

Compositional heterogeneities of cell membranes are thought to play an important role in many physiological processes. We study how variations in the membrane composition can be driven by nonthermal fluctuating forces and therefore show how these can occur relatively far from any critical point for the membrane. We show that the membrane steady state is not only controlled by the strength of the forces and how they couple to the membrane, but also by their dynamics: In a simple class of models this is captured by a single force correlation time. We conclude that the coupling of membrane composition to normal mechanical forces, such as might be exerted by polymerizing cytoskeleton filaments, could play an important role in controlling the steady state of a cell membrane that exhibits transient lateral modulations of its composition on length scales in the 10-100 nm regime.

Hydrodynamic coupling between two fluid membranes

The coupled in-plane diffusion dynamics between point-particles embedded in stacked fluid membranes is investigated. We calculate the contributions to the coupling longitudinal and transverse diffusion coefficients due to particle motion within the different as well as the same membranes. The stacked geometry leads to a hydrodynamic coupling between the two membranes.

Correlated dynamics of inclusions in a supported membrane

The hydrodynamic theory of heterogeneous fluid membranes is extended to the case of a membrane adjacent to a solid substrate. We derive the coupling diffusion coefficients of pairs of membrane inclusions in the limit of large separation compared to the inclusion size. Two-dimensional compressive stresses in the membrane make the coupling coefficients decay asymptotically as 1/r(2) with interparticle distance r. For the common case, where the distance to the substrate is of submicrometer scale, we present expressions for the coupling between distant disklike inclusions, which are valid for arbitrary inclusion size. We calculate the effect of inclusions on the response of the membrane and the associated corrections to the coupling diffusion coefficients to leading order in the concentration of inclusions. While at short distances the response is modified as if the membrane were a two-dimensional suspension, the large-distance response is not renormalized by the inclusions.