Protein diffusion on charged membranes: a dynamic mean-field model describes time evolution and lipid reorganization

As charged macromolecules adsorb and diffuse on cell membranes in a large variety of cell signaling processes, they can attract or repel oppositely charged lipids. This results in lateral membrane rearrangement and affects the dynamics of protein function. To address such processes quantitatively we introduce a dynamic mean-field scheme that allows self-consistent calculations of the equilibrium state of membrane-protein complexes after such lateral reorganization of the membrane components, and serves to probe kinetic details of the process. Applicable to membranes with heterogeneous compositions containing several types of lipids, this comprehensive method accounts for mobile salt ions and charged macromolecules in three dimensions, as well as for lateral demixing of charged and net-neutral lipids in the membrane plane. In our model, the mobility of membrane components is governed by the diffusion-like Cahn-Hilliard equation, while the local electrochemical potential is based on nonlinear Poisson-Boltzmann theory. We illustrate the method by applying it to the adsorption of the anionic polypeptide poly-Lysine on negatively charged lipid membranes composed of binary mixtures of neutral and monovalent lipids, or onto ternary mixtures of neutral, monovalent, and multivalent lipids. Consistent with previous calculations and experiments, our results show that at steady-state multivalent lipids (such as PIP₂), but not monovalent lipid (such as phosphatidylserine), will segregate near the adsorbing macromolecules. To address the corresponding diffusion of the adsorbing protein in the membrane plane, we couple lipid mobility with the propagation of the adsorbing protein through a dynamic Monte Carlo scheme. We find that due to their higher mobility dictated by the electrochemical potential, multivalent lipids such as PIP₂ more quickly segregate near oppositely charged proteins than do monovalent lipids, even though their diffusion constants may be similar. The segregation, in turn, slows protein diffusion, as lipids introduce an effective drag on the motion of the adsorbate. In contrast, monovalent lipids such as phosphatidylserine only weakly segregate, and the diffusions of protein and lipid remain largely uncorrelated.     

Two distinct anionic phospholipid-dependent events involved in SecA-mediated protein translocation

Protein translocation across the bacterial cytoplasmic membrane is an essential process catalyzed by the Sec translocase, which in its minimal form consists of the protein-conducting channel SecYEG, and the motor ATPase SecA. SecA binds via its positively charged N-terminus to membranes containing anionic phospholipids, leading to a lipid-bound intermediate. This interaction induces a conformational change in SecA, resulting in a high-affinity association with SecYEG, which initiates protein translocation. Here, we examined the effect of anionic lipids on the SecA-SecYEG interaction in more detail, and discovered a second, yet unknown, anionic lipid-dependent event that stimulates protein translocation. Based on molecular dynamics simulations we identified an anionic lipid-enriched region in vicinity of the lateral gate of SecY. Here, the anionic lipid headgroup accesses the lateral gate, thereby stabilizing the pre-open state of the channel. The simulations suggest flip-flop movement of phospholipid along the lateral gate. Electrostatic contribution of the anionic phospholipids at the lateral gate may directly stabilize positively charged residues of the signal sequence of an incoming preprotein. Such a mechanism allows for the correct positioning of the entrant peptide, thereby providing a long-sought explanation for the role of anionic lipids in signal sequence folding during protein translocation.     

Effect of polyelectrolyte adsorption on lateral distribution and dynamics of anionic lipids: a Monte Carlo study of a coarse-grain model

We employ Monte Carlo simulations to investigate the interaction between an adsorbing linear flexible cationic polyelectrolyte and a ternary mixed fluid membrane containing neutral (phosphatidylcholine, PC), monovalent (phosphatidylserine, PS), and multivalent (phosphatidylinositol, PIP₂) anionic lipids. We systematically explore the influences of polyelectrolyte chain length, polyelectrolyte charge density, polyelectrolyte total charge amount, and salt solution ionic strength on the static and dynamic properties of different anionic lipid species. Our results show that the multivalent PIP₂ lipids dominate the polyelectrolyte–membrane interaction and competitively inhibit polyelectrolyte–PS binding. When the total charge amount of the polyelectrolyte is less than that of the local oppositely charged PIP₂ lipids, the polyelectrolyte can drag the bound multivalent lipids to diffuse on the membrane, but cannot interact with the PS lipids. Under this condition, the diffusion behaviors of the polyelectrolyte closely follow the prediction of the Rouse model, and the polyelectrolyte chain properties determine the adsorption amount, concentration gradients, and hierarchical mobility of the bound PIP₂ lipids. However, when the total charge amount of the polyelectrolyte is larger than that of the local PIP₂ lipids, the polyelectrolyte further binds the PS lipids around the polyelectrolyte–PIP₂ complex to achieve local electrical neutrality. In this condition, parts of the polyelectrolyte desorb from the membrane and show faster mobility, and the bound PS presents much faster mobility than the segregated PIP₂. This work provides an explanation for heterogeneity formation in different anionic lipids induced by polyelectrolyte adsorption.     

Anionic phospholipid gradients: an uncharacterized frontier of the plant endomembrane network

Anionic phospholipids include phosphatidic acid (PA), phosphatidylserine (PS), phosphatidylinositol (PI), and its phosphorylated derivatives the phosphoinositides (e.g. phosphatidylinositol-4-phosphate [PI4P] and phosphatidylinositol-4,5-bisphosphate [PI(4,5)P₂]). Although anionic phospholipids are low-abundant lipids, they are particularly important for membrane functions. In particular, anionic lipids act as biochemical and biophysical landmarks that contribute to the establishment of membrane identity, signaling activities, and compartment morphodynamics. Each anionic lipid accumulates in different endomembranes according to a unique subcellular pattern, where they locally provide docking platforms for proteins. As such, they are mostly believed to act in the compartments in which they accumulate. However, mounting evidence throughout eukaryotes suggests that anionic lipids are not as compartment-specific as initially thought and that they are instead organized as concentration gradients across different organelles. In this update, we review the evidence for the existence of anionic lipid gradients in plants. We then discuss the possible implication of these gradients in lipid dynamics and homeostasis, and also in coordinating subcellular activities. Finally, we introduce the notion that anionic lipid gradients at the cellular scale may translate into gradients at the tissue level, which could have implications for plant development.

One-sentence  summaryMounting evidence reveals that anionic lipids exist in gradients within the plant endomembrane network that may have implications on lipid homeostasis and membrane dynamics.     

Structural investigation on the adsorption of the MARCKS peptide on anionic lipid monolayers - effects beyond electrostatic

The presence of charged lipids in the cell membrane constitutes the background for the interaction with numerous membrane proteins. As a result, the valence of the lipids plays an important role concerning their lateral organization in the membrane and therefore the very manner of this interaction. This present study examines this aspect, particularly regarding to the interaction of the anionic lipid DPPS with the highly basic charged effector domain of the MARCKS protein, examined in monolayer model systems. Film balance, fluorescence microscopy and X-ray reflection/diffraction measurements were used to study the behavior of DPPS in a mixture with DPPC for its dependance on the presence of MARCKS (151-175). In the mixed monolayer, both lipids are completely miscible therefore DPPS is incorporated in the ordered crystalline DPPC domains as well. The interaction of MARCKS peptide with the mixed monolayer leads to the formation of lipid/peptide clusters causing an elongation of the serine group of the DPPS up to 7? in direction to surface normal into the subphase. The large cationic charge of the peptide pulls out the serine group of the interface which simultaneously causes an elongation of the phosphodiester group of the lipid fraction too. The obtained results were used to compare the interaction of MARCKS peptide with the polyvalent PIP(2) in mixed monolayers. On this way we surprisingly find out, that the relative small charge difference of the anionic lipids causes a significant different interaction with MARCKS (151-175). The lateral arrangement of the anionic lipids depends on their charge values and determines the diffusion of the electrostatic binding clusters within the membrane.     

Role of GAP-43 in sequestering phosphatidylinositol 4,5-bisphosphate to Raft bilayers

The lipid phosphatidylinositol 4,5-bisphosphate (PIP₂) is critical for a number of physiological functions, and its presence in membrane microdomains (rafts) appears to be important for several of these spatially localized events. However, lipids like PIP₂ that contain polyunsaturated hydrocarbon chains are usually excluded from rafts, which are enriched in phospholipids (such as sphingomyelin) containing saturated or monounsaturated chains. Here we tested a mechanism by which multivalent PIP₂ molecules could be transferred into rafts through electrostatic interactions with polybasic cytoplasmic proteins, such as GAP-43, which bind to rafts via their acylated N-termini. We analyzed the interactions between lipid membranes containing raft microdomains and a peptide (GAP-43P) containing the linked N-terminus and the basic effector domain of GAP-43. In the absence or presence of nonacylated GAP-43P, PIP₂ was found primarily in detergent-soluble membranes thought to correspond to nonraft microdomains. However, when GAP-43P was acylated by palmitoyl coenzyme A, both the peptide and PIP₂ were greatly enriched in detergent-resistant membranes that correspond to rafts; acylation of GAP-43P changed the free energy of transfer of PIP₂ from detergent-soluble membranes to detergent-resistant membranes by −1.3 kcal/mol. Confocal microscopy of intact giant unilamellar vesicles verified that in the absence of GAP-43P PIP₂ was in nonraft microdomains, whereas acylated GAP-43P laterally sequestered PIP₂ into rafts. These data indicate that sequestration of PIP₂ to raft microdomains could involve interactions with acylated basic proteins such as GAP-43.     

Minimal Effect of Lipid Charge on Membrane Miscibility Phase Behavior in Three Ternary Systems

Giant unilamellar vesicles composed of a ternary mixture of phospholipids and cholesterol exhibit coexisting liquid phases over a range of temperatures and compositions. A significant fraction of lipids in biological membranes are charged. Here, we present phase diagrams of vesicles composed of phosphatidylcholine (PC) lipids, which are zwitterionic; phosphatidylglycerol (PG) lipids, which are anionic; and cholesterol (Chol). Specifically, we use DiPhyPG-DPPC-Chol and DiPhyPC-DPPG-Chol. We show that miscibility in membranes containing charged PG lipids occurs over similarly high temperatures and broad lipid compositions as in corresponding membranes containing only uncharged lipids, and that the presence of salt has a minimal effect. We verified our results in two ways. First, we used mass spectrometry to ensure that charged PC/PG/Chol vesicles formed by gentle hydration have the same composition as the lipid stocks from which they are made. Second, we repeated the experiments by substituting phosphatidylserine for PG as the charged lipid and observed similar phenomena. Our results consistently support the view that monovalent charged lipids have only a minimal effect on lipid miscibility phase behavior in our system.     

Analytic Approaches to Stochastic Gene Expression in Multicellular Systems

Giant unilamellar vesicles composed of a ternary mixture of phospholipids and cholesterol exhibit coexisting liquid phases over a range of temperatures and compositions. A significant fraction of lipids in biological membranes are charged. Here, we present phase diagrams of vesicles composed of phosphatidylcholine (PC) lipids, which are zwitterionic; phosphatidylglycerol (PG) lipids, which are anionic; and cholesterol (Chol). Specifically, we use DiPhyPG-DPPC-Chol and DiPhyPC-DPPG-Chol. We show that miscibility in membranes containing charged PG lipids occurs over similarly high temperatures and broad lipid compositions as in corresponding membranes containing only uncharged lipids, and that the presence of salt has a minimal effect. We verified our results in two ways. First, we used mass spectrometry to ensure that charged PC/PG/Chol vesicles formed by gentle hydration have the same composition as the lipid stocks from which they are made. Second, we repeated the experiments by substituting phosphatidylserine for PG as the charged lipid and observed similar phenomena. Our results consistently support the view that monovalent charged lipids have only a minimal effect on lipid miscibility phase behavior in our system.     

The effect of charged lipids on bacteriorhodopsin membrane reconstitution and its photochemical activities

Bacteriorhodopsin (BR) was reconstituted into artificial lipid membrane containing various charged lipid compositions. The proton pumping activity of BR under flash and continuous illumination, proton permeability across membrane, as well as the decay kinetics of the photocycle intermediate M₄₁₂ were studied. The results showed that lipid charges would significantly affect the orientation of BR inserted into lipid membranes. In liposomes containing anionic lipids, BRs were more likely to take natural orientation as in living cells. In neutral or positively charged liposomes, most BRs were reversely assembled, assuming an inside out orientation. Moreover, the lipids charges also affect BR’s M intermediate kinetics, especially the slow component in M intermediate decay. The half-life M_412s increased significantly in BRs in liposomes containing cationic lipids, while decreased in those in anionic liposomes.     

Residue specific partitioning of KL4 into phospholipid bilayers

KL₄, which has demonstrated success in the treatment of respiratory distress, is a synthetic helical, amphipathic peptide mimetic of lung surfactant protein B. The unusual periodicity of charged residues within KL₄ and its relatively high hydrophobicity distinguish it from canonical amphipathic helical peptides. Here we utilized site specific spin labeling of both lipids and the peptide coupled with EPR spectroscopy to discern the effects of KL₄ on lipid dynamics, the residue specific dynamics of hydrophobic regions within KL₄, and the partitioning depths of specific KL₄ residues into the DPPC/POPG and POPC/POPG lipid bilayers under physiologically relevant conditions. KL₄ induces alterations in acyl chain dynamics in a lipid-dependent manner, with the peptide partitioning more deeply into DPPC-rich bilayers. Combined with an earlier NMR study of changes in lipid dynamics on addition of KL₄ (V.C. Antharam et al., 2009), we are able to distinguish how KL₄ affects both collective bilayer motions and intramolecular acyl chain dynamics in a lipid-dependent manner. EPR power saturation results for spin labeled lipids demonstrate that KL₄ also alters the accessibility profiles of paramagnetic colliders in a lipid-dependent manner. Measurements of dynamics and depth parameters for individual spin-labeled residues within KL₄ are consistent with a model where the peptide partitions deeply into the lipid bilayers but lies parallel to the bilayer interface in both lipid environments; the depth of partitioning is dependent on the degree of lipid acyl chain saturation within the bilayer.     

Membrane interactions in small fast-tumbling bicelles as studied by 31P NMR

Small fast-tumbling bicelles are ideal for studies of membrane interactions at molecular level; they allow analysis of lipid properties using solution-state NMR. In the present study we used ³¹P NMR relaxation to obtain detailed information on lipid head-group dynamics. We explored the effect of two topologically different membrane-interacting peptides on bicelles containing either dimyristoylphosphocholine (DMPC), or a mixture of DMPC and dimyristoylphosphoglycerol (DMPG), and dihexanoylphosphocholine (DHPC). KALP21 is a model transmembrane peptide, designed to span a DMPC bilayer and dynorphin B is a membrane surface active neuropeptide. KALP21 causes significant increase in bicelle size, as evidenced by both dynamic light scattering and ³¹P T₂ relaxation measurements. The effect of dynorphin B on bicelle size is more modest, although significant effects on T₂ relaxation are observed at higher temperatures. A comparison of ³¹P T₁ values for the lipids with and without the peptides showed that dynorphin B has a greater effect on lipid head-group dynamics than KALP21, especially at elevated temperatures. From the field-dependence of T₁ relaxation data, a correlation time describing the overall lipid motion was derived. Results indicate that the positively charged dynorphin B decreases the mobility of the lipid molecules  – in particular for the negatively charged DMPG – while KALP21 has a more modest influence. Our results demonstrate that while a transmembrane peptide has severe effects on overall bilayer properties, the surface bound peptide has a more dramatic effect in reducing lipid head-group mobility. These observations may be of general importance for understanding peptide–membrane interactions.     

A Molecular Dynamics Investigation of Lipid Bilayer Perturbation by PIP2

Phosphoinositides like phosphatidylinositol 4,5-bisphosphate (PIP₂) are negatively charged lipids that play a pivotal role in membrane trafficking, signal transduction, and protein anchoring. We have designed a force field for the PIP₂ headgroup using quantum mechanical methods and characterized its properties inside a lipid bilayer using molecular dynamics simulations. Macroscopic properties such as area/headgroup, density profiles, and lipid order parameters calculated from these simulations agree well with the experimental values. However, microscopically, the PIP₂ introduces a local perturbation of the lipid bilayer. The average PIP₂ headgroup orientation of 45° relative to the bilayer normal induces a unique, distance-dependent organization of the lipids that surround PIP₂. The headgroups of these lipids preferentially orient closer to the bilayer normal. This perturbation creates a PIP₂ lipid microdomain with the neighboring lipids. We propose that the PIP₂ lipid microdomain enables the PIP₂ to function as a membrane-bound anchoring molecule.     

Binding of β-Amyloid (1-42) Peptide to Negatively Charged Phospholipid Membranes in the Liquid-Ordered State: Modeling and Experimental Studies

To explore the initial stages of amyloid β peptide (Aβ42) deposition on membranes, we have studied the interaction of Aβ42 in the monomeric form with lipid monolayers and with bilayers in either the liquid-disordered or the liquid-ordered (L_o) state, containing negatively charged phospholipids. Molecular dynamics (MD) simulations of the system have been performed, as well as experimental measurements. For bilayers in the L_o state, in the absence of the negatively charged lipids, interaction is weak and it cannot be detected by isothermal calorimetry. However, in the presence of phosphatidic acid, or of cardiolipin, interaction is detected by different methods and in all cases interaction is strongest with lower (2.5–5 mol %) than higher (10–20 mol %) proportions of negatively charged phospholipids. Liquid-disordered bilayers consistently allowed a higher Aβ42 binding than L_o ones. Thioflavin T assays and infrared spectroscopy confirmed a higher proportion of β-sheet formation under conditions when higher peptide binding was measured. The experimental results were supported by MD simulations. We used 100 ns MD to examine interactions between Aβ42 and three different 512 lipid bilayers consisting of palmitoylsphingomyelin, dimyristoyl phosphatidic acid, and cholesterol in three different proportions. MD pictures are different for the low- and high-charge bilayers, in the former case the peptide is bound through many contact points to the bilayer, whereas for the bilayer containing 20 mol % anionic phospholipid only a small fragment of the peptide appears to be bound. The MD results indicate that the binding and fibril formation on the membrane surface depends on the composition of the bilayer, and is the result of a subtle balance of many inter- and intramolecular interactions between the Aβ42 and membrane.     

Flexible charged macromolecules on mixed fluid lipid membranes: theory and Monte Carlo simulations

Fluid membranes containing charged lipids enhance binding of oppositely charged proteins by mobilizing these lipids into the interaction zone, overcoming the concomitant entropic losses due to lipid segregation and lower conformational freedom upon macromolecule adsorption. We study this energetic-entropic interplay using Monte Carlo simulations and theory. Our model system consists of a flexible cationic polyelectrolyte, interacting, via Debye-Hückel and short-ranged repulsive potentials, with membranes containing neutral lipids, 1% tetravalent, and 10% (or 1%) monovalent anionic lipids. Adsorption onto a fluid membrane is invariably stronger than to an equally charged frozen or uniform membrane. Although monovalent lipids may suffice for binding rigid macromolecules, polyvalent counter-lipids (e.g., phosphatidylinositol 4,5 bisphosphate), whose entropy loss upon localization is negligible, are crucial for binding flexible macromolecules, which lose conformational entropy upon adsorption. Extending Rosenbluth's Monte Carlo scheme we directly simulate polymer adsorption on fluid membranes. Yet, we argue that similar information could be derived from a biased superposition of quenched membrane simulations. Using a simple cell model we account for surface concentration effects, and show that the average adsorption probabilities on annealed and quenched membranes coincide at vanishing surface concentrations. We discuss the relevance of our model to the electrostatic-switch mechanism of, e.g., the myristoylated alanine-rich C kinase substrate protein.     

Calcium-Lipid Interactions Observed with Isotope-Edited Infrared Spectroscopy

Calcium ions bind to lipid membranes containing anionic lipids; however, characterizing the specific ion-lipid interactions in multicomponent membranes has remained challenging because it requires nonperturbative lipid-specific probes. Here, using a combination of isotope-edited infrared spectroscopy and molecular dynamics simulations, we characterize the effects of a physiologically relevant (2 mM) Ca²⁺ concentration on zwitterionic phosphatidylcholine and anionic phosphatidylserine lipids in mixed lipid membranes. We show that Ca²⁺ alters hydrogen bonding between water and lipid headgroups by forming a coordination complex involving the lipid headgroups and water. These interactions distort interfacial water orientations and prevent hydrogen bonding with lipid ester carbonyls. We demonstrate, experimentally, that these effects are more pronounced for the anionic phosphatidylserine lipids than for zwitterionic phosphatidylcholine lipids in the same membrane.     

Interaction of the cholesterol reducing agent simvastatin with zwitterionic DPPC and charged DPPG phospholipid membranes

Simvastatin is a lipid-lowering drug in the pharmaceutical group statins. Interaction of a drug with lipids may define its role in the system and be critical for its pharmacological activity. We examined the interactions of simvastatin with zwitterionic dipalmitoyl phosphatidylcholine (DPPC) and anionic dipalmitoyl phosphatidylglycerol (DPPG) multilamellar vesicles (MLVs) as a function of temperature at different simvastatin concentrations using Fourier transform infrared (FTIR) spectroscopy and differential scanning calorimetry (DSC). The FTIR results indicate that the effect of simvastatin on membrane structure and dynamics depends on the type of membrane lipids. In anionic DPPG MLVs, high simvastatin concentrations (12, 18, 24 mol%) change the position of the CH₂ antisymmetric stretching mode to lower wavenumber values, implying an ordering effect. However, in zwitterionic DPPC MLVs, high concentrations of simvastatin disorder systems both in the gel and liquid crystalline phases. Moreover, in DPPG and DPPC MLVs, simvastatin has opposite dual effects on membrane dynamics. The bandwidth of the CH₂ antisymmetric stretching modes increases in DPPG MLVs, implying an increase in the dynamics, whereas it decreases in DPPC MLVs. Simvastatin caused broadening of the phase transition peaks and formation of shoulders on the phase transition peaks in DSC curves, indicating multi-domain formations in the phospholipid membranes. Because physical features of membranes such as lipid order and fluidity may be changed with the bioactivity of drugs, opposing effects of simvastatin on the order and dynamics of neutral and charged phospholipids may be critical to deduce the action mechanism of the drug and estimate drug-membrane interactions.     

Dual Regulation of Diacylglycerol Kinase (DGK)-θ: POLYBASIC PROTEINS PROMOTE ACTIVATION BY PHOSPHOLIPIDS AND INCREASE SUBSTRATE AFFINITY*

Diacylglycerol kinases are important mediators of lipid signaling cascades, and insight into their regulation is of increasing interest. Using purified DGK-θ, we show that this isoform is subject to dual regulation and that the previously characterized stimulation by acidic phospholipids is dependent on the presence of a positively charged protein or peptide. Polybasic cofactors lowered the K_m for diacylglycerol at the membrane surface (K_m(surf)), and worked synergistically with acidic phospholipids to increase activity 10- to 30-fold, suggesting that the purified enzyme is autoinhibited. Vesicle pulldown studies showed that acidic phospholipids recruit polybasic cofactors to the vesicle surface but have little effect on the membrane association of DGK-θ, suggesting that a triad of enzyme, acidic lipid and basic protein are necessary for interfacial activity. Importantly, these data demonstrate that the interfacial association and catalytic activity of DGK-θ are independently regulated. Finally, we show that DGK-θ directly interacts with, and is activated by, basic proteins such as histone H1 and Tau with nm affinity, consistent with a potential role for a polybasic protein or protein domain in the activation of this enzyme.     

Molecular dynamics simulation of cation–phospholipid clustering in phospholipid bilayers: Possible role in stalk formation during membrane fusion

In this study, we performed all-atom long-timescale molecular dynamics simulations of phospholipid bilayers incorporating three different proportions of negatively charged lipids in the presence of K⁺, Mg^2 +, and Ca^2 + ions to systemically determine how membrane properties are affected by cations and lipid compositions. Our simulations revealed that the binding affinity of Ca^2 + ions with lipids is significantly stronger than that of K⁺ and Mg^2 + ions, regardless of the composition of the lipid bilayer. The binding of Ca^2 + ions to the lipids resulted in bilayers having smaller lateral areas, greater thicknesses, greater order, and slower rotation of their lipid head groups, relative to those of corresponding K⁺- and Mg^2 +-containing systems. The Ca^2 + ions bind preferentially to the phosphate groups of the lipids. The complexes formed between the cations and the lipids further assembled to form various multiple-cation-centered clusters in the presence of anionic lipids and at higher ionic strength—most notably for Ca^2 +. The formation of cation–lipid complexes and clusters dehydrated and neutralized the anionic lipids, creating a more-hydrophobic environment suitable for membrane aggregation. We propose that the formation of Ca^2 +–phospholipid clusters across apposed lipid bilayers can work as a “cation glue” to adhere apposed membranes together, providing an adequate configuration for stalk formation during membrane fusion.     

Molecular dynamics simulation of cation-phospholipid clustering in phospholipid bilayers: possible role in stalk formation during membrane fusion

In this study, we performed all-atom long-timescale molecular dynamics simulations of phospholipid bilayers incorporating three different proportions of negatively charged lipids in the presence of K(+), Mg(2+), and Ca(2+) ions to systemically determine how membrane properties are affected by cations and lipid compositions. Our simulations revealed that the binding affinity of Ca(2+) ions with lipids is significantly stronger than that of K(+) and Mg(2+) ions, regardless of the composition of the lipid bilayer. The binding of Ca(2+) ions to the lipids resulted in bilayers having smaller lateral areas, greater thicknesses, greater order, and slower rotation of their lipid head groups, relative to those of corresponding K(+)- and Mg(2+)-containing systems. The Ca(2+) ions bind preferentially to the phosphate groups of the lipids. The complexes formed between the cations and the lipids further assembled to form various multiple-cation-centered clusters in the presence of anionic lipids and at higher ionic strength-most notably for Ca(2+). The formation of cation-lipid complexes and clusters dehydrated and neutralized the anionic lipids, creating a more-hydrophobic environment suitable for membrane aggregation. We propose that the formation of Ca(2+)-phospholipid clusters across apposed lipid bilayers can work as a "cation glue" to adhere apposed membranes together, providing an adequate configuration for stalk formation during membrane fusion.     

Alzheimer Aβ peptide interactions with lipid membranes: Fibrils,oligomers and polymorphic amyloid channels

Fibrillar aggregates of misfolded amyloid proteins are involved in a variety of diseases such as Alzheimer disease (AD), type 2 diabetes, Parkinson, Huntington and prion-related diseases. In the case of AD amyloid β (Aβ) peptides, the toxicity of amyloid oligomers and larger fibrillar aggregates is related to perturbing the biological function of the adjacent cellular membrane. We used atomistic molecular dynamics (MD) simulations of Aβ_9–40 fibrillar oligomers modeled as protofilament segments, including lipid bilayers and explicit water molecules, to probe the first steps in the mechanism of Aβ-membrane interactions. Our study identified the electrostatic interaction between charged peptide residues and the lipid headgroups as the principal driving force that can modulate the further penetration of the C-termini of amyloid fibrils or fibrillar oligomers into the hydrophobic region of lipid membranes. These findings advance our understanding of the detailed molecular mechanisms and the effects related to Aβ-membrane interactions, and suggest a polymorphic structural character of amyloid ion channels embedded in lipid bilayers. While inter-peptide hydrogen bonds leading to the formation of β-strands may still play a stabilizing role in amyloid channel structures, these may also present a significant helical content in peptide regions (e.g., termini) that are subject to direct interactions with lipids rather than with neighboring Aβ peptides.