Bilayer Asymmetry Influences Integrin Sequestering in Raft-Mimicking Lipid Mixtures

There is growing recognition that lipid heterogeneities in cellular membranes play an important role in the distribution and functionality of membrane proteins. However, the detection and characterization of such heterogeneities at the cellular level remains challenging. Here we report on the poorly understood relationship between lipid bilayer asymmetry and membrane protein sequestering in raft-mimicking model membrane mixtures using a powerful experimental platform comprised of confocal spectroscopy XY-scan and photon-counting histogram analyses. This experimental approach is utilized to probe the domain-specific sequestering and oligomerization state of α_vβ₃ and α₅β₁ integrins in bilayers, which contain coexisting liquid-disordered/liquid-ordered (l_d/l_o) phase regions exclusively in the top leaflet of the bilayer (bottom leaflet contains l_d phase). Comparison with previously reported integrin sequestering data in bilayer-spanning l_o-l_d phase separations demonstrates that bilayer asymmetry has a profound influence on α_vβ₃ and α₅β₁ sequestering behavior. For example, both integrins sequester preferentially to the l_o phase in asymmetric bilayers, but to the l_d phase in their symmetric counterparts. Furthermore, our data show that bilayer asymmetry significantly influences the role of native ligands in integrin sequestering.     

Phospholipid Chain Interactions with Cholesterol Drive Domain Formation in Lipid Membranes

Cholesterol is a key component of eukaryotic membranes, but its role in cellular biology in general and in lipid rafts in particular remains controversial. Model membranes are used extensively to determine the phase behavior of ternary mixtures of cholesterol, a saturated lipid, and an unsaturated lipid with liquid-ordered and liquid-disordered phase coexistence. Despite many different experiments that determine lipid-phase diagrams, we lack an understanding of the molecular-level driving forces for liquid phase coexistence in bilayers with cholesterol. Here, we use atomistic molecular dynamics computer simulations to address the driving forces for phase coexistence in ternary lipid mixtures. Domain formation is directly observed in a long-timescale simulation of a mixture of 1,2-distearoyl-sn-glycero-3-phosphocholine, unsaturated 1,2-dilinoleoyl-sn-glycero-3-phosphocholine, and cholesterol. Free-energy calculations for the exchange of the saturated and unsaturated lipids between the ordered and disordered phases give insight into the mixing behavior. We show that a large energetic contribution to domain formation is favorable enthalpic interactions of the saturated lipid in the ordered phase. This favorable energy for forming an ordered, cholesterol-rich phase is opposed by a large unfavorable entropy. Martini coarse-grained simulations capture the unfavorable free energy of mixing but do not reproduce the entropic contribution because of the reduced representation of the phospholipid tails. Phospholipid tails and their degree of unsaturation are key energetic contributors to lipid phase separation.     

Drunken Membranes: Short-Chain Alcohols Alter Fusion of Liposomes to Planar Lipid Bilayers

Although the effects of ethanol on protein receptors and lipid membranes have been studied extensively, ethanol’s effect on vesicles fusing to lipid bilayers is not known. To determine the effect of alcohols on fusion rates, we utilized the nystatin/ergosterol fusion assay to measure fusion of liposomes to a planar lipid bilayer (BLM). The addition of ethanol excited fusion when applied on the cis (vesicle) side, and inhibited fusion on the trans side. Other short-chain alcohols followed a similar pattern. In general, the inhibitory effect of alcohols (trans) occurs at lower doses than the excitatory (cis) effect, with a decrease of 29% in fusion rates at the legal driving limit of 0.08% (w/v) ethanol (IC₅₀ = 0.2% v/v, 34 mM). Similar inhibitory effects were observed with methanol, propanol, and butanol, with ethanol being the most potent. Significant variability was observed with different alcohols when applied to the cis side. Ethanol and propanol enhanced fusion, butanol also enhanced fusion but was less potent, and low doses of methanol mildly inhibited fusion. The inhibition by trans addition of alcohols implies that they alter the planar membrane structure and thereby increase the activation energy required for fusion, likely through an increase in membrane fluidity. The cis data are likely a combination of the above effect and a proportionally greater lowering of the vesicle lysis tension and hydration repulsive pressure that combine to enhance fusion. Alternate hypotheses are also discussed. The inhibitory effect of ethanol on liposome-membrane fusion is large enough to provide a possible biophysical explanation of compromised neuronal behavior.     

Ionomycin-Induced Changes in Membrane Potential Alter Electroporation Outcomes in HL-60 Cells

Previous studies have shown greater fluorophore uptake during electroporation on the anode-facing side of the cell than on the cathode-facing side. Based on these observations, we hypothesized that hyperpolarizing a cell before electroporation would decrease the requisite pulsed electric field intensity for electroporation outcomes, thereby yielding a higher probability of reversible electroporation at lower electric field strengths and a higher probability of irreversible electroporation (IRE) at higher electric field strengths. In this study, we tested this hypothesis by hyperpolarizing HL-60 cells using ionomycin before electroporation. These cells were then electroporated in a solution containing propidium iodide, a membrane integrity indicator. After 20 min, we added trypan blue to identify IRE cells. Our results showed that hyperpolarizing cells before electroporation alters the pulsed electric field intensity thresholds for reversible electroporation and IRE, allowing for greater control and selectivity of electroporation outcomes.     

Probing Conformational Changes during the Gating Cycle of a Potassium Channel in Lipid Bilayers

Ion conduction across the cellular membrane requires the simultaneous opening of activation and inactivation gates of the K⁺ channel pore. The bacterial KcsA channel has served as a powerful system for dissecting the structural changes that are related to four major functional states associated with K⁺ gating. Yet, the direct observation of the full gating cycle of KcsA has remained structurally elusive, and crystal structures mimicking these gating events require mutations in or stabilization of functionally relevant channel segments. Here, we found that changes in lipid composition strongly increased the KcsA open probability. This enabled us to probe all four major gating states in native-like membranes by combining electrophysiological and solid-state NMR experiments. In contrast to previous crystallographic views, we found that the selectivity filter and turret region, coupled to the surrounding bilayer, were actively involved in channel gating. The increase in overall steady-state open probability was accompanied by a reduction in activation-gate opening, underscoring the important role of the surrounding lipid bilayer in the delicate conformational coupling of the inactivation and activation gates.     

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.     

Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions

Rhodopsin, a prototypical G protein-coupled receptor, is a membrane protein that can sense dim light. This highly effective photoreceptor is known to be sensitive to the composition of its lipidic environment, but the molecular mechanisms underlying this fine-tuned modulation of the receptor’s function and structural stability are not fully understood. There are two competing hypotheses to explain how this occurs: 1) lipid modulation occurs via solvent-like interactions, where lipid composition controls membrane properties like hydrophobic thickness, which in turn modulate the protein’s conformational equilibrium; or 2) protein-lipid interactions are ligand-like, with specific hot spots and long-lived binding events. By analyzing an ensemble of all-atom molecular dynamics simulations of five different states of rhodopsin, we show that a local ordering effect takes place in the membrane upon receptor activation. Likewise, docosahexaenoic acid acyl tails and phosphatidylethanolamine headgroups behave like weak ligands, preferentially binding to the receptor in inactive-like conformations and inducing subtle but significant structural changes.     

Ligand Binding Alters Dimerization and Sequestering of Urokinase Receptors in Raft-Mimicking Lipid Mixtures

Lipid heterogeneities, such as lipid rafts, are widely considered to be important for the sequestering of membrane proteins in plasma membranes, thereby influencing membrane protein functionality. However, the underlying mechanisms of such sequestration processes remain elusive, in part, due to the small size and often transient nature of these functional membrane heterogeneities in cellular membranes. To overcome these challenges, here we report the sequestration behavior of urokinase receptor (uPAR), a glycosylphosphatidylinositol-anchored protein, in a planar model membrane platform with raft-mimicking lipid mixtures of well-defined compositions using a powerful optical imaging platform consisting of confocal spectroscopy XY-scans, photon counting histogram, and fluorescence correlation spectroscopy analyses. This methodology provides parallel information about receptor sequestration, oligomerization state, and lateral mobility with single molecule sensitivity. Most notably, our experiments demonstrate that moderate changes in uPAR sequestration are not only associated with modifications in uPAR dimerization levels, but may also be linked to ligand-mediated allosteric changes of these membrane receptors. Our data show that these modifications in uPAR sequestration can be induced by exposure to specific ligands (urokinase plasminogen activator, vitronectin), but not via adjustment of the cholesterol level in the planar model membrane system. Good agreement of our key findings with published results on cell membranes confirms the validity of our model membrane approach. We hypothesize that the observed mechanism of receptor translocation in the presence of raft-mimicking lipid mixtures is also applicable to other glycosylphosphatidylinositol-anchored proteins.     

Gramicidin Increases Lipid Flip-Flop in Symmetric and Asymmetric Lipid Vesicles

Unlike most transmembrane proteins, phospholipids can migrate from one leaflet of the membrane to the other. Because this spontaneous lipid translocation (flip-flop) tends to be very slow, cells facilitate the process with enzymes that catalyze the transmembrane movement and thereby regulate the transbilayer lipid distribution. Nonenzymatic membrane-spanning proteins with unrelated primary functions have also been found to accelerate lipid flip-flop in a nonspecific manner and by various hypothesized mechanisms. Using deuterated phospholipids, we examined the acceleration of flip-flop by gramicidin channels, which have well-defined structures and known functions, features that make them ideal candidates for probing the protein-membrane interactions underlying lipid flip-flop. To study compositionally and isotopically asymmetric proteoliposomes containing gramicidin, we expanded a recently developed protocol for the preparation and characterization of lipid-only asymmetric vesicles. Channel incorporation, conformation, and function were examined with small angle x-ray scattering, circular dichroism, and a stopped-flow spectrofluorometric assay, respectively. As a measure of lipid scrambling, we used differential scanning calorimetry to monitor the effect of gramicidin on the melting transition temperatures of the two bilayer leaflets. The two calorimetric peaks of the individual leaflets merged into a single peak over time, suggestive of scrambling, and the effect of the channel on the transbilayer lipid distribution in both symmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and asymmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles was quantified from proton NMR measurements. Our results show that gramicidin increases lipid flip-flop in a complex, concentration-dependent manner. To determine the molecular mechanism of the process, we used molecular dynamics simulations and further computational analysis of the trajectories to estimate the extent of membrane deformation. Together, the experimental and computational approaches were found to constitute an effective means for studying the effects of transmembrane proteins on lipid distribution in both symmetric and asymmetric model membranes.     

Ligand Binding Alters Dimerization and Sequestering of Urokinase Receptors in Raft-Mimicking Lipid Mixtures

Lipid heterogeneities, such as lipid rafts, are widely considered to be important for the sequestering of membrane proteins in plasma membranes, thereby influencing membrane protein functionality. However, the underlying mechanisms of such sequestration processes remain elusive, in part, due to the small size and often transient nature of these functional membrane heterogeneities in cellular membranes. To overcome these challenges, here we report the sequestration behavior of urokinase receptor (uPAR), a glycosylphosphatidylinositol-anchored protein, in a planar model membrane platform with raft-mimicking lipid mixtures of well-defined compositions using a powerful optical imaging platform consisting of confocal spectroscopy XY-scans, photon counting histogram, and fluorescence correlation spectroscopy analyses. This methodology provides parallel information about receptor sequestration, oligomerization state, and lateral mobility with single molecule sensitivity. Most notably, our experiments demonstrate that moderate changes in uPAR sequestration are not only associated with modifications in uPAR dimerization levels, but may also be linked to ligand-mediated allosteric changes of these membrane receptors. Our data show that these modifications in uPAR sequestration can be induced by exposure to specific ligands (urokinase plasminogen activator, vitronectin), but not via adjustment of the cholesterol level in the planar model membrane system. Good agreement of our key findings with published results on cell membranes confirms the validity of our model membrane approach. We hypothesize that the observed mechanism of receptor translocation in the presence of raft-mimicking lipid mixtures is also applicable to other glycosylphosphatidylinositol-anchored proteins.     

Controlling Anomalous Diffusion in Lipid Membranes

Diffusion in cell membranes is not just simple two-dimensional Brownian motion but typically depends on the timescale of the observation. The physical origins of this anomalous subdiffusion are unresolved, and model systems capable of quantitative and reproducible control of membrane diffusion have been recognized as a key experimental bottleneck. Here, we control anomalous diffusion using supported lipid bilayers containing lipids derivatized with polyethylene glycol (PEG) headgroups. Bilayers with specific excluded area fractions are formed by control of PEG lipid mole fraction. These bilayers exhibit a switch in diffusive behavior, becoming anomalous as bilayer continuity is disrupted. Using a combination of single-molecule fluorescence and interferometric imaging, we measure the anomalous behavior in this model over four orders of magnitude in time. Diffusion in these bilayers is well described by a power-law dependence of the mean-square displacement with observation time. Anomaleity in this system can be tailored by simply controlling the mole fraction of PEG lipid, producing bilayers with diffusion parameters similar to those observed for anomalous diffusion in biological membranes.     

Coexistence and Pattern Formation in Bacterial Mixtures with Contact-Dependent Killing

Multistrain microbial communities often exhibit complex spatial organization that emerges because of the interplay of various cooperative and competitive interaction mechanisms. One strong competitive mechanism is contact-dependent neighbor killing enabled by the type VI secretion system. It has been previously shown that contact-dependent killing can result in bistability of bacterial mixtures so that only one strain survives and displaces the other. However, it remains unclear whether stable coexistence is possible in such mixtures. Using a population dynamics model for two interacting bacterial strains, we found that coexistence can be made possible by the interplay of contact-dependent killing and long-range growth inhibition, leading to the formation of various cellular patterns. These patterns emerge in a much broader parameter range than that required for the linear Turing-like instability, suggesting this may be a robust mechanism for pattern formation.     

Intrinsic Curvature-Mediated Transbilayer Coupling in Asymmetric Lipid Vesicles

We measured the effect of intrinsic lipid curvature, $J_0$, on structural properties of asymmetric vesicles made of palmitoyl-oleoyl-phosphatidylethanolamine (POPE; $J_0<0$) and palmitoyl-oleoyl-phosphatidylcholine (POPC; $J_0～0$). Electron microscopy and dynamic light scattering were used to determine vesicle size and morphology, and x-ray and neutron scattering, combined with calorimetric experiments and solution NMR, yielded insights into leaflet-specific lipid packing and melting processes. Below the lipid melting temperature we observed strong interleaflet coupling in asymmetric vesicles with POPE inner bilayer leaflets and outer leaflets enriched in POPC. This lipid arrangement manifested itself by lipids melting cooperatively in both leaflets, and a rearrangement of lipid packing in both monolayers. On the other hand, no coupling was observed in vesicles with POPC inner bilayer leaflets and outer leaflets enriched in POPE. In this case, the leaflets melted independently and did not affect each other’s acyl chain packing. Furthermore, we found no evidence for transbilayer structural coupling above the melting temperature of either sample preparation. Our results are consistent with the energetically preferred location of POPE residing in the inner leaflet, where it also resides in natural membranes, most likely causing the coupling of both leaflets. The loss of this coupling in the fluid bilayers is most likely the result of entropic contributions.     

A Cholesterol Analog Induces an Oligomeric Reorganization of VDAC

The oligomeric organization of the voltage-dependent anion-selective channel (VDAC) and its interactions with hexokinase play integral roles in mitochondrially mediated apoptotic signaling. Various small to large assemblies of VDAC are observed in mitochondrial outer membranes, but they do not predominate in detergent-solubilized VDAC samples. In this study, a cholesterol analog, cholesteryl-hemisuccinate (CHS), was shown to induce the formation of detergent-soluble VDAC multimers. The various oligomeric states of VDAC induced by the addition of CHS were deciphered through an integrated biophysics approach using microscale thermophoresis, analytical ultracentrifugation, and size-exclusion chromatography small angle x-ray scattering. Furthermore, CHS stabilizes the interaction between VDAC and hexokinase (K_d of 27 ± 6 μM), confirming the biological relevance of oligomers generated. Thus, sterols such as cholesterol in higher eukaryotes or ergosterol in fungi may regulate the VDAC oligomeric state and may provide a potential target for the modulation of apoptotic signaling by effecting VDAC-VDAC and VDAC-hexokinase interactions. In addition, the integrated biophysical approach described provides a powerful platform for the study of membrane protein complexes in solution.     

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na⁺ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the ³¹P chemical shift anisotropy and ¹⁴N and CD₃ NMR quadrupole splittings.     

An Asymmetry in Monolayer Tension Regulates Lipid Droplet Budding Direction

Cells store excess energy in the form of neutral lipids that are synthesized and encapsulated within the endoplasmic reticulum intermonolayer space. The lipids next demix to form lipid droplets (LDs), which, surprisingly, bud off mostly toward the cytosol. This directional LD formation is critical to energy metabolism, but its mechanism remains poorly understood. Here, we reconstituted the LD formation topology by embedding artificial LDs into the intermonolayer space of bilayer vesicles. We provide experimental evidence that the droplet behavior in the membrane is recapitulated by the physics of three-phase wetting systems, dictated by the equilibrium of surface tensions. We thereupon determined that slight tension asymmetries between the membrane monolayers regulate the droplet budding side. A differential regulation of lipid or protein composition around a forming LD can generate a monolayer tension asymmetry that will determine the LD budding side. Our results offer, to our knowledge, new insights on how the proteins might regulate LD formation side by generating a monolayer tension asymmetry.