CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

CGDB: a database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

Defining the Role of the Tension Sensor in the Mechanosensitive Channel of Small Conductance

Mutations that alter the phenotypic behavior of the Escherichia coli mechanosensitive channel of small conductance (MscS) have been identified; however, most of these residues play critical roles in the transition between the closed and open states of the channel and are not directly involved in lipid interactions that transduce the tension response. In this study, we use molecular dynamic simulations to predict critical lipid interacting residues in the closed state of MscS. The physiological role of these residues was then investigated by performing osmotic downshock assays on MscS mutants where the lipid interacting residues were mutated to alanine. These experiments identified seven residues in the first and second transmembrane helices as lipid-sensing residues. The majority of these residues are hydrophobic amino acids located near the extracellular interface of the membrane. All of these residues interact strongly with the lipid bilayer in the closed state of MscS, but do not face the bilayer directly in structures associated with the open and desensitized states of the channel. Thus, the position of these residues relative to the lipid membrane appears related to the ability of the channel to sense tension in its different physiological states.     

Structure of the Mechanosensitive Channel MscS Embedded in the Membrane Bilayer

Since life has emerged, gradients of osmolytes over the cell membrane cause pressure changes in the cell and require tight regulation to prevent cell rupture. The mechanosensitive channel of small conductance (MscS) releases solutes and water when a hypo-osmotic shock raises the pressure in the cell. It is a member of a large family of MscS-like channels found in bacteria, archaea, fungi and plants and model for mechanosensation. MscS senses the increase of tension in the membrane directly by the force from the lipids, but the molecular mechanism is still elusive. We determined the lipid interactions of MscS by resolving the structure of Escherichia coli MscS embedded in membrane discs to 2.9-Å resolution using cryo-electron microscopy. The membrane is attached only to parts of the sensor paddles of MscS, but phospholipid molecules move through grooves into remote pockets on the cytosolic side. On the periplasmic side, a lipid bound by R88 at the pore entrance is separated from the membrane by TM1 helices. The N-terminus interacts with the periplasmic membrane surface. We demonstrate that the unique membrane domain of MscS promotes deep penetration of lipid molecules and shows multimodal interaction with the membrane to fine-tune tension sensing.     

Lipid-protein interaction of the MscS mechanosensitive channel examined by scanning mutagenesis

The mechanosensitive channel of small conductance (MscS) is a bacterial mechanosensitive channel that opens in response to rapid hypoosmotic stress. Since MscS can be opened solely by membrane stretch without help from any accessory protein, the lipid-protein interface must play a crucial role in sensing membrane tension. In this study, the hydrophobic residues in the lipid-protein interface were substituted one by one with a hydrophilic amino acid, asparagine, to modify the interaction between the protein and the lipid. Function of the mutant MscSs was examined by patch-clamp and hypoosmotic shock experiments. An increase in the gating threshold and a decrease in the viability on hypoosmotic shock were observed when the hydrophobic residues near either end of the first or the second transmembrane helix (TM1 or TM2) were replaced with asparagine. This observation indicates that the lipid-protein interaction at the ends of both helices (TM1 and TM2) is essential to MscS function.     

Three-dimensional architecture of membrane-embedded MscS in the closed conformation

The mechanosensitive channel of small conductance (MscS) is part of a coordinated response to osmotic challenges in Escherichia coli. MscS opens as a result of membrane tension changes, thereby releasing small solutes and effectively acting as an osmotic safety valve. Both the functional state depicted by its crystal structure and its gating mechanism remain unclear. Here, we combine site-directed spin labeling, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations with novel energy restraints based on experimental electron paramagnetic resonance data to investigate the native transmembrane (TM) and periplasmic molecular architecture of closed MscS in a lipid bilayer. In the closed conformation, MscS shows a more compact TM domain than in the crystal structure, characterized by a realignment of the TM segments towards the normal of the membrane. The previously unresolved NH₂-terminus forms a short helical hairpin capping the extracellular ends of TM1 and TM2 and is in close interaction with the bilayer interface. The present three-dimensional model of membrane-embedded MscS in the closed state represents a key step in determining the molecular mechanism of MscS gating.     

Gap junction hemichannel interactions with zwitterionic lipid, anionic lipid, and cholesterol: molecular simulation studies

Interactions with membrane lipids can exert dramatic functional consequences on gap junction proteins. Recent experimental work has highlighted the importance of anionic lipids and cholesterol in facilitating channel activity. In this work, we have employed a coarse-grained molecular model in conjunction with molecular dynamics (MD) simulations to study the interactions between a connexin 26 (Cx26) hemichannel and a number of lipid species, including palmitoyloleoylphosphatidylcholine (POPC), anionic palmitoyloleoylphosphatidic acid (POPA), and cholesterol, in order to identify sites at the protein interface which may exhibit preferential, specific binding to these lipids, as well as determine the characteristics of these interactions. We have also employed an atomistic model of Cx26 embedded in a mixed PA/PC bilayer as a comparison and to elucidate further lipid-protein interactions. Our simulation results suggest enrichment of interfacial PA at the intracellular leaflet at high bulk PA concentrations. PC can form tight binding interactions with the hemichannel, particularly at intersubunit crevices (classical nonannular sites). In mixed bilayers, however, POPA competes with POPC for these sites, displacing the latter in some cases. While the residues responsible for interactions with PC and PA are similar, the latter exhibits a unique property of being capable of forming stable hydrophilic contacts with multiple residues spanning two different adjacent subunits at both leaflets of the bilayer, as opposed to POPC which can only do so at the extracellular side. These results suggest that POPA may be essential to channel function by acting as an intersubunit lipid bridge. Additionally, we propose that the compositional enrichment of POPA at the Cx26 interface may serve important roles in voltage gating. Simulation of a mixed POPC:cholesterol bilayer suggests that the hemichannel enhances the transbilayer mobility of vicinal cholesterols, increasing the likelihood of site-hopping and interleaflet flip-flop transitions.     

Lipid-protein interactions as determinants of membrane protein structure and function

To determine how the lipid environment affects membrane protein structure and function, strains of Escherichia coli were developed in which normal phospholipid composition can be altered or foreign lipids can be introduced. The properties of LacY (lactose permease) were investigated as a function of lipid environment. Assembly of LacY in membranes lacking PE (phosphatidylethanolamine) results in misorientation of the N-terminal six-TM (transmembrane domain) helical bundle with loss of energy-dependent uphill transport and retention of energy-independent downhill transport. Post-assembly introduction of PE results in nearly native orientation of TMs and restoration of uphill transport. Foreign lipids with no net charge can substitute for PE in supporting native LacY topology, but restoration of uphill transport is dependent on native topology and the proper folding of a solvent-exposed domain. Increasing the positive charge density of the cytoplasmically exposed surface of LacY counters TM misorientation in the absence of neutral lipids, demonstrating that charge interactions between these domains and the surface of the membrane bilayer are determinants of TM orientation. Therefore membrane protein organization or reorganization is determined either during initial assembly or post-insertionally through direct interactions between the protein and the lipid environment, which affects the topogenic potency of opposing charged residues as topological signals independent of the translocon.     

The role of tryptophan residues in the function and stability of the mechanosensitive channel MscS from Escherichia coli

Tryptophan (Trp) residues play important roles in many proteins. In particular they are enriched in protein surfaces involved in protein docking and are often found in membrane proteins close to the lipid head groups. However, they are usually absent from the membrane domains of mechanosensitive channels. Three Trp residues occur naturally in the Escherichia coli MscS (MscS-Ec) protein: W16 lies in the periplasm, immediately before the first transmembrane span (TM1), whereas W240 and W251 lie at the subunit interfaces that create the cytoplasmic vestibule portals. The role of these residues in MscS function and stability were investigated using site-directed mutagenesis. Functional channels with altered properties were created when any of the Trp residues were replaced by another amino acid, with the greatest retention of function associated with phenylalanine (Phe) substitutions. Analysis of the fluorescence properties of purified mutant MscS proteins containing single Trp residues revealed that W16 and W251 are relatively inaccessible, whereas W240 is accessible to quenching agents. The data point to a significant role for W16 in the gating of MscS, and an essential role for W240 in MscS oligomer stability.     

Structural Interactions between Lipids, Water and S1-S4 Voltage-Sensing Domains

Membrane proteins serve crucial signaling and transport functions, yet relatively little is known about their structures in membrane environments or how lipids interact with these proteins. For voltage-activated ion channels, X-ray structures suggest that the mobile voltage-sensing S4 helix would be exposed to the membrane, and functional studies reveal that lipid modification can profoundly alter channel activity. Here, we use solid-state NMR to investigate structural interactions of lipids and water with S1-S4 voltage-sensing domains and to explore whether lipids influence the structure of the protein. Our results demonstrate that S1-S4 domains exhibit extensive interactions with lipids and that these domains are heavily hydrated when embedded in a membrane. We also find evidence for preferential interactions of anionic lipids with S1-S4 domains and that these interactions have lifetimes on the timescale of ≤ 10^− 3 s. Arg residues within S1-S4 domains are well hydrated and are positioned in close proximity to lipids, exhibiting local interactions with both lipid headgroups and acyl chains. Comparative studies with a positively charged lipid lacking a phosphodiester group reveal that this lipid modification has only modest effects on the structure and hydration of S1-S4 domains. Taken together, our results demonstrate that Arg residues in S1-S4 voltage-sensing domains reside in close proximity to the hydrophobic interior of the membrane yet are well hydrated, a requirement for carrying charge and driving protein motions in response to changes in membrane voltage.     

Molecular dynamics study of gating in the mechanosensitive channel of small conductance MscS

Mechanosensitive channels are a class of ubiquitous membrane proteins gated by mechanical strain in the cellular membrane. MscS, the mechanosensitive channel of small conductance, is found in the inner membrane of Escherichia coli and its crystallographic structure in an open form has been recently solved. By means of molecular dynamics simulations we studied the stability of the channel conformation suggested by crystallography in a fully solvated lipid (POPC) bilayer, the combined system encompassing 224,340 atoms. When restraining the backbone of the protein, the channel remained in the open form and the simulation revealed intermittent permeation of water molecules through the channel. Abolishing the restraints under constant pressure conditions led to spontaneous closure of the transmembrane channel, whereas abolishing the restraints when surface tension (20 dyn/cm) was applied led to channel widening. The large balloon-shaped cytoplasmic domain of MscS exhibited spontaneous diffusion of ions through its side openings. Interaction between the transmembrane domain and the cytoplasmic domain of MscS was observed and involved formation of salt bridges between residues Asp62 and Arg128; this interaction may be essential for the gating of MscS. K+ and Cl- ions showed distinctively different distributions in and around the channel.     

Identification of the hydrophobic thickness of a membrane protein using fluorescence spectroscopy: studies with the mechanosensitive channel MscL

The hydrophobic thickness of a membrane protein is an important parameter, defining how the protein sits within the hydrocarbon core of the lipid bilayer that surrounds it in a membrane. Here we show that Trp scanning mutagenesis combined with fluorescence spectroscopy can be used to define the hydrophobic thickness of a membrane protein. The mechanosensitive channel of large conductance (MscL) contains two transmembrane alpha-helices, of which the second (TM2) is lipid-exposed. The region of TM2 that spans the hydrocarbon core of the bilayer when MscL is reconstituted into bilayers of dioleoylphosphatidylcholine runs from Leu-69 to Leu-92, giving a hydrophobic thickness of ca. 25 A. The results obtained using Trp scanning mutagenesis were confirmed using Cys residues labeled with the N-methyl-amino-7-nitroben-2-oxa-1,3-diazole [NBD] group; both fluorescence emission maxima and fluorescence lifetimes for the NBD group are sensitive to solvent dielectric constant over the range (2-40) thought to span the lipid headgroup region of a lipid bilayer. Changing phospholipid fatty acyl chain lengths from C14 and C24 results in no significant change for the fluorescence of the interfacial residues, suggesting very efficient hydrophobic matching between the protein and the surrounding lipid bilayer.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

Manipulating the folding of membrane proteins: using the bilayer to our advantage

The folding mechanisms of integral membrane proteins have largely eluded detailed study. This is owing to the inherent difficulties in folding these hydrophobic proteins in vitro, which, in turn, reflects the often apparently insurmountable problem of mimicking the natural membrane bilayer with lipid or detergent mixtures. There is, however, a large body of information on lipid properties and, in particular, on phosphatidylcholine and phosphatidylethanolamine lipids, which are common to many biological membranes. We have exploited this knowledge to develop efficient in vitro lipid-bilayer folding systems for the membrane protein, bacteriorhodopsin. Furthermore, we have shown that a rate-limiting apoprotein folding step and the overall folding efficiency appear to be controlled by particular properties of the lipid bilayer. The properties of interest are the stored curvature elastic energy within the bilayer, and the lateral pressure that the lipid chains exert on the their neighbouring folding proteins. These are generic properties of the bilayer that can be achieved with simple mixtures of biological lipids, and are not specific to the lipids studied here. These bilayer properties also seem to be important in modulating the function of several membrane proteins, as well as the function of membranes in vivo. Thus, it seems likely that careful manipulations of lipid properties will shed light on the forces that drive membrane protein folding, and will aid the development of bilayer folding systems for other membrane proteins.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

Lipids in the structure,folding, and function of the KcsA K+ channel

Lipid molecules surround an ion channel in its native environment of cellular membranes. The importance of the lipid bilayer and the role of lipid protein interactions in ion channel structure and function are not well understood. Here we demonstrate that the bacterial potassium channel KcsA binds a negatively charged lipid molecule. We have defined the potential binding site of the lipid molecule on KcsA by X-ray crystallographic analysis of a complex of KcsA with a monoclonal antibody Fab fragment. We also demonstrate that lipids are required for the in vitro refolding of the KcsA tetramer from the unfolded monomeric state. The correct refolding of the KcsA tetramer requires lipids, but it is not dependent on negatively charged lipids as refolding takes place in the absence of such lipids. We confirm that the presence of negatively charged lipids is required for ion conduction through the KcsA potassium channel, suggesting that the lipid bound to KcsA is important for ion channel function.     

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.     

Lipid-protein interactions

Intrinsic membrane proteins are solvated by a shell of lipid molecules interacting with the membrane-penetrating surface of the protein; these lipid molecules are referred to as annular lipids. Lipid molecules are also found bound between transmembrane α-helices; these are referred to as non-annular lipids. Annular lipid binding constants depend on fatty acyl chain length, but the dependence is less than expected from models based on distortion of the lipid bilayer alone. This suggests that hydrophobic matching between a membrane protein and the surrounding lipid bilayer involves some distortion of the transmembrane α-helical bundle found in most membrane proteins, explaining the importance of bilayer thickness for membrane protein function. Annular lipid binding constants also depend on the structure of the polar headgroup region of the lipid, and hotspots for binding anionic lipids have been detected on some membrane proteins; binding of anionic lipid molecules to these hotspots can be functionally important. Binding of anionic lipids to non-annular sites on membrane proteins such as the potassium channel KcsA can also be important for function. It is argued that the packing preferences of the membrane-spanning α-helices in a membrane protein result in a structure that matches nicely with that of the surrounding lipid bilayer, so that lipid and protein can meet without either having to change very much.     

Structural and Functional Profiling of the Lateral Gate of the Sec61 Translocon

The evolutionarily conserved Sec61 translocon mediates the translocation and membrane insertion of proteins. For the integration of proteins into the membrane, the Sec61 translocon opens laterally to the lipid bilayer. Previous studies suggest that the lateral opening of the channel is mediated by the helices TM2b and TM7 of a pore-forming subunit of the Sec61 translocon. To map key residues in TM2b and TM7 in yeast Sec61 that modulate lateral gating activity, we performed alanine scanning and in vivo site-directed photocross-linking experiments. Alanine scanning identified two groups of critical residues in the lateral gate, one group that leads to defects in the translocation and membrane insertion of proteins and the other group that causes faster translocation and facilitates membrane insertion. Photocross-linking data show that the former group of residues is located at the interface of the lateral gate. Furthermore, different degrees of defects for the membrane insertion of single- and double-spanning membrane proteins were observed depending on whether the mutations were located in TM2b or TM7. These results demonstrate subtle differences in the molecular mechanism of the signal sequence binding/opening of the lateral gate and membrane insertion of a succeeding transmembrane segment in a polytopic membrane protein.