Electrohydrodynamic model of vesicle deformation in alternating electric fields

Compound giant unilamellar vesicles as a bio-mimetic model for electrohydrodynamics of a nucleate cell

The understanding obtained by studies on the electrohydrodynamics (EHD) of single giant unilamellar vesicles (sGUVs) has contributed significantly towards a better comprehension of the response of biological cells to electric fields. This has subsequently helped in developing technologies such as cell dielectrophoresis and cell electroporation. For nucleate eukaryotic cells though, a vesicle-in-vesicle compound giant unilamellar vesicle (cGUV) is a more appropriate bio-mimic than a sGUV. In this work, we present an improvised method for the formation of cGUVs, wherein the electrical conductivities of the inner, annular and outer regions of the cGUVs can be modified. A comprehensive experimental study is presented on the EHD of these cGUVs under weak AC fields over a wide range of frequencies, and an encouraging agreement is observed between the experiments and earlier published theoretical studies on concentric cGUVs. The spherical, prolate or oblate spheroidal deformations of a cGUV under AC electric fields depend upon the membrane electromechanical properties as well as the magnitude and direction of the electric traction at the membrane produced by the Maxwell stress that varies with the relative timescales associated with the frequency of the applied AC electric field and that of the membrane charging time and the Maxwell-Wagner relaxation time. This work establishes cGUVs as appropriate bio-mimics for conducting EHD studies relevant to eukaryotic cells.

Giant vesicles form in physiological saline and encapsulate pDNA by the modified electroformation method

Giant vesicles (GVs) are used to study the structures and functions of cells and cell membranes. Electroformation is the most commonly used method for GV preparation. However, the electroformation of GVs is hindered in highly concentrated ionic solutions, limiting their application as cell models for research under physiological conditions. In this study, giant multilayer vesicles were successfully generated in physiological saline using a modified electroformation device by adding an insulating layer between the two electrode plates. The influence of the electric frequency and strength on the electroformation of GVs in physiological saline was explored, and a possible mechanism for this improvement was assessed. It has been shown that an insulating layer between the two electrodes can improve the electroformation of GVs in physiological saline by increasing the electrical impedance, which is weakened by the saline solution, thereby restoring the reduced effective electric field strength. Furthermore, macromolecular plasmid DNA (pDNA) was successfully encapsulated in the electroformed GVs of the modified device. This modified electroformation method may be useful for generating eukaryotic cell models under physiological conditions.

Stationary shapes of axisymmetric vesicles beyond lowest-energy configurations

We conduct a systematic exploration of the energy landscape of vesicle morphologies within the framework of the Helfrich model. Vesicle shapes are determined by minimizing the elastic energy subject to constraints of constant area and volume. The results show that pressurized vesicles can adopt higher-energy spindle-like configurations that require the action of point forces at the poles. If the internal pressure is lower than the external one, multilobed shapes are predicted. We utilize our results to rationalize experimentally observed spindle shapes of giant vesicles in a uniform AC electric field.

Electromechanical interactions between cell membrane and nuclear envelope: Beyond the standard Schwan's model of biological cells

We investigate little-appreciated features of the hierarchical core-shell (CS) models of the electrical, mechanical, and electromechanical interactions between the cell membrane (CM) and nuclear envelope (NE). We first consider a simple model of an individual cell based on a coupled resistor-capacitor (Schwan model (SM)) network and show that the CM, when exposed to ac electric fields, acts as a low pass filter while the NE acts as a wide and asymmetric bandpass filter. We provide a simplified calculation for characteristic time associated with the capacitive charging of the NE and parameterize its range of behavior. We furthermore observe several new features dealing with mechanical analogs of the SM based on elementary spring-damper combinations. The chief merit of these models is that they can predict creep compliance responses of an individual cell under static stress and their effective retardation time constants. Next, we use an alternative and a more accurate CS physical model solved by finite element simulations for which geometrical cell reshaping under electromechanical stress (electrodeformation (ED)) is included in a continuum approach with spatial resolution. We show that under an electric field excitation, the elongated nucleus scales differently compared to the electrodeformed cell.

Stationary shapes of axisymmetric vesicles beyond lowest-energy configurations

We conduct a systematic exploration of the energy landscape of vesicle morphologies within the framework of the Helfrich model. Vesicle shapes are determined by minimizing the elastic energy subject to constraints of constant area and volume. The results show that pressurized vesicles can adopt higher-energy spindle-like configurations that require the action of point forces at the poles. If the internal pressure is lower than the external one, multilobed shapes are predicted. We utilize our results to rationalize the experimentally observed spindle shapes of giant vesicles in a uniform AC field.

Photomanipulation of Minimal Synthetic Cells: Area Increase, Softening, and Interleaflet Coupling of Membrane Models Doped with Azobenzene-Lipid Photoswitches

Light can effectively interrogate biological systems in a reversible and physiologically compatible manner with high spatiotemporal precision. Understanding the biophysics of photo-induced processes in bio-systems is crucial for achieving relevant clinical applications. Employing membranes doped with the photolipid azobenzene-phosphatidylcholine (azo-PC), a holistic picture of light-triggered changes in membrane kinetics, morphology, and material properties obtained from correlative studies on cell-sized vesicles, Langmuir monolayers, supported lipid bilayers, and molecular dynamics simulations is provided. Light-induced membrane area increases as high as ≈25% and a ten-fold decrease in the membrane bending rigidity is observed upon trans-to-cis azo-PC isomerization associated with membrane leaflet coupling and molecular curvature changes. Vesicle electrodeformation measurements and atomic force microscopy reveal that trans azo-PC bilayers are thicker than palmitoyl-oleoyl phosphatidylcholine (POPC) bilayers but have higher specific membrane capacitance and dielectric constant suggesting an increased ability to store electric charges across the membrane. Lastly, incubating POPC vesicles with azo-PC solutions results in the insertion of azo-PC in the membrane enabling them to become photoresponsive. All these results demonstrate that light can be used to finely manipulate the shape, mechanical and electric properties of photolipid-doped minimal cell models, and liposomal drug carriers, thus, presenting a promising therapeutic alternative for the repair of cellular disorders.

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress-tensor autocorrelation function modeled by a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation, and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly 20 times that of DPPC. Furthermore, we find that inclusion of the entire range of the dispersion interaction increases viscosity by up to 140%. The simulated viscosities are similar to experimental values obtained from the rotational dynamics of small chromophores and from the diffusion of integral membrane proteins but significantly lower than recent measurements based on the deformation of giant vesicles.

Influence of surface viscosities on the electrodeformation of a prolate viscous drop

Contaminants and other agents are often present at the interface between two fluids, giving rise to rheological properties such as surface shear and dilatational viscosities. The dynamics of viscous drops with interfacial viscosities has attracted greater interest in recent years, due to the influence of surface rheology on deformation and the surrounding flows. We investigate the effects of shear and dilatational viscosities on the electro-deformation of a viscous drop using the Taylor-Melcher leaky dielectric model. We use a large deformation analysis to derive an ordinary differential equation for the drop shape. Our model elucidates the contributions of each force to the overall deformation of the drop and reveals a rich range of dynamic behaviors that show the effects of surface viscosities and their dependence on rheological and electrical properties of the system. We also examine the physical mechanisms underlying the observed behaviors by analyzing the surface dilatation and surface deformation.

Assessing membrane material properties from the response of giant unilamellar vesicles to electric fields

Knowledge of the material properties of membranes is crucial to understanding cell viability and physiology. A number of methods have been developed to probe membranes in vitro, utilizing the response of minimal biomimetic membrane models to an external perturbation. In this review, we focus on techniques employing giant unilamellar vesicles (GUVs), model membrane systems, often referred to as minimal artificial cells because of the potential they offer to mimick certain cellular features. When exposed to electric fields, GUV deformation, dynamic response and poration can be used to deduce properties such as bending rigidity, pore edge tension, membrane capacitance, surface shear viscosity, excess area and membrane stability. We present a succinct overview of these techniques, which require only simple instrumentation, available in many labs, as well as reasonably facile experimental implementation and analysis.

A lateral electric field inhibits gel-to-fluid transition in lipid bilayers

We report evidence of lateral electric field-induced changes in the phase transition temperatures of lipid bilayers. Our atomic scale molecular dynamics simulations show that a lateral electric field increases the melting temperatures of DPPC, POPC and POPE bilayers. Remarkably, these shifts in the melting temperatures are only induced by lateral electric fields, and not normal electric fields. This mechanism could provide new mechanistic insights into lipid-lipid and lipid-protein interactions in the presence of endogenous and exogenous electric fields.

Bending Rigidity, Capacitance, and Shear Viscosity of Giant Vesicle Membranes Prepared by Spontaneous Swelling, Electroformation, Gel-Assisted, and Phase Transfer Methods: A Comparative Study

Closed lipid bilayers in the form of giant unilamellar vesicles (GUVs) are commonly used membrane models. Various methods have been developed to prepare GUVs, however it is unknown if all approaches yield membranes with the same elastic, electric, and rheological properties. Here, we combine flickering spectroscopy and electrodefomation of GUVs to measure, at identical conditions, membrane capacitance, bending rigidity and shear surface viscosity of palmitoyloleoylphosphatidylcholine (POPC) membranes formed by several commonly used preparation methods: thin film hydration (spontaneous swelling), electroformation, gel-assisted swelling using poly(vinyl alcohol) (PVA) or agarose, and phase-transfer. We find relatively similar bending rigidity value across all the methods except for the agarose hydration method. In addition, the capacitance values are similar except for vesicles prepared via PVA gel hydration. Intriguingly, membranes prepared by the gel-assisted and phase-transfer methods exhibit much higher shear viscosity compared to electroformation and spontaneous swelling, likely due to remnants of polymers (PVA and agarose) and oils (hexadecane and mineral) in the lipid bilayer structure.

Nonlinear large deformation of a spherical red blood cell induced by ultrasonic standing wave

A computational model is developed to investigate the nonlinear static deformation of a spherical (osmotically swollen) red blood cell (RBC) induced by ultrasonic standing wave. The ultrasonic standing wave can generate steady acoustic radiation stress to deform the cell, and in turn, the deformed cell reshapes the acoustic field. This is a real-time coupling problem between the acoustic field and the mechanical field. In the computational model, the acoustic radiation stress acting on the RBC membrane is modeled by adopting the nonviscous momentum flux theory. The RBC membrane is modeled as a hyperelastic shell considering the in-plane elasticity, bending elasticity, and surface tension of the membrane. The volume conservation constraint of the membrane sealing fluid is applied to ensure the osmotic balance of the membrane. To address this real-time coupling problem, the computational model is implemented by a finite element method algorithm. The numerical results are compared with the existing theoretical model and experimental data, and the strain hardening trend of the experimental data is successfully predicted, which verifies the accuracy and effectiveness of the computational model. The computational model can accurately extract the mechanical properties of cells from acoustic deformation experiments, which is helpful for the diagnosis of some human diseases.

Nonlinear Transient and Steady State Stretching of Deflated Vesicles in Flow

Membrane-bound vesicles and organelles exhibit a wide array of nonspherical shapes at equilibrium, including biconcave and tubular morphologies. Despite recent progress, the stretching dynamics of deflated vesicles is not fully understood, particularly far from equilibrium where complex nonspherical shapes undergo large deformations in flow. Here, we directly observe the transient and steady-state nonlinear stretching dynamics of deflated vesicles in extensional flow using a Stokes trap. Automated flow control is used to observe vesicle dynamics over a wide range of flow rates, shape anisotropy, and viscosity contrast. Our results show that deflated vesicle membranes stretch into highly deformed shapes in flow above a critical capillary number Ca_c1. We further identify a second critical capillary number Ca_c2, above which vesicle stretch diverges in flow. Vesicles are robust to multiple nonlinear stretch-relax cycles, evidenced by relaxation of dumbbell-shaped vesicles containing thin lipid tethers following flow cessation. An analytical model is developed for vesicle deformation in flow, which enables comparison of nonlinear steady-state stretching results with theories for different reduced volumes. Our results show that the model captures the steady-state stretching of moderately deflated vesicles; however, it underpredicts the steady-state nonlinear stretching of highly deflated vesicles. Overall, these results provide a new understanding of the nonlinear stretching dynamics and membrane mechanics of deflated vesicles in flow.

Dielectric Properties of Phosphatidylcholine Membranes and the Effect of Sugars

Simple carbohydrates are associated with the enhanced risk of cardiovascular disease and adverse changes in lipoproteins in the organism. Conversely, sugars are known to exert a stabilizing effect on biological membranes, and this effect is widely exploited in medicine and industry for cryopreservation of tissues and materials. In view of elucidating molecular mechanisms involved in the interaction of mono- and disaccharides with biomimetic lipid systems, we study the alteration of dielectric properties, the degree of hydration, and the rotational order parameter and dipole potential of lipid bilayers in the presence of sugars. Frequency-dependent deformation of cell-size unilamellar lipid vesicles in alternating electric fields and fast Fourier transform electrochemical impedance spectroscopy are applied to measure the specific capacitance of phosphatidylcholine lipid bilayers in sucrose, glucose and fructose aqueous solutions. Alteration of membrane specific capacitance is reported in sucrose solutions, while preservation of membrane dielectric properties is established in the presence of glucose and fructose. We address the effect of sugars on the hydration and the rotational order parameter for 1-palmitoyl-2-oleoyl-sn-glycero-3- phosphocholine (POPC) and 1-stearoyl-2-oleoyl-sn-glycero-3- phosphocholine (SOPC). An increased degree of lipid packing is reported in sucrose solutions. The obtained results provide evidence that some small carbohydrates are able to change membrane dielectric properties, structure, and order related to membrane homeostasis. The reported data are also relevant to future developments based on the response of lipid bilayers to external physical stimuli such as electric fields and temperature changes.

Modeling cell membrane electrodeformation by alternating electric fields

With the aim of characterizing and gaining insight into the frequency response of cells suspended in a fluid medium and deformed with a controlled alternating electric field, a continuum-based analysis is presented for modeling electrodeformation (ED) via Maxwell stress tensor (MST) calculation. Our purpose here is to apply this approach to explain the fact that the electric field anisotropy and electrical conductivity ratio Λ of the cytoplasm and the extracellular medium significantly impact the MST exerted on the cytoplasm-membrane interface. One important finding is that the modulation of electrical cues and MST force by the frequency of the applied electric field provides an extremely rich tool kit for manipulating cells. We show the extreme sensitivity of proximity-induced capacitive coupling arising concomitantly when the magnitude of the MST increases as the distance between cells is decreased and the spatial anisotropy becomes important. Moreover, our model highlights the strongly localized character of the electrostatic field effect emanating from neighboring cells and suggests the possibility of exploiting cell distribution as a powerful tool to engineer the functional performance of cell assemblies by controlling ED and capacitive coupling. We furthermore show that frequency has a significant impact on the attenuation-amplification transition of MST, suggesting that shape anisotropy has a much weaker influence on ED of the cell membrane compared to the anisotropy induced by the orientation angle itself.

Electromechanical characterization of biomimetic membranes using electrodeformation of vesicles

We describe a facile method to simultaneously measure the bending rigidity and capacitance of biomimetic lipid bilayers. Our approach utilizes the ellipsoidal deformation of quasi-spherical giant unilamellar vesicles induced by a uniform AC electric field. Vesicle shape depends on the electric field frequency and amplitude. Membrane bending rigidity can be obtained from the variation of the vesicle elongation on either field amplitude at fixed frequency or frequency at fixed field amplitude. Membrane capacitance is determined from the frequency at which the vesicle shape changes from prolate to oblate ellipsoid as the frequency is increased at a given field amplitude.

Single-cell mechanical analysis and tension quantification via electrodeformation relaxation

The mechanical behavior and cortical tension of single cells are analyzed using electrodeformation relaxation. Four types of cells, namely, MCF-10A, MCF-7, MDA-MB-231, and GBM, are studied, with pulse durations ranging from 0.01 to 10 s. Mechanical response in the long-pulse regime is characterized by a power-law behavior, consistent with soft glassy rheology resulting from unbinding events within the cortex network. In the subsecond short-pulse regime, a single timescale well describes the process and indicates the naive tensioned (prestressed) state of the cortex with minimal force-induced alteration. A mathematical model is employed and the simple ellipsoidal geometry allows for use of an analytical solution to extract the cortical tension. At the shortest pulse of 0.01 s, tensions for all four cell types are on the order of 10^{-2} N/m.

Fluctuation spectroscopy of giant unilamellar vesicles using confocal and phase contrast microscopy

A widely used method to measure the bending rigidity of bilayer membranes is fluctuation spectroscopy, which analyses the thermally-driven membrane undulations of giant unilamellar vesicles recorded with either phase-contrast or confocal microscopy. Here, we analyze the fluctuations of the same vesicle using both techniques and obtain consistent values for the bending modulus. We discuss the factors that may lead to discrepancies.

Development of transmembrane potential in concentric spherical, confocal spheroidal, and bispherical vesicles subjected to nanosecond-pulse electric field

Electroporation of concentric compound spherical and confocal spheroidal as well as eccentric compound spherical vesicles, considered to be good models for corresponding nucleate cells, are investigated with an emphasis on their response to nanosecond pulse electric field (nsPEF). Analytical models are developed for the estimation of the transmembrane potential (TMP) across the bilayers of the inner and the outer vesicles and finite-element simulations are also carried out for the eccentric case. Our calculations show that with an increase in the aspect ratio, while the TMP decreases when nsPEF is used, it increases for confocal spheroids when the pulse width is greater than the membrane charging time, leading to fully charged vesicles. Bipolar pulses are shown to effectively control the TMP for a desired time period in the nsPEF regime, and a fast decay of the TMP to zero can be achieved by judicious use of pulse polarity. The external conductivity is found to significantly influence the TMP in nsPEF, unlike millisecond pulses where its effect is insignificant. Additionally the critical electric field required to induce a TMP of 1 V at the inner vesicle is presented for different pulse widths, rise time, as well as membrane capacitance, and the TMP of the outer vesicle is found to be within limits of reversible poration. It is found that the maximum TMP has a roughly linear dependence on the outer aspect ratio of the vesicle. We also introduce a new method to obtain the particular solution to the Laplace equation for bispherical system, and it is validated with finite-element simulations. Our study on nsPEF electroporation of bispherical vesicles shows that the north pole TMP is typically greater than the south pole, thereby suggesting the typical pathway a charged species might take inside an eccentric nucleate cell under electroporation.

Electrohydrodynamics of Vesicles and Capsules

Giant unilamellar vesicles (GUVs) made up of phospholipid bilayer membranes (liposomes) and elastic capsules with a cross-linked, polymerized membrane, have emerged as biomimetic alternatives to investigating biological cells such as leukocytes and erythrocytes. This feature article looks at the similarities and differences in the electrohydrodynamics (EHD) of vesicles and capsules under electric fields that determines their electromechanical response. The physics of EHD is illustrated through several examples such as the electrodeformation of single and compound, spherical and cylindrical, and charged and uncharged vesicles in uniform and nonuniform electric fields, and the relevance and challenges are discussed. Both small and large deformation results are discussed. The use of EHD in understanding complex interfacial kinetics in capsules and the synthesis of nonspherical capsules using electric fields are also presented. Finally, the review looks at the large electrodeformation of water-in-water capsules and the relevance of constitutive laws in their response.

Vesicle formation mechanisms: an overview

Vesicle structures primarily embody spherical capsules composed of a single or multiple bilayers, entrapping a pool of aqueous solution in their interior. The bilayers can be synthesised by phospholipids or other amphiphiles (surfactants, block copolymers, etc.). Vesicles with broad-spectrum applications in numerous scientific disciplines, including biochemistry, biophysics, biology, and various pharmaceutical industries, have attracted widespread attention. Consequently, a multitude of protocols have been devised and proposed for their fabrication. In this review, with a motivation to derive the basic conditions for the formation of vesicles, the associated thermodynamic and kinetic aspects are comprehensively appraised. Contextually, an all-purpose overview of the underlying thermodynamics of bilayer/membrane generation and deformation, including the chemical potential of aggregates, geometric packing and the concept of elastic properties, is presented. Additionally, the current review highlights the probable, inherent mechanisms of vesicle formation under distinct modes of manufacturing. We lay focus on vesicle formation from pre-existing bilayers, as well as from bilayers, which form when lipids from an organic solvent are transferred into an aqueous medium. Furthermore, we outline the kinetic effects on vesicle formation from the lamellar phase, with and without the presence of shearing force. Wherever required, the experimental and/or theoretical outcomes, the driving forces for vesicle size selection, and various scaling laws are also reviewed, all of which facilitate an overall improved understanding of the vesicle formation mechanisms.

Movement of giant lipid vesicles induced by millimeter wave radiation change when they contain magnetic nanoparticles

Superparamagnetic iron oxide nanoparticles are used in a rapidly expanding number of research and practical applications in biotechnology and biomedicine. Recent developments in iron oxide nanoparticle design and understanding of nanoparticle membrane interactions have led to applications in magnetically triggered, liposome delivery vehicles with controlled structure. Here we study the effect of external physical stimuli-such as millimeter wave radiation-on the induced movement of giant lipid vesicles in suspension containing or not containing iron oxide maghemite (γ-Fe₂O₃) nanoparticles (MNPs). To increase our understanding of this phenomenon, we used a new microscope image-based analysis to reveal millimeter wave (MMW)-induced effects on the movement of the vesicles. We found that in the lipid vesicles not containing MNPs, an exposure to MMW induced collective reorientation of vesicle motion occurring at the onset of MMW switch "on." Instead, no marked changes in the movements of lipid vesicles containing MNPs were observed at the onset of first MMW switch on, but, importantly, by examining the course followed; once the vesicles are already irradiated, a directional motion of vesicles was induced. The latter vesicles were characterized by a planar motion, absence of gravitational effects, and having trajectories spanning a range of deflection angles narrower than vesicles not containing MNPs. An explanation for this observed delayed response could be attributed to the possible interaction of MNPs with components of lipid membrane that, influencing, e.g., phospholipids density and membrane stiffening, ultimately leads to change vesicle movement.

Giant Vesicles and Their Use in Assays for Assessing Membrane Phase State, Curvature, Mechanics, and Electrical Properties

Giant unilamellar vesicles represent a promising and extremely useful model biomembrane system for systematic measurements of mechanical, thermodynamic, electrical, and rheological properties of lipid bilayers as a function of membrane composition, surrounding media, and temperature. The most important advantage of giant vesicles over other model membrane systems is that the membrane responses to external factors such as ions, (macro)molecules, hydrodynamic flows, or electromagnetic fields can be directly observed under the microscope. Here, we briefly review approaches for giant vesicle preparation and describe several assays used for deducing the membrane phase state and measuring a number of material properties, with further emphasis on membrane reshaping and curvature.

Shape deformation of a vesicle under an axisymmetric non-uniform alternating electric field

We suggest that non-uniform electric fields that are commonly used to study vesicle dielectrophoresis can be employed in hitherto relatively unexplored areas of vesicle deformation (for electromechanical characterization) and electroporation. Conventionally, the tension generated in vesicles is commonly modeled to be entropic or enthalpic in origin. A comparison of the configuration of a vesicle in the enthalpic and entropic regimes as well as the cross over between the two regimes during vesicle deformation has eluded understanding. A lucid demonstration of this concept is provided by the study of vesicle deformation under axisymmetric quadrupole electric field and the shapes of the vesicles obtained using the entropic and the enthalpic approaches, show significant differences. A strong dependence of the final vesicle shapes on the ratio of electrical conductivities of the fluids inside and outside the vesicle as well as on the frequency of the applied quadrupole electric field is observed. A comparison with experimental data from the literature is also made. Moreover, an excess area dependent transition between the entropic and enthalpic regimes is observed. The method could be used to estimate electromechanical properties of the vesicle.

Electrodeformation of Vesicles Suspended in a Liquid Medium

Deformation of flexible vesicles suspended in a fluid medium due to an applied electric field can provide valuable insight into deformation dynamics at a very small scale. In an electric field, the response of the vesicle membrane is strongly influenced by the conductivity of surrounding fluid, vesicle size and shape, and the magnitude of applied field. We studied the electrodeformation of vesicles immersed in a fluid media under a DC electric field. An immersed interface method is used to solve the electric field over the domain with conductive or non-conductive vesicles while an immersed boundary method is employed to solve fluid flow, fluid-solid interaction, membrane mechanics and vesicle deformation. Initial force analysis on the membrane surface reveals almost linear influence of vesicle size, but the vesicle size does not affect the long-term deformation which is consistent with experimental evidence. Highly nonlinear effect of the applied field as well as the conductivity ratios inside and outside of the vesicle are observed. Results also point towards an early linear deformation regime followed by an equilibrium stage for the membranes. Modeling results suggest that electrodeforming vesicles can create unique external flows for different conductivity ratios. Moreover, significant influence of the initial aspect ratio of the vesicle on the force distribution is observed across a range of conductivity ratios.

Three-dimensional multicomponent vesicles: dynamics and influence of material properties

In this work, the nonlinear hydrodynamics of a three-dimensional multicomponent vesicle in shear flow are explored. Using a volume- and area-conserving projection method coupled to a gradient-augmented level set and surface phase field method, the dynamics are systematically studied as a function of the membrane bending rigidity difference between the components, the speed of diffusion compared to the underlying shear flow, and the strength of the phase domain energy compared to the bending energy. Using a pre-segregated vesicle, three dynamics are observed: stationary phase, phase-treading, and a new dynamic called vertical banding. These regimes are very sensitive to the strength of the domain line energy, as the vertical banding regime is not observed when the line energy is larger than the bending energy. The findings demonstrate that a complete understanding of multicomponent vesicle dynamics requires that the full three-dimensional system be modeled, and show the complexity obtained when considering heterogeneous material properties.

Electric Field Deformation of Protein-Coated Droplets in Thin Channels

High-strength droplet interfaces are attractive for many applications, specifically in cases where droplets are channeled through fluidic devices and manipulated by electromagnetic fields. Using models and experiments, we study the deformation of droplets and capsules with protein interfaces in an electric field in thin and wide electrode gaps. Proteins are chosen from candidates expected to display qualitatively different interfacial interactions and strengths: a globular protein (bovine serum albumin), a reversible cross-linking peptide (AFD4), and a hydrophobin (cerato ulmin). Dilute protein additives can lead to over 1 order of magnitude stronger oil-water interfaces than those stabilized by small surfactants. We develop small deformation models to evaluate a protein membrane's interfacial elasticity, notably accounting for the electric field perturbation encountered in a gap and a careful treatment of a generalized elastic interface with both surface tension and interfacial elasticity. Results indicate that globular proteins, which typically have comparable surface tension and interfacial elasticity, can be modeled well by this generalized elastic interface. We further find that when in a gap, droplets and capsules migrate toward one electrode, deform asymmetrically, exhibit polar spreading on the electrode, and predictably stretch more than in the infinite gap scenario at constant field strength.

Large deformation electrohydrodynamics of a Skalak elastic capsule in AC electric field

The axisymmetric electrohydrodynamic deformation of an elastic capsule with a capacitive membrane obeying the Skalak law under a uniform AC electric field is investigated using analytical and boundary integral theory. The low capillary number (the ratio of destabilizing shear or electric force to the stabilizing elastic force) regime shows that time-averaged prolate and oblate spheroid deformations, and the time-periodic prolate-sphere, oblate-sphere breathing modes are commensurate with the time averaged-deformation. A novel prolate-oblate breathing mode is observed due to an interplay of finite membrane charging time and the field reversal of the AC field. The study, when extended to high capillary numbers, shows new breathing modes of cylinder-prolate, cylinder-oblate, and biconcave-prolate deformation. These are the results of highly compressive normal Maxwell stress at the poles and are aided by a weak compressive equatorial stress, characteristic of a capacitive membrane. The findings of this work should form the basis for the understanding of more complex biological cells and synthetic capsules for industrial applications.

Lipid vesicles in pulsed electric fields: Fundamental principles of the membrane response and its biomedical applications

The present review focuses on the effects of pulsed electric fields on lipid vesicles ranging from giant unilamellar vesicles (GUVs) to small unilamellar vesicles (SUVs), from both fundamental and applicative perspectives. Lipid vesicles are the most popular model membrane systems for studying biophysical and biological processes in living cells. Furthermore, as vesicles are made from biocompatible and biodegradable materials, they provide a strategy to create safe and functionalized drug delivery systems in health-care applications. Exposure of lipid vesicles to pulsed electric fields is a common physical method to transiently increase the permeability of the lipid membrane. This method, termed electroporation, has shown many advantages for delivering exogenous molecules including drugs and genetic material into vesicles and living cells. In addition, electroporation can be applied to induce fusion between vesicles and/or cells. First, we discuss in detail how research on cell-size GUVs as model cell systems has provided novel insight into the basic mechanisms of cell electroporation and associated phenomena. Afterwards, we continue with a thorough overview how electroporation and electrofusion have been used as versatile methods to manipulate vesicles of all sizes in different biomedical applications. We conclude by summarizing the open questions in the field of electroporation and possible future directions for vesicles in the biomedical field.

Electrohydrodynamics of a compound vesicle under an AC electric field

Compound vesicles are relevant as simplified models for biological cells as well as in technological applications such as drug delivery. Characterization of these compound vesicles, especially the inner vesicle, remains a challenge. Similarly their response to electric field assumes importance in light of biomedical applications such as electroporation. Fields lower than that required for electroporation cause electrodeformation in vesicles and can be used to characterize their mechanical and electrical properties. A theoretical analysis of the electrohydrodynamics of a compound vesicle with outer vesicle of radius R _o and an inner vesicle of radius [Formula: see text], is presented. A phase diagram for the compound vesicle is presented and elucidated using detailed plots of electric fields, free charges and electric stresses. The electrohydrodynamics of the outer vesicle in a compound vesicle shows a prolate-sphere and prolate-oblate-sphere shape transitions when the conductivity of the annular fluid is greater than the outer fluid, and vice-versa respectively, akin to single vesicle electrohydrodynamics reported in the literature. The inner vesicle in contrast shows sphere-prolate-sphere and sphere-prolate-oblate-sphere transitions when the inner fluid conductivity is greater and smaller than the annular fluid, respectively. Equations and methodology are provided to determine the bending modulus and capacitance of the outer as well as the inner membrane, thereby providing an easy way to characterize compound vesicles and possibly biological cells.

The Influence of Vesicle Shape and Medium Conductivity on Possible Electrofusion under a Pulsed Electric Field

The effects of electric field on lipid membrane and cells have been extensively studied in the last decades. The phenomena of electroporation and electrofusion are of particular interest due to their wide use in cell biology and biotechnology. However, numerical studies on the electrofusion of cells (or vesicles) with different deformed shapes are still rare. Vesicle, being of cell size, can be treated as a simple model of cell to investigate the behaviors of cell in electric field. Based on the finite element method, we investigate the effect of vesicle shape on electrofusion of contact vesicles in various medium conditions. The transmembrane voltage (TMV) and pore density induced by a pulsed field are examined to analyze the possibility of vesicle fusion. In two different medium conditions, the prolate shape is observed to have selective electroporation at the contact area of vesicles when the exterior conductivity is smaller than the interior one; selective electroporation is more inclined to be found at the poles of the oblate vesicles when the exterior conductivity is larger than the interior one. Furthermore, we find that when the exterior conductivity is lower than the internal conductivity, the pulse can induce a selective electroporation at the contact area between two vesicles regardless of the vesicle shape. Both of these two findings have important practical applications in guiding electrofusion experiments.

Wrinkling dynamics of fluctuating vesicles in time-dependent viscous flow

We study the fully nonlinear, nonlocal dynamics of two-dimensional vesicles in a time-dependent, incompressible viscous flow at finite temperature. We focus on a transient instability that can be observed when the direction of applied flow is suddenly reversed, which induces compressive forces on the vesicle interface, and small-scale interface perturbations known as wrinkles develop. These wrinkles are driven by regions of negative elastic tension on the membrane. Using a stochastic immersed boundary method with a biophysically motivated choice of thermal fluctuations, we investigate the wrinkling dynamics numerically. Different from deterministic wrinkling dynamics, thermal fluctuations lead to symmetry-breaking wrinkling patterns by exciting higher order modes. This leads to more rapid and more realistic wrinkling dynamics. Our results are in excellent agreement with the experimental data by Kantsler et al. [Kantsler et al., Phys. Rev. Lett., 2007, 99, 17802]. We compare the nonlinear simulation results with perturbation theory, modified to account for thermal fluctuations. The strength of the applied flow strongly influences the most unstable wavelength characterizing the wrinkles, and there are significant differences between the results from perturbation theory and the fully nonlinear simulations, which suggests that the perturbation theory misses important nonlinear interactions. Strikingly, we find that thermal fluctuations actually have the ability to attenuate variability of the characteristic wavelength of wrinkling by exciting a wider range of modes than the deterministic case, which makes the evolution less constrained and enables the most unstable wavelength to emerge more readily. We further find that thermal noise helps prevent the vesicle from rotating if it is misaligned with the direction of the applied extensional flow.

Membrane tubulation from giant lipid vesicles in alternating electric fields

We report on the formation of tubular membrane protrusions from giant unilamellar vesicles in alternating electric fields. The construction of the experimental chamber permitted the application of external AC fields with strength of dozens of V/mm and kHz frequency during relatively long time periods (several minutes). Besides the vesicle electrodeformation from quasispherical to prolate ellipsoidal shape, the formation of long tubular membrane protrusions with length of up to several vesicle diameters, arising from the vesicular surface in the field direction, was registered and analyzed. The threshold electric field at which the electro-induced protrusions appeared was lower than the field strengths inducing membrane electroporation.

Voltage-morphology coupling in biomimetic membranes: dynamics of giant vesicles in applied electric fields

An electric potential difference across the plasma membrane is common to all living cells and is essential to physiological functions such as the generation of action potentials for cell-to-cell communication. While the basics of cell electrical activity are well established (e.g. the Hodgkin-Huxley model of the action potential), the reciprocal coupling of voltage and membrane deformation has received limited attention. In recent years, studies of biomimetic membranes in externally applied electric fields have revealed a plethora of intriguing dynamics (formation of edges, pearling, and phase separation) that challenge the current understanding of membrane electromechanics.

Ellipsoidal Relaxation of Deformed Vesicles

Theoretical analysis and experimental quantification on the ellipsoidal relaxation of vesicles are presented. The current work reveals the simplicity and universal aspects of this process. The Helfrich formula is shown to apply to the dynamic relaxation of moderate-to-high tension membranes, and a closed-form solution is derived which predicts the vesicle aspect ratio as a function of time. Scattered data are unified by a time scale, which leads to a similarity behavior, governed by a distinctive solution for each vesicle type. Two separate regimes in the relaxation are identified, namely, the "entropic" and the "constant-tension" regimes. The bending rigidity and the initial membrane tension can be simultaneously extracted from the data analysis, posing the current approach as an effective means for the mechanical analysis of biomembranes.

Deformation of phospholipid vesicles in an optical stretcher

Phospholipid vesicles are common model systems for cell membranes. Important aspects of the membrane function relate to its mechanical properties. Here we have investigated the deformation behaviour of phospholipid vesicles in a dual-beam laser trap, also called an optical stretcher. This study explicitly makes use of the inherent heating present in such traps to investigate the dependence of vesicle deformation on temperature. By using lasers with different wavelengths, optically induced mechanical stresses and temperature increase can be tuned fairly independently with a single setup. The phase transition temperature of vesicles can be clearly identified by an increase in deformation. In the case of no heating effects, a minimal model for drop deformation in an optical stretcher and a more specific model for vesicle deformation that takes explicitly into account the angular dependence of the optical stress are presented to account for the experimental results. Elastic constants are extracted from the fitting procedures, which agree with literature data. This study demonstrates the utility of optical stretching, which is easily combined with microfluidic delivery, for the future serial, high-throughput study of the mechanical and thermodynamic properties of phospholipid vesicles.

Dynamics of three-dimensional vesicles in dc electric fields

A numerical and systematic parameter study of three-dimensional vesicle electrohydrodynamics is presented to investigate the effects of varying electric field strength and different fluid and membrane properties. The dynamics of vesicles in the presence of dc electric fields is considered, in both the presence and absence of linear shear flow. For suspended vesicles it is shown that the conductivity ratio and viscosity ratio between the interior and exterior fluids, as well as the vesicle membrane capacitance, substantially affect the minimum electric field strength required to induce a full prolate-oblate-prolate transition. In addition, there exists a critical electric field strength above which a vesicle will no longer tumble when exposed to linear shear flow.

Vesicle dynamics in uniform electric fields: squaring and breathing

We computationally investigate the dynamics of a vesicle exposed to uniform DC or AC electric fields. We employ the two-dimensional boundary integral method in order to simulate vesicle deformation under experimental conditions where peculiar drum-like ("squared") shapes have been observed. The vesicle membrane is modeled as an infinitely thin, capacitive, area-incompressible interface, with the surrounding fluids acting as leaky dielectrics. Our simulations capture the "squaring" phenomenon, in which vesicles deform into rectangular profiles with corner-like regions of high curvature, as vesicles undergo dynamic transitions between oblate and prolate ellipsoidal shapes.

An engineered membrane to measure electroporation: effect of tethers and bioelectronic interface

This article reports on the construction and predictive models for a platform comprised of an engineered tethered membrane. The platform provides a controllable and physiologically relevant environment for the study of the electroporation process. The mixed self-assembled membrane is formed via a rapid solvent exchange technique. The membrane is tethered to the gold electrode and includes an ionic reservoir separating the membrane and gold surface. Above the membrane, there is an electrolyte solution, and a gold counterelectrode. A voltage is applied between the gold electrodes and the current measured. The current is dependent on the energy required to form aqueous pores and the conductance of each pore. A two-level predictive model, consisting of a macroscopic and a continuum model, is developed to relate the pore dynamics to the measured current. The macroscopic model consists of an equivalent circuit model of the tethered membrane, and asymptotic approximations to the Smoluchowski-Einstein equation of electroporation that is dependent on the pore conductance and the energy required to form aqueous pores. The continuum model is a generalized Poisson-Nernst-Planck (GPNP) system where an activity coefficient to account for steric effects of ions is added to the standard PNP system. The GPNP is used to evaluate the conductance of aqueous pores, and the electrical energy required to form the pores. As an outcome of the setup of the device and the two-level model, biologically important variables can be estimated from experimental measurements. To validate the accuracy of the two-level model, the predicted current is compared with experimentally measured current for different tethering densities.

Vesicle biomechanics in a time-varying magnetic field

Background

Cells exhibit distortion when exposed to a strong electric field, suggesting that the field imposes control over cellular biomechanics. Closed pure lipid bilayer membranes (vesicles) have been widely used for the experimental and theoretical studies of cellular biomechanics under this electrodeformation. An alternative method used to generate an electric field is by electromagnetic induction with a time-varying magnetic field. References reporting the magnetic control of cellular mechanics have recently emerged. However, theoretical analysis of the cellular mechanics under a time-varying magnetic field is inadequate.

We developed an analytical theory to investigate the biomechanics of a modeled vesicle under a time-varying magnetic field. Following previous publications and to simplify the calculation, this model treated the inner and suspending media as lossy dielectrics, the membrane thickness set at zero, and the electric resistance of the membrane assumed to be negligible. This work provided the first analytical solutions for the surface charges, electric field, radial pressure, overall translational forces, and rotational torques introduced on a vesicle by the time-varying magnetic field. Frequency responses of these measures were analyzed, particularly the frequency used clinically by transcranial magnetic stimulation (TMS).

Results

The induced surface charges interacted with the electric field to produce a biomechanical impact upon the vesicle. The distribution of the induced surface charges depended on the orientation of the coil and field frequency. The densities of these charges were trivial at low frequency ranges, but significant at high frequency ranges. The direction of the radial force on the vesicle was dependent on the conductivity ratio between the vesicle and the medium. At relatively low frequencies (<200 KHz), including the frequency used in TMS, the computed radial pressure and translational forces on the vesicle were both negligible.

Conclusions

This work provides an analytical framework and insight into factors affecting cellular biomechanics under a time-varying magnetic field. Biological effects of clinical TMS are not likely to occur via alteration of the biomechanics of brain cells.

Recent developments in the field of bending rigidity measurements on membranes

This review gives a brief overview of experimental approaches used to assess the bending rigidity of membranes. Emphasis is placed on techniques based on the use of giant unilamellar vesicles. We summarize the effect on the bending rigidity of membranes as a function of membrane composition, presence of various inclusions in the bilayer and molecules and ions in the bathing solutions. Examples for the impact of temperature, cholesterol, some peptides and proteins, sugars and salts are provided and the literature data are discussed critically. Future directions, open questions and possible developments in this research field are also included.

Vesicle deformation in DC electric pulses

The transient deformation of giant vesicles in square DC electric pulses is investigated. We experimentally observe the theoretically predicted transition from an oblate to prolate ellipsoidal shape in the case of a quasi-spherical vesicle encapsulating solution less conducting than the suspending medium. The transition is detected by utilizing a two-step pulse in order to avoid electroporation and vesicle collapse. We develop a theoretical model to describe both the deformation under an electric field and relaxation after the field is turned off. Agreement between experiment and theory demonstrates that the time-dependent vesicle shape can be used to measure membrane properties such as viscosity and capacitance.