Molecular structure of membrane tethers

Molecular Mechanism of Autophagosome-Lysosome Fusion in Mammalian Cells

In eukaryotes, targeting intracellular components for lysosomal degradation by autophagy represents a catabolic process that evolutionarily regulates cellular homeostasis. The successful completion of autophagy initiates the engulfment of cytoplasmic materials within double-membrane autophagosomes and subsequent delivery to autolysosomes for degradation by acidic proteases. The formation of autolysosomes relies on the precise fusion of autophagosomes with lysosomes. In recent decades, numerous studies have provided insights into the molecular regulation of autophagosome-lysosome fusion. In this review, an overview of the molecules that function in the fusion of autophagosomes with lysosomes is provided. Moreover, the molecular mechanism underlying how these functional molecules regulate autophagosome-lysosome fusion is summarized.

Biophysical Interaction Landscape of Mycobacterial Mycolic Acids and Phenolic Glycolipids with Host Macrophage Membranes

Lipidic adjuvant formulations consisting of immunomodulatory mycobacterial cell wall lipids interact with host cells following administration. The impact of this cross-talk on the host membrane's structure and function is rarely given enough consideration but is imperative to rule out nonspecific perturbation underlying the adjuvant. In this work, we investigated changes in the plasma membranes of live mammalian cells after exposure to mycobacterial mycolic acid (MA) and phenolic glycolipids, two strong candidates for lipidic adjuvant therapy. We found that phenolic glycolipid 1 softened the plasma membrane, lowering membrane tension and stiffness, but MA did not significantly change the membrane characteristics. Further, phenolic glycolipid 1 had a fluidizing impact on the host plasma membrane, increasing the fluidity and the abundance of fluid-ordered-disordered coexisting lipid domains. Notably, lipid diffusion was not impacted. Overall, MA and, to a lesser extent, phenolic glycolipid 1, due to minor disruption of host cell membranes, may serve as appropriate lipids in adjuvant formulations.

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

From reversible to irreversible: When the membrane nanotube pearling is coupled with phase separation

Membrane nanotubes, which are ubiquitous in biology and act as channels maintaining transport between different cells and organelles, readily undergo pearling in response to external stimuli. Membrane nanotube pearling involves generation of heterogeneous curvature coupled with redistribution of membrane components that may interfere with the shape recovery of pearled nanotubes. However, the mechanism underlying such delicate process remains unclear and difficult to study at the molecular scale in vivo. By means of molecular dynamics simulation, here we investigate pearling of multi-component membrane nanotubes and reversibility through manipulating system temperature and osmotic pressure. With the equilibrium shape of membrane nanotubes controlled by the osmotic pressure, our results demonstrate that the process of membrane nanotube pearling can be reversible or irreversible, depending on the phase segregation state. For the pearled nanotube releasing high surface energy, different lipid components redistribute along the tube axial direction. Lipids with unsaturated tails prefer gathering at the high-curvature shrinking region, whereas the swelling region is constituted by saturated lipids forming the liquid-ordered phase of a higher bending rigidity. Such curvature sensitive phase segregation minimizes the system free energy by reducing both the membrane bending energy and line tension at the phase boundary. As such, the pearled nanotube fails to recover its shape upon retracting stimuli, suggesting irreversibility of the membrane nanotube pearling coupled with phase separation. Given importance of membrane nanotube pearling in various cellular activities, these results provide a new mechanism of controlling equilibrium shapes of membrane nanotubes in complex cellular environment.

FAM134B-RHD Protein Clustering Drives Spontaneous Budding of Asymmetric Membranes

Living cells constantly remodel the shape of their lipid membranes. In the endoplasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attractions. At a critical size, as in a nucleation process, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetic barrier for membrane remodeling.

Influence of membrane-cortex linkers on the extrusion of membrane tubes

The cell membrane is an inhomogeneous system composed of phospholipids, sterols, carbohydrates, and proteins that can be directly attached to underlying cytoskeleton. The protein linkers between the membrane and the cytoskeleton are believed to have a profound effect on the mechanical properties of the cell membrane and its ability to reshape. Here, we investigate the role of membrane-cortex linkers on the extrusion of membrane tubes using computer simulations and experiments. In simulations, we find that the force for tube extrusion has a nonlinear dependence on the density of membrane-cortex attachments: at a range of low and intermediate linker densities, the force is not significantly influenced by the presence of the membrane-cortex attachments and resembles that of the bare membrane. For large concentrations of linkers, however, the force substantially increases compared with the bare membrane. In both cases, the linkers provided membrane tubes with increased stability against coalescence. We then pulled tubes from HEK cells using optical tweezers for varying expression levels of the membrane-cortex attachment protein Ezrin. In line with simulations, we observed that overexpression of Ezrin led to an increased extrusion force, while Ezrin depletion had a negligible effect on the force. Our results shed light on the importance of local protein rearrangements for membrane reshaping at nanoscopic scales.

Superspreading - Has the mystery been unraveled?

Superspreading is a fascinating phenomenon first observed about 30 years ago with dilute solutions of trisiloxane surfactants on hydrophobic substrates. Although many groups all over the world have contributed considerably to solve the scientific challenges involved, the reasons why only some trisiloxane surfactants promote superspreading, whereas others of similar chemical structure behave more like ordinary surfactants, has remained a mystery up to now. A number of original papers and reviews on superspreading have been published in recent years. The driving force still proposed today is most often Marangoni flow. This is, however, in contradiction with recent results showing that superspreading only starts after a surface tension gradient between apex and leading edge has been eliminated. From foam film experiments unrelated to wetting, there is evidence for "dangling" bilayers attached to the air/water interface only in case of the superspreading trisiloxane surfactants. By combining this and other published experimental findings, a new hypothesis of the mode of action is put forward: Advancing by "rolling action" at the leading edge, and the supply of surfactant by "unzippering" of the dangling bilayers all over the surface of the drop; this hypothesis even fulfills basic thermodynamic requirements.

Determination of the Mechanical Properties of Model Lipid Bilayers Using Atomic Force Microscopy Indentation

By conducting a systematic study of model lipid membranes using the atomic force microscopy (AFM) indentation, we demonstrate the importance of an experimental protocol on the determination of their mechanical parameters. We refine the experimental approach by analyzing the influence of the contact mechanics models used to process the data, substrate preparation, and indenter geometry. We show that both bending rigidity and area compressibility can be determined from a single AFM indentation measurement.

The transbilayer distribution of polyunsaturated phospholipids determines their facilitating effect on membrane deformation

In the cell, membrane deformation and fission (collectively referred to as 'budding') is driven by specific protein machineries but is also influenced by lipid composition. We previously reported that phospholipids with polyunsaturated acyl chains facilitate membrane budding because they adapt their shape to membrane curvature, thereby decreasing membrane bending rigidity. The facilitating effect of polyunsaturated lipids was observed in experiments and simulations performed on membranes where the two bilayer leaflets had the same lipid composition. However, biological membranes are generally asymmetric. Here, we present coarse-grained molecular dynamics simulations on asymmetric phospholipid bilayers undergoing deformation via a pulling force along the bilayer normal. One leaflet contains monounsaturated C18:0-C18:1-phospholipids, whereas the opposite leaflet contains polyunsaturated C18:0-C22:6-phospholipids. When present in the monolayer orientated towards the pulling force and thereby in the convex face of the forming tube, C18:0-C22:6-phospholipids facilitate membrane tubulation. In contrast, C18:0-C22:6-phospholipids in the concave face of the tube have no effect. Analysis of lipid shape indicates that these contrasting effects arise from the superior ability of polyunsaturated phospholipids to swell in the convex leaflet, whereas mono and polyunsaturated phospholipids behave similarly in the concave leaflet. The leaflet-dependent effect of polyunsaturated phospholipids matches well their asymmetric distribution in biological membranes, notably in synaptic vesicles, which are produced by the fastest budding event in the body.

How a lipid bilayer membrane responds to an oscillating nanoparticle: Promoted membrane undulation and directional wave propagation

Mechanical forces acting on a plasma membrane are of essential importance to cellular functioning via inducing delicate change of the membrane shape with the underlying mechanism yet to be elucidated. Here, we introduce an oscillating nanoparticle (NP) interaction with a lipid bilayer membrane, using the coarse-grained simulation to investigate the dynamic membrane response to constrained mechanical stimulation, which is ubiquitous in biology. Our results demonstrate that, the membrane responds to an oscillating NP by generating nanoscale undulation waves, which immediately propagate through the membrane. In dynamics, propagation of the generated membrane undulation waves always starts from flattening of the region where the NP locates, thus producing a lateral force to propel the waves away from the point of stimulation. The speed of membrane undulation wave propagation is proportional to that of NP oscillation and accelerated by increasing the integral membrane surface tension, suggesting that both the membrane bending and stretching contribute to the energy driving the unique response of membrane undulation wave propagation.

Phospholipid Scramblases Remodel the Shape of Asymmetric Membranes

The cell membrane and many organellar membranes are asymmetric and highly curved. In experiments, it is challenging to reconstitute and characterize membranes that differ in the lipid composition of their leaflets. Here we use molecular dynamics simulations to study the large-scale membrane shape changes associated with lipid shuttling between asymmetric leaflets. We exploit leaflet asymmetry to create a stable, near-spherical vesicle bud connected to a flat bilayer under periodic boundary conditions. Then we demonstrate how the lipid scramblase nhTMEM16 relaxes the lipid-number asymmetry. By mediating the flipping of lipids, this transmembrane protein dissipates the mechanochemical gradient between the leaflets and drives a large-scale membrane reorganization, converting the vesicle bud into a flat membrane. Our procedure to exploit bilayer asymmetry for simulations of highly curved membranes can be used to study the function of other lipid transporters and membrane-shaping proteins.

Surface Shear Viscosity and Interleaflet Friction from Nonequilibrium Simulations of Lipid Bilayers

Nonequilibrium simulation protocols based on shear deformations are applied to determine the surface viscosity and interleaflet friction of lipid bilayers. At high shear rates, a non-Newtonian shear thinning regime is observed, but lower shear rates yield a Newtonian plateau and results that are consistent with equilibrium measurements based on fluctuation-dissipation theorems. Application to all-atom bilayers modeled with the CHARMM36 parameter set yields values for the surface viscosity that are consistent with microscopic measurements based on membrane protein diffusion but are approximately 10 times lower than more macroscopic experimental measurements. The interleaflet friction is about 10 times lower than experimental measurements. Trends across different lipids, temperatures, and ternary liquid-disordered phase mixtures produce results that are consistent with experimental diffusion constants. Application of the protocol to the liquid-ordered phase fails to yield a Newtonian plateau, suggesting more complex rheology.

The late stage of COPI vesicle fission requires shorter forms of phosphatidic acid and diacylglycerol

Studies on vesicle formation by the Coat Protein I (COPI) complex have contributed to a basic understanding of how vesicular transport is initiated. Phosphatidic acid (PA) and diacylglycerol (DAG) have been found previously to be required for the fission stage of COPI vesicle formation. Here, we find that PA with varying lipid geometry can all promote early fission, but only PA with shortened acyl chains promotes late fission. Moreover, diacylglycerol (DAG) acts after PA in late fission, with this role of DAG also requiring shorter acyl chains. Further highlighting the importance of the short-chain lipid geometry for late fission, we find that shorter forms of PA and DAG promote the vesiculation ability of COPI fission factors. These findings advance a general understanding of how lipid geometry contributes to membrane deformation for vesicle fission, and also how proteins and lipids coordinate their actions in driving this process.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.

Mechanics of the Formation, Interaction, and Evolution of Membrane Tubular Structures

Membrane nanotubes, also known as membrane tethers, play important functional roles in many cellular processes, such as trafficking and signaling. Although considerable progresses have been made in understanding the physics regulating the mechanical behaviors of individual membrane nanotubes, relatively little is known about the formation of multiple membrane nanotubes due to the rapid occurring process involving strong cooperative effects and complex configurational transitions. By exerting a pair of external extraction upon two separate membrane regions, here, we combine molecular dynamics simulations and theoretical analysis to investigate how the membrane nanotube formation and pulling behaviors are regulated by the separation between the pulling forces and how the membrane protrusions interact with each other. As the force separation increases, different membrane configurations are observed, including an individual tubular protrusion, a relatively less deformed protrusion with two nanotubes on its top forming a V shape, a Y-shaped configuration through nanotube coalescence via a zipper-like mechanism, and two weakly interacting tubular protrusions. The energy profile as a function of the separation is determined. Moreover, the directional flow of lipid molecules accompanying the membrane shape transition is analyzed. Our results provide new, to our knowledge, insights at a molecular level into the interaction between membrane protrusions and help in understanding the formation and evolution of intra- and intercellular membrane tubular networks involved in numerous cell activities.

BUMPy: A Model-Independent Tool for Constructing Lipid Bilayers of Varying Curvature and Composition

Molecular dynamics is a powerful tool to investigate atomistic and mesoscopic phenomena in lipid bilayer systems. These studies have progressed with the advent of increased computational power, and efforts are now increasingly being directed toward investigating the role of curvature and bilayer morphology, as these are critical features of biological processes. Computational studies of lipid bilayers benefit from tools that can create starting configurations for molecular dynamics simulations, but the majority of such tools are restricted to generating flat bilayers. Generating curved bilayer configurations comes with practical complications and potential ramifications on physical properties in the simulated system if the bilayer is initiated in a high-strain state. We present a new tool for creating curved lipid bilayers that combines flexibility of shape, force field, model resolution, and bilayer composition. A key aspect of our approach is the use of the monolayer pivotal plane location to accurately estimate interleaflet area differences in a curved bilayer. Our tool is named BUMPy (Building Unique Membranes in Python), is written in Python, is fast, and has a simple command line interface.

Implicit-solvent dissipative particle dynamics force field based on a four-to-one coarse-grained mapping scheme

A new set of efficient solvent-free dissipative particle dynamics (DPD) force fields was developed for phospholipids and peptides. To enhance transferability, this model maps around four heavy atoms and their connected hydrogen atoms into a coarse-grained elementary bead based on functional group. The effective hybrid potential between any pair of beads is composed of a short-range repulsive soft-core potential that directly adopts the form of an explicit-solvent DPD model and a long-range attractive hydrophobic potential. The parameters of the attractive potentials for lipid molecules were obtained by fitting the explicit-solvent DPD simulation of one bead of any type in a water box, then finely tuning it until the bilayer membrane properties obtained in the explicit-solvent model were matched. These parameters were further extended to amino acids according to bead type. The structural and elastic properties of bilayer membranes, free energy profiles for a lipid flip-flop and amino acid analogues translocating across the membrane, and membrane pore formation induced by antimicrobial peptides obtained from this solvent-free DPD force field considerably agreed with the explicit-solvent DPD results. Importantly, the efficiency of this method is guaranteed to accelerate the assembly of vesicles composed of several thousand lipids by up to 50-fold, rendering the experimental liposome dynamics as well as membrane-peptide interactions feasible at accessible computational expense.

Incorporation of DDR2 clusters into collagen matrix via integrin-dependent posterior remnant tethering

Cell-matrix interactions play critical roles in cell adhesion, tissue remodeling and cancer metastasis. Discoidin domain receptor 2 (DDR2) is a collagen receptor belonging to receptor tyrosine kinase (RTK) family. It is a powerful regulator of collagen deposition in the extracellular matrix (ECM). Although the oligomerization of DDR extracellular domain (ECD) proteins can affect matrix remodeling by inhibiting fibrillogenesis, it is still unknown how cellular DDR2 is incorporated into collagen matrix. Using 3-dimentional (3D) imaging for migrating cells, we identified a novel mechanism that explains how DDR2 incorporating into collagen matrix, which we named as posterior remnant tethering. We followed the de novo formation of these remnants and identified that DDR2 clusters formed at the retracting phase of a pseudopodium, then these clusters were tethered to fibrillar collagen and peeled off from the cell body to generate DDR2 containing posterior remnants. Inhibition of β1-integrin or Rac1 activity abrogated the remnant formation. Thus, our findings unveil a special cellular mechanism for DDR2 clusters incorporating into collagen matrix in an integrin-dependent manner.

Acyl chain asymmetry and polyunsaturation of brain phospholipids facilitate membrane vesiculation without leakage

Phospholipid membranes form cellular barriers but need to be flexible enough to divide by fission. Phospholipids generally contain a saturated fatty acid (FA) at position sn1 whereas the sn2-FA is saturated, monounsaturated or polyunsaturated. Our understanding of the impact of phospholipid unsaturation on membrane flexibility and fission is fragmentary. Here, we provide a comprehensive view of the effects of the FA profile of phospholipids on membrane vesiculation by dynamin and endophilin. Coupled to simulations, this analysis indicates that: (i) phospholipids with two polyunsaturated FAs make membranes prone to vesiculation but highly permeable; (ii) asymmetric sn1-saturated-sn2-polyunsaturated phospholipids provide a tradeoff between efficient membrane vesiculation and low membrane permeability; (iii) When incorporated into phospholipids, docosahexaenoic acid (DHA; omega-3) makes membranes more deformable than arachidonic acid (omega-6). These results suggest an explanation for the abundance of sn1-saturated-sn2-DHA phospholipids in synaptic membranes and for the importance of the omega-6/omega-3 ratio on neuronal functions.

Surrounding each living cell is a membrane that is mainly made of fat molecules called phospholipids. Similar membranes also surround many of the structures inside cells. It is important for life that these membranes are impermeable to many molecules; for example, they do not allow ions to cross them freely. The membranes also need to be flexible and allow cells to form different shapes. Flexible membranes also allow cells to move molecules around and to divide to produce new cells.

Each phospholipid includes two long chains of atoms called fatty acids. There are many fatty acids but they are typically grouped into saturated and unsaturated based on their chemical structures. The omega-3 and omega-6 fats are both groups of unsaturated fatty acids that are found in brain cells. Many phospholipids in cell membranes contain one saturated and one unsaturated fatty acid but it is not clear why.

By studying fat molecules in the laboratory and combining this with simulations, Manni et al. have now examined the effects of fatty acids on membranes. The investigation showed that phospholipids with both saturated and unsaturated fatty acids strike a balance between impermeable and flexible membranes. More unsaturated fatty acids make more flexible membranes but they are too permeable to be used in cells. The experiments also revealed that omega-3 unsaturated fats aid flexibility more than omega-6. This finding may help to explain why the relative amounts of omega-3 and -6 are so important in the membranes of brain cells.

The connection between the fats we eat and the fatty acids in our cells is complex. Yet, findings like these serve to remind us that we need a balanced diet of different fats to keep all our cells healthy.

Size-dependent formation of membrane nanotubes: continuum modeling and molecular dynamics simulations

With continuum theory and molecular dynamics simulations we demonstrated that the lipid membrane upon extraction exhibits size- and tension-dependent mechanical behaviors, and different structural lipid rearrangements in different leaflets.

High aspect ratio silicon nanowires control fibroblast adhesion and cytoskeleton organization

Cell-cell and cell-matrix interactions are essential to the survival and proliferation of most cells, and are responsible for triggering a wide range of biochemical pathways. More recently, the biomechanical role of those interactions was highlighted, showing, for instance, that adhesion forces are essential for cytoskeleton organization. Silicon nanowires (Si NWs) with their small size, high aspect ratio and anisotropic mechanical response represent a useful model to investigate the forces involved in the adhesion processes and their role in cellular development. In this work we explored and quantified, by single cell force spectroscopy (SCFS), the interaction of mouse embryonic fibroblasts with a flexible forest of Si NWs. We observed that the cell adhesion forces are comparable to those found on collagen and bare glass coverslip, analogously the membrane tether extraction forces are similar to that on collagen but stronger than that on bare flat glass. Cell survival did not depend significantly on the substrate, although a reduced proliferation after 36 h was observed. On the contrary both cell morphology and cytoskeleton organization revealed striking differences. The cell morphology on Si-NW was characterized by a large number of filopodia and a significant decrease of the cell mobility. The cytoskeleton organization was characterized by the absence of actin fibers, which were instead dominant on collagen and flat glass support. Such findings suggest that the mechanical properties of disordered Si NWs, and in particular their strong asymmetry, play a major role in the adhesion, morphology and cytoskeleton organization processes. Indeed, while adhesion measurements by SCFS provide out-of-plane forces values consistent with those measured on conventional substrates, weaker in-plane forces hinder proper cytoskeleton organization and migration processes.

Exploiting Lipid Permutation Symmetry to Compute Membrane Remodeling Free Energies

A complete physical description of membrane remodeling processes, such as fusion or fission, requires knowledge of the underlying free energy landscapes, particularly in barrier regions involving collective shape changes, topological transitions, and high curvature, where Canham-Helfrich (CH) continuum descriptions may fail. To calculate these free energies using atomistic simulations, one must address not only the sampling problem due to high free energy barriers, but also an orthogonal sampling problem of combinatorial complexity stemming from the permutation symmetry of identical lipids. Here, we solve the combinatorial problem with a permutation reduction scheme to map a structural ensemble into a compact, nondegenerate subregion of configuration space, thereby permitting straightforward free energy calculations via umbrella sampling. We applied this approach, using a coarse-grained lipid model, to test the CH description of bending and found sharp increases in the bending modulus for curvature radii below 10 nm. These deviations suggest that an anharmonic bending term may be required for CH models to give quantitative energetics of highly curved states.

Combined optical micromanipulation and interferometric topography (COMMIT)

Optical tweezers have emerged as a prominent light-based tool for pico-Newton (pN) force microscopy in mechanobiological studies. However, the efficacy of optical tweezers are limited in applications where concurrent metrology of the nano-sized structures under interrogation is essential to the quantitative analysis of its mechanical properties and various mechanotransduction events. We have developed an all-optical platform delivering pN force resolution in parallel with nano-scale structural imaging of the biological sample by combining optical tweezers with interferometric quantitative phase microscopy. These capabilities allow real-time micromanipulation and label-free measurement of sample's nanostructures and nanomechanical responses, opening avenues to a wide range of new research possibilities and applications in biology.

Binding of Lipopolysaccharide and Cholesterol-Modified Gelatin on Supported Lipid Bilayers: Effect of Bilayer Area Confinement and Bilayer Edge Tension

Binding of amphiphilic molecules to supported lipid bilayers (SLBs) often results in lipid fibril extension from the SLBs. Previous studies proposed that amphiphiles with large and flexible hydrophilic regions trigger lipid fibril formation in SLBs by inducing membrane curvature via their hydrophilic regions. However, no experimental studies have verified this mechanism of fibril formation. In this work, we investigated the binding of lipopolysaccharide (LPS) and cholesterol-modified gelatin to SLBs using fluorescence microscopy. SLBs with restricted and unrestricted bilayer areas were employed to identify the mechanism of fibril generation. We show that the main cause of lipid fibril formation is an approximately 20% expansion in the bilayer area rather than increased membrane curvature. The data indicate that bilayer area confinement plays a critical role in morphological changes of SLBs even when bound amphiphilic molecules have a large hydrophilic domain. We also show that bilayer area change after LPS insertion is dependent on the patch shape of the SLB. When an SLB patch consists of a broad bilayer segment connected to a long thin streak, bilayer area expansion mainly occurs within the bilayer streak. The results indicate that LPS insertion causes net lipid flow from the broad bilayer region to the streak area. The differential increase in area is explained by the instability of planar bilayer streaks that originate from the large energetic contribution of line tension arising along the bilayer edge.

Inter-tube adhesion mediates a new pearling mechanism

A common mechanism for intracellular transport is the controlled shape transformation, also known as pearling, of membrane tubes. Exploring how tube pearling takes place is thus of quite importance to not only understand the bio-functions of tubes, but also promote their potential biomedical applications. While the pearling mechanism of one single tube is well understood, both the pathway and the mechanism of pearling of multiple tubes still remain unclear. Herein, by means of computer simulations we show that the tube pearling can be mediated by the inter-tube adhesion. By increasing the inter-tube adhesion strength, each tube undergoes a discontinuous transition from no pearling to thorough pearling. The discontinuous pearling transition is ascribed to the competitive variation between tube surface tension and the extent of inter-tube adhesion. Besides, the final pearling instability is also affected by tube diameter and inter-tube orientation. Thinner tubes undergo inter-tube lipid diffusion before completion of pearling. The early lipid diffusion reduces the extent of inter-tube adhesion and thus restrains the subsequent pearling. Therefore, only partial or no pearling can take place for two thinner tubes. For two perpendicular tubes, the pearling is also observed, but with different pathways and higher efficiency. The finite size effect is discussed by comparing the pearling of tubes with different lengths. It is expected that this work will not only provide new insights into the mechanism of membrane tube pearling, but also shed light on the potential applications in biomaterials science and nanomedicine.

Lipid extraction mediates aggregation of carbon nanospheres in pulmonary surfactant monolayers

Our MD simulations demonstrate that the aggregation of carbon nanospheres in PSM is in fact size-dependent and mediated by lipid extractions.

Quantifying Molecular-Level Cell Adhesion on Electroactive Conducting Polymers using Electrochemical-Single Cell Force Spectroscopy

Single Cell Force Spectroscopy was combined with Electrochemical-AFM to quantify the adhesion between live single cells and conducting polymers whilst simultaneously applying a voltage to electrically switch the polymer from oxidized to reduced states. The cell-conducting polymer adhesion represents the non-specific interaction between cell surface glycocalyx molecules and polymer groups such as sulfonate and dodecylbenzene groups, which rearrange their orientation during electrical switching. Single cell adhesion significantly increases as the polymer is switched from an oxidized to fully reduced state, indicating stronger cell binding to sulfonate groups as opposed to hydrophobic groups. This increase in single cell adhesion is concomitant with an increase in surface hydrophilicity and uptake of cell media, driven by cation movement, into the polymer film during electrochemical reduction. Binding forces between the glycocalyx and polymer surface are indicative of molecular-level interactions and during electrical stimulation there is a decrease in both the binding force and stiffness of the adhesive bonds. The study provides insight into the effects of electrochemical switching on cell adhesion at the cell-conducting polymer interface and is more broadly applicable to elucidating the binding of cell adhesion molecules in the presence of electrical fields and directly at electrode interfaces.

Alchembed: A Computational Method for Incorporating Multiple Proteins into Complex Lipid Geometries

A necessary step prior to starting any membrane protein computer simulation is the creation of a well-packed configuration of protein(s) and lipids. Here, we demonstrate a method, alchembed, that can simultaneously and rapidly embed multiple proteins into arrangements of lipids described using either atomistic or coarse-grained force fields. During a short simulation, the interactions between the protein(s) and lipids are gradually switched on using a soft-core van der Waals potential. We validate the method on a range of membrane proteins and determine the optimal soft-core parameters required to insert membrane proteins. Since all of the major biomolecular codes include soft-core van der Waals potentials, no additional code is required to apply this method. A tutorial is included in the Supporting Information.

Molecular modeling of interaction between lipid monolayer and graphene nanosheets: implications for pulmonary nanotoxicity and pulmonary drug delivery

Understanding how nanoparticles interact with the pulmonary surfactant monolayer (PSM) is of great importance for safe applications in biomedicine and for evaluation of both health and environment impacts.

Coarse-grained force fields for molecular simulations

Molecular dynamics (MD) simulations at the atomic scale are a powerful tool to study the structure and dynamics of model biological systems. However, because of their high computational cost, the time and length scales of atomistic simulations are limited. Biologically important processes, such as protein folding, ion channel gating, signal transduction, and membrane remodeling, are difficult to investigate using atomistic simulations. Coarse-graining reduces the computational cost of calculations by reducing the number of degrees of freedom in the model, allowing simulations of larger systems for longer times. In the first part of this chapter we review briefly some of the coarse-grained models available for proteins, focusing on the specific scope of each model. Then we describe in more detail the MARTINI coarse-grained force field, and we illustrate how to set up and run a simulation of a membrane protein using the Gromacs software package. We explain step-by-step the preparation of the protein and the membrane, the insertion of the protein in the membrane, the equilibration of the system, the simulation itself, and the analysis of the trajectory.

Computer investigations of influences of molar fraction and acyl chain length of lipids on the nanoparticle–biomembrane interactions

Structural variations of the heterogeneous membrane: (a) a water defect, (b) the membrane buckling.

Interleaflet sliding in lipidic bilayers under shear flow: comparison of the gel and fluid phases using reversed non-equilibrium molecular dynamics simulations

The friction between two rubbing surfaces lubricated by water can be diminished if they are coated with phospholipidic bilayers or brushes of polyelectrolytes. In the case of a coating by lipid membranes, the friction is lower when the lipids are in the gel phase rather than in the liquid phase. We investigated the response of fluid or gel bilayers to a mechanical load or under shear using non-equilibrium molecular dynamics simulations (NEMD) to understand whether this difference could come from intermonolayer sliding. The system is composed of a single fully hydrated bilayer of coarse grained phospholipids under a parallel shear with vorticity parallel to the bilayer. In both the liquid and the gel phases, an intermonolayer slip was measured in the velocity profile. In the liquid phase this slip is proportional to the shear stress. In the tilted gel phase of our model the stress is not systematically linear and relaxes differently when the shear is in the direction of the tilt or perpendicular to it. The impact of surface tension (or load) on the friction is different for the liquid and gel phases, but grossly the slip remains of the same order of magnitude.

Lipid cell biology. Polyunsaturated phospholipids facilitate membrane deformation and fission by endocytic proteins

Phospholipids (PLs) with polyunsaturated acyl chains are extremely abundant in a few specialized cellular organelles such as synaptic vesicles and photoreceptor discs, but their effect on membrane properties is poorly understood. Here, we found that polyunsaturated PLs increased the ability of dynamin and endophilin to deform and vesiculate synthetic membranes. When cells incorporated polyunsaturated fatty acids into PLs, the plasma membrane became more amenable to deformation by a pulling force and the rate of endocytosis was accelerated, in particular, under conditions in which cholesterol was limiting. Molecular dynamics simulations and biochemical measurements indicated that polyunsaturated PLs adapted their conformation to membrane curvature. Thus, by reducing the energetic cost of membrane bending and fission, polyunsaturated PLs may help to support rapid endocytosis.

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Tubular membrane formation of binary giant unilamellar vesicles composed of cylinder and inverse-cone-shaped lipids

We have succeeded in controlling tubular membrane formations in binary giant unilamellar vesicles (GUVs) using a simple temperature changing between the homogeneous one-phase region and the two-phase coexistence region. The binary GUV is composed of inverse-cone (bulky hydrocarbon chains and a small headgroup) and cylinder-shaped lipids. When the temperature was set in the two-phase coexistence region, the binary GUV had a spherical shape with solidlike domains. By increasing the temperature to the homogeneous one-phase region, the excess area created by the chain melting of the lipid produced tubes inside the GUV. The tubes had a radius on the micrometer scale and were stable in the one-phase region. When we again decreased the temperature to the two-phase coexisting region, the tubes regressed and the GUVs recovered their phase-separated spherical shape. We infer that the tubular formation was based on the mechanical balance of the vesicle membrane (spontaneous tension) coupled with the asymmetric distribution of the inverse-cone-shaped lipids between the inner and outer leaflets of the vesicle (lipid sorting).

The power of coarse graining in biomolecular simulations

Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling.

Molecular modeling of membrane tube pearling and the effect of nanoparticle adsorption

DPD simulations suggest that the membrane tube pearling can be regulated by the inner water pressure and NP adsorption.

Computer simulation studies on passive recruitment dynamics of lipids induced by the adsorption of charged nanoparticles

The recruitment dynamics of lipids in the biomembrane is believed to play an important role in a variety of cellular processes. (a) flip-flops, (b) lateral diffusion, (c) nanoscale domain, (d) protrusions.

Gaussian curvature elasticity determined from global shape transformations and local stress distributions: a comparative study using the MARTINI model

We calculate the Gaussian curvature modulus kappa of a systematically coarse-grained (CG) one-component lipid membrane by applying the method recently proposed by Hu et al. [Biophys. J., 2012, 102, 1403] to the MARTINI representation of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). We find the value kappa/kappa = -1.04 +/- 0.03 for the elastic ratio between the Gaussian and the mean curvature modulus and deduce kappa(m)/kappa(m) = -0.98 +/- 0.09 for the monolayer elastic ratio, where the latter is based on plausible assumptions for the distance z0 of the monolayer neutral surface from the bilayer midplane and the spontaneous lipid curvature K(0m). By also analyzing the lateral stress profile sigma0(z) of our system, two other lipid types and pertinent data from the literature, we show that determining K(0m) and kappa through the first and second moment of sigma0(z) gives rise to physically implausible values for these observables. This discrepancy, which we previously observed for a much simpler CG model, suggests that the moment conditions derived from simple continuum assumptions miss the effect of physically important correlations in the lipid bilayer.

Constant pH simulations with the coarse-grained MARTINI model — Application to oleic acid aggregates

Long chain fatty acids are biologically important molecules with complex and pH sensitive aggregation behavior. The carboxylic head group of oleic acid is ionizable, with the pKa shifting to larger values, even above a value of 7, in certain aggregate states. While experiments have determined the macroscopic phase behavior, we have yet to understand the molecular level details for this complex behavior. This level of detail is likely required to fully appreciate the role of fatty acids in biology and for nanoscale biotechnological and industrial applications. Here, we introduce the use of constant pH molecular dynamics (MD) simulations with the coarse-grained MARTINI model and apply the method to oleic acid aggregates and a model lipid bilayer. By running simulations at different constant pH values, we determined titration curves and the resulting pKa for oleic acid in different environments. The coarse-grained model predicts positive pKa shifts, with a shift from 4.8 in water to 6.5 in a small micelle, and 6.6 in a dioleoylphosphatidylcholine (DOPC) bilayer, similar to experimental estimates. The size of the micelles increased as the pH increased, and correlated with the fraction of deprotonated oleic acid. We show this combination of constant pH MD and the coarse-grained MARTINI model can be used to model pH-dependent surfactant phase behavior. This suggests a large number of potential new applications of large-scale MARTINI simulations in other biological systems with ionizable molecules.

Mechanical model of vertical nanowire cell penetration

Direct access into cells' interiors is essential for biomolecular delivery, gene transfection, and electrical recordings yet is challenging due to the cell membrane barrier. Recently, molecular delivery using vertical nanowires (NWs) has been demonstrated for introducing biomolecules into a large number of cells in parallel. However, the microscopic understanding of how and when the nanowires penetrate cell membranes is still lacking, and the degree to which actual membrane penetration occurs is controversial. Here we present results from a mechanical continuum model of elastic cell membrane penetration through two mechanisms, namely through "impaling" as cells land onto a bed of nanowires, and through "adhesion-mediated" penetration, which occurs as cells spread on the substrate and generate adhesion force. Our results reveal that penetration is much more effective through the adhesion mechanism, with NW geometry and cell stiffness being critically important. Stiffer cells have higher penetration efficiency, but are more sensitive to NW geometry. These results provide a guide to designing nanowires for applications in cell membrane penetration.