State-dependent protein-lipid interactions of a pentameric ligand-gated ion channel in a neuronal membrane

Pentameric ligand-gated ion channels (pLGICs) are receptor proteins that are sensitive to their membrane environment, but the mechanism for how lipids modulate function under physiological conditions in a state dependent manner is not known. The glycine receptor is a pLGIC whose structure has been resolved in different functional states. Using a realistic model of a neuronal membrane coupled with coarse-grained molecular dynamics simulations, we demonstrate that some key lipid-protein interactions are dependent on the receptor state, suggesting that lipids may regulate the receptor’s conformational dynamics. Comparison with existing structural data confirms known lipid binding sites, but we also predict further protein-lipid interactions including a site at the communication interface between the extracellular and transmembrane domain. Moreover, in the active state, cholesterol can bind to the binding site of the positive allosteric modulator ivermectin. These protein-lipid interaction sites could in future be exploited for the rational design of lipid-like allosteric drugs.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.     

Transbilayer movement of phospholipids in biogenic membranes

Biogenic membranes contain the enzymes that synthesize the cell's membrane lipids, of which the phospholipids are the most widespread throughout nature. Being synthesized at and inserted into the cytoplasmic leaflet of biogenic membranes, the phospholipids must migrate to the opposite leaflet to ensure balanced growth of the membrane. In this review, the current knowledge of transbilayer movement of phospholipids in biogenic membranes is summarized and the available data are compared to what is known about lipid translocation in other membranes. On the basis of this, a mechanism is proposed, in which phospholipid translocation in biogenic membranes is mediated via membrane-spanning segments of a subset of proteins, characterized by a small number of transmembrane helices. We speculate that proteins of this subset facilitate lipid translocation via the protein-lipid interface, because they display more dynamic behavior and engage in less stable protein-lipid interactions than larger membrane proteins.     

Coarse grain lipid-protein molecular interactions and diffusion with MsbA flippase

Coarse-grained (CG) modeling has proven effective for simulating lipid bilayer dynamics on scales of biological interest. Modeling the dynamics of flexible membrane proteins within the bilayer, on the other hand, poses a considerable challenge due to the complexity of the folding or conformational landscape. In the present work, the multiscale coarse-graining method is applied to atomistic peptide-lipid "soup" simulations to develop a general set of CG protein-lipid interaction potentials. The reduced model was constructed to be compatible with recent solvent-free CG models developed for protein-protein folding and lipid-lipid model bilayer interactions. The utility of the force field was demonstrated by molecular dynamics simulation of the MsbA ABC transporter in a mixed DOPC/DOPE bilayer. An elastic network was parameterized to restrain the MsbA dimer in its open, closed and hydrolysis intermediate conformations and its impact on domain flexibility was examined. Conformational stability enabled long-time dynamics simulation of MsbA freely diffusing in a 25 nm membrane patch. Three-dimensional density analysis revealed that a shell of weakly bound "annular lipids" solvate the membrane accessible surface of MsbA and its internal substrate-binding chamber. The annular lipid binding modes, along with local perturbations in head group structure, are a function of the orientation of grooves formed between transmembrane helices and may influence the alternating access mechanism of substrate entry and translocation.     

Stoichiometry of lipid interactions with transmembrane proteins--Deduced from the 3D structures

The stoichiometry of the first shell of lipids interacting with a transmembrane protein is defined operationally by the population of spin-labeled lipid chains whose motion is restricted directly by the protein. Interaction stoichiometries have been determined experimentally for a wide range of alpha-helical integral membrane proteins by using spin-label ESR spectroscopy. Here, we determine the spatially defined number of first-shell lipids at the hydrophobic perimeter of integral membrane proteins whose 3D structure has been determined by X-ray crystallography and lipid-protein interactions characterized by spin-labeling. Molecular modeling is used to build a single shell of lipids surrounding transmembrane structures derived from the PDB. Constrained energy optimization of the protein-lipid assemblies is performed by molecular mechanics. For relatively small proteins (up to 7-12 transmembrane helices), the geometrical first shell corresponds to that defined experimentally by perturbation of the lipid-chain dynamics. For larger, multi-subunit alpha-helical proteins, the lipids perturbed directly by the protein may either exceed or be less in number than those that can be accommodated at the intramembranous perimeter. In these latter cases, the motionally restricted spin-labeled lipids can be augmented by intercalation, or can correspond to a specific subpopulation at the protein interface, respectively. For monomeric beta-barrel proteins, the geometrical lipid stoichiometry corresponds to that determined from lipid mobility for a 22-stranded barrel, but fewer lipids are motionally restricted than can be accommodated around an eight-stranded barrel. Deviations from the geometrical first shell, in the beta-barrel case, are for the smaller protein with a highly curved barrel.     

Assembly of the M2 Tetramer Is Strongly Modulated by Lipid Chain Length

The influenza virus matrix protein 2 (M2) assembles into a tetramer in the host membrane during viral uncoating and maturation. It has been used as a model system to understand the relative contributions of protein-lipid and protein-protein interactions to membrane protein structure and association. Here we investigate the effect of lipid chain length on the association of the M2 transmembrane domain into tetramers using Förster resonance energy transfer. We observe that the interactions between the M2 helices are much stronger in 1,2-dilauroyl-sn-glycero-3-phosphocholine than in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. Thus, lipid chain length and bilayer thickness not only modulate peptide interactions, but could also be a major determinant of the association of transmembrane helices into functional membrane protein oligomers.     

Protein-lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring

Biological membranes are characterized by a heterogeneous composition, which is not only manifested in the wide variety of their components, but also in aspects like the lateral organization, topology, and conformation of proteins and lipids. In bringing about the correct membrane structure, protein-lipid interactions can be expected to play a prominent role. The extent of hydrophobic matching between transmembrane protein segments and lipids potentially constitutes a versatile director of membrane organization, because a tendency to avoid hydrophobic mismatch could result in compensating adaptations such as tilt of the transmembrane segment or segregation into distinct domains. Also, interfacial interactions between lipid headgroups and the aromatic and charged residues that typically flank transmembrane domains may act as an organizing element. In this review, we discuss the numerous model studies that have systematically explored the influence of hydrophobic matching and interfacial anchoring on membrane structure. Designed peptides consisting of a polyleucine or polyleucine/alanine hydrophobic stretch, which is flanked on both sides by tryptophan or lysine residues, reflect the general layout of transmembrane protein segments. It is shown for phosphatidylcholine bilayers and for other model membranes that these peptides adapt a transmembrane topology without extensive peptide or lipid adaptations under conditions of hydrophobic matching, but that significant rearrangements can result from hydrophobic mismatch. Moreover, these effects depend on the nature of the flanking residues, implying a modulation of the mismatch response by interfacial interactions of the flanking residues. The implications of these model studies for the organization of biomembranes are discussed in the context of recent experiments with more complex systems.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Molecular dynamics simulations of the ErbB-2 transmembrane domain within an explicit membrane environment: comparison with vacuum simulations

Two 500-ps molecular dynamics simulations performed on the single transmembrane domain of the ErbB-2 tyrosine kinase receptor immersed in a fully solvated dilauroylphosphatidyl-ethanolamine bilayer (DLPE) are compared to vacuum simulations. One membrane simulation shows that the initial alpha helix undergoes a local pi helix conversion in the peptide part embedded in the membrane core similar to that found in simulation vacuum. Lipid/water/peptide interaction analysis shows that in the helix core, the intramolecular peptide interactions are largely dominant compared to the interactions with water and lipids whereas the helix extremities are much more sensitive to these interactions at the membrane interfaces. Our results suggest that simulations in a lipid environment are required to understand the dynamics of transmembrane helices, but can be reasonably supplemented by in vacuo simulations to explore rapidly its conformational space and to describe the internal deformation of the hydrophobic core.     

Influence of the extracellular domain size on the dynamic behavior of membrane proteins

The dynamic behavior of plasma membrane proteins mediates various cellular processes such as cellular motility, communication, and signaling. It is widely accepted that the dynamics of the membrane proteins is determined either by the interactions of the transmembrane domain with the surrounding lipids or by the interactions of the intracellular domain with cytosolic components such as cortical actin. Although initiation of different cellular signaling events at the plasma membrane has been attributed to the extracellular domain (ECD) properties recently, the impact of ECDs on the dynamic behavior of membrane proteins is rather unexplored. Here, we investigate how ECD properties influence protein dynamics in the lipid bilayer by reconstituting ECDs of different sizes or glycosylation in model membrane systems and analyzing ECD-driven protein sorting in lipid domains as well as protein mobility. Our data show that increasing the ECD mass or glycosylation leads to a decrease in ordered domain partitioning and diffusivity. Our data reconcile different mechanisms proposed for the initiation of cellular signaling by linking the ECD size of membrane proteins with their localization and diffusion dynamics in the plasma membrane.     

Functional competition within a membrane: Lipid recognition vs. transmembrane helix oligomerization

Binding of specific lipids to large, polytopic membrane proteins is well described, and it is clear that such lipids are crucial for protein stability and activity. In contrast, binding of defined lipid species to individual transmembrane helices and regulation of transmembrane helix monomer–oligomer equilibria by binding of distinct lipids is a concept, which has emerged only lately. Lipids bind to single-span membrane proteins, both in the juxta-membrane region as well as in the hydrophobic membrane core. While some interactions counteract transmembrane helix oligomerization, in other cases lipid binding appears to enhance oligomerization. As reversible oligomerization is involved in activation of many membrane proteins, binding of defined lipids to single-span transmembrane proteins might be a mechanism to regulate and/or fine-tune the protein activity. But how could lipid binding trigger the activity of a protein? How can binding of a single lipid molecule to a transmembrane helix affect the structure of a transmembrane helix oligomer, and consequently its signaling state? These questions are discussed in the present article based on recent results obtained with simple, single-span transmembrane proteins. This article is part of a Special Issue entitled: Lipid–protein interactions.     

Structural determinants of purple membrane assembly

The purple membrane is a two-dimensional crystalline lattice formed by bacteriorhodopsin and lipid molecules in the cytoplasmic membrane of Halobacterium salinarum. High-resolution structural studies, in conjunction with detailed knowledge of the lipid composition, make the purple membrane one of the best models for elucidating the forces that are responsible for the assembly and stability of integral membrane protein complexes. In this review, recent mutational efforts to identify the structural features of bacteriorhodopsin that determine its assembly in the purple membrane are discussed in the context of structural, calorimetric and reconstitution studies. Quantitative evidence is presented that interactions between transmembrane helices of neighboring bacteriorhodopsin molecules contribute to purple membrane assembly. However, other specific interactions, particularly between bacteriorhodopsin and lipid molecules, may provide the major driving force for assembly. Elucidating the molecular basis of protein-protein and protein-lipid interactions in the purple membrane may provide insights into the formation of integral membrane protein complexes in other systems.     

Self-association of transmembrane alpha-helices in model membranes: importance of helix orientation and role of hydrophobic mismatch

Interactions between transmembrane helices play a key role in almost all cellular processes involving membrane proteins. We have investigated helix-helix interactions in lipid bilayers with synthetic tryptophan-flanked peptides that mimic the membrane spanning parts of membrane proteins. The peptides were functionalized with pyrene to allow the self-association of the helices to be monitored by pyrene fluorescence and Trp-pyrene fluorescence resonance energy transfer (FRET). Specific labeling of peptides at either their N or C terminus has shown that helix-helix association occurs almost exclusively between antiparallel helices. Furthermore, computer modeling suggested that antiparallel association arises primarily from the electrostatic interactions between alpha-helix backbone atoms. We propose that such interactions may provide a force for the preferentially antiparallel association of helices in polytopic membrane proteins. Helix-helix association was also found to depend on the lipid environment. In bilayers of dioleoylphosphatidylcholine, in which the hydrophobic length of the peptides approximately matched the bilayer thickness, association between the helices was found to require peptide/lipid ratios exceeding 1/25. Self-association of the helices was promoted by either increasing or decreasing the bilayer thickness, and by adding cholesterol. These results indicate that helix-helix association in membrane proteins can be promoted by unfavorable protein-lipid interactions.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

The Molecular Basis of Polyunsaturated Fatty Acid Interactions with the Shaker Voltage-Gated Potassium Channel

Voltage-gated potassium (K_V) channels are membrane proteins that respond to changes in membrane potential by enabling K⁺ ion flux across the membrane. Polyunsaturated fatty acids (PUFAs) induce channel opening by modulating the voltage-sensitivity, which can provide effective treatment against refractory epilepsy by means of a ketogenic diet. While PUFAs have been reported to influence the gating mechanism by electrostatic interactions to the voltage-sensor domain (VSD), the exact PUFA-protein interactions are still elusive. In this study, we report on the interactions between the Shaker K_V channel in open and closed states and a PUFA-enriched lipid bilayer using microsecond molecular dynamics simulations. We determined a putative PUFA binding site in the open state of the channel located at the protein-lipid interface in the vicinity of the extracellular halves of the S3 and S4 helices of the VSD. In particular, the lipophilic PUFA tail covered a wide range of non-specific hydrophobic interactions in the hydrophobic central core of the protein-lipid interface, while the carboxylic head group displayed more specific interactions to polar/charged residues at the extracellular regions of the S3 and S4 helices, encompassing the S3-S4 linker. Moreover, by studying the interactions between saturated fatty acids (SFA) and the Shaker K_V channel, our study confirmed an increased conformational flexibility in the polyunsaturated carbon tails compared to saturated carbon chains, which may explain the specificity of PUFA action on channel proteins.     

Molecular dynamics simulation of human immunodeficiency virus protein U (Vpu) in lipid/water Langmuir monolayer

Virus protein U (Vpu) is an accessory membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation and that the cytoplasmic domain adopts the helix-loop-helix-turn motif. This 3.5-ns molecular dynamics (MD) simulation of Vpu in a lipid/membrane environment has fully reproduced these structural characteristics. Membrane propensities of two amphipathic helices in the cytoplasmic domain are further compared here to understand better their complicated orientational behavior known from experiment. This study first reveals that the highly conserved loop region in the cytoplasmic domain can be closely associated with the membrane surface. It is known from the simulation that Vpu is associated with 34 lipids in this Langmuir monolayer. The lipids that are located between the Vpu transmembrane helix and the first helix in the cytoplasmic domain are pushed up by Vpu. These elevated lipids have increased P-N tilt angles for the head groups but unchanged acyl-chain tilt angles compared with lipids that do not interact with Vpu. This study verifies the significance of applying MD simulation in refining protein structure and revealing detailed protein-lipid interaction in membrane/water environment. Figure XZ view of a snapshot of Vpu/DLGPC/water system after 3.5 ns NP(N)gamma T MD simulation. Coloring scheme: Vpu, red; C, green; H, pink; N, blue; O, orange; P, magenta; water, light blue     

How protein transmembrane segments sense the lipid environment

Integral membrane proteins have central roles in a vast number of vital cellular processes. A structural feature that most membrane proteins have in common is the presence of one or more alpha-helices with which they traverse the lipid bilayer. Because of the interaction with the surrounding lipids, the organization of these transmembrane helices will be sensitive to lipid properties like lateral packing, hydrophobic thickness, and headgroup charge. The helices may adapt to the lipids in different ways, which in turn can influence the structure and function of the intact membrane protein. In this review, we will focus on how the lipid environment influences two specific properties of transmembrane segments: their lateral association and their tilt with respect to the bilayer normal.     

跨膜蛋白CD3ε与磷脂酰肌醇二磷酸在脂筏域体系中蛋白质-脂质相互作用的粗粒化模拟研究

细胞膜局部区域可形成富含饱和脂质、胆固醇、鞘脂的脂筏域作为其信号转导调控平台。传统实验手段在研究脂筏及其功能时受到系统复杂度高及区域结构瞬时性强等制约。近年来,分子动力学模拟技术为细胞膜的组织原则提供了重要的理论支撑,从简单的单一组分模型到多组分系统转变,最终形成了越来越多的细胞膜仿真模型。其中,粗粒化模拟由于其简化模型,可大副拓展模拟体系的复杂程度与模拟时间,在细胞膜以及蛋白质-脂质相互作用相关研究中得到了广泛应用。本文采用MARTINI粗粒化力场模拟,构建了一种含有阴离子脂质磷脂酰肌醇二磷酸(phosphatidylinositol diphosphate, PIP2)的混合膜体系。模拟结果表明,该体系在适当温度及饱和度条件下,能自发分层形成脂筏域;膜厚度、膜组分分布、膜组分流动性等多种参数均表明,脂筏结构形成且符合其结构特征;少量PIP2添加不影响分层特性且PIP2对脂筏具有显著亲和性。此外,利用该模型以跨膜信号蛋白CD3ε为例研究了脂筏域体系中蛋白质-脂质相互作用。结果表明,PIP2-CD3ε胞内区相互作用可能是脂筏招募CD3ε的驱动力,且该过程可受钙离子调控。本工作体现了粗粒化模拟在仿真膜相关研究中的巨大优势及良好应用前景。