Spontaneous adsorption of coiled-coil model peptides K and E to a mixed lipid bilayer

Enhanced Liposomal Drug Delivery Via Membrane Fusion Triggered by Dimeric Coiled-Coil Peptides

An ideal nanomedicine system improves the therapeutic efficacy of drugs. However, most nanomedicines enter cells via endosomal/lysosomal pathways and only a small fraction of the cargo enters the cytosol inducing therapeutic effects. To circumvent this inefficiency, alternative approaches are desired. Inspired by fusion machinery found in nature, synthetic lipidated peptide pair E4/K4 is used to induce membrane fusion previously. Peptide K4 interacts specifically with E4, and it has a lipid membrane affinity and resulting in membrane remodeling. To design efficient fusogens with multiple interactions, dimeric K4 variants are synthesized to improve fusion with E4-modified liposomes and cells. The secondary structure and self-assembly of dimers are studied; the parallel PK4 dimer forms temperature-dependent higher-order assemblies, while linear K4 dimers form tetramer-like homodimers. The structures and membrane interactions of PK4 are supported by molecular dynamics simulations. Upon addition of E4, PK4 induced the strongest coiled-coil interaction resulting in a higher liposomal delivery compared to linear dimers and monomer. Using a wide spectrum of endocytosis inhibitors, membrane fusion is found to be the main cellular uptake pathway. Doxorubicin delivery results in efficient cellular uptake and concomitant antitumor efficacy. These findings aid the development of efficient delivery systems of drugs into cells using liposome-cell fusion strategies.

Not sorcery after all: Roles of multiple charged residues in membrane insertion of gasdermin-A3

Gasdermins execute programmatory cell death, known as pyroptosis, by forming medium-sized membrane pores. Recently, the molecular structure of those pores as well as the diversity in their shape and size have been revealed by cryoTEM and atomic force microscopy, respectively. Even though a growth of smaller to larger oligomers and reshaping from slits to rings could be documented, the initiation of the gasdermin pore formation remains a mystery. In one hypothesis, gasdermin monomers insert into membranes before associating into oligomeric pores. In the other hypothesis, gasdermin oligomers preassemble on the membrane surface prior to membrane insertion. Here, by studying the behavior of monomeric membrane-inserted gasdermin-A3 (GSDMA3), we unveil that a monomeric gasdermin prefers the membrane-adsorbed over the membrane-inserted state. Our results thus support the hypothesis of oligomers preassembling on the membrane surface before membrane penetration. At the same time, our simulations of small membrane-inserted arcs of GSDMA3 suggest that the inserting oligomer can be small and does not have to comprise a full ring of approximately 26-30 subunits. Moreover, our simulations have revealed an astonishingly large impact of salt-bridge formation and protein surroundings on the transmembrane passage of charged residues, reducing the energetic cost by up to 53% as compared to their free forms. The here observed free energy barrier of mere 5.6 kcal/mol for the membrane insertion of monomeric GSDMA3 explains the surprising ability of gasdermins to spontaneously self-insert into cellular membranes.

Theater in the Self-Cleaning Cell: Intrinsically Disordered Proteins or Protein Regions Acting with Membranes in Autophagy

Intrinsically disordered proteins and protein regions (IDPs/IDPRs) are mainly involved in signaling pathways, where fast regulation, temporal interactions, promiscuous interactions, and assemblies of structurally diverse components including membranes are essential. The autophagy pathway builds, de novo, a membrane organelle, the autophagosome, using carefully orchestrated interactions between proteins and lipid bilayers. Here, we discuss molecular mechanisms related to the protein disorder-based interactions of the autophagy machinery with membranes. We describe not only membrane binding phenomenon, but also examples of membrane remodeling processes including membrane tethering, bending, curvature sensing, and/or fragmentation of membrane organelles such as the endoplasmic reticulum, which is an important membrane source as well as cargo for autophagy. Summary of the current state of knowledge presented here will hopefully inspire new studies. A profound understanding of the autophagic protein-membrane interface is essential for advancements in therapeutic interventions against major human diseases, in which autophagy is involved including neurodegeneration, cancer as well as cardiovascular, metabolic, infectious, musculoskeletal, and other disorders.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of β2-Adrenergic Receptor and β-arrestin2

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β₂-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

Monitoring the binding and insertion of a single transmembrane protein by an insertase

Cells employ highly conserved families of insertases and translocases to insert and fold proteins into membranes. How insertases insert and fold membrane proteins is not fully known. To investigate how the bacterial insertase YidC facilitates this process, we here combine single-molecule force spectroscopy and fluorescence spectroscopy approaches, and molecular dynamics simulations. We observe that within 2 ms, the cytoplasmic α-helical hairpin of YidC binds the polypeptide of the membrane protein Pf3 at high conformational variability and kinetic stability. Within 52 ms, YidC strengthens its binding to the substrate and uses the cytoplasmic α-helical hairpin domain and hydrophilic groove to transfer Pf3 to the membrane-inserted, folded state. In this inserted state, Pf3 exposes low conformational variability such as typical for transmembrane α-helical proteins. The presence of YidC homologues in all domains of life gives our mechanistic insight into insertase-mediated membrane protein binding and insertion general relevance for membrane protein biogenesis.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity

BACKGROUND

Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments.

RESULTS

To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly.

CONCLUSIONS

The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

Coarse-Grained Molecular Model for the Glycosylphosphatidylinositol Anchor with and without Protein

Glycosylphosphatidylinositol (GPI) anchors are a unique class of complex glycolipids that anchor a great variety of proteins to the extracellular leaflet of plasma membranes of eukaryotic cells. These anchors can exist either with or without an attached protein called GPI-anchored protein (GPI-AP) both in vitro and in vivo. Although GPIs are known to participate in a broad range of cellular functions, it is to a large extent unknown how these are related to GPI structure and composition. Their conformational flexibility and microheterogeneity make it difficult to study them experimentally. Simplified atomistic models are amenable to all-atom computer simulations in small lipid bilayer patches but not suitable for studying their partitioning and trafficking in complex and heterogeneous membranes. Here, we present a coarse-grained model of the GPI anchor constructed with a modified version of the MARTINI force field that is suited for modeling carbohydrates, proteins, and lipids in an aqueous environment using MARTINI's polarizable water. The nonbonded interactions for sugars were reparametrized by calculating their partitioning free energies between polar and apolar phases. In addition, sugar-sugar interactions were optimized by adjusting the second virial coefficients of osmotic pressures for solutions of glucose, sucrose, and trehalose to match with experimental data. With respect to the conformational dynamics of GPI-anchored green fluorescent protein, the accessible time scales are now at least an order of magnitude larger than for the all-atom system. This is particularly important for fine-tuning the mutual interactions of lipids, carbohydrates, and amino acids when comparing to experimental results. We discuss the prospective use of the coarse-grained GPI model for studying protein-sorting and trafficking in membrane models.

Modulation of Coiled-Coil Binding Strength and Fusogenicity through Peptide Stapling

Peptide stapling is a technique which has been widely employed to constrain the conformation of peptides. One of the effects of such a constraint can be to modulate the interaction of the peptide with a binding partner. Here, a cysteine bis-alkylation stapling technique was applied to generate structurally isomeric peptide variants of a heterodimeric coiled-coil forming peptide. These stapled variants differed in the position and size of the formed macrocycle. C-terminal stapling showed the most significant changes in peptide structure and stability, with calorimetric binding analysis showing a significant reduction of binding entropy for stapled variants. This entropy reduction was dependent on cross-linker size and was accompanied by a change in binding enthalpy, illustrating the effects of preorganization. The stapled peptide, along with its binding partner, were subsequently employed as fusogens in a liposome model system. An increase in both lipid- and content-mixing was observed for one of the stapled peptide variants: this increased fusogenicity was attributed to increased coiled-coil binding but not to membrane affinity, an interaction theorized to be a primary driving force in this fusion system.

Interaction of SNARE Mimetic Peptides with Lipid bilayers: Effects of Secondary Structure, Bilayer Composition and Lipid Anchoring

The coiled-coil forming peptides 'K' enriched in lysine and 'E' enriched in glutamic acid have been used as a minimal SNARE mimetic system for membrane fusion. Here we describe atomistic molecular dynamics simulations to characterize the interactions of these peptides with lipid bilayers for two different compositions. For neutral phosphatidylcholine (PC)/phosphatidylethanolamine (PE) bilayers the peptides experience a strong repulsive barrier against adsorption, also observed in potential of mean force (PMF) profiles calculated with umbrella sampling. For peptide K, a minimum of -12 k_BT in the PMF provides an upper bound for the binding free energy whereas no stable membrane bound state could be observed for peptide E. In contrast, the electrostatic interactions with negatively charged phosphatidylglycerol (PG) lipids lead to fast adsorption of both peptides at the head-water interface. Experimental data using fluorescently labeled peptides confirm the stronger binding to PG containing bilayers. Lipid anchors have little effect on the peptide-bilayer interactions or peptide structure, when the peptide also binds to the bilayer in the absence of a lipid anchor. For peptide E, which does not bind to the PC bilayer without a lipid anchor, the presence of such an anchor strengthens the electrostatic interactions between the charged side chains and the zwitterionic head-groups and leads to a stabilization of the peptide's helical fold by the membrane.

Intrinsic lipid binding activity of ATG16L1 supports efficient membrane anchoring and autophagy

Membrane targeting of autophagy-related complexes is an important step that regulates their activities and prevents their aberrant engagement on non-autophagic membranes. ATG16L1 is a core autophagy protein implicated at distinct phases of autophagosome biogenesis. In this study, we dissected the recruitment of ATG16L1 to the pre-autophagosomal structure (PAS) and showed that it requires sequences within its coiled-coil domain (CCD) dispensable for homodimerisation. Structural and mutational analyses identified conserved residues within the CCD of ATG16L1 that mediate direct binding to phosphoinositides, including phosphatidylinositol 3-phosphate (PI3P). Mutating putative lipid binding residues abrogated the localisation of ATG16L1 to the PAS and inhibited LC3 lipidation. On the other hand, enhancing lipid binding of ATG16L1 by mutating negatively charged residues adjacent to the lipid binding motif also resulted in autophagy inhibition, suggesting that regulated recruitment of ATG16L1 to the PAS is required for its autophagic activity. Overall, our findings indicate that ATG16L1 harbours an intrinsic ability to bind lipids that plays an essential role during LC3 lipidation and autophagosome maturation.

Structural Model of the mIgM B-Cell Receptor Transmembrane Domain From Self-Association Molecular Dynamics Simulations

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-β heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-β heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes.

Folding and Lipid Composition Determine Membrane Interaction of the Disordered Protein COR15A

Plants from temperate climates, such as the model plant Arabidopsis thaliana, are challenged with seasonal low temperatures that lead to increased freezing tolerance in fall in a process termed cold acclimation. Among other adaptations, this involves the accumulation of cold-regulated (COR) proteins, such as the intrinsically disordered chloroplast-localized protein COR15A. Together with its close homolog COR15B, it stabilizes chloroplast membranes during freezing. COR15A folds into amphipathic α-helices in the presence of high concentrations of low-molecular-mass crowders or upon dehydration. Under these conditions, the (partially) folded protein binds peripherally to membranes. In our study, we have used coarse-grained molecular dynamics simulations to elucidate the details of COR15A-membrane binding and its effects on membrane structure and dynamics. Simulation results indicate that at least partial folding of COR15A and the presence of highly unsaturated galactolipids in the membranes are necessary for efficient membrane binding. The bound protein is stabilized on the membrane by interactions of charged and polar amino acids with galactolipid headgroups and by interactions of hydrophobic amino acids with the upper part of the fatty acyl chains. Experimentally, the presence of liposomes made from a mixture of lipids mimicking chloroplast membranes induces additional folding in COR15A under conditions of partial dehydration, in agreement with the simulation results.

Membrane-Fusogen Distance Is Critical for Efficient Coiled-Coil-Peptide-Mediated Liposome Fusion

We have developed a model system for membrane fusion that utilizes lipidated derivatives of a heterodimeric coiled-coil pair dubbed E3 (EIAALEK)3 and K3 (KIAALKE)3. In this system, peptides are conjugated to a lipid anchor via a poly(ethylene glycol) (PEG) spacer, and this contribution studies the influence of the PEG spacer length, coupled with the type of lipid anchor, on liposome-liposome fusion. The effects of these modifications on peptide secondary structure, their interactions with liposomes, and their ability to mediate fusion were studied using a variety of different content mixing experiments and CD spectroscopy. Our results demonstrate the asymmetric role of the peptides in the fusion process because alterations to the PEG spacer length affect E3 and K3 differently. We conclude that negatively charged E3 acts as a "handle" for positively charged K3 and facilitates liposome docking, the first stage of the fusion process, through coiled-coil formation. The efficacy of this E3 handle is enhanced by longer spacer lengths. K3 directs the fusion process via peptide-membrane interactions, but the length of the PEG spacer plays two competing roles: a PEG4/PEG8 spacer length is optimal for membrane destabilization; however, a PEG12 spacer increases the fusion efficiency over time by improving the peptide accessibility for successive fusion events. Both the anchor type and spacer length affect the peptide structure; a cholesterol anchor appears to enhance K3-membrane interactions and thus mediates fusion more efficiently.

Geometry of the Contact Zone between Fused Membrane-Coated Beads Mimicking Cell-Cell Fusion

The fusion of lipid membranes is a key process in biology. It enables cells and organelles to exchange molecules with their surroundings, which otherwise could not cross the membrane barrier. To study such complex processes we use simplified artificial model systems, i.e., an optical fusion assay based on membrane-coated glass spheres. We present a technique to analyze membrane-membrane interactions in a large ensemble of particles. Detailed information on the geometry of the fusion stalk of fully fused membranes is obtained by studying the diffusional lipid dynamics with fluorescence recovery after photobleaching experiments. A small contact zone is a strong obstruction for the particle exchange across the fusion spot. With the aid of computer simulations, fluorescence-recovery-after-photobleaching recovery times of both fused and single-membrane-coated beads allow us to estimate the size of the contact zones between two membrane-coated beads. Minimizing delamination and bending energy leads to minimal angles close to those geometrically allowed.

A Coiled-Coil Peptide Shaping Lipid Bilayers upon Fusion

A system based on two designed peptides, namely the cationic peptide K, (KIAALKE)3, and its complementary anionic counterpart called peptide E, (EIAALEK)3, has been used as a minimal model for membrane fusion, inspired by SNARE proteins. Although the fact that docking of separate vesicle populations via the formation of a dimeric E/K coiled-coil complex can be rationalized, the reasons for the peptides promoting fusion of vesicles cannot be fully explained. Therefore it is of significant interest to determine how the peptides aid in overcoming energetic barriers during lipid rearrangements leading to fusion. In this study, investigations of the peptides' interactions with neutral PC/PE/cholesterol membranes by fluorescence spectroscopy show that tryptophan-labeled K∗ binds to the membrane (KK∗ ∼6.2 103 M-1), whereas E∗ remains fully water-solvated. 15N-NMR spectroscopy, depth-dependent fluorescence quenching, CD-spectroscopy experiments, and MD simulations indicate a helix orientation of K∗ parallel to the membrane surface. Solid-state 31P-NMR of oriented lipid membranes was used to study the impact of peptide incorporation on lipid headgroup alignment. The membrane-immersed K∗ is found to locally alter the bilayer curvature, accompanied by a change of headgroup orientation relative to the membrane normal and of the lipid composition in the vicinity of the bound peptide. The NMR results were supported by molecular dynamics simulations, which showed that K reorganizes the membrane composition in its vicinity, induces positive membrane curvature, and enhances the lipid tail protrusion probability. These effects are known to be fusion relevant. The combined results support the hypothesis for a twofold role of K in the mechanism of membrane fusion: 1) to bring opposing membranes into close proximity via coiled-coil formation and 2) to destabilize both membranes thereby promoting fusion.

Membrane-Mediated Oligomerization of G Protein Coupled Receptors and Its Implications for GPCR Function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signaling complexes. Recent findings further suggest that the surrounding membrane, i.e., its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function.

The function of the two-pore channel TPC1 depends on dimerization of its carboxy-terminal helix

Two-pore channels (TPCs) constitute a family of intracellular cation channels with diverse permeation properties and functions in animals and plants. In the model plant Arabidopsis, the vacuolar cation channel TPC1 is involved in propagation of calcium waves and in cation homeostasis. Here, we discovered that the dimerization of a predicted helix within the carboxyl-terminus (CTH) is essential for the activity of TPC1. Bimolecular fluorescence complementation and co-immunoprecipitation demonstrated the interaction of the two C-termini and pointed towards the involvement of the CTH in this process. Synthetic CTH peptides dimerized with a dissociation constant of 3.9 µM. Disruption of this domain in TPC1 either by deletion or point mutations impeded the dimerization and cation transport. The homo-dimerization of the CTH was analyzed in silico using coarse-grained molecular dynamics (MD) simulations for the study of aggregation, followed by atomistic MD simulations. The simulations revealed that the helical region of the wild type, but not a mutated CTH forms a highly stable, antiparallel dimer with characteristics of a coiled-coil. We propose that the voltage- and Ca(2+)-sensitive conformation of TPC1 depends on C-terminal dimerization, adding an additional layer to the complex regulation of two-pore cation channels.

Synaptobrevin Transmembrane Domain Dimerization Studied by Multiscale Molecular Dynamics Simulations

Synaptic vesicle fusion requires assembly of the SNARE complex composed of SNAP-25, syntaxin-1, and synaptobrevin-2 (sybII) proteins. The SNARE proteins found in vesicle membranes have previously been shown to dimerize via transmembrane (TM) domain interactions. While syntaxin homodimerization is supposed to promote the transition from hemifusion to complete fusion, the role of synaptobrevin's TM domain association in the fusion process remains poorly understood. Here, we combined coarse-grained and atomistic simulations to model the homodimerization of the sybII transmembrane domain and of selected TM mutants. The wild-type helix is shown to form a stable, right-handed dimer with the most populated helix-helix interface, including key residues predicted in a previous mutagenesis study. In addition, two alternative binding interfaces were discovered, which are essential to explain the experimentally observed higher-order oligomerization of sybII. In contrast, only one dimerization interface was found for a fusion-inactive poly-Leu mutant. Moreover, the association kinetics found for this mutant is lower as compared to the wild-type. These differences in dimerization between the wild-type and the poly-Leu mutant are suggested to be responsible for the reported differences in fusogenic activity between these peptides. This study provides molecular insight into the role of TM sequence specificity for peptide aggregation in membranes.

All-atom simulations and free-energy calculations of coiled-coil peptides with lipid bilayers: binding strength, structural transition, and effect on lipid dynamics

Peptides E and K, which are synthetic coiled-coil peptides for membrane fusion, were simulated with lipid bilayers composed of lipids and cholesterols at different ratios using all-atom models. We first calculated free energies of binding from umbrella sampling simulations, showing that both E and K peptides tend to adsorb onto the bilayer surface, which occurs more strongly in the bilayer composed of smaller lipid headgroups. Then, unrestrained simulations show that K peptides more deeply insert into the bilayer with partially retaining the helical structure, while E peptides less insert and predominantly become random coils, indicating the structural transition from helices to random coils, in quantitative agreement with experiments. This is because K peptides electrostatically interact with lipid phosphates, as well as because hydrocarbons of lysines of K peptide are longer than those of glutamic acids of E peptide and thus form stronger hydrophobic interactions with lipid tails. This deeper insertion of K peptide increases the bilayer dynamics and a vacancy below the peptide, leading to the rearrangement of smaller lipids. These findings help explain the experimentally observed or proposed differences in the insertion depth, binding strength, and structural transition of E and K peptides, and support the snorkeling effect.

Interplay between Lipid Interaction and Homo-coiling of Membrane-Tethered Coiled-Coil Peptides

The designed coiled-coil-forming peptides E [(EIAALEK)3] and K [(KIAALKE)3] are known to trigger efficient membrane fusion when they are tethered to lipid vesicles in the form of lipopeptides. Knowledge of their secondary structure is a key element in understanding their role in membrane fusion. Special conditions can be found at the interface of the membrane, where the peptides are confined in close proximity to other peptide molecules as well as to the lipid interface. Consequently, different structural states were proposed for the peptides when tethered to this interface. Due to the multitude of possible states, determining the structure solely on the basis of circular dichroism (CD) spectra at a single temperature can be misleading. In addition, it has not yet been possible to unambiguously distinguish between the membrane-bound and the coiled-coil states of these peptides by means of infrared (IR) spectroscopy due to their very similar amide I' bands. Here, the molecular basis of this similarity is investigated by means of site-specific (13)C-labeled FTIR spectroscopy. Structural similarities between the membrane-interacting helix of K and the homo-coiled-coil-forming helix of E are shown to cause the similar spectroscopic properties. Furthermore, the peptide structure is investigated using temperature-dependent CD and IR spectroscopy, and it is shown that the different states can be distinguished on the basis of their thermal behavior. It is shown that the two peptides behave fundamentaly differently when tethered to the lipid membrane, which implies that their role during membrane fusion is different and the mechanism of this process is asymmetric.

Phase Transition of Glycolipid Membranes Studied by Coarse-Grained Simulations

Glycolipids are important components of biological membranes. High concentrations of glycolipids are particularly found in lipid rafts, which take part in many physiological phenomena. This different partitioning and interaction pattern of glycolipids in the membrane as compared to those of phospholipids are likely due to their different chemical structures: the polar regions of glycosphingolipids can be even larger than for their hydrophobic moieties, giving rise to a rich conformational landscape. Here we study the influence of glycosphingolipids galactosylceramide (GCER) and monosialotetrahexosylganglioside (GM1) on the structural and thermodynamic properties of a phospholipid (DPPC) bilayer. Using the method of coarse-grained molecular dynamics simulation we show that both glycolipids increase the phase-transition temperature of phospholipid membranes and that the extent of this increase depends on the headgroup size and structure. GM1 shows a strong tendency to form mixed clusters with phospholipids, thereby stabilizing the membrane. In contrast, GCER is dispersed in the membrane. By occupying the interstitial space between phospholipids it causes a tighter packing of the lipids in the membrane.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.