Nonequilibrium membrane fluctuations driven by active proteins

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.

Non-Markovian Modeling of Nonequilibrium Fluctuations and Dissipation in Active Viscoelastic Biomatter

Based on a Hamiltonian that incorporates the elastic coupling between a tracer particle and the embedding active viscoelastic biomatter, we derive a generalized non-Markovian Langevin model for the nonequilibrium mechanical tracer response. Our analytical expressions for the frequency-dependent tracer response function and the tracer positional autocorrelation function agree quantitatively with experimental data for red blood cells and actomyosin networks with and without adenosine triphosphate over the entire frequency range and in particular reproduce the low-frequency violation of the fluctuation-dissipation theorem. The viscoelastic power laws, the elastic constants and effective friction coefficients extracted from the experimental data allow straightforward physical interpretation.

A new method for the calculation of functional and path integrals

This paper addresses a disconnect between the pivotal role of functional (path) integrals in modern theories, such as quantum mechanics and statistical thermodynamics, and the currently limited ability to perform the actual calculation. We present a new method for calculating functional integrals, based on a finite-element formulation, which solves all limitations of existing methods. This approach is far more robust, versatile, and powerful than the prevailing methods, thus allowing for more sophisticated computations and the study of problems that could not previously be tackled. Importantly, existing procedures, element libraries and shape functions, which have been developed throughout the years in the context of engineering analysis and partial differential equations, may be directly employed for this purpose.

Receptor-Ligand Binding: Effect of Mechanical Factors

Gaining insight into the in situ receptor-ligand binding is pivotal for revealing the molecular mechanisms underlying the physiological and pathological processes and will contribute to drug discovery and biomedical application. An important issue involved is how the receptor-ligand binding responds to mechanical stimuli. This review aims to provide an overview of the current understanding of the effect of several representative mechanical factors, such as tension, shear stress, stretch, compression, and substrate stiffness on receptor-ligand binding, wherein the biomedical implications are focused. In addition, we highlight the importance of synergistic development of experimental and computational methods for fully understanding the in situ receptor-ligand binding, and further studies should focus on the coupling effects of these mechanical factors.

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.

A continuum membrane model can predict curvature sensing by helix insertion

Protein domains, such as ENTH (epsin N-terminal homology) and BAR (bin/amphiphysin/rvs), contain amphipathic helices that drive preferential binding to curved membranes. However, predicting how the physical parameters of these domains control this 'curvature sensing' behavior is challenging due to the local membrane deformations generated by the nanoscopic helix on the surface of a large sphere. We here use a deformable continuum model that accounts for the physical properties of the membrane and the helix insertion to predict curvature sensing behavior, with direct validation against multiple experimental datasets. We show that the insertion can be modeled as a local change to the membrane's spontaneous curvature, cins0, producing excellent agreement with the energetics extracted from experiments on ENTH binding to vesicles and cylinders, and of ArfGAP helices to vesicles. For small vesicles with high curvature, the insertion lowers the membrane energy by relieving strain on a membrane that is far from its preferred curvature of zero. For larger vesicles, however, the insertion has the inverse effect, de-stabilizing the membrane by introducing more strain. We formulate here an empirical expression that accurately captures numerically calculated membrane energies as a function of both basic membrane properties (bending modulus κ and radius R) as well as stresses applied by the inserted helix (cins0 and area A_ins). We therefore predict how these physical parameters will alter the energetics of helix binding to curved vesicles, which is an essential step in understanding their localization dynamics during membrane remodeling processes.

More than just a barrier: using physical models to couple membrane shape to cell function

The correct execution of many cellular processes, such as division and motility, requires the cell to adopt a specific shape. Physically, these shapes are determined by the interplay of the plasma membrane and internal cellular driving factors. While the plasma membrane defines the boundary of the cell, processes inside the cell can result in the generation of forces that deform the membrane. These processes include protein binding, the assembly of protein superstructures, and the growth and contraction of cytoskeletal networks. Due to the complexity of the cell, relating observed membrane deformations back to internal processes is a challenging problem. Here, we review cell shape changes in endocytosis, cell adhesion, cell migration and cell division and discuss how by modeling membrane deformations we can investigate the inner working principles of the cell.

Correlated flickering of erythrocytes membrane observed with dual time resolved membrane fluctuation spectroscopy under different D-glucose concentrations

A correlated human red blood cell membrane fluctuation dependent on D-glucose concentration was found with dual time resolved membrane fluctuation spectroscopy (D-TRMFS). This new technique is a modified version of the dual optical tweezers method that has been adapted to measure the mechanical properties of red blood cells (RBCs) at distant membrane points simultaneously, enabling correlation analysis. Mechanical parameters under different D-glucose concentrations were obtained from direct membrane flickering measurements, complemented with membrane fluidity measurements using Laurdan Generalized Polarization (GP) Microscopy. Our results show an increase in the fluctuation amplitude of the lipid bilayer, and a decline in tension value, bending modulus and fluidity as D-glucose concentration increases. Metabolic mechanisms are proposed as explanations for the results.

Large-scale simulation of biomembranes incorporating realistic kinetics into coarse-grained models

Biomembranes are two-dimensional assemblies of phospholipids that are only a few nanometres thick, but form micrometre-sized structures vital to cellular function. Explicit molecular modelling of biologically relevant membrane systems is computationally expensive due to the large number of solvent particles and slow membrane kinetics. Coarse-grained solvent-free membrane models offer efficient sampling but sacrifice realistic kinetics, thereby limiting the ability to predict pathways and mechanisms of membrane processes. Here, we present a framework for integrating coarse-grained membrane models with continuum-based hydrodynamics. This framework facilitates efficient simulation of large biomembrane systems with large timesteps, while achieving realistic equilibrium and non-equilibrium kinetics. It helps to bridge between the nanometer/nanosecond spatiotemporal resolutions of coarse-grained models and biologically relevant time- and lengthscales. As a demonstration, we investigate fluctuations of red blood cells, with varying cytoplasmic viscosities, in 150-milliseconds-long trajectories, and compare kinetic properties against single-cell experimental observations.

Active Contact Forces Drive Nonequilibrium Fluctuations in Membrane Vesicles

We analyze the nonequilibrium shape fluctuations of giant unilamellar vesicles encapsulating motile bacteria. Owing to bacteria-membrane collisions, we experimentally observe a significant increase in the magnitude of membrane fluctuations at low wave numbers, compared to the well-known thermal fluctuation spectrum. We interrogate these results by numerically simulating membrane height fluctuations via a modified Langevin equation, which includes bacteria-membrane contact forces. Taking advantage of the lengthscale and timescale separation of these contact forces and thermal noise, we further corroborate our results with an approximate theoretical solution to the dynamical membrane equations. Our theory and simulations demonstrate excellent agreement with nonequilibrium fluctuations observed in experiments. Moreover, our theory reveals that the fluctuation-dissipation theorem is not broken by the bacteria; rather, membrane fluctuations can be decomposed into thermal and active components.

Spontaneous cellular vibratory motions of osteocytes are regulated by ATP and spectrin network

Vibration at high frequency has been demonstrated to be anabolic for bone and embedded osteocytes. The response of osteocytes to vibration is frequency-dependent, but the mechanism remains unclear. Our previous computational study using an osteocyte finite element model has predicted a resonance effect involving in the frequency-dependent response of osteocytes to vibration. However, the cellular spontaneous vibratory motion of osteocytes has not been confirmed. In the present study, the cellular vibratory motions (CVM) of osteocytes were recorded by a custom-built digital holographic microscopy and quantitatively analyzed. The roles of ATP and spectrin network in the CVM of osteocytes were studied. Results showed the MLO-Y4 osteocytes displayed dynamic vibratory motions with an amplitude of ~80 nm, which is relied both on the ATP content and spectrin network. Spectrum analysis showed several frequency peaks in CVM of MLO-Y4 osteocytes at 30 Hz, 39 Hz, 83 Hz and 89 Hz. These peak frequencies are close to the commonly used effective frequencies in animal training and in-vitro cell experiments, and show a correlation with the computational predictions of the osteocyte finite element model. These results implicate that osteocytes are dynamic and the cellular dynamic motion is involved in the cellular mechanotransduction of vibration.

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.

Equilibrium Modeling of the Mechanics and Structure of the Cancer Glycocalyx

The glycocalyx is a thick coat of proteins and carbohydrates on the outer surface of all eukaryotic cells. Overproduction of large, flexible or rod-like biopolymers, including hyaluronic acid and mucins, in the glycocalyx strongly correlates with the aggression of many cancer types. However, theoretical frameworks to predict the effects of these changes on cancer cell adhesion and other biophysical processes remain limited. Here, we propose a detailed modeling framework for the glycocalyx incorporating important physical effects of biopolymer flexibility, excluded volume, counterion mobility, and coupled membrane deformations. Because mucin and hyaluronic biopolymers are proposed to extend and rigidify depending on the extent of their decoration with side chains, we propose and consider two limiting cases for the structural elements of the glycocalyx: stiff beams and flexible chains. Simulations predict the mechanical response of the glycocalyx to compressive loads, which are imposed on cells residing in the highly confined spaces of the solid tumor or invaded tissues. Notably, the shape of the mechanical response transitions from hyperbolic to sigmoidal for more rod-like glycocalyx elements. These mechanical responses, along with the corresponding equilibrium protein organizations and membrane topographies, are summarized to aid in hypothesis generation and the evaluation of future experimental measurements. Overall, the modeling framework developed provides a theoretical basis for understanding the physical biology of the glycocalyx in human health.

Statistical Mechanics of an Elastically Pinned Membrane: Static Profile and Correlations

The relation between thermal fluctuations and the mechanical response of a free membrane has been explored in great detail, both theoretically and experimentally. However, understanding this relationship for membranes locally pinned by proteins is significantly more challenging. Given that the coupling of the membrane to the cell cytoskeleton, to the extracellular matrix, and to other internal structures is crucial for the regulation of a number of cellular processes, understanding the role of the pinning is of great interest. In this manuscript, we consider a single protein (elastic spring of a finite rest length) pinning a membrane modeled in the Monge gauge. First, we determine the Green's function for the system and complement this approach by the calculation of the mode-coupling coefficients for the plane wave expansion and the orthonormal fluctuation modes, in turn building a set of tools for numerical and analytic studies of a pinned membrane. Furthermore, we explore static correlations of the free and the pinned membrane, as well as the membrane shape, showing that all three are mutually interdependent and have an identical long-range behavior characterized by the correlation length. Interestingly, the latter displays a nonmonotonic behavior as a function of membrane tension. Importantly, exploiting these relations allows for the experimental determination of the elastic parameters of the pinning. Last but not least, we calculate the interaction potential between two pinning sites and show that even in the absence of the membrane deformation, the pinnings will be subject to an attractive force because of changes in membrane fluctuations.

Living Matter: Mesoscopic Active Materials

An introduction to the physical properties of living active matter at the mesoscopic scale (tens of nanometers to micrometers) and their unique features compared with "dead," nonactive matter is presented. This field of research is increasingly denoted as "biological physics" where physics includes chemical physics, soft matter physics, hydrodynamics, mechanics, and the related engineering sciences. The focus is on the emergent properties of these systems and their collective behavior, which results in active self-organization and how they relate to cellular-level biological function. These include locomotion (cell motility and migration) forces that give rise to cell division, the growth and form of cellular assemblies in development, the beating of heart cells, and the effects of mechanical perturbations such as shear flow (in the bloodstream) or adhesion to other cells or tissues. An introduction to the fundamental concepts and theory with selected experimental examples related to the authors' own research is presented, including red-blood-cell membrane fluctuations, motion of the nucleus within an egg cell, self-contracting acto-myosin gels, and structure and beating of heart cells (cardiomyocytes), including how they can be driven by an oscillating, mechanical probe.

Broken detailed balance and non-equilibrium dynamics in living systems: a review

Living systems operate far from thermodynamic equilibrium. Enzymatic activity can induce broken detailed balance at the molecular scale. This molecular scale breaking of detailed balance is crucial to achieve biological functions such as high-fidelity transcription and translation, sensing, adaptation, biochemical patterning, and force generation. While biological systems such as motor enzymes violate detailed balance at the molecular scale, it remains unclear how non-equilibrium dynamics manifests at the mesoscale in systems that are driven through the collective activity of many motors. Indeed, in several cellular systems the presence of non-equilibrium dynamics is not always evident at large scales. For example, in the cytoskeleton or in chromosomes one can observe stationary stochastic processes that appear at first glance thermally driven. This raises the question how non-equilibrium fluctuations can be discerned from thermal noise. We discuss approaches that have recently been developed to address this question, including methods based on measuring the extent to which the system violates the fluctuation-dissipation theorem. We also review applications of this approach to reconstituted cytoskeletal networks, the cytoplasm of living cells, and cell membranes. Furthermore, we discuss a more recent approach to detect actively driven dynamics, which is based on inferring broken detailed balance. This constitutes a non-invasive method that uses time-lapse microscopy data, and can be applied to a broad range of systems in cells and tissue. We discuss the ideas underlying this method and its application to several examples including flagella, primary cilia, and cytoskeletal networks. Finally, we briefly discuss recent developments in stochastic thermodynamics and non-equilibrium statistical mechanics, which offer new perspectives to understand the physics of living systems.

Tracking fast cellular membrane dynamics with sub-nm accuracy in the normal direction

Cellular membranes are important biomaterials with highly dynamic structures. Membrane dynamics plays an important role in numerous cellular processes, but precise tracking it is challenging due to the lack of tools with a highly sensitive and fast detection capability. Here we demonstrate a broad bandwidth optical imaging technique to measure cellular membrane displacements in the normal direction at sub-nm level detection limits and 20 μs temporal resolution (1 Hz-50 kHz). This capability allows us to study the intrinsic cellular membrane dynamics over a broad temporal and spatial spectrum. We measured the nanometer-scale stochastic fluctuations of the plasma membrane of HEK-293 cells, and found them to be highly dependent on the cytoskeletal structure of the cells. By analyzing the fluctuations, we further determine the mechanical properties of the cellular membranes. We anticipate that the method will contribute to the understanding of the basic cellular processes, and applications, such as mechanical phenotyping of cells at the single-cell level.

Nonequilibrium dynamics of probe filaments in actin-myosin networks

Active dynamic processes of cells are largely driven by the cytoskeleton, a complex and adaptable semiflexible polymer network, motorized by mechanoenzymes. Small dimensions, confined geometries, and hierarchical structures make it challenging to probe dynamics and mechanical response of such networks. Embedded semiflexible probe polymers can serve as nonperturbing multiscale probes to detect force distributions in active polymer networks. We show here that motor-induced forces transmitted to the probe polymers are reflected in nonequilibrium bending dynamics, which we analyze in terms of spatial eigenmodes of an elastic beam under steady-state conditions. We demonstrate how these active forces induce correlations among the mode amplitudes, which furthermore break time-reversal symmetry. This leads to a breaking of detailed balance in this mode space. We derive analytical predictions for the magnitude of resulting probability currents in mode space in the white-noise limit of motor activity. We relate the structure of these currents to the spatial profile of motor-induced forces along the probe polymers and provide a general relation for observable currents on two-dimensional hyperplanes.

Direct Cytoskeleton Forces Cause Membrane Softening in Red Blood Cells

Erythrocytes are flexible cells specialized in the systemic transport of oxygen in vertebrates. This physiological function is connected to their outstanding ability to deform in passing through narrow capillaries. In recent years, there has been an influx of experimental evidence of enhanced cell-shape fluctuations related to metabolically driven activity of the erythroid membrane skeleton. However, no direct observation of the active cytoskeleton forces has yet been reported to our knowledge. Here, we show experimental evidence of the presence of temporally correlated forces superposed over the thermal fluctuations of the erythrocyte membrane. These forces are ATP-dependent and drive enhanced flickering motions in human erythrocytes. Theoretical analyses provide support for a direct force exerted on the membrane by the cytoskeleton nodes as pulses of well-defined average duration. In addition, such metabolically regulated active forces cause global membrane softening, a mechanical attribute related to the functional erythroid deformability.

Minimal model of stochastic athermal systems: origin of non-Gaussian noise

For a wide class of stochastic athermal systems, we derive Langevin-like equations driven by non-Gaussian noise, starting from master equations and developing a new asymptotic expansion. We found an explicit condition whereby the non-Gaussian properties of the athermal noise become dominant for tracer particles associated with both thermal and athermal environments. Furthermore, we derive an inverse formula to infer microscopic properties of the athermal bath from the statistics of the tracer particle. We apply our formulation to a granular motor under viscous friction and analytically obtain the angular velocity distribution function. Our theory demonstrates that the non-Gaussian Langevin equation is the minimal model of athermal systems.

Nonequilibrium collective dynamics in photoexcited lipid multilayers by time resolved diffuse x-ray scattering

We study the nonequilibrium shape fluctuations in fluorescence labeled phospholipid multibilayers composed of the model lipid DOPC and the well-known lipid dye Texas red, driven out of equilibrium by short laser pulses. The temporal evolution of the lipid bilayer undulations after excitation was recorded by time resolved x-ray diffraction. Already at moderate peak intensities (Pp≤10(5)  W/cm2), pulsed laser illumination leads to significant changes of the undulation modes in a well-defined lateral wavelength band. The observed phenomena evolve on nano- to microsecond time scales after optical excitation, and can be described in terms of a modulation instability in the lipid multilamellar stack.

Coarse-grain simulations of active molecular machines in lipid bilayers

A coarse-grain method for simulations of the dynamics of active protein inclusions in lipid bilayers is described. It combines the previously proposed hybrid simulations of bilayers [M.-J. Huang, R. Kapral, A. S. Mikhailov, and H.-Y. Chen, J. Chem. Phys. 137, 055101 (2012)], based on molecular dynamics for the lipids and multi-particle collision dynamics for the solvent, with an elastic-network description of active proteins. The method is implemented for a model molecular machine which performs active conformational motions induced by ligand binding and its release after reaction. The situation characteristic for peripheral membrane proteins is considered. Statistical investigations of the effects of single active or passive inclusions on the shape of the membrane are carried out. The results show that the peripheral machine produces asymmetric perturbations of the thickness of two leaflets of the membrane. It also produces a local saddle in the midplane height of the bilayer. Analysis of the power spectrum of the fluctuations of the membrane midplane shows that the conformational motion of the machine perturbs these membrane fluctuations. The hydrodynamic lipid flows induced by cyclic conformational changes in the machine are analyzed. It is shown that such flows are long-ranged and should provide an additional important mechanism for interactions between active inclusions in biological membranes.

Optical phase nanoscopy in red blood cells using low-coherence spectroscopy

We propose a low-coherence spectral-domain phase microscopy (SDPM) system for accurate quantitative phase measurements in red blood cells (RBCs) for the prognosis and monitoring of disease conditions that affect the visco-elastic properties of RBCs. Using the system, we performed time-recordings of cell membrane fluctuations, and compared the nano-scale fluctuation dynamics of healthy and glutaraldehyde-treated RBCs. Glutaraldehyde-treated RBCs possess lower amplitudes of fluctuations, reflecting an increased membrane stiffness. To demonstrate the ability of our system to measure fluctuations of lower amplitudes than those measured by the commonly used holographic phase microscopy techniques, we also constructed wide-field digital interferometry (WFDI) system and compared the performances of both systems. Due to its common-path geometry, the optical-path-delay stability of SDPM was found to be less than 0.3 nm in liquid environment, at least three times better than WFDI under the same conditions. In addition, due to the compactness of SDPM and its inexpensive and robust design, the system possesses a high potential for clinical applications.

Effective tension and fluctuations in active membranes

We calculate the fluctuation spectrum of the shape of a lipid vesicle or cell exposed to a nonthermal source of noise. In particular, we take constraints on the membrane area and the volume of fluid that it encapsulates into account when obtaining expressions for the dependency of the membrane tension on the noise. We then investigate three possible origins of the nonthermal noise taken from the literature: A direct force, which models an external medium pushing on the membrane, a curvature force, which models a fluctuating spontaneous curvature, and a permeation force coming from an active transport of fluid through the membrane. For the direct force and curvature force cases, we compare our results to existing experiments on active membranes.

Holographic tissue dynamics spectroscopy

Tissue dynamics spectroscopy uses digital holography as a coherence gate to extract depth-resolved quasi-elastic dynamic light scattering from inside multicellular tumor spheroids. The temporal speckle contrast provides endogenous dynamical images of proliferating and hypoxic or necrotic tissues. Fluctuation spectroscopy similar to diffusing wave spectroscopy is performed on the dynamic speckle to generate tissue-response spectrograms that track time-resolved changes in intracellular motility in response to environmental perturbations. The spectrograms consist of several frequency bands that range from 0.005 to 5 Hz. The fluctuation spectral density and temporal autocorrelations show the signature of constrained anomalous diffusion, but with large fluctuation amplitudes caused by active processes far from equilibrium. Differences in the tissue-response spectrograms between the proliferating outer shell and the hypoxic inner core differentiate normal from starved conditions. The differential spectrograms provide an initial library of tissue-response signatures to environmental conditions of temperature, osmolarity, pH, and serum growth factors.

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.

Biophysics of malarial parasite exit from infected erythrocytes

Upon infection and development within human erythrocytes, P. falciparum induces alterations to the infected RBC morphology and bio-mechanical properties to eventually rupture the host cells through parasitic and host derived proteases of cysteine and serine families. We used previously reported broad-spectrum inhibitors (E64d, EGTA-AM and chymostatin) to inhibit these proteases and impede rupture to analyze mechanical signatures associated with parasite escape. Treatment of late-stage iRBCs with E64d and EGTA-AM prevented rupture, resulted in no major RBC cytoskeletal reconfiguration but altered schizont morphology followed by dramatic re-distribution of three-dimensional refractive index (3D-RI) within the iRBC. These phenotypes demonstrated several-fold increased iRBC membrane flickering. In contrast, chymostatin treatment showed no 3D-RI changes and caused elevated fluctuations solely within the parasitophorous vacuole. We show that E64d and EGTA-AM supported PV breakdown and the resulting elevated fluctuations followed non-Gaussian pattern that resulted from direct merozoite impingement against the iRBC membrane. Optical trapping experiments highlighted reduced deformability of the iRBC membranes upon rupture-arrest, more specifically in the treatments that facilitated PV breakdown. Taken together, our experiments provide novel mechanistic interpretations on the role of parasitophorous vacuole in maintaining the spherical schizont morphology, the impact of PV breakdown on iRBC membrane fluctuations leading to eventual parasite escape and the evolution of membrane stiffness properties of host cells in which merozoites were irreversibly trapped, recourse to protease inhibitors. These findings provide a comprehensive, previously unavailable, body of information on the combined effects of biochemical and biophysical factors on parasite egress from iRBCs.

Effective temperature of red-blood-cell membrane fluctuations

Biologically driven nonequilibrium fluctuations are often characterized by their non-Gaussianity or by an "effective temperature", which is frequency dependent and higher than the ambient temperature. We address these two measures theoretically by examining a randomly kicked particle, with a variable number of kicking motors, and show how these two indicators of nonequilibrium behavior can contradict. Our results are compared with new experiments on shape fluctuations of red-blood cell membranes, and demonstrate how the physical nature of the motors in this system can be revealed using these global measures of nonequilibrium.

Effective 3D viscoelasticity of red blood cells measured by diffraction phase microscopy

We present optical measurements of nanoscale red blood cell fluctuations obtained by highly sensitive quantitative phase imaging. These spatio-temporal fluctuations are modeled in terms of the bulk viscoelastic response of the cell. Relating the displacement distribution to the storage and loss moduli of the bulk has the advantage of incorporating all geometric and cortical effects into a single effective medium behavior. The results on normal cells indicate that the viscous modulus is much larger than the elastic one throughout the entire frequency range covered by the measurement, indicating fluid behavior.

Electromagnetic cellular interactions

Chemical and electrical interaction within and between cells is well established. Just the opposite is true about cellular interactions via other physical fields. The most probable candidate for an other form of cellular interaction is the electromagnetic field. We review theories and experiments on how cells can generate and detect electromagnetic fields generally, and if the cell-generated electromagnetic field can mediate cellular interactions. We do not limit here ourselves to specialized electro-excitable cells. Rather we describe physical processes that are of a more general nature and probably present in almost every type of living cell. The spectral range included is broad; from kHz to the visible part of the electromagnetic spectrum. We show that there is a rather large number of theories on how cells can generate and detect electromagnetic fields and discuss experimental evidence on electromagnetic cellular interactions in the modern scientific literature. Although small, it is continuously accumulating.

Metabolic remodeling of the human red blood cell membrane

The remarkable deformability of the human red blood cell (RBC) results from the coupled dynamic response of the phospholipid bilayer and the spectrin molecular network. Here we present quantitative connections between spectrin morphology and membrane fluctuations of human RBCs by using dynamic full-field laser interferometry techniques. We present conclusive evidence that the presence of adenosine 5'-triphosphate (ATP) facilitates non-equilibrium dynamic fluctuations in the RBC membrane that are highly correlated with the biconcave shape of RBCs. Spatial analysis of the fluctuations reveals that these non-equilibrium membrane vibrations are enhanced at the scale of spectrin mesh size. Our results indicate that the dynamic remodeling of the coupled membranes powered by ATP results in non-equilibrium membrane fluctuations manifesting from both metabolic and thermal energies and also maintains the biconcave shape of RBCs.

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.

Systematic multiscale simulation of membrane protein systems

Current multiscale simulation approaches for membrane protein systems vary depending on their degree of connection to the underlying molecular scale interactions. Various approaches have been developed that include such information into coarse-grained models of both the membrane and the proteins. By contrast, other approaches employ parameterizations obtained from experimental data. Mesoscopic models operate at larger scales and have also been employed to examine membrane remodeling, protein inclusions, and ion channel gating. When bridged together such that molecular-level information is propagated between the different scales, a systematic multiscale methodology for membrane protein systems can be achieved.

Refractive index maps and membrane dynamics of human red blood cells parasitized by Plasmodium falciparum

Parasitization by malaria-inducing Plasmodium falciparum leads to structural, biochemical, and mechanical modifications to the host red blood cells (RBCs). To study these modifications, we investigate two intrinsic indicators: the refractive index and membrane fluctuations in P. falciparum-invaded human RBCs (Pf-RBCs). We report experimental connections between these intrinsic indicators and pathological states. By employing two noninvasive optical techniques, tomographic phase microscopy and diffraction phase microscopy, we extract three-dimensional maps of refractive index and nanoscale cell membrane fluctuations in isolated RBCs. Our systematic experiments cover all intraerythrocytic stages of parasite development under physiological and febrile temperatures. These findings offer potential, and sufficiently general, avenues for identifying, through cell membrane dynamics, pathological states that cause or accompany human diseases.

The stochastic dynamics of filopodial growth

A filopodium is a cytoplasmic projection, exquisitely built and regulated, which extends from the leading edge of the migrating cell, exploring the cell's neighborhood. Commonly, filopodia grow and retract after their initiation, exhibiting rich dynamical behaviors. We model the growth of a filopodium based on a stochastic description which incorporates mechanical, physical, and biochemical components. Our model provides a full stochastic treatment of the actin monomer diffusion and polymerization of each individual actin filament under stress of the fluctuating membrane. We investigated the length distribution of individual filaments in a growing filopodium and studied how it depends on various physical parameters. The distribution of filament lengths turned out to be narrow, which we explained by the negative feedback created by the membrane load and monomeric G-actin gradient. We also discovered that filopodial growth is strongly diminished upon increasing retrograde flow, suggesting that regulating the retrograde flow rate would be a highly efficient way to control filopodial extension dynamics. The filopodial length increases as the membrane fluctuations decrease, which we attributed to the unequal loading of the membrane force among individual filaments, which, in turn, results in larger average polymerization rates. We also observed significant diffusional noise of G-actin monomers, which leads to smaller G-actin flux along the filopodial tube compared with the prediction using the diffusion equation. Overall, partial cancellation of these two fluctuation effects allows a simple mean field model to rationalize most of our simulation results. However, fast fluctuations significantly renormalize the mean field model parameters. The biological significance of our filopodial model and avenues for future development are also discussed.

Fluctuations of coupled fluid and solid membranes with application to red blood cells

The fluctuation spectra and the intermembrane interaction of two membranes at a fixed average distance are investigated. Each membrane can either be a fluid or a solid membrane, and in isolation, its fluctuations are described by a bare or a wave-vector-dependent bending modulus, respectively. The membranes interact via their excluded-volume interaction; the average distance is maintained by an external, homogeneous pressure. For strong coupling, the fluctuations can be described by a single, effective membrane that combines the elastic properties. For weak coupling, the fluctuations of the individual, noninteracting membranes are recovered. The case of a composite membrane consisting of one fluid and one solid membrane can serve as a microscopic model for the plasma membrane and cytoskeleton of the red blood cell. We find that, despite the complex microstructure of bilayers and cytoskeletons in a real cell, the fluctuations with wavelengths lambda greater, similar 400 nm are well described by the fluctuations of a single, polymerized membrane (provided that there are no inhomogeneities of the microstructure). The model is applied to the fluctuation data of discocytes ("normal" red blood cells), a stomatocyte, and an echinocyte. The elastic parameters of the membrane and an effective temperature that quantifies active, metabolically driven fluctuations are extracted from the experiments.

Coherence properties of red blood cell membrane motions

We use a highly sensitive, noncontact, optical interferometric technique to quantify the red blood cell membrane fluctuations at the nanometer and millisecond scales. The results reveal significant properties of both temporal and spatial coherence associated with the membrane dynamics. We show that these correlations can be accounted for by the viscoelastic properties of the cell membrane. From this measurement, we extract the loss and storage moduli associated with the membrane and find a crossover frequency at which the buffer viscosity seems to become dominant.

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.

Active membranes studied by X-ray scattering

In view of recent theories of "active" membranes, we have studied multilamellar phospholipid membrane stacks with reconstituted transmembrane protein bacteriorhodopsin (BR) under different illumination conditions by X-ray scattering. The light-active protein is considered as an active constituent which drives the system out of equilibrium and is predicted to change the collective fluctuation properties of the membranes. Using X-ray reflectivity, X-ray non-specular (diffuse) scattering, and grazing incidence scattering, we find no detectable change in the scattering curves when changing the illumination condition. In particular the intermembrane spacing d remains constant, after eliminating hydration-related artifacts by design of a suitable sample environment. The absence of any observable non-equilibrium effects in the experimental window is discussed in view of the relevant parameters and recent theories.

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.

Mapping correlated membrane pulsations and fluctuations in human cells

The cell membrane and cytoskeleton are dynamic structures that are strongly influenced by the thermo-mechanical background in addition to biologically driven mechanical processes. We used atomic force microscopy (AFM) to measure the local membrane motion of human foreskin fibroblasts (HFFs) which were found to be governed by random and non-random correlated mechanical processes. Interphase cells displayed distinct membrane pulsations in which the membrane was observed to slowly contract upwards followed by a recovery to its initial position. These pulsations occurred one to three times per minute with variable amplitudes (20-100 pN) separated by periods of random baseline fluctuations with amplitudes of <20 pN. Cells were exposed to actin and microtubule (MT) destabilizing drugs and induced into early apoptosis. Mechanical pulsations (20-80 pN) were not prevented by actin or MT depolymerization but were prevented in early apoptotic cells which only displayed small amplitude baseline fluctuations (<20 pN). Correlation analysis revealed that the cell membrane motion is largely random; however several non-random processes, with time constants varying between approximately 2 and 35 s are present. Results were compared to measured cardiomyocyte motion which was well defined and highly correlated. Employing automated positioning of the AFM tip, interphase HFF correlation time constants were also mapped over a 10 microm2 area above the nucleus providing some insights into the spatial variability of membrane correlations. Here, we are able to show that membrane pulsations and fluctuations can be linked to physiological state and cytoskeletal dynamics through distinct sets of correlation time constants in human cells.

Analysis of shape, fluctuations, and dynamics in intermembrane junctions

A dynamic-elastic model for weakly adhered intermembrane junctions is presented. Helfrich membrane energetics coupled to hydrodynamic modes of the surrounding solvent reproduce the average shape, fluctuations, and dynamics of these junctions as measured experimentally. Comparison between numerical results and experimental data provides the first direct measure of surface tension in these systems (0.01-0.06 dyn/cm). The measurements suggest bilayer-bilayer adhesion energetics as the dominant source of surface tension in the experimental systems.

Fluctuation spectrum of quasispherical membranes with force-dipole activity

The fluctuation spectrum of a quasispherical vesicle with active membrane proteins is calculated. The activity of the proteins is modeled as the proteins pushing on their surroundings giving rise to nonlocal force distributions. Both the contributions from the thermal fluctuations of the active protein densities and the temporal noise in the individual active force distributions of the proteins are taken into account. The noise in the individual force distributions is found to become significant at short wavelengths.

Adhesion of membranes via switchable molecules

Biomimetic membranes that adhere to a solid substrate or another interface via switchable crosslinker molecules are studied theoretically using analytical methods and Monte Carlo simulations. The flexible crosslinkers exhibit two conformations which have a different end-to-end distance and, thus, lead to different local separations of the membrane from the substrate surface. Transitions between the molecular conformations can be induced by light, electric potential, or changes in pH and lead to active shape fluctuations of the membrane and, thus, to an increased membrane roughness. The forward and backward transitions are characterized by two transition rates, omega(+) and omega(-), respectively, which define the average fraction X=omega(+) /(omega(+) + omega(-)) of + (or on) states and the mean switching rate omega = (omega(+) + omega(-)) / 2. The membrane roughness is explicitly calculated as a function of X and omega. It is shown that the interplay of active and thermal fluctuations is subtle and that it is, in general, not possible to describe the active fluctuations in terms of an effective temperature.