Molecular dynamics simulations of pore formation in stretched phospholipid/cholesterol bilayers

Changes in free energy barrier for water permeation by stretch-induced phase transitions in phospholipid/cholesterol bilayers

Water permeation through phospholipid/cholesterol bilayers is the key to understanding tension-induced rupture of biological cell membranes. We performed molecular dynamics simulations of stretched phospholipid/cholesterol bilayers to investigate changes in the free energy profile of water molecules across the bilayer and the lipid structure responsible for water permeation. We modeled stretching of the bilayer by applying areal strain. In stretched phospholipid/cholesterol bilayers, the hydrophobic tail of the phospholipids became disordered and the free energy barrier to water permeation decreased. Upon exceeding the critical areal strain, a phase transition to an interdigitated gel phase occurred before rupture, and the hydrophobic tail ordering as well as the free energy barrier were restored. In pure phospholipid bilayers, we did not observe such recoveries. These transient recoveries in the phospholipid/cholesterol bilayer suppressed water permeation and membrane rupture, followed by an increase in the critical areal strain at which the bilayer ruptured. This result agrees with experimental results and provides a reasonable molecular mechanism for the toughness of phospholipid/cholesterol bilayers under tension.Communicated by Ramaswamy H. Sarma.

Influence of phospholipid head and tail molecular structures on cell membrane mechanical response under tension

Biological cell membranes are primarily comprised of a diverse lipid bilayer with multiple phospholipid (lipid) types, each of which is comprised of a hydrophilic headgroup and two hydrophobic hydrocarbon tails. The lipid type determines the molecular structure of head and tail groups, which can affect membrane mechanics at nanoscale and subsequently cell viability under mechanical loading. Hence, using molecular dynamics simulations, the current study investigated seven membrane phospholipids and the effect of their structural differences on physical deformation, mechanoporation damage, and mechanical failure of the membranes under tension. The inspected phospholipids showed similar yield stresses and strains, as well as pore evolution and damage, but significantly different failure strains. In general, failure occurred at a lower strain for lipids with a larger equilibrium area per lipid. The obtained results suggest that larger headgroup structure, greater degree of unsaturation, and tail-length asymmetry influenced the phospholipids' ability to pack against each other, increased the fluidity and equilibrium area per lipid of the membrane, and resulted in lower failure strain. Overall, this study provides insights on how different phospholipid structures affect membrane physical responses at the molecular level and serves as a reference for future studies of more complex membrane systems with intricate biophysical properties.

Feasibility Study for the Use of Gene Electrotransfer and Cell Electrofusion as a Single-Step Technique for the Generation of Activated Cancer Cell Vaccines

Cell-based therapies hold great potential for cancer immunotherapy. This approach is based on manipulation of dendritic cells to activate immune system against specific cancer antigens. For the development of an effective cell vaccine platform, gene transfer, and cell fusion have been used for modification of dendritic or tumor cells to express immune (co)stimulatory signals and to load dendritic cells with tumor antigens. Both, gene transfer and cell fusion can be achieved by single technique, a cell membrane electroporation. The cell membrane exposed to external electric field becomes temporarily permeable, enabling introduction of genetic material, and also fusogenic, enabling the fusion of cells in the close contact. We tested the feasability of combining gene electrotransfer and electrofusion into a single-step technique and evaluated the effects of electroporation buffer, pulse parameters, and cell membrane fluidity for single or combined method of gene delivery or cell fusdion. We determined the percentage of fused cells expressing green fluorescence protein (GFP) in a murine cell model of melanoma B16F1, cell line used in our previous studies. Our results suggest that gene electrotransfer and cell electrofusion can be applied in a single step. The percentage of viable hybrid cells expressing GFP depends on electric pulse parameters and the composition of the electroporation buffer. Furthermore, our results suggest that cell membrane fluidity is not related to the efficiency of the gene electrotransfer and electrofusion. The protocol is compatible with microfluidic devices, however further optimization of electric pulse parameters and buffers is still needed.

Effect of Force Field Resolution on Membrane Mechanical Response and Mechanoporation Damage under Deformation Simulations

Damage induced by transient disruption and mechanoporation in an intact cell membrane is a vital nanoscale biomechanical mechanism that critically affects cell viability. To complement experimental studies of mechanical membrane damage and disruption, molecular dynamics (MD) simulations have been performed at different force field resolutions, each of which follows different parameterization strategies and thus may influence the properties and dynamics of membrane systems. Therefore, the current study performed tensile deformation MD simulations of bilayer membranes using all-atom (AA), united-atom (UA), and coarse-grained Martini (CG-M) models to investigate how the damage biomechanics differs across atomistic and coarse-grained (CG) simulations. The mechanical response and mechanoporation damage were qualitatively similar but quantitatively different in the three models, including some progressive changes based on the coarse-graining level. The membranes yielded and reached ultimate strength at similar strains; however, the coarser systems exhibited lower average yield stresses and failure strains. The average failure strain in the UA model was approximately 7% lower than the AA, and the CG-M was 20% lower than UA and 27% lower than AA. The CG systems also nucleated a higher number of pores and larger pores, which resulted in higher damage during the deformation process. Overall, the study provides insight on the impact of force field-a critical factor in modeling biomolecular systems and their interactions-in inspecting membrane mechanosensitive responses and serves as a reference for justifying the appropriate force field for future studies of more complex membranes and more diverse biomolecular assemblies.

Studying the roles of salt ions in the pore initiation and closure stages in the biomembrane electroporation

Electroporation shows great potential in biology and biomedical applications. However, there is still a lack of reliable protocol for cell electroporation to achieve a high perforation efficiency due to the unclear influence mechanism of various factors, especially the salt ions in buffer solution. The tiny membrane structure of a cell and the electroporation scale make it difficult to monitor the electroporation process. In this study, we used both molecular dynamics (MD) simulation and experimental methods to explore the influence of salt ions on the electroporation process. Giant unilamellar vesicles (GUVs) were constructed as the model, and sodium chloride (NaCl) was selected as the representative salt ion in this study. The results show that the electroporation process follows lag-burst kinetics, where the lag period first appears after applying the electric field, followed by a rapid pore expansion. For the first time, we find that the salt ion plays opposite roles in different stages of the electroporation process. The accumulation of salt ions near the membrane surface provides an extra potential to promote the pore initiation, while the charge screening effect of the ions within the pore increases the line tension of the pore to induce the instability of the pore and lead to the closure. The GUV electroporation experiments obtain qualitatively consistent results with MD simulations. This work can provide guidance for the selection of parameters for cell electroporation process.

Effect of the cholesterol on electroporation of planar lipid bilayer

Electroporation threshold depends on the membrane composition, with cholesterol being one of its key components already studied in the past, but the results were inconclusive. The aim of our study was to determine behaviour of planar lipid bilayers with varying cholesterol concentrations under electric field. This would give us a better insight into cholesterol effect on membrane properties during electroporation process, since cholesterol is one of the major components of biological membranes and plays a crucial role in membrane organisation, dynamics, and function. Planar lipid bilayers were prepared from phosphatidylcholine lipids with 0, 20, 30, 50 and 80 mol% cholesterol. Capacitance was measured using the discharge method. Results show no statistical difference of cBLM between the cholesterol concentrations. Breakdown voltage Ubr of planar lipid bilayers was measured by means of linear rising voltage with seven different slopes. Obtained results were fitted to a strength-duration curve, where parameter Ubrmin represents minimal breakdown voltage, and parameter τRC represents the inclination of the strength-duration curve. Adding cholesterol to planar lipid bilayer gradually increased its Ubrmin until 50 mol% cholesterol concentration. Afterwards at 80 mol% Ubrmin does not further increase, in fact it reduces by 20% of the Ubrmin at 50 mol% cholesterol concentration.

Impact of Lipid Peroxidation on the Response of Cell Membranes to High-Speed Equibiaxial Stretching: A Computational Study

The difference between diseased and healthy cellular membranes in response to mechanical stresses is crucial for biology, as well as in the development of medical devices. However, the biomolecular mechanisms by which mechanical stresses interact with diseased cellular components remain largely unknown. In this work, we focus on the response of diseased cellular membranes with lipid peroxidation to high-speed tensile loadings. We find that the critical areal strain (ξ_c, when the pore forms) is highly sensitive to lipid peroxidation. For example, ξ_c of a fully oxidized bilayer is only 64 and 69% of the nonoxidized one at the stretching speed of 0.1 and 0.6 m/s, respectively. ξ_c decreases with the increase in the oxidized lipid ratio, regardless of the speeds. Also, the critical rupture tension of membranes exhibits a similar change. It is obvious that the oxidized membranes are more easily damaged than normal ones by high-speed stretching, which coincides with experimental findings. The reason is that peroxidation introduces a polar group to the tail of lipids, increases the hydrophilicity of tails, and warps the tails to the membrane-water interface, which causes loose accumulation and disorder of lipid tails. This can be deduced from the variation in the area per lipid and order parameter. In addition, the lowering stretching modulus and line tension of membranes (i.e., softening) after lipid peroxidation is also a significant factor. We reveal the difference between the peroxidized (diseased) and normal membrane in response to high-speed stretching, give the ξ_c value in the pore formation of membranes and analyze the influence of the stretching speed, peroxidation ratio, and molecular structure of phospholipids. We hope that the molecular-level information will be useful for the development of biological and medical devices in the future.

Elastic property of sickle cell anemia and sickle cell trait red blood cells

SIGNIFICANCE

We introduce a model for better calibration of the trapping force using an equal but oppositely directed drag force acting on a trapped red blood cell (RBC). We demonstrate this approach by studying RBCs' elastic properties from deidentified sickle cell anemia (SCA) and sickle cell trait (SCT) blood samples.

AIM

A laser trapping (LT) force was formulated and analytically calculated in a cylindrical model. Using this trapping force relative percent difference, the maximum (longitudinal) and minimum (transverse) radius rate and stiffness were used to study the elasticity.

APPROACH

The elastic property of SCA and SCT RBCs was analyzed using LT technique with computer controlled piezo-driven stage, in order to trap and stretch the RBCs.

RESULTS

For all parameters, the results show that the SCT RBC samples have higher elastic property than the SCA RBCs. The higher rigidity in the SCA cell may be due to the lipid composition of the membrane, which was affected by the cholesterol concentration.

CONCLUSIONS

By developing a theoretical model for different trapping forces, we have also studied the elasticity of RBCs in SCT (with hemoglobin type HbAS) and in SCA (with hemoglobin type HbSS). The results for the quantities describing the elasticity of the cells consistently showed that the RBCs in the SCT display lower rigidity and higher deformability than the RBCs with SCA.

Effects of cholesterol on bilayers with various degrees of unsaturation of their phospholipid tails under mechanical stress

Cholesterol is one of the essential components of the cell membrane. It has a significant influence on various mechanical properties of biomembranes, such as fluidity and elasticity, which have attracted much attention. It is also well known that the concentration of cholesterol affects the mechanical strength of cell membranes. In this paper, we aim to explore the influence of the degree of unsaturation of phospholipid tails on the concentration-effect of cholesterol. Three different phospholipids (DPPC, DIPC and DAPC) were selected as the respective main components of the bilayers and several concentrations of cholesterol were also added to the systems. Our coarse-grained molecular dynamics simulations show that as the cholesterol concentration increases, the saturated phospholipid bilayer is first strengthened, by increasing the rupture tension from 68.9 to 110 mN m-1, and then weakened. The non-monotonic concentration-effect gradually decreases as the degree of unsaturation of the phospholipid tails increases, and in particular, the mechanical strength of the DAPC bilayer hardly changes. The results suggest that cholesterol does not influence a bilayer composed of highly unsaturated phospholipids. Furthermore, lateral density distributions reveal that the distribution of cholesterol in the bilayer is related to the carbon tail unsaturation of the phospholipids.

Hierarchy of interactions dictating the thermodynamics of real cell membranes: Following the insulin secretory granules paradigm up to fifteen-components vesicles

A fifteen-components model membrane that reflected the 80 % of phospholipids present in Insulin Secretory Granules was obtained and thermodynamic exploitation was performed, through micro-DSC, in order to assess the synergic contributions to the stability of a mixed complex system very close to real membranes. Simpler systems were also stepwise investigated, to complete a previous preliminary study and to highlight a hierarchy of interactions that can be now summarized as phospholipid tail unsaturation > phospholipid tail length > phospholipid headgroup > membrane curvature. In particular, Small Unilamellar Vesicles (SUVs) that consisted in phospholipids with different headgroups (choline, ethanolamine and serine), was step by step considered, following inclusion of sphingomyelins and lysophosphatidylcholines together with a more complete fatty acids distribution characterizing the phospholipid bilayer of the Insulin Secretory Granules. The inclusion of cholesterol was finally considered and the influence of three FFAs (stearic, oleic and elaidic acids) was investigated in comparison with simpler systems, highlighting the magnitude of the effects on such a detailed membrane in the frame of Type 2 Diabetes Mellitus alterations.

Role of lipid composition on the structural and mechanical features of axonal membranes: a molecular simulation study

The integrity of cellular membranes is critical for the functionality of axons. Failure of the axonal membranes (plasma membrane and/or myelin sheath) can be the origin of neurological diseases. The two membranes differ in the content of sphingomyelin and galactosylceramide lipids. We investigate the relation between lipid content and bilayer structural-mechanical properties, to better understand the dependency of membrane properties on lipid composition. A sphingomyelin/phospholipid/cholesterol bilayer is used to mimic a plasma membrane and a galactosylceramide/phospholipid/cholesterol bilayer to mimic a myelin sheath. Molecular dynamics simulations are performed at atomistic and coarse-grained levels to characterize the bilayers at equilibrium and under deformation. For comparison, simulations of phospholipid and phospholipid/cholesterol bilayers are also performed. The results clearly show that the bilayer biomechanical and structural features depend on the lipid composition, independent of the molecular models. Both galactosylceramide or sphingomyelin lipids increase the order of aliphatic tails and resistance to water penetration. Having 30% galactosylceramide increases the bilayers stiffness. Galactosylceramide lipids pack together via sugar-sugar interactions and hydrogen-bond phosphocholine with a correlated increase of bilayer thickness. Our findings provide a molecular insight on role of lipid content in natural membranes.

Mechanisms underlying interactions between PAMAM dendron-grafted surfaces with DPPC membranes

Biofouling is a pervasive problem which demands the creation of smart, antifouling surfaces. Towards this end, we examine the interactions between a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer and a polyamidoamine (PAMAM) dendron-grafted surface. In addition, we investigate the impact of dendron generation on the system behavior. To resolve the multiscale dynamical processes occurring over a large spatial scale, we employ Molecular Dynamics simulations with a coarse-grained implicit solvent force field. Our results demonstrate the transient and equilibrium system dynamics to be determined by the PAMAM dendron generation along with the underlying mechanisms. Higher generation dendrons are observed to favor penetration of the DPPC molecules into the dendron branches, thereby enabling sustained interactions between the membrane and the dendron-grafted surface. Under equilibrium, the membrane adopts a bowl-shaped morphology whose dimensions are determined by the dendron generation and density of interactions. The results from our study can be used to guide the design of novel surfaces with selective antifouling properties which can prevent the adsorption of microorganisms onto lipid membranes.

The interactions between a DPPC lipid membrane and a PAMAM dendron-grafted surface.

Branching pattern effect and co-assembly with lipids of amphiphilic Janus dendrimersomes

The influence of the branching patterns on the membrane properties of Janus dendrimers in water has been investigated by dissipative particle dynamics simulations.

Molecular dynamics simulations of heterogeneous cell membranes in response to uniaxial membrane stretches at high loading rates

The chemobiomechanical signatures of diseased cells are often distinctively different from that of healthy cells. This mainly arises from cellular structural/compositional alterations induced by disease development or therapeutic molecules. Therapeutic shock waves have the potential to mechanically destroy diseased cells and/or increase cell membrane permeability for drug delivery. However, the biomolecular mechanisms by which shock waves interact with diseased and healthy cellular components remain largely unknown. By integrating atomistic simulations with a novel multiscale numerical framework, this work provides new biomolecular mechanistic perspectives through which many mechanosensitive cellular processes could be quantitatively characterised. Here we examine the biomechanical responses of the chosen representative membrane complexes under rapid mechanical loadings pertinent to therapeutic shock wave conditions. We find that their rupture characteristics do not exhibit significant sensitivity to the applied strain rates. Furthermore, we show that the embedded rigid inclusions markedly facilitate stretch-induced membrane disruptions while mechanically stiffening the associated complexes under the applied membrane stretches. Our results suggest that the presence of rigid molecules in cellular membranes could serve as “mechanical catalysts” to promote the mechanical destructions of the associated complexes, which, in concert with other biochemical/medical considerations, should provide beneficial information for future biomechanical-mediated therapeutics.

Caveolins and cavins in the trafficking, maturation, and degradation of caveolae: implications for cell physiology

Caveolins (Cavs) are ~20 kDa scaffolding proteins that assemble as oligomeric complexes in lipid raft domains to form caveolae, flask-shaped plasma membrane (PM) invaginations. Caveolae ("little caves") require lipid-lipid, protein-lipid, and protein-protein interactions that can modulate the localization, conformational stability, ligand affinity, effector specificity, and other functions of proteins that are partners of Cavs. Cavs are assembled into small oligomers in the endoplasmic reticulum (ER), transported to the Golgi for assembly with cholesterol and other oligomers, and then exported to the PM as an intact coat complex. At the PM, cavins, ~50 kDa adapter proteins, oligomerize into an outer coat complex that remodels the membrane into caveolae. The structure of caveolae protects their contents (i.e., lipids and proteins) from degradation. Cellular changes, including signal transduction effects, can destabilize caveolae and produce cavin dissociation, restructuring of Cav oligomers, ubiquitination, internalization, and degradation. In this review, we provide a perspective of the life cycle (biogenesis, degradation), composition, and physiologic roles of Cavs and caveolae and identify unanswered questions regarding the roles of Cavs and cavins in caveolae and in regulating cell physiology.¹.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

The Interaction between Influenza HA Fusion Peptide and Transmembrane Domain Affects Membrane Structure

Viral glycoproteins, such as influenza hemagglutinin (HA) and human immunodeficiency virus gp41, are anchored by a single helical segment transmembrane domain (TMD) on the viral envelope membrane. The fusion peptides (FP) of the glycoproteins insert into the host membrane and initiate membrane fusion. Our previous study showed that the FP or TMD alone perturbs membrane structure. Interaction between the influenza HA FP and TMD has previously been shown, but its role is unclear. We used PC spin labels dipalmitoylphospatidyl-tempo-choline (on the headgroup), 5PC and 14PC (5-C and 14-C positions on the acyl chain) to detect the combined effect of FP-TMD interaction by titrating HA FP to TMD-reconstituted 1,2-dimyristoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol)/cholesterol lipid bilayers using electron spin resonance. We found that the FP-TMD increases the lipid order at all positions, which has a greater lipid ordering effect than the sum of the FP or TMD alone, and this effect reaches deeper into the membranes. Although HA-mediated membrane fusion is pH dependent, this combined effect is observed at both pH 5 and pH 7. In addition to increasing lipid order, multiple components are found for 5PC at increased concentration of FP-TMD, indicating that distinct domains are induced. However, the mutation of Gly1 in the FP and L187 in the TMD eliminates the perturbations, consistent with their fusogenic phenotypes. Electron spin resonance on spin-labeled peptides confirms these observations. We suggest that this interaction may provide a driving force in different stages of membrane fusion: initialization, transition from hemifusion stalk to transmembrane contact, and fusion pore formation.

Effects of Stretching Speed on Mechanical Rupture of Phospholipid/Cholesterol Bilayers: Molecular Dynamics Simulation

Rupture of biological cell membrane under mechanical stresses is critical for cell viability. It is triggered by local rearrangements of membrane molecules. We investigated the effects of stretching speed on mechanical rupture of phospholipid/cholesterol bilayers using unsteady molecular dynamics simulations. We focused on pore formation, the trigger of rupture, in a 40 mol% cholesterol-including bilayer. The unsteady stretching was modeled by proportional and temporal scaling of atom positions at stretching speeds from 0.025 to 30 m/s. The effects of the stretching speed on the critical areal strain, where the pore forms, is composed of two regimes. At low speeds (<1.0 m/s), the critical areal strain is insensitive to speed, whereas it significantly increases at higher speeds. Also, the strain is larger than that of a pure bilayer, regardless of the stretching speeds, which qualitatively agrees with available experimental data. Transient recovery of the cholesterol and phospholipid molecular orientations was evident at lower speeds, suggesting the formation of a stretch-induced interdigitated gel-like phase. However, this recovery was not confirmed at higher speeds or for the pure bilayer. The different responses of the molecular orientations may help explain the two regimes for the effect of stretching speed on pore formation.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.