Giant vesicles in electric fields

Considering Joule heating in coupled electroporation and electrodeformation modeling of glioblastoma cells

Cells exposed to a pulsed electric field undergo electroporation(EP) and electrodeformation(ED) under electric field stress, and a coupled model of EP and ED of glioblastoma(GBM) taking into account Joule heating is proposed. The model geometry is extracted from real cell boundaries, and the effects of Joule heating-induced temperature rise on the EP and ED processes are considered. The results show that the temperature rise will increase the cell's local conductivity, leading to a decrease in the transmembrane potential(TMP). The temperature rise also causes a decrease in the dynamic Young's modulus of the cell membrane, making the cell less resistant to deformation. In addition, GBM cells are more susceptible to EP in the middle portion of the cell and ED in the three tentacle portions under pulsed electric fields, and the cells undergo significant positional shifts. The ED of the nucleus is similar to spherical cells, but the degree of ED is smaller.

Probability and kinetics of rupture and electrofusion in giant unilamellar vesicles under various frequencies of direct current pulses

Irreversible electroporation induces permanent permeabilization of lipid membranes of vesicles, resulting in vesicle rupture upon the application of a pulsed electric field. Electrofusion is a phenomenon wherein neighboring vesicles can be induced to fuse by exposing them to a pulsed electric field. We focus how the frequency of direct current (DC) pulses of electric field impacts rupture and electrofusion in cell-sized giant unilamellar vesicles (GUVs) prepared in a physiological buffer. The average time, probability, and kinetics of rupture and electrofusion in GUVs have been explored at frequency 500, 800, 1050, and 1250 Hz. The average time of rupture of many 'single GUVs' decreases with the increase in frequency, whereas electrofusion shows the opposite trend. At 500 Hz, the rupture probability stands at 0.45 ± 0.02, while the electrofusion probability is 0.71 ± 0.01. However, at 1250 Hz, the rupture probability increases to 0.69 ± 0.03, whereas the electrofusion probability decreases to 0.46 ± 0.03. Furthermore, when considering kinetics, at 500 Hz, the rate constant of rupture is (0.8 ± 0.1)×10-2 s-1, and the rate constant of fusion is (2.4 ± 0.1)×10-2 s-1. In contrast, at 1250 Hz, the rate constant of rupture is (2.3 ± 0.8)×10-2 s-1, and the rate constant of electrofusion is (1.0 ± 0.1)×10-2 s-1. These results are discussed by considering the electrical model of the lipid bilayer and the energy barrier of a prepore.

Photoswitchable Endocytosis of Biomolecular Condensates in Giant Vesicles

Interactions between membranes and biomolecular condensates can give rise to complex phenomena such as wetting transitions, mutual remodeling, and endocytosis. In this study, light-triggered manipulation of condensate engulfment is demonstrated using giant vesicles containing photoswitchable lipids. UV irradiation increases the membrane area, which can be stored in nanotubes. When in contact with a condensate droplet, the UV light triggers rapid condensate endocytosis, which can be reverted by blue light. The affinity of the protein-rich condensates to the membrane and the reversibility of the engulfment processes is quantified from confocal microscopy images. The degree of photo-induced engulfment, whether partial or complete, depends on the vesicle excess area and the relative sizes of vesicles and condensates. Theoretical estimates suggest that utilizing the light-induced excess area to increase the vesicle-condensate adhesion interface is energetically more favorable than the energy gain from folding the membrane into invaginations and tubes. The overall findings demonstrate that membrane-condensate interactions can be easily and quickly modulated via light, providing a versatile system for building platforms to control cellular events and design intelligent drug delivery systems for cell repair.

Dynamics of giant vesicle assembly from thin lipid films

MOTIVATION

Giant unilamellar vesicles (GUVs), cell-like synthetic micrometer size structures, assemble when thin lipid films are hydrated in aqueous solutions. Quantitative measurements of static yields and distribution of sizes of GUVs obtained from thin film hydration methods were recently reported. Dynamic data such as the time evolution of yields and distribution of sizes, however, is not known. Dynamic data can provide insights into the assembly pathway of GUVs and guidelines for choosing conditions to obtain populations with desired size distributions.

APPROACH

We develop the 'stopped-time' technique to characterize the time evolution of the distribution of sizes and molar yields of populations of free-floating GUVs. We additionally capture high resolution time-lapse images of surface-attached GUV buds on the lipid films. We systematically study the dynamics of assembly of GUVs from three widely used thin film hydration methods, PAPYRUS (Paper-Abetted amPhiphile hYdRation in aqUeous Solutions), gentle hydration, and electroformation.

FINDINGS

We find that the molar yield versus time curves of GUVs demonstrate a characteristic sigmoidal shape, with an initial yield, a transient, and then a steady state plateau for all three methods. The population of GUVs showed a right-skewed distribution of diameters. The variance of the distributions increased with time. The systems reached steady state within 120 min. We rationalize the dynamics using the thermodynamically motivated budding and merging (BNM) model. These results further the understanding of lipid dynamics and provide for the first-time practical parameters to tailor the production of GUVs of specific sizes for applications.

(2+δ)-dimensional theory of the electromechanics of lipid membranes: Electrostatics

The coupling of electric fields to the mechanics of lipid membranes gives rise to intriguing electromechanical behavior, as, for example, evidenced by the deformation of lipid vesicles in external electric fields. Electromechanical effects are relevant for many biological processes, such as the propagation of action potentials in axons and the activation of mechanically gated ion channels. Currently, a theoretical framework describing the electromechanical behavior of arbitrarily curved and deforming lipid membranes does not exist. Purely mechanical models commonly treat lipid membranes as two-dimensional surfaces, ignoring their finite thickness. While holding analytical and numerical merit, this approach cannot describe the coupling of lipid membranes to electric fields and is thus unsuitable for electromechanical models. In a sequence of articles, we derive an effective surface theory of the electromechanics of lipid membranes, called the (2+δ)-dimensional theory, which has the advantages of surface descriptions while accounting for finite thickness effects. The present article proposes a generic dimension reduction procedure relying on low-order spectral expansions. This procedure is applied to the electrostatics of lipid membranes to obtain the (2+δ)-dimensional theory that captures potential differences across and electric fields within lipid membranes. This model is tested on different geometries relevant for lipid membranes, showing good agreement with the corresponding three-dimensional electrostatics theory.

Light-Activated Synthetic Rotary Motors in Lipid Membranes Induce Shape Changes Through Membrane Expansion

Membranes are the key structures to separate and spatially organize cellular systems. Their rich dynamics and transformations during the cell cycle are orchestrated by specific membrane-targeted molecular machineries, many of which operate through energy dissipation. Likewise, man-made light-activated molecular rotary motors have previously shown drastic effects on cellular systems, but their physical roles on and within lipid membranes remain largely unexplored. Here, the impact of rotary motors on well-defined biological membranes is systematically investigated. Notably, dramatic mechanical transformations are observed in these systems upon motor irradiation, indicative of motor-induced membrane expansion. The influence of several factors on this phenomenon is systematically explored, such as motor concentration and membrane composition., Membrane fluidity is found to play a crucial role in motor-induced deformations, while only minor contributions from local heating and singlet oxygen generation are observed. Most remarkably, the membrane area expansion under the influence of the motors continues as long as irradiation is maintained, and the system stays out-of-equilibrium. Overall, this research contributes to a comprehensive understanding of molecular motors interacting with biological membranes, elucidating the multifaceted factors that govern membrane responses and shape transitions in the presence of these remarkable molecular machines, thereby supporting their future applications in chemical biology.

Determination of pore edge tension from the kinetics of rupture of giant unilamellar vesicles using the Arrhenius equation: effects of sugar concentration, surface charge and cholesterol

The pore edge tension (Γ) of a membrane closely intertwines with membrane stability and plays a vital role in the mechanisms that facilitate membrane resealing following pore formation caused by electrical and mechanical tensions. We have explored a straightforward procedure to determine Γ by fitting the inverse of the tension-dependent logarithm of the rate constant of rupture of giant unilamellar vesicles (GUVs) using the Arrhenius equation. The GUVs were prepared using a combination of 1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DOPG) and 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in a physiological environment. The effects of sugar concentration, membrane surface charge density, and membrane cholesterol concentration on Γ have been investigated. The values of Γ increase with sugar concentration in the physiological buffer, measuring 9.6 ± 0.3, 10.4 ± 0.1, and 16.2 ± 0.1 pN for 40, 100, and 300 mM, respectively. A higher concentration of anionic lipids (70 mol% of DOPG) significantly reduces Γ. An increasing trend of Γ with cholesterol content was observed; specifically, the values of Γ were 11.9 ± 0.9, 13.9 ± 0.7, and 16.2 ± 0.4 pN for 15, 29, and 40 mol% cholesterol, respectively. Thus, the presence of higher anionic lipids in the bilayer led to a decrease in membrane stability. In contrast, the presence of higher sugar concentrations in the buffer and increased cholesterol concentration in the membranes enhanced membrane stability.

Dynamics of viscous droplet coalescence in the confined geometry of optical cells

The dynamics of quasi-two-dimensional coalescence of isotropic droplets in nematic liquid crystal environment was studied. Investigations were made in confined geometry of a Hele-Shaw optical cell with different transverse droplet sizes. The existence of three distinct dynamic regimes was found for coalescence, namely, short-, middle-, and long-time regimes. The fast dynamics of bridge transformation was visualized. At short time the dynamics of droplet transformation is similar to the transformation of free (three-dimensional) droplets. At later stages, two regimes of the coalescence at different timescales are determined by Poiseuille flow. Experimental data are discussed on the basis of existing theories.

Electrochemistry in sensing of molecular interactions of proteins and their behavior in an electric field

Electrochemical methods can be used not only for the sensitive analysis of proteins but also for deeper research into their structure, transport functions (transfer of electrons and protons), and sensing their interactions with soft and solid surfaces. Last but not least, electrochemical tools are useful for investigating the effect of an electric field on protein structure, the direct application of electrochemical methods for controlling protein function, or the micromanipulation of supramolecular protein structures. There are many experimental arrangements (modalities), from the classic configuration that works with an electrochemical cell to miniaturized electrochemical sensors and microchip platforms. The support of computational chemistry methods which appropriately complement the interpretation framework of experimental results is also important. This text describes recent directions in electrochemical methods for the determination of proteins and briefly summarizes available methodologies for the selective labeling of proteins using redox-active probes. Attention is also paid to the theoretical aspects of electron transport and the effect of an external electric field on the structure of selected proteins. Instead of providing a comprehensive overview, we aim to highlight areas of interest that have not been summarized recently, but, at the same time, represent current trends in the field.

Identification of Proteins Involved in Cell Membrane Permeabilization by Nanosecond Electric Pulses (nsEP)

The study was aimed at identifying endogenous proteins which assist or impede the permeabilized state in the cell membrane disrupted by nsEP (20 or 40 pulses, 300 ns width, 7 kV/cm). We employed a LentiArray CRISPR library to generate knockouts (KOs) of 316 genes encoding for membrane proteins in U937 human monocytes stably expressing Cas9 nuclease. The extent of membrane permeabilization by nsEP was measured by the uptake of Yo-Pro-1 (YP) dye and compared to sham-exposed KOs and control cells transduced with a non-targeting (scrambled) gRNA. Only two KOs, for SCNN1A and CLCA1 genes, showed a statistically significant reduction in YP uptake. The respective proteins could be part of electropermeabilization lesions or increase their lifespan. In contrast, as many as 39 genes were identified as likely hits for the increased YP uptake, meaning that the respective proteins contributed to membrane stability or repair after nsEP. The expression level of eight genes in different types of human cells showed strong correlation (R > 0.9, p < 0.02) with their LD₅₀ for lethal nsEP treatments, and could potentially be used as a criterion for the selectivity and efficiency of hyperplasia ablations with nsEP.

Temperature-Promoted Giant Unilamellar Vesicle (GUV) Aggregation: A Way of Multicellular Formation

The evolution of unicellular to multicellular life is considered to be an important step in the origin of life, and it is crucial to study the influence of environmental factors on this process through cell models in the laboratory. In this paper, we used giant unilamellar vesicles (GUVs) as a cell model to investigate the relationship between environmental temperature changes and the evolution of unicellular to multicellular life. The zeta potential of GUVs and the conformation of the headgroup of phospholipid molecules at different temperatures were examined using phase analysis light scattering (PALS) and attenuated total reflection-Fourier transform infrared spectroscopy (ATR-FTIR), respectively. In addition, the effect of increasing temperature on the aggregation of GUVs was further investigated in ionic solutions, and the possible mechanisms involved were explored. The results showed that increasing temperature reduced the repulsive forces between cells models and promoted their aggregation. This study could effectively contribute to our understanding of the evolution of primitive unicellular to multicellular life.

Microfluidic techniques for mechanical measurements of biological samples

The use of microfluidics to make mechanical property measurements is increasingly common. Fabrication of microfluidic devices has enabled various types of flow control and sensor integration at micrometer length scales to interrogate biological materials. For rheological measurements of biofluids, the small length scales are well suited to reach high rates, and measurements can be made on droplet-sized samples. The control of flow fields, constrictions, and external fields can be used in microfluidics to make mechanical measurements of individual bioparticle properties, often at high sampling rates for high-throughput measurements. Microfluidics also enables the measurement of bio-surfaces, such as the elasticity and permeability properties of layers of cells cultured in microfluidic devices. Recent progress on these topics is reviewed, and future directions are discussed.

Domain dynamics of phase-separated lipid membranes under shear flow

The dynamical behaviour of lateral domains on phase-separated lipid vesicles under external flow is reported. A microfluidic chamber was used for the immobilization of vesicles and the application of shear. Microscopic observation revealed that domains tended to be localized at the vortex center and to exhibit a stripe morphology as the flow speed increased. We clarified the dependency of domain behaviors on the flow speed and lipid mixing fraction. The cholesterol ratio in the membrane affected these domain behaviors. Next, we investigated the growth of domains under flow. We discuss the mechanism of these trends by considering the free energy of phase separation, and reproduce the experimental results by numerical simulations. These findings may lead to a better understanding of the dynamical properties of the membrane under nonequilibrium situations and the biophysical mechanism of cellular mechanotransduction.

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.

Effects of cholesterol on the anionic magnetite nanoparticle-induced deformation and poration of giant lipid vesicles

We have investigated the effects of cholesterol on the deformation and poration of giant unilamellar vesicles (GUVs) induced by anionic magnetite nanoparticles (NPs). Negatively charged lipid, neutral lipid, and cholesterol were used to prepare the charged GUVs (surface charge density of membranes - 0.16 C m^-2), while only neutral lipid and cholesterol were used to prepare the neutral GUVs. Cholesterol content varied from 0 to 40 mole% for preparing the biologically relevant membranes. The degree of deformation has been characterized by compactness, the value of which remains at 1.0 for spherical GUVs. The value of compactness increases with time for both membranes, but this increase depends on cholesterol content. The average compactness decreases with cholesterol content, and at 60 min, the values are 1.280 ± 0.002 and 1.131 ± 0.010 for 0 and 40 mole% cholesterol containing charged GUVs. The average compactness is relatively lower for neutral GUVs for the corresponding cholesterol. Membrane poration has been investigated by the leakage of calcein, which indicates a two-state transition model. The fraction of deformation is higher for charged GUVs than for neutral ones, while the fraction of poration shows the opposite result. Both the fractions decrease with cholesterol content.

Influence of sugar concentration on the vesicle compactness, deformation and membrane poration induced by anionic nanoparticles

Sugar plays a vital role in the structural and functional characteristics of cells. Hence, the interaction of NPs with cell membranes in the presence of sugar concentrations is important for medicinal and pharmacological innovations. This study integrated three tools: giant unilamellar vesicles (GUVs), anionic magnetite nanoparticles (NPs), and sugar concentrations, to understand a simplified mechanism for interactions between the vesicle membranes and NPs under various sugar concentrations. We focused on changing the sugar concentration in aqueous solution; more precisely, sucrose inside the GUVs and glucose outside with equal osmolarity. 1,2-dioleoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (sodium salt) (DOPG) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) were used to prepare the charged membranes of 40mole%DOPG/60mole%DOPC-GUVs, whereas only DOPC was used to prepare the neutral membranes. Phase contrast fluorescence microscopy shows that the adherence of 18 nm magnetite NPs with anionic charge depends on the sugar concentration. The alterations of GUVs induced by the NPs are characterized in terms of i) vesicle compactness, ii) deformation, and iii) membrane poration. The presence of sugar provides additional structural stability to the GUVs and reduces the effects of the NPs with respect to these parameters; more precisely, the higher the sugar concentration, the smaller the alteration induced by the NPs. The differences in NPs effects are explained by the change in the type of interaction between sugar molecules and lipid membranes, namely enthalpy and entropy-driven interaction, respectively. In addition, such alterations are influenced by the surface charge density of the lipid bilayer. The surface pressure of membranes due to the adsorption of NPs is responsible for inducing the poration in membranes. The differences in deformation and poration in charged and neutral GUVs under various sugar concentrations are discussed based on the structure of the head of lipid molecules.

GM1 asymmetry in the membrane stabilizes pores

Cell membranes are highly asymmetric and their stability against poration is crucial for survival. We investigated the influence of membrane asymmetry on electroporation of giant unilamellar vesicles with membranes doped with GM1, a ganglioside asymmetrically enriched in the outer leaflet of neuronal cell membranes. Compared with symmetric membranes, the lifetimes of micronsized pores are about an order of magnitude longer suggesting that pores are stabilized by GM1. Internal membrane nanotubes caused by the GM1 asymmetry, obstruct and additionally slow down pore closure, effectively reducing pore edge tension and leading to leaky membranes. Our results point to the drastic effects this ganglioside can have on pore resealing in biotechnology applications based on poration as well as on membrane repair processes.

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.

The effect of cholesterol on highly curved membranes measured by nanosecond Fluorescence Correlation Spectroscopy

Small lipid vesicles (with diameter ≤ 100nm) with their highly curved membranes comprise a special class of biological lipid bilayers. The mechanical properties of such membranes are critical for their function, e.g. exocytosis. Cholesterol is a near-universal regulator of membrane properties in animal cells. Yet measurements of the effect of cholesterol on the mechanical properties of membranes have remained challenging, and the interpretation of such measurements has remained a matter of debate. Here we show that nanosecond fluorescence correlation spectroscopy (FCS) can directly measure the ns-microsecond rotational correlation time (τr) of a lipid probe in high curvature vesicles with extraordinary sensitivity. Using a home-built 4-Pi fluorescence cross-correlation spectrometer containing polarization-modulating elements, we measure the rotational correlation time (τr) of Nile Red in neurotransmitter vesicle mimics. As the cholesterol mole fraction increases from 0 to 50 %, τr increases from 17 ± 1 to 112 ± 12 ns, indicating a viscosity change of nearly a factor of 7. These measurements are corroborated by solid-state NMR results, which show that the order parameter of the lipid acyl chains increases by about 50% for the same change in cholesterol concentration. Additionally, we measured the spectral parameters of polarity-sensitive fluorescence dyes, which provide an indirect measure of viscosity. The green/red ratio of Nile Red and the generalized polarization of Laurdan show consistent increases of 1.3x and 2.6x, respectively. Our results demonstrate that rotational FCS can directly measure the viscosity of highly curved membranes with higher sensitivity and wider dynamic range compared to other conventional techniques. Significantly, we observe that the viscosity of neurotransmitter vesicle mimics is remarkably sensitive to their cholesterol content.

Effects of sugar concentration on the electroporation, size distribution and average size of charged giant unilamellar vesicles

We investigated the effects of sugar concentration on the electroporation, size distribution and average size of giant unilamellar vesicles (GUVs). GUVs were prepared from 40 mol% of 1,2-dioleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DOPG) and 60 mol% of 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipids. Pulsed electric field was applied to the 40%DOPG/60%DOPC-GUVs and it induced lateral electric tension (σ_c) in the membranes of vesicles. The σ_c-induced probability of rupture (P_pore) and the rate constant of rupture (k_p) of GUVs under the sugar concentration, c = 40, 100 and 300 mM, were determined. Both the P_pore and k_p increased with the increase of σ_c, but higher tension was required to generate the same values of P_pore and k_p with increasing c. We also investigated average sizes of GUVs from the size distribution of vesicles under various sugar concentrations. With the increase of c, the peak of the size distribution histograms shifted to the region of smaller vesicles. The average size decreased 1.6-fold when c increased from 10 to 300 mM. These investigations help to understand various biomedical, biophysical, and biochemical processes in vesicles and cells. Electroporation, size distribution and average size of charged GUVs were investigated under various sugar concentrations. The sugar concentration influences the electroporation of vesicles and the average size of GUVs.

Shear stress induced lipid order and permeability changes of giant unilamellar vesicles

BACKGROUND

The permeability of a lipid bilayer is a function of its phase state and depends non-linearly on thermodynamic variables such as temperature, pressure or pH. We investigated how shear forces influence the phase state of giant unilamellar vesicles and their membrane permeability.

METHODS

We determined the permeability of giant unilamellar vesicles composed of different phospholipid species under shear flow in a tube at various temperatures around and far off the melting point by analyzing the release of fluorescently labelled dextran. Furthermore, we quantified phase state changes of these vesicles under shear forces using spectral decomposition of the membrane embedded fluorescent dye Laurdan.

RESULTS

We observed that the membrane permeability follows a step function with increasing permeability at the transition from the gel to the fluid phase and vice versa. Second, there was an all-or-nothing permeabilization near the main phase transition temperature and a gradual dye release far off the melting transition. Third, the Laurdan phase state analysis suggests that shear forces induce a reversible melting temperature shift in giant unilamellar vesicle membranes.

MAJOR CONCLUSIONS

The observed effects can be explained best in a scenario in which shear forces directly induce membrane pores that possess relatively long pore lifetimes in proximity to the phase transition.

GENERAL SIGNIFICANCE

Our study elucidates the release mechanism of thermo-responsive drug carriers as we found that liposome permeabilization is not continuous but quantized. Furthermore, the shear force induced melting temperature shift must be taken into consideration when thermo-responsive liposomes are designed.

Actin protein inside DMPC GUVs and its mechanical response to AC electric fields

Cells are dynamic systems with complex mechanical properties, regulated by the presence of different species of proteins capable to assemble (and disassemble) into filamentous forms as required by different cells functions. Giant unilamellar vesicles (GUVs) of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) are systems frequently used as a simplified model of cells because they offer the possibility of assaying separately different stimuli, which is no possible in living cells. Here we present a study of the effect of acting protein on mechanical properties of GUVs, when the protein is inside the vesicles in either monomeric G-actin or filamentous F-actin. For this, rabbit skeletal muscle G-actin is introduced inside GUVs by the electroformation method. Protein polymerization inside the GUVs is promoted by adding to the solution MgCl₂ and the ion carrier A23187 to allow the transport of Mg⁺² ions into the GUVs. To determine how the presence of actin changes the mechanical properties of GUVs, the vesicles are deformed by the application of an AC electric field in both cases with G-actin and with polymerized F-actin. The changes in shape of the vesicles are characterized by optical microscopy and from them the bending stiffness of the membrane are determined. It is found that G-actin has no appreciable effect on the bending stiffness of DMPC GUVs, but the polymerized actin makes the vesicles more rigid and therefore more resistant to deformations. This result is supported by evidence that actin filaments tend to accumulate near the membrane.

PoET: automated approach for measuring pore edge tension in giant unilamellar vesicles

Motivation

A reliable characterization of the membrane pore edge tension of single giant unilamellar vesicles (GUVs) requires the measurement of micrometer sized pores in hundreds to thousands of images. When manually performed, this procedure has shown to be extremely time-consuming and to generate inconsistent results among different users and imaging systems. A user-friendly software for such analysis allowing quick processing and generation of reproducible data had not yet been reported.

Results

We have developed a software (PoET) for automatic pore edge tension measurements on GUVs. The required image processing steps and the characterization of the pore dynamics are performed automatically within the software and its use allowed for a 30-fold reduction in the analysis time. We demonstrate the applicability of the software by comparing the pore edge tension of GUVs of different membrane compositions and surface charges. The approach was applied to electroporated GUVs but is applicable to other means of pore formation.

Availability and implementation

The complete software is implemented in Python and available for Windows at https://dx.doi.org/10.17617/3.7h.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Pressure-Biased Nanopores for Excluded Volume Metrology, Lipid Biomechanics, and Cell-Adhesion Rupturing

Nanopore sensing has been widely used in applications ranging from DNA sequencing to disease diagnosis. To improve these capabilities, pressure-biased nanopores have been explored in the past to-primarily-increase the residence time of the analyte inside the pore. Here, we studied the effect of pressure on the ability to accurately quantify the excluded volume which depends on the current drop magnitude produced by a single entity. Using the calibration standard, the inverse current drop (1/ΔI) decreases linearly with increasing pressure, while the dwell drop reduces exponentially. We therefore had to derive a pressure-corrected excluded volume equation to accurately assess the volume of translocating species under applied pressure. Moreover, a method to probe deformation in nanoliposomes and a single cell is developed as a result. We show that the soft nanoliposomes and even cells deform significantly under applied pressure which can be probed in terms of the shape factor which was introduced in the excluded volume equation. The proposed work has practical applications in mechanobiology, namely, assessing the stiffness and mechanical rigidity of liposomal drug carriers. Pressure-biased pores also enabled multiple observations of cell-cell aggregates as well as their subsequent rupture, potentially allowing for the study of microbial symbioses or pathogen recognition by the human immune system.

Recent developments in the kinetics of ruptures of giant vesicles under constant tension

External tension in membranes plays a vital role in numerous physiological and physicochemical phenomena. In this review, recent developments in the constant electric- and mechanical-tension-induced rupture of giant unilamellar vesicles (GUVs) are considered. We summarize the results relating to the kinetics of GUV rupture as a function of membrane surface charge, ions in the bathing solution, lipid composition, cholesterol content in the membrane, and osmotic pressure. The mechanical stability and line tension of the membrane under these conditions are discussed. The membrane tension due to osmotic pressure and the critical tension of rupture for various membrane compositions are also discussed. The results and their analysis provide a biophysical description of the kinetics of rupture, along with insight into biological processes. Future directions and possible developments in this research area are included.

On the Coupling between Mechanical Properties and Electrostatics in Biological Membranes

Cell membrane structure is proposed as a lipid matrix with embedded proteins, and thus, their emerging mechanical and electrostatic properties are commanded by lipid behavior and their interconnection with the included and absorbed proteins, cytoskeleton, extracellular matrix and ionic media. Structures formed by lipids are soft, dynamic and viscoelastic, and their properties depend on the lipid composition and on the general conditions, such as temperature, pH, ionic strength and electrostatic potentials. The dielectric constant of the apolar region of the lipid bilayer contrasts with that of the polar region, which also differs from the aqueous milieu, and these changes happen in the nanometer scale. Besides, an important percentage of the lipids are anionic, and the rest are dipoles or higher multipoles, and the polar regions are highly hydrated, with these water molecules forming an active part of the membrane. Therefore, electric fields (both, internal and external) affects membrane thickness, density, tension and curvature, and conversely, mechanical deformations modify membrane electrostatics. As a consequence, interfacial electrostatics appears as a highly important parameter, affecting the membrane properties in general and mechanical features in particular. In this review we focus on the electromechanical behavior of lipid and cell membranes, the physicochemical origin and the biological implications, with emphasis in signal propagation in nerve cells.

To Close or to Collapse: The Role of Charges on Membrane Stability upon Pore Formation

Resealing of membrane pores is crucial for cell survival. Membrane surface charge and medium composition are studied as defining regulators of membrane stability. Pores are generated by electric field or detergents. Giant vesicles composed of zwitterionic and negatively charged lipids mixed at varying ratios are subjected to a strong electric pulse. Interestingly, charged vesicles appear prone to catastrophic collapse transforming them into tubular structures. The spectrum of destabilization responses includes the generation of long-living submicroscopic pores and partial vesicle bursting. The origin of these phenomena is related to the membrane edge tension, which governs pore closure. This edge tension significantly decreases as a function of the fraction of charged lipids. Destabilization of charged vesicles upon pore formation is universal-it is also observed with other poration stimuli. Disruption propensity is enhanced for membranes made of lipids with higher degree of unsaturation. It can be reversed by screening membrane charge in the presence of calcium ions. The observed findings in light of theories of stability and curvature generation are interpreted and mechanisms acting in cells to prevent total membrane collapse upon poration are discussed. Enhanced membrane stability is crucial for the success of electroporation-based technologies for cancer treatment and gene transfer.

Effects of osmotic pressure on the irreversible electroporation in giant lipid vesicles

Irreversible electroporation (IRE) is a nonthermal tumor/cell ablation technique in which a series of high-voltage short pulses are used. As a new approach, we aimed to investigate the rupture of giant unilamellar vesicles (GUVs) using the IRE technique under different osmotic pressures (Π), and estimated the membrane tension due to Π. Two categories of GUVs were used in this study. One was prepared with a mixture of dioleoylphosphatidylglycerol (DOPG), dioleoylphosphatidylcholine (DOPC) and cholesterol (chol) for obtaining more biological relevance while other with a mixture of DOPG and DOPC, with specific molar ratios. We determined the rate constant (kp) of rupture of DOPG/DOPC/chol (46/39/15)-GUVs and DOPG/DOPC (40/60)-GUVs induced by constant electric tension (σc) under different Π. The σc dependent kp values were fitted with a theoretical equation, and the corresponding membrane tension (σoseq) at swelling equilibrium under Π was estimated. The estimated membrane tension agreed well with the theoretical calculation within the experimental error. Interestingly, the values of σoseq were almost same for both types of synthesized GUVs under same osmotic pressure. We also examined the sucrose leakage, due to large osmotic pressure-induced pore formation, from the inside of DOPG/DOPC/chol(46/39/15)-GUVs. The estimated membrane tension due to large Π at which sucrose leaked out was very similar to the electric tension at which GUVs were ruptured without Π. We explained the σc and Π induced pore formation in the lipid membranes of GUVs.

Membrane Tension Modifies Redox Loading and Release in Single Liposome Electroanalysis

Here, we present a study of how liposomes are loaded and release their contents during their electrochemical detection. We loaded 200 nm liposomes with a redox mediator, ferrocyanide, and used amperometry to detect their collision on a carbon-fiber microelectrode (CFE). We found that we could control the favorability of their electroporation process and the amount of ferrocyanide released by modifying the osmolarity of the buffer in which the liposomes were suspended. Interestingly, we observed that the quantity of the released ferrocyanide varied significantly with buffer osmolarity in a nonmonotonic fashion. Using stimulated Raman scattering (SRS), we confirmed that this behavior was partly explained by fluctuations in the intravesicular redox concentration in response to osmotic pressure. To our surprise, the redox concentration obtained from SRS was much greater than that obtained from amperometry, implying that liposomes may release only a fraction of their contents during electroporation. Consistent with this hypothesis, we observed barrages of electrochemical signals that far exceeded the frequency predicted by Poisson statistics, suggesting that single liposomes can collide with the CFE and electroporate multiple times. With this study, we have resolved some outstanding questions surrounding electrochemical detection of liposomes while extending observations from giant unilamellar vesicles to 200 nm liposomes with high temporal resolution and sensitivity.

Pulsed Electric Fields Promote Liposome Buddings

Background: Liposomes have been a useful tool to analyze membrane behavior. Various studies have attempted to induce biological activities, for example, buddings, divisions, and endocytosis, on liposomes, focusing on lipid rafts that move along electric fields. Materials and Methods: Liposomes consisting of soybean lecithin, phosphatidylcholine, and cholesterol were prepared, with inner and outer liquid conductivities of 0.595 and 1.564 S/m, respectively. Results: We tried to induce buddings by pulsed electric fields (PEFs) on liposomes. Results demonstrated that 1.248 kV/cm, 400 μs PEF promoted postpulse liposome buddings, which were preceded by a membrane relaxation. Although a transient thick area (a lipid raft-like area) on the membrane just after PEF application preceded buddings, it was not the sufficient factor for buddings. Conclusion: We established a brief model as follows: 1.248 kV/cm, 400 μs PEF induced the lipid membrane relaxation without electroporation to trigger buddings. The current results could be a new frontier in bioelectrics.

Modification of lipid membrane compressibility induced by an electric field

Changes in membrane deformation and compressibility, induced by an external electric field, are investigated using coarse-grained martini force field simulations in a salt-free environment. We observe changes in the area of the membrane above a critical electric field. Below this value, the membrane compressibility modulus is found to decrease monotonically. For higher electric fields, the membrane projected area remains constant while the net interfacial area increases, with the corresponding compressibility moduli, show the opposite behavior. We find that the mechanical parameters, surface tension and bending modulus, of a freely floating membrane in the absence of explicit ions, are unaffected by the presence of the electric field. We believe these results have a bearing on our understanding of the electroformation of uncharged lipids in a salt-free environment.

An investigation into the critical tension of electroporation in anionic lipid vesicles

Irreversible electroporation (IRE) is a technique for the disruption of localized cells or vesicles by a series of short and high-frequency electric pulses which has been used for tissue ablation and treatment in certain diseases. It is well reported that IRE induces lateral tension in the membranes of giant unilamellar vesicles (GUVs). The GUVs are prepared by a mixture of anionic lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. Here the influence of DOPG mole fraction, XDOPG, on the critical tension of electroporation in GUVs has been investigated in sodium chloride-containing PIPES buffer. The critical tension decreases from 9.0 ± 0.3 to 6.0 ± 0.2 mN/m with the increase of XDOPG from 0.0 to 0.60 in the membranes of GUVs. Hence an increase in XDOPG greatly decreases the mechanical stability of membranes. We develop a theoretical equation that fits the XDOPG dependent normalized critical tension, and obtain a binding constant for the lipid-ion interaction of 0.75 M-1. The decrease in the energy barrier for formation of the nano-size nascent or prepore state, due to the increase in XDOPG, is the main factor explaining the decrease in critical tension of electroporation in vesicles.

Electrostatic effects on the electrical tension-induced irreversible pore formation in giant unilamellar vesicles

Irreversible electroporation (IRE) is a new technique in which a series of short pulses with high frequency electrical energy is applied on the targeted regions of cells or vesicles for their destruction or rupture formation. IRE induces lateral tension in the membranes of vesicles. We have investigated the electrostatic interaction effects on the constant electrical tension-induced rate constant of irreversible pore formation in the membranes of giant unilamellar vesicles (GUVs). The electrostatic interaction has been varied by changing the salt concentration in buffer and the surface charge density of membranes. The membranes of GUVs are synthesized by a mixture of negatively charged lipid dioleoylphosphatidylglycerol (DOPG) and neutral lipid dioleoylphosphatidylcholine (DOPC) using the natural swelling method. The rate constant of pore formation increases with the decrease of salt concentration in buffer along with the increase of surface charge density of membranes. The tension dependent probability of pore formation and the rate constant of pore formation are fitted to the theoretical equation, and obtained the line tension of membranes. The decrease in energy barrier of a prepore due to electrostatic interaction is the key factor causing an increase of rate constant of pore formation.

Influence of cholesterol on electroporation in lipid membranes of giant vesicles

Irreversible electroporation (IRE) is primarily a nonthermal ablative technology that uses a series of high-voltage and ultra-short pulses with high-frequency electrical energy to induce cell death. This paper presents the influence of cholesterol on the IRE-induced probability of pore formation and the rate constant of pore formation in giant unilamellar vesicles (GUVs). The GUVs are prepared by a mixture of dioleoylphosphatidylglycerol (DOPG), dioleoylphosphatidylcholine (DOPC) and cholesterol using the natural swelling method. An IRE signal of frequency 1.1 kHz is applied to the membranes of GUVs. The probability of pore formation and the rate constant of pore formation events are obtained using statistical analysis from several single GUVs. The time-dependent fraction of intact GUVs among all those examined is fitted to a single exponential decay function from where the rate constant of pore formation is calculated. The probability of pore formation and the rate constant of pore formation decreases with an increase in cholesterol content in the membranes of GUVs. Theoretical equations are fitted to the tension-dependent rate constant of pore formation and to the probability of pore formation, which allows us to obtain the line tension of membranes. The obtained line tension increases with an increase in cholesterol in the membranes. The increase in the energy barrier of the prepore state, due to the increase of cholesterol in membranes, is the main factor explaining the decrease in the rate constant of pore formation.

Kinetics of irreversible pore formation under constant electrical tension in giant unilamellar vesicles

Stretching in the plasma membranes of cells and lipid membranes of vesicles plays important roles in various physiological and physicochemical phenomena. Irreversible electroporation (IRE) is a minimally invasive non-thermal tumor ablation technique where a series of short electrical energy pulses with high frequency is applied to destabilize the cell membranes. IRE also induces lateral tension due to stretching in the membranes of giant unilamellar vesicles (GUVs). Here, the kinetics of irreversible pore formation under constant electrical tension in GUVs has been investigated. The GUVs are prepared by a mixture of dioleoylphosphatidylglycerol and dioleoylphosphatidylcholine using the natural swelling method. An IRE signal of frequency 1.1 kHz is applied to the GUVs through a gold-coated electrode system. Stochastic pore formation is observed for several 'single GUVs' at a particular constant tension. The time course of the fraction of intact GUVs among all the examined GUVs is fitted with a single-exponential decay function from which the rate constant of pore formation in the vesicle, k_p, is calculated. The value of k_p increases with an increase of membrane tension. An increase in the proportion of negatively charged lipids in a membrane gives a higher k_p. Theoretical equations are fitted to the tension-dependent k_p and to the probability of pore formation, which allows us to obtain the line tension of the membranes. The decrease in the energy barrier for formation of the nano-size nascent or prepore state, due to the increase in electrical tension, is the main factor explaining the increase of k_p.

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.

Mechanical characterization of vesicles and cells: A review

Vesicles perform many essential functions in all living organisms. They respond like a transducer to mechanical stress in converting the applied force into mechanical and biological responses. At the same time, both biochemical and biophysical signals influence the vesicular response in bearing mechanical loads. In recent years, liposomes, artificial lipid vesicles, have gained substantial attention from the pharmaceutical industry as a prospective drug carrier which can also serve as an artificial cell-mimetic system. The ability of these vesicles to enter through pores of even smaller size makes them ideal candidates for therapeutic agents to reach the infected sites effectively. Engineering of vesicles with desired mechanical properties that can encapsulate drugs and release as required is the prime challenge in this field. This requirement has led to the modifications of the composition of the bilayer membrane by adding cholesterol, sphingomyelin, etc. In this article, we review the manufacturing and characterization techniques of various artificial/synthetic vesicles. We particularly focus on the electric field-driven characterization techniques to determine different properties of vesicle and its membranes, such as bending rigidity, viscosity, capacitance, conductance, etc., which are indicators of their content and mobility. Similarities and differences between artificial vesicles, natural vesicles, and cells are highlighted throughout the manuscript since most of these artificial vesicles are intended for cell mimetic functions.

Polymerization-Induced Polymersome Fusion

The dynamic interactions of membranes, particularly their fusion and fission, are critical for the transmission of chemical information between cells. Fusion is primarily driven by membrane tension built up through membrane deformation. For artificial polymersomes, fusion is commonly induced via the external application of a force field. Herein, fusion-promoted development of anisotropic tubular polymersomes (tubesomes) was achieved in the absence of an external force by exploiting the unique features of aqueous ring-opening metathesis polymerization-induced self-assembly (ROMPISA). The out-of-equilibrium tubesome morphology was found to arise spontaneously during polymerization, and the composition of each tubesome sample (purity and length distribution) could be manipulated simply by targeting different core-block degrees of polymerization (DPs). The evolution of tubesomes was shown to occur via fusion of "monomeric" spherical polymersomes, evidenced most notably by a step-growth-like relationship between the fraction of tubular to spherical nano-objects and the average number of fused particles per tubesome (analogous to monomer conversion and DP, respectively). Fusion was also confirmed by Förster resonance energy transfer (FRET) studies to show membrane blending and confocal microscopy imaging to show mixing of the polymersome lumens. We term this unique phenomenon polymerization-induced polymersome fusion, which operates via the buildup of membrane tension exerted by the growing polymer chains. Given the growing body of evidence demonstrating the importance of nanoparticle shape on biological activity, our methodology provides a facile route to reproducibly obtain samples containing mixtures of spherical and tubular polymersomes, or pure samples of tubesomes, of programmed length. Moreover, the capability to mix the interior aqueous compartments of polymersomes during polymerization-induced fusion also presents opportunities for its application in catalysis, small molecule trafficking, and drug delivery.

Physicochemical considerations for bottom-up synthetic biology

The bottom-up construction of synthetic cells from molecular components is arguably one of the most challenging areas of research in the life sciences. We review the impact of confining biological systems in synthetic vesicles. Complex cell-like systems require control of the internal pH, ionic strength, (macro)molecular crowding, redox state and metabolic energy conservation. These physicochemical parameters influence protein activity and need to be maintained within limits to ensure the system remains in steady-state. We present the physicochemical considerations for building synthetic cells with dimensions ranging from the smallest prokaryotes to eukaryotic cells.

Compartments for Synthetic Cells: Osmotically Assisted Separation of Oil from Double Emulsions in a Microfluidic Chip

Liposomes are used in synthetic biology as cell-like compartments and their microfluidic production through double emulsions allows for efficient encapsulation of various components. However, residual oil in the membrane remains a critical bottleneck for creating pristine phospholipid bilayers. It has been discovered that osmotically driven shrinking leads to detachment of the oil drop. Separation inside a microfluidic chip has been realized to automate the procedure, which allows for controlled continuous production of monodisperse liposomes.

Ejection of Large Particulate Materials from Giant Unilamellar Vesicles Induced by Electropulsation

Electroporation or electropermealization is a technique to open pores in the lipid bilayer membrane of cells and vesicles transiently to increase its permeability to otherwise impermeable molecules. However, the upper size limit of the materials permeable through this operation has not been studied in the past. Here, we investigate the size of the material that can be released (ejected) from giant unilamellar vesicles (GUVs) upon electrical pulsation. We confirm that the volume of GUV shrinks in a stepwise manner upon periodical pulsation, in accordance with previous studies. When the same operation is applied to GUVs that encapsulate microbeads, we find that beads as large as 20 μm can be ejected across the membrane without rupturing the whole GUV structure. We also demonstrate that functional bioactive particulate materials, such as gel balls, vesicles, and cells can be encapsulated in and ejected from GUVs. We foresee that this phenomenon can be applied to precisely regulate the time and location of release of these particulate materials in the microenvironment.

Effects of electrically-induced constant tension on giant unilamellar vesicles using irreversible electroporation

Stretching in membranes of cells and vesicles plays important roles in various physiological and physicochemical phenomena. Irreversible electroporation (IRE) is the irreversible permeabilization of the membrane through the application of a series of electrical field pulses of micro- to millisecond duration. IRE induces lateral tension due to stretching in the membranes of giant unilamellar vesicles (GUVs). However, the effects of electrically induced (i.e., IRE) constant tension in the membranes of GUVs have not been investigated yet in detail. To explore the effects of electrically induced tension on GUVs, firstly a microcontroller-based IRE technique is developed which produces electric field pulses (332 V/cm) with pulse width 200 µs. Then the electrodeformation, electrofusion and membrane rupture of GUVs are investigated at various constant tensions in which the membranes of GUVs are composed of dioleoylphosphatidylglycerol (DOPG) and dioleoylphosphatidylcholine (DOPC). Stochastic electropore formation is observed in the membranes at an electrically induced constant tension in which the probability of pore formation is increased with the increase of tension from 2.5 to 7.0 mN/m. The results of pore formation at different electrically-induced constant tensions are in agreement with those reported for mechanically-induced constant tension. The decrease in the energy barrier of the pre-pore state due to the increase of electrically-induced tension is the main factor increasing the probability of electropore formation. These investigations help to provide an understanding of the complex behavior of cells/vesicles in electric field pulses and can form the basis for practical applications in biomedical technology.

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.

A finite difference method with subsampling for immersed boundary simulations of the capsule dynamics with viscoelastic membranes

The membrane or interfacial viscosity is an important property in many multiphase and biofluidic situations, such as the red blood cell dynamics and emulsion stability. The immersed boundary method (IBM), which incorporates the dynamic flow-membrane interaction via force distribution and velocity interpolation, has been extensively employed in simulations of such systems. Unfortunately, direct implementation of membrane viscosity in IBM suffers severe numerical instability, which causes an IBM calculation to break down before generating any useful results. Few attempts have been recently reported; however, several concerns exist in these attempts, such as the inconsistency to the classical definition of membrane viscosity, the inability to model the shear and dilatational viscosities separately, the unjustified mathematical formulations, and the complicated algorithms and computation. To overcome these concerns, in this paper, we propose a finite difference approach for implementing membrane viscosity in immersed boundary simulations. The viscous stress is obtained via finite difference approximations to the differential strain-stress relationship, with the help of a subsampling scheme to reduce the numerical noise in the calculated strain rates. This simple method has also avoided the complicated matrix calculations in previous attempts, and hence, a better computational efficiency is expected. Detailed mathematical description of the method and key steps for its implementation in immersed boundary programs are provided. Validation and illustration calculations are performed, and our results are compared with analytical solutions and previous publications with satisfactory agreement. The influences of membrane mesh resolution and simulation time step are also examined; and the results show no indication that our finite difference method has downgraded the general IBM accuracy. Based on these simulations and analysis, we believe that our method would be a better choice for future IBM simulations of capsule dynamics with viscoelastic membranes.

Response of an actin network in vesicles under electric pulses

We study the role of a biomimetic actin network during the application of electric pulses that induce electroporation or electropermeabilization, using giant unilamellar vesicles (GUVs) as a model system. The actin cortex, a subjacently attached interconnected network of actin filaments, regulates the shape and mechanical properties of the plasma membrane of mammalian cells, and is a major factor influencing the mechanical response of the cell to external physical cues. We demonstrate that the presence of an actin shell inhibits the formation of macropores in the electroporated GUVs. Additionally, experiments on the uptake of dye molecules after electroporation show that the actin network slows down the resealing process of the permeabilized membrane. We further analyze the stability of the actin network inside the GUVs exposed to high electric pulses. We find disruption of the actin layer that is likely due to the electrophoretic forces acting on the actin filaments during the permeabilization of the GUVs. Our findings on the GUVs containing a biomimetic network provide a step towards understanding the discrepancies between the electroporation mechanism of a living cell and its simplified model of the empty GUV.