Dynamic simulations of membranes with cytoskeletal interactions

Iron Overload-Related Oxidative Stress Leads to Hyperphosphorylation and Altered Anion Exchanger 1 (Band 3) Function in Erythrocytes from Subjects with β-Thalassemia Minor

β-thalassemia, a hereditary hemoglobinopathy, is caused by reduced or absent synthesis of the β-globin chains of hemoglobin. Three clinical conditions are recognized: β-thalassemia major, β-thalassemia intermedia, and β-thalassemia minor (β-Thal+). This latter condition occurs when an individual inherits a mutated β-globin gene from one parent. In erythrocytes from β-Thal+ subjects, the excess α-globin chains produce unstable α-tetramers, which can induce substantial oxidative stress leading to plasma membrane and cytoskeleton damage, as well as deranged cellular function. In the present study, we hypothesized that increased oxidative stress might lead to structural rearrangements in erythrocytes from β-Thal+ volunteers and functional alterations of ion transport proteins, including band 3 protein. The data obtained showed significant modifications of the cellular shape in erythrocytes from β-Thal+ subjects. In particular, a significantly increased number of elliptocytes was observed. Interestingly, iron overload, detected in erythrocytes from β-Thal+ subjects, provoked a significant production of reactive oxygen species (ROS), overactivation of the endogenous antioxidant enzymes catalase and superoxide dismutase, and glutathione depletion, resulting in (a) increased lipid peroxidation, (b) protein sulfhydryl group (-SH) oxidation. Iron overload-related oxidative stress affected Na+/K+-ATPase activity, which in turn may have contributed to impaired β-Thal+ erythrocyte deformability. As a result, alterations in the distribution of cytoskeletal proteins, including α/β-spectrin, protein 4.1, and α-actin, in erythrocytes from β-Thal+ subjects have been detected. Significantly, oxidative stress was also associated with increased phosphorylation and altered band 3 ion transport activity, as well as increased oxidized hemoglobin, which led to abnormal clustering and redistribution of band 3 on the plasma membrane. Taken together, these findings contribute to elucidating potential oxidative stress-related perturbations of ion transporters and associated cytoskeletal proteins, which may affect erythrocyte and systemic homeostasis in β-Thal+ subjects.

Measuring sub-nanometer undulations at microsecond temporal resolution with metal- and graphene-induced energy transfer spectroscopy

Out-of-plane fluctuations, also known as stochastic displacements, of biological membranes play a crucial role in regulating many essential life processes within cells and organelles. Despite the availability of various methods for quantifying membrane dynamics, accurately quantifying complex membrane systems with rapid and tiny fluctuations, such as mitochondria, remains a challenge. In this work, we present a methodology that combines metal/graphene-induced energy transfer (MIET/GIET) with fluorescence correlation spectroscopy (FCS) to quantify out-of-plane fluctuations of membranes with simultaneous spatiotemporal resolution of approximately one nanometer and one microsecond. To validate the technique and spatiotemporal resolution, we measure bending undulations of model membranes. Furthermore, we demonstrate the versatility and applicability of MIET/GIET-FCS for studying diverse membrane systems, including the widely studied fluctuating membrane system of human red blood cells, as well as two unexplored membrane systems with tiny fluctuations, a pore-spanning membrane, and mitochondrial inner/outer membranes.

AAPH-induced oxidative damage reduced anion exchanger 1 (SLC4A1/AE1) activity in human red blood cells: protective effect of an anthocyanin-rich extract

Introduction: During their lifespan in the bloodstream, red blood cells (RBCs) are exposed to multiple stressors, including increased oxidative stress, which can affect their morphology and function, thereby contributing to disease. Aim: This investigation aimed to explore the cellular and molecular mechanisms related to oxidative stress underlying anion exchanger 1 activity (band 3, SLC4A1/AE1) in human RBCs. To achieve this aim, the relationship between RBC morphology and functional and metabolic activity has been explored. Moreover, the potential protective effect of an anthocyanin-enriched fraction extracted from Callistemon citrinus flowers was studied. Methods: Cellular morphology, parameters of oxidative stress, as well as the anion exchange capability of band 3 have been analyzed in RBCs treated for 1 h with 50 mM of the pro-oxidant 2,2'-azobis (2-methylpropionamide)-dihydrochloride (AAPH). Before or after the oxidative insult, subsets of cells were exposed to 0.01 μg/mL of an anthocyanin-enriched fraction for 1 h. Results: Exposure to AAPH caused oxidative stress, exhaustion of reduced glutathione, and over-activation of the endogenous antioxidant machinery, resulting in morphological alterations of RBCs, specifically the formation of acanthocytes, increased lipid peroxidation and oxidation of proteins, as well as abnormal distribution and hyper-phosphorylation of band 3. Expected, oxidative stress was also associated with a decreased band 3 ion transport activity and an increase of oxidized haemoglobin, which led to abnormal clustering of band 3. Exposure of cells to the anthocyanin-enriched fraction prior to, but not after, oxidative stress efficiently counteracted oxidative stress-related alterations. Importantly, protection of band3 function from oxidative stress could only be achieved in intact cells and not in RBC ghosts. Conclusion: These findings contribute a) to clarify oxidative stress-related physiological and biochemical alterations in human RBCs, b) propose anthocyanins as natural antioxidants to neutralize oxidative stress-related modifications, and 3) suggest that cell integrity, and therefore a cytosolic component, is required to reverse oxidative stress-related pathophysiological derangements in human mature RBCs.

Aging Injury Impairs Structural Properties and Cell Signaling in Human Red Blood Cells; Açaì Berry Is a Keystone

Red blood cell (RBC) deformability is the ability of cells to modulate their shape to ensure transit through narrow capillaries of the microcirculation. A loss of deformability can occur in several pathological conditions, during natural RBC aging through an increase in membrane protein phosphorylation, and/or through the structural rearrangements of cytoskeletal proteins due to oxidative conditions, with a key role played by band 3. Due to the close relationship between aging and oxidative stress, flavonoid-rich foods are good candidates to counteract age-related alterations. This study aims to verify the beneficial role of Açaì extract in a d-Galactose (d-Gal)-induced model of aging in human RBCs. To this end, band 3 phosphorylation and structural rearrangements in membrane cytoskeleton-associated proteins, namely spectrin, ankyrin, and/or protein 4.1, are analyzed in RBCs treated with 100 mM d-Gal for 24 h, with or without pre-incubation with 10 μg/mL Açaì extract for 1 h. Furthermore, RBC deformability is also measured. Tyrosine phosphorylation of band 3, membrane cytoskeleton-associated proteins, and RBC deformability (elongation index) are analyzed using western blotting analysis, FACScan flow cytometry, and ektacytometry, respectively. The present data show that: (i) Açaì berry extract restores the increase in band 3 tyrosine phosphorylation and Syk kinase levels after exposure to 100 mM d-Gal treatment; and (ii) Açaì berry extract partially restores alterations in the distribution of spectrin, ankyrin, and protein 4.1. Interestingly, the significant decrease in membrane RBC deformability associated with d-Gal treatment is alleviated by pre-treatment with Açaì extract. These findings further contribute to clarify mechanisms of natural aging in human RBCs, and propose flavonoid substances as potential natural antioxidants for the treatment and/or prevention of oxidative-stress-related disease risk.

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.

Interplay between surface and bending energy helps membrane protrusion formation

We consider a one-dimensional elastic membrane, which is pushed by growing filaments. The filaments tend to grow by creating local protrusions in the membrane and this process has surface energy and bending energy costs. Although it is expected that with increasing surface tension and bending rigidity, it should become more difficult to create a protrusion, we find that for a fixed bending rigidity, as the surface tension increases, protrusions are more easily formed. This effect also gives rise to nontrivial dependence of membrane velocity on the surface tension, characterized by a dip and a peak. We explain this unusual phenomenon by studying in detail the interplay of the surface and the bending energy and show that this interplay is responsible for a qualitative shape change of the membrane, which gives rise to the above effect.

Statistical Mechanics of an Elastically Pinned Membrane: Equilibrium Dynamics and Power Spectrum

In biological settings, membranes typically interact locally with other membranes: the extracellular matrix in the exterior or internal cellular structures such as the cytoskeleton, locally pinning the membrane. Characterizing the dynamical properties of such interactions presents a difficult task. Significant progress has been achieved through simulations and experiments, yet analytical progress in modeling pinned membranes has been impeded by the complexity of governing equations. Here, we circumvent these difficulties by calculating analytically the time-dependent Green's function of the operator governing the dynamics of an elastically pinned membrane in a hydrodynamic surrounding and subject to external forces. This enables us to calculate the equilibrium power spectral density for an overdamped membrane pinned by an elastic, permanently attached spring subject to thermal excitations. By considering the effects of the finite experimental resolution on the measured spectra, we show that the elasticity of the pinning can be extracted from the experimentally measured spectrum. Membrane fluctuations can thus be used as a tool to probe mechanical properties of the underlying structures. Such a tool may be particularly relevant in the context of cell mechanics, in which the elasticity of the membrane's attachment to the cytoskeleton could be measured.

Simulation Best Practices for Lipid Membranes [Article v1.0]

We establish a reliable and robust standardization of settings for practical molecular dynamics (MD) simulations of pure and mixed (single- and multi-component) lipid bilayer membranes. In lipid membranes research, particle-based molecular simulations are a powerful tool alongside continuum theory, lipidomics, and model, in vitro, and in vivo experiments. Molecular simulations can provide precise and reproducible spatiotemporal (atomic- and femtosecond-level) information about membrane structure, mechanics, thermodynamics, kinetics, and dynamics. Yet the simulation of lipid membranes can be a daunting task, given the uniqueness of lipid membranes relative to conventional liquid-liquid and solid-liquid interfaces, the immense and complex thermodynamic and statistical mechanical theory, the diversity of multiscale lipid models, limitations of modern computing power, the difficulty and ambiguity of simulation controls, finite size effects, competitive continuum simulation alternatives, and the desired application, including vesicle experiments and biological membranes. These issues can complicate an essential understanding of the field of lipid membranes, and create major bottlenecks to simulation advancement. In this article, we clarify these issues and present a consistent, thorough, and user-friendly framework for the design of state-of-the-art lipid membrane MD simulations. We hope to allow early-career researchers to quickly overcome common obstacles in the field of lipid membranes and reach maximal impact in their simulations.

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.

Model for probing membrane-cortex adhesion by micropipette aspiration and fluctuation spectroscopy

We propose a model for membrane-cortex adhesion that couples membrane deformations, hydrodynamics, and kinetics of membrane-cortex ligands. In its simplest form, the model gives explicit predictions for the critical pressure for membrane detachment and for the value of adhesion energy. We show that these quantities exhibit a significant dependence on the active acto-myosin stresses. The model provides a simple framework to access quantitative information on cortical activity by means of micropipette experiments. We also extend the model to incorporate fluctuations and show that detailed information on the stability of membrane-cortex coupling can be obtained by a combination of micropipette aspiration and fluctuation spectroscopy measurements.

Simulating Membrane Dynamics in Nonhomogeneous Hydrodynamic Environments

Two previously introduced simulation algorithms for the dynamics of elastic membrane sheets embedded in a fluid medium are extended to account for inhomogeneous hydrodynamic environments. We calculate the height autocorrelation function for a lipid bilayer randomly pinned to a flat substrate and the influence of fluid confinement by the spectrin cytoskeleton on short wavelength membrane undulations of the human red blood cell. Altering the hydrodynamic environment of the membrane leads to significant changes in dynamics, and we discuss these effects in the context of recent experiments.

Fluctuating hydrodynamics of multicomponent membranes with embedded proteins

A simulation method for the dynamics of inhomogeneous lipid bilayer membranes is presented. The membrane is treated using stochastic Saffman-Delbrück hydrodynamics, coupled to a phase-field description of lipid composition and discrete membrane proteins. Multiple applications are considered to validate and parameterize the model. The dynamics of membrane composition fluctuations above the critical point and phase separation dynamics below the critical point are studied in some detail, including the effects of adding proteins to the mixture.

Toroidal membrane vesicles in spherical confinement

We investigate the morphology of a toroidal fluid membrane vesicle confined inside a spherical container. The equilibrium shapes are assembled in a geometrical phase diagram as a function of scaled area and reduced volume of the membrane. For small area the vesicle can adopt its free form. When increasing the area, the membrane cannot avoid contact and touches the confining sphere along a circular contact line, which extends to a zone of contact for higher area. The elastic energies of the equilibrium shapes are compared to those of their confined counterparts of spherical topology to predict under which conditions a topology change is favored energetically.

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.

The ratio of the lateral correlation length and particle radius determines the density profile of spherical molecules near a fluctuating membrane

Interactions between membranes and molecules are important for many biological processes, e.g., transport of molecules across cell membranes. However, the detailed physical description of the membrane-biomolecule system remains a challenge and simplified schemes allow capturing its main intrinsic features. In this work, by means of Monte Carlo computer simulations, we systematically study the distribution of uncharged spherical molecules in contact with a flexible surface. Our results show that the distribution for finite size particles has the same simple functional form as the one obtained for point-like particles and depends only on the ratio of the lateral correlation length of the membrane and the radius of the molecules.

Membrane sorting via the extracellular matrix

We consider the coupling between a membrane and the extracellular matrix. Computer simulations demonstrate that the latter coupling is able to sort lipids. It is assumed that membranes are elastic manifolds, and that this manifold is disrupted by the extracellular matrix. For a solid-supported membrane with an actin network on top, regions of positive curvature are induced below the actin fibers. A similar mechanism is conceivable by assuming that the proteins which connect the cytoskeleton to the membrane induce local membrane curvature. The regions of non-zero curvature exist irrespective of any phase transition the lipids themselves may undergo. For lipids that prefer certain curvature, the extracellular matrix thus provides a spatial template for the resulting lateral domain structure of the membrane.

Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins

Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.

Dynamic scaling in phase separation kinetics for quasi-two-dimensional membranes

We consider the dynamics of phase separation in lipid bilayer membranes, modeled as flat two-dimensional liquid sheets within a bulk fluid, both in the creeping flow approximation. We present scaling arguments that suggest asymptotic coarsening in these systems is characterized by a length scale R(t) ~ t(1/2) for critical (bicontinuous) phase separation and R(t) ~t(1/3) for off-critical concentrations (droplet morphology). In this limit, the bulk fluid is the primary source of dissipation. We also address these questions with continuum stochastic hydrodynamic simulations. We see evidence of scaling violation in critical phase separation, where isolated circular domains coarsen slower than elongated ones. However, we also find a region of apparent scaling where R(t) ~ t(1/2) is observed. This appears to be due to the competition of thermal and hydrodynamic effects. We argue that the diversity of scaling exponents measured in experiment and prior simulations can in part be attributed to certain measurements lying outside the asymptotic long-length-scale regime, and provide a framework to help understand these results. We also discuss a few simple generalizations to confined membranes and membranes in which inertia is relevant.

Beyond the creeping viscous flow limit for lipid bilayer membranes: theory of single-particle microrheology, domain flicker spectroscopy, and long-time tails

Recent experiments suggest that lipid bilayer membranes may be viscoelastic. We present a generalized "Saffman-Einstein" relation that may be used to determine the linear viscoelastic shear modulus from single-bead microrheology experiments on membranes. We show that viscoelastic parameters can also be extracted from membrane domain flicker spectroscopy experiments. Contributions from fluid inertia are expected to be negligible in both microrheology and domain flicker spectroscopy experiments, but can create a "long-time tail" in the membrane velocity autocorrelation function. In a viscous membrane, this tail crosses over from t(-1) at intermediate times, as in a two-dimensional fluid, to t(-3/2) at long times, as in a three-dimensional fluid. If the membrane is viscoelastic, the velocity autocorrelation function may be negative at intermediate times.

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.

Dynamic simulations of multicomponent lipid membranes over long length and time scales

We present a stochastic phase-field model for multicomponent lipid bilayers that explicitly accounts for the quasi-two-dimensional hydrodynamic environment unique to a thin fluid membrane immersed in aqueous solution. Dynamics over a wide range of length scales (from nanometers to microns) for durations up to seconds and longer are readily accessed and provide a direct comparison to fluorescence microscopy measurements in ternary lipid-cholesterol mixtures. Simulations of phase separation kinetics agree with experiment and elucidate the importance of hydrodynamics in the coarsening process.

Nanomechanics of multiple units in the erythrocyte membrane skeletal network

Erythrocytes undergo deformations when they transport O₂ and CO₂ across the membrane, yet the 3D nanomechanics of the skeletal network remains poorly understood. Expanding from our previous single isolated unit, we now simulate networks consisting of 1–10 concentric rings of repeating units in equibiaxial deformation. The networks are organized with (1) a 3D model for a single unit, (2) a wrap-around mode between Sp and actin protofilament in the intra-unit interaction, and (3) a random inter-unit connectivity. These assumptions permit efficient five-degrees-of-freedom (5DOF) simulations when up to 30 pN of radial forces are applied to the boundary spectrin (Sp) and the center and other units are analyzed. As 6 Sp balance their tensions, hexagonal units become irregular. While actin protofilaments remain tangent to the network, their yaw (Φ) angles change drastically with addition of neighboring units or an Sp unfolding. It is anticipated that during deformation, transmembrane complexes associated with the network move laterally through the lipid bilayer and increase the diffusion of molecules across the membrane. When protofilament/Sp sweeps under the lipid bilayer, they mix up the submembrane concentration gradient. Thus, the nanomechanics of actin protofilaments and Sp may enhance the transport of molecules during erythrocyte deformation.

Out-of-equilibrium behavior of Casimir-type fluctuation-induced forces for free classical fields

We present a general method to study the nonequilibrium behavior of Casimir-type fluctuation-induced forces for classical free scalar field theories. In particular, we analyze the temporal evolution of the force toward its equilibrium value when the field dynamics is given by a general class of overdamped stochastic dynamics (including the model A and model B classes). The steady-state force is also analyzed for systems which have nonequilibrium steady states, for instance, where they are driven by colored noise. The key to the method is that the out of equilibrium force is computed by specifying an energy of interaction between the field and the surfaces in the problem. In general, we find that there is a mapping of the dynamical problem onto a corresponding static one and in the case where the latter can be solved, the full dynamical behavior of the force can be extracted. The method is used to compute the nonequilibrium Casimir force induced between two parallel plates by a fluctuating field, in the cases of Dirichlet, Neumann, and mixed boundary conditions. Various other examples, such as the fluctuation-induced force between inclusions in fluctuating media, are discussed.

Adhesive dynamics of lubricated films

Membrane waves have been observed near the leading edge of a motile cell. Such phenomenon is the result of the interplay between hydrodynamics and adhesive dynamics. Here we consider membrane dynamics on a thin fluid gap supported by adhesive bonds. Using coupled lubrication theory and adhesive dynamics, we derive an evolution equation to account for membrane tension, bending, adhesion, and viscous lubrication. Four adhesion scenarios are examined: no adhesion, uniform adhesion, clustered adhesion, and focal adhesion. Two contrasting traveling wave types are found, namely, tension and adhesion waves. Tension waves disperse with time and space, whereas adhesion waves show increased amplitudes and are highly persistent. We show that the transition from tension to adhesion waves depends on a necessary, but insufficient, criterion that the wave amplitude must exceed a critical gap height, which is dependent on adhesion binding probability. We also show that strong adhesion results in sharp tension-to-adhesion wave transitions. The present work could explain the strong persistence of the waves observed in adhered cells using differential inference contrast (DIC) microscopy and the observation that the wavelengths decrease shortly after leading edge retraction.

Multiscale simulation of erythrocyte membranes

To quantitatively predict the mechanical response and mechanically induced remodeling of red blood cells, we developed a multiscale method to correlate distributions of internal stress with overall cell deformation. This method consists of three models at different length scales: in the complete cell level the membrane is modeled as two distinct layers of continuum shells using finite element method (Level III), in which the skeleton-bilayer interactions are depicted as a slide in the lateral (i.e., in-plane) direction (caused by the mobility of the skeleton-bilayer pinning points) and a normal contact force; the constitutive laws of the inner layer (the protein skeleton) are obtained from a molecular-based model (Level II); the mechanical properties of the spectrin (Sp, a key component of the skeleton), including its folding/unfolding reactions, are obtained with a stress-strain model (Level I). Model verification is achieved through comparisons with existing numerical and experimental studies in terms of the resting shape of the cell as well as cell deformations induced by micropipettes and optical tweezers. Detailed distributions of the interaction force between the lipid bilayer and the skeleton that may cause their dissociation and lead to phenomena such as vesiculation are predicted. Specifically, our model predicts correlation between the occurrence of Sp unfolding and increase in the mechanical load upon individual skeleton-bilayer pinning points. Finally a simulation of the necking process after skeleton-bilayer dissociation, a precursor of vesiculation, is conducted.

Diffusing proteins on a fluctuating membrane: analytical theory and simulations

Using analytical calculations and computer simulations, we consider both the lateral diffusion of a membrane protein and the fluctuation spectrum of the membrane in which the protein is embedded. The membrane protein interacts with the membrane shape through its spontaneous curvature and bending rigidity. The lateral motion of the protein may be viewed as diffusion in an effective potential, hence, the effective mobility is always reduced compared to the case of free diffusion. Using a rigorous path-integral approach, we derive an analytical expression for the effective diffusion coefficient for small ratios of temperature and bending rigidity, which is the biologically relevant limit. Simulations show very good quantitative agreement with our analytical result. The analysis of the correlation functions contributing to the diffusion coefficient shows that the correlations between the stochastic force of the protein and the response in the membrane shape are responsible for the reduction. Our quantitative analysis of the membrane height correlation spectrum shows an influence of the protein-membrane interaction causing a distinctly altered wave-vector dependence compared to a free membrane. Furthermore, the time correlations exhibit the two relevant time scales of the system: that of membrane fluctuations and that of lateral protein diffusion with the latter typically much longer than the former. We argue that the analysis of the long-time decay of membrane height correlations can thus provide a new means to determine the effective diffusion coefficient of proteins in the membrane.

Multiscale simulation of protein mediated membrane remodeling

Proteins interacting with membranes can result in substantial membrane deformations and curvatures. This effect is known in its broadest terms as membrane remodeling. This review article will survey current multiscale simulation methodologies that have been employed to examine protein mediated membrane remodeling.

Calculation of free energies in fluid membranes subject to heterogeneous curvature fields

We present a computational methodology for incorporating thermal effects and calculating relative free energies for elastic fluid membranes subject to spatially dependent intrinsic curvature fields using the method of thermodynamic integration. Based on a simple model for the intrinsic curvature imposed only in a localized region of the membrane, we employ thermodynamic integration to calculate the free-energy change as a function of increasing strength of the intrinsic curvature field and a thermodynamic cycle to compute free-energy changes for different sizes of the localized region. By explicitly computing the free-energy changes and by quantifying the loss of entropy accompanied with increasing membrane deformation, we show that the membrane stiffness increases with increasing intrinsic field, thereby, renormalizing the membrane bending rigidity. The second main conclusion of this work is that the entropy of the membrane decreases with increasing size of the localized region subject to the curvature field. Our results help to quantify the free-energy change when a planar membrane deforms under the influence of curvature-inducing proteins at a finite temperature.

Hybrid elastic and discrete-particle approach to biomembrane dynamics with application to the mobility of curved integral membrane proteins

We introduce a simulation strategy to consistently couple continuum biomembrane dynamics to the motion of discrete biological macromolecules residing within or on the membrane. The methodology is used to study the diffusion of integral membrane proteins that impart a curvature on the bilayer surrounding them. Such proteins exhibit a substantial reduction in diffusion coefficient relative to "flat" proteins; this effect is explained by elementary hydrodynamic considerations.

Hybrid coarse-graining approach for lipid bilayers at large length and time scales

A hybrid analytic-systematic (HAS) coarse-grained (CG) lipid model is developed and employed in a large-scale simulation of a liposome. The methodology is termed hybrid analytic-systematic because one component of the interaction between CG sites is variationally determined from the multiscale coarse-graining (MS-CG) methodology, whereas the remaining component utilizes an analytic potential. The systematic component models the in-plane center-of-mass interaction of the lipids as determined from an atomistic-level MD simulation of a bilayer. The analytic component is based on the well-known Gay-Berne ellipsoid-of-revolution liquid-crystal model and is designed to model the highly anisotropic interactions at a highly coarse-grained level. The HAS CG approach is the first step in an "aggressive" CG methodology designed to model multicomponent biological membranes at very large length and time scales.

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.

Cytoskeleton mediated effective elastic properties of model red blood cell membranes

The plasma membrane of human red blood cells consists of a lipid bilayer attached to a regular network of underlying cytoskeletal polymers. We model this system at a dynamic coarse-grained level, treating the bilayer as an elastic sheet and the cytoskeletal network as a series of phantom entropic springs. In contrast to prior simulation efforts, we explicitly account for dynamics of the cytoskeletal network, both via motion of the protein anchors that attach the cytoskeleton to the bilayer and through breaking and reconnection of individual cytoskeletal filaments. Simulation results are explained in the context of a simple mean field percolation model and comparison is made to experimental measurements of red blood cell fluctuation amplitudes.

Landscape of finite-temperature equilibrium behaviour of curvature-inducing proteins on a bilayer membrane explored using a linearized elastic free energy model

Using a recently developed multiscale simulation methodology, we describe the equilibrium behaviour of bilayer membranes under the influence of curvature-inducing proteins using a linearized elastic free energy model. In particular, we describe how the cooperativity associated with a multitude of protein-membrane interactions and protein diffusion on a membrane-mediated energy landscape elicits emergent behaviour in the membrane phase. Based on our model simulations, we predict that, depending on the density of membrane-bound proteins and the degree to which a single protein molecule can induce intrinsic mean curvature in the membrane, a range of membrane phase behaviour can be observed including two different modes of vesicle-bud nucleation and repressed membrane undulations. A state diagram as a function of experimentally tunable parameters to classify the underlying states is proposed.

Analysis and interpretation of two-dimensional single-particle tracking microscopy measurements: effect of local surface roughness

Methodological advances in light microscopy have made it possible to record the motions of individual lipid and protein molecules resident in the membrane of living cells down to the nanometer level of precision in the x, y plane. Such measurement of a single molecule's trajectory for a sufficiently long period of time or the measurement of multiple molecules' trajectories for a shorter period of time can in principle provide the necessary information to derive the particle's macroscopic two-dimensional-diffusion coefficient-a quantity of vital biological interest. However, one drawback of the light microscopy procedures used in such experiments is their relatively poor discriminatory capability for determining spatial differences along the z axis in comparison to those in the x, y plane. In this study we used computer simulation to examine the likely effect of local surface roughness over the nanometer to micrometer scale on the determination of diffusion constants in the membrane bilayer by the use of such optical-microscope-based single-particle tracking (SPT) procedures. We specifically examined motion of a single molecule along (i) a locally planar and (ii) a locally rough surface. Our results indicate a need for caution in applying overly simplistic analytical strategies to the analysis of data from SPT measurements and provide upper and lower bounds for the likely degree of error introduced on the basis of surface roughness effects alone. Additionally we present an empirical method based on an autocorrelation function approach that may prove useful in identifying the existence of surface roughness and give some idea of its extent.

Spectrin folding versus unfolding reactions and RBC membrane stiffness

Spectrin (Sp), a key component of the erythrocyte membrane, is routinely stretched to near its fully folded contour length during cell deformations. Such dynamic loading may induce domain unfolding as suggested by recent experiments. Herein we develop a model to describe the folding/unfolding of spectrin during equilibrium or nonequilibrium extensions. In both cases, our model indicates that there exists a critical extension beyond which unfolding occurs. We further deploy this model, together with a three-dimensional model of the junctional complex in the erythrocyte membrane, to explore the effect of Sp unfolding on the membrane's mechanical properties, and on the thermal fluctuation of membrane-attached beads. At large deformations our results show a distinctive strain-induced unstiffening behavior, manifested in the slow decrease of the shear modulus, and accompanied by an increase in bead fluctuation. Bead fluctuation is also found to be influenced by mode switching, a phenomenon predicted by our three-dimensional model. The amount of stiffness reduction, however, is modest compared with that reported in experiments. A possible explanation for the discrepancy is the occurrence of spectrin head-to-head disassociation which is also included within our modeling framework and used to analyze bead motion as observed via experiment.

Diffusion on ruffled membrane surfaces

We present a position Langevin equation for overdamped particle motion on rough two-dimensional surfaces. A Brownian dynamics algorithm is suggested to evolve this equation numerically, allowing for the prediction of effective (projected) diffusion coefficients over corrugated surfaces. In the case of static surface roughness, we find that a simple area-scaling prediction for the projected diffusion coefficient leads to seemingly quantitative agreement with numerical results. To study the effect of dynamic surface evolution on the diffusive process, we consider particle diffusion over a thermally fluctuating elastic membrane. Surface fluctuation has the effect of increasing the effective diffusivity toward a limiting annealed-surface value discussed previously. We argue that protein motion over cell surfaces spans a variety of physical regimes, making it impossible to identify a single approximation scheme appropriate to all measurements of interest.

A hybrid model for erythrocyte membrane: a single unit of protein network coupled with lipid bilayer

To investigate the nanomechanics of the erythrocyte membrane we developed a hybrid model that couples the actin-spectrin network to the lipid bilayer. This model features a Fourier space Brownian dynamics model of the bilayer, a Brownian dynamics model of the actin protofilament, and a modified wormlike-chain model of the spectrin (including a cable-dynamics model to predict the oscillation in tension). This model enables us to predict the nanomechanics of single or multiple units of the protein network, the lipid bilayer, and the effect of their interactions. The present work is focused on the attitude of the actin protofilament at the equilibrium states coupled with the elevations of the lipid bilayer through their primary linkage at the suspension complex in deformations. Two different actin-spectrin junctions are considered at the junctional complex. With a point-attachment junction, large pitch angles and bifurcation of yaw angles are predicted. Thermal fluctuations at bifurcation may lead to mode-switching, which may affect the network and the physiological performance of the membrane. In contrast, with a wrap-around junction, pitch angles remain small, and the occurrence of bifurcation is greatly reduced. These simulations suggest the importance of three-dimensional molecular junctions and the lipid bilayer/protein network coupling on cell membrane mechanics.

Membrane remodeling from N-BAR domain interactions: insights from multi-scale simulation

Liposome remodeling processes (e.g., vesiculation and tubulation) due to N-BAR domain interactions with the lipid bilayer are explored with a multi-scale simulation approach. Results from atomistic-level molecular dynamics simulations of membrane binding to the concave face of N-BAR domains are used along with discretized mesoscopic field-theoretic simulations to examine how the spontaneous curvature fields generated by N-BAR domains result in membrane remodeling. It is found that tubulation can be generated by anisotropic N-BAR spontaneous curvature fields, whereas vesiculation is only observed with isotropic N-BAR spontaneous curvature fields at high density. The results of the multi-scale simulations provide insight into recent experimental observations.

Hybrid simulations of lateral diffusion in fluctuating membranes

In this paper we introduce a method to simulate lateral diffusion of inclusions in a fluctuating membrane. The regarded systems are governed by two dynamic processes: the height fluctuations of the membrane and the diffusion of the inclusion along the membrane. While membrane fluctuations can be expressed in terms of a dynamic equation which follows from the Helfrich Hamiltonian, the dynamics of the diffusing particle is described by a Langevin or Smoluchowski equation. In the latter equations, the curvature of the surface needs to be accounted for, which makes particle diffusion a function of membrane fluctuations. In our scheme these coupled dynamic equations, the membrane equation and the Langevin equation for the particle, are numerically integrated to simulate diffusion in a membrane. The simulations are used to study the ratio of the diffusion coefficient projected on a flat plane and the intramembrane diffusion coefficient for the case of free diffusion. We compare our results with recent analytical results that employ a preaveraging approximation and analyze the validity of this approximation. A detailed simulation study of the relevant correlation functions reveals a surprisingly large range where the approximation is applicable.

Thermal undulations of lipid bilayers relax by intermonolayer friction at submicrometer length scales

The time correlation functions of the thermal undulations of a lipid membrane have been studied by molecular dynamics simulations of a coarse-grained bilayer model. We observe a double-exponential decay, with relaxation rates in good agreement with the theory by Seifert and Langer, [Europhys. Lett. 23, 71 (1993)]. Intermonolayer friction resulting from local velocity differences between the two monolayers is shown to be the dominant dissipative mechanism for fluctuations with wave lengths below approximately -0.1 microm.

Nonequilibrium membrane fluctuations driven by active proteins

We extend a model for nonthermal membrane undulations driven by active (adenosine triphosphate-dependent or light-harvesting) membrane proteins [N. Gov, Phys. Rev. Lett. 93, 268104 (2004)]. The present model accounts for the fact that proteins can diffuse laterally across the membrane surface and that individual proteins are expected to exert forces preferentially in one normal direction over the other (due to their orientation within the bilayer). The addition of these effects alters the scaling of fluctuation amplitudes with system size. Additionally, theoretical arguments and dynamic simulations both suggest that, in certain regimes, the probability distribution of fluctuation amplitudes is expected to be non-Gaussian (in contrast to thermal systems).

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.