Lateral pressure profile, spontaneous curvature frustration, and the incorporation and conformation of proteins in membranes

Lipid-protein interactions are an important determinant of the stability and function of integral and transmembrane proteins. In addition to local interactions at the lipid-protein interface, global interactions such as the distribution of internal lateral pressure may also influence protein conformation. It is shown here that the effects of the membrane lateral pressure profile on the conformation or insertion of proteins in membranes are equivalent to the elastic response to the frustrated spontaneous curvature, c_o, of the component lipid monolayer leaflets. The chemical potential of the protein in the membrane is predicted to depend linearly on the spontaneous curvature of the lipid leaflets, just as does the contribution of the protein to the elastic bending energy of the lipid, and to be independent of the hydrophobic tension, γ_phob, at the lipid-water interface. Analysis of the dependence of protein partitioning or conformational transitions on spontaneous curvature of the constituent lipids gives an experimental estimate for the cross-sectional intramembrane shape of the protein or its difference between conformations. Values in the region of 50-110 Å² are estimated for the effective cross-sectional shape changes on the insertion and conductance transitions of alamethicin, and on the activation of CTP:phosphocholine cytidylyltransferase or rhodopsin in lipid membranes. Much larger values are estimated for the mechanosensitive channel, MscL. Values for the change in intramembrane shape may also be used, together with determinations of lipid relative association constants, to estimate contributions of direct lipid-protein interactions to the lateral pressure experienced by the protein. Changes in chemical potential ∼12 kJ mol⁻¹ can be estimated for radial changes of 1 Å in a protein of diameter 40 Å.     

Protein shape change has a major effect on the gating energy of a mechanosensitive channel

Increasing experimental evidence has shown that membrane protein functionality depends on molecular composition of cell membranes. However, the origin of this dependence is not fully understood. It is reasonable to assume that specific lipid-protein interactions are important, yet more generic effects due to mechanical properties of lipid bilayers likely play a significant role too. Previously it has been demonstrated using models for elastic properties of membranes and lateral pressure profiles of lipid bilayers that the mechanical properties of a lipid bilayer can contribute as much as ∼10 k_BT to the free energy difference associated with a change in protein conformational state. Here, we extend those previous approaches to a more realistic model for a large mechanosensitive channel (MscL). We use molecular dynamics together with the MARTINI model to simulate the open and closed states of MscL embedded in a DOPC bilayer. We introduce a procedure to calculate the mechanical energy change in the channel gating using a three-dimensional pressure distribution inside a membrane, computed from the molecular dynamics simulations. We decompose the mechanical energy to terms associated with area dilation and shape contribution. Our results highlight that the lateral pressure profile of a lipid bilayer together with the shape change in gating can induce a contribution of ∼30 k_BT on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.     

Spring constants for channel-induced lipid bilayer deformations. Estimates using gramicidin channels.

Hydrophobic interactions between a bilayer and its embedded membrane proteins couple protein conformational changes to changes in the packing of the surrounding lipids. The energetic cost of a protein conformational change therefore includes a contribution from the associated bilayer deformation energy (DeltaGdef0), which provides a mechanism for how membrane protein function depends on the bilayer material properties. Theoretical studies based on an elastic liquid-crystal model of the bilayer deformation show that DeltaGdef0 should be quantifiable by a phenomenological linear spring model, in which the bilayer mechanical characteristics are lumped into a single spring constant. The spring constant scales with the protein radius, meaning that one can use suitable reporter proteins for in situ measurements of the spring constant and thereby evaluate quantitatively the DeltaGdef0 associated with protein conformational changes. Gramicidin channels can be used as such reporter proteins because the channels form by the transmembrane assembly of two nonconducting monomers. The monomerleft arrow over right arrow dimer reaction thus constitutes a well characterized conformational transition, and it should be possible to determine the phenomenological spring constant describing the channel-induced bilayer deformation by examining how DeltaGdef0 varies as a function of a mismatch between the hydrophobic channel length and the unperturbed bilayer thickness. We show this is possible by analyzing experimental studies on the relation between bilayer thickness and gramicidin channel duration. The spring constant in nominally hydrocarbon-free bilayers agrees well with estimates based on a continuum analysis of inclusion-induced bilayer deformations using independently measured material constants.     

Regulation of sodium channel function by bilayer elasticity: the importance of hydrophobic coupling. Effects of Micelle-forming amphiphiles and cholesterol

Membrane proteins are regulated by the lipid bilayer composition. Specific lipid-protein interactions rarely are involved, which suggests that the regulation is due to changes in some general bilayer property (or properties). The hydrophobic coupling between a membrane-spanning protein and the surrounding bilayer means that protein conformational changes may be associated with a reversible, local bilayer deformation. Lipid bilayers are elastic bodies, and the energetic cost of the bilayer deformation contributes to the total energetic cost of the protein conformational change. The energetics and kinetics of the protein conformational changes therefore will be regulated by the bilayer elasticity, which is determined by the lipid composition. This hydrophobic coupling mechanism has been studied extensively in gramicidin channels, where the channel-bilayer hydrophobic interactions link a "conformational" change (the monomer<-->dimer transition) to an elastic bilayer deformation. Gramicidin channels thus are regulated by the lipid bilayer elastic properties (thickness, monolayer equilibrium curvature, and compression and bending moduli). To investigate whether this hydrophobic coupling mechanism could be a general mechanism regulating membrane protein function, we examined whether voltage-dependent skeletal-muscle sodium channels, expressed in HEK293 cells, are regulated by bilayer elasticity, as monitored using gramicidin A (gA) channels. Nonphysiological amphiphiles (beta-octyl-glucoside, Genapol X-100, Triton X-100, and reduced Triton X-100) that make lipid bilayers less "stiff", as measured using gA channels, shift the voltage dependence of sodium channel inactivation toward more hyperpolarized potentials. At low amphiphile concentration, the magnitude of the shift is linearly correlated to the change in gA channel lifetime. Cholesterol-depletion, which also reduces bilayer stiffness, causes a similar shift in sodium channel inactivation. These results provide strong support for the notion that bilayer-protein hydrophobic coupling allows the bilayer elastic properties to regulate membrane protein function.     

Conformational exchange of fatty acid binding protein induced by protein-nanodisc interactions

Fatty acid binding proteins (FABPs) can facilitate the transfer of long-chain fatty acids between intracellular membranes across considerable distances. The transfer process involves fatty acids, their donor membrane and acceptor membrane, and FABPs, implying that potential protein-membrane interactions exist. Despite intensive studies on FABP-membrane interactions, the interaction mode remains elusive, and the protein-membrane association and dissociation rates are inconsistent. In this study, we used nanodiscs (NDs) as mimetic membranes to investigate FABP-membrane interactions. Our NMR experiments showed that human intestinal FABP interacts weakly with both negatively charged and neutral membranes, but it prefers the negatively charged one. Through simultaneous analysis of NMR relaxation in the rotating-frame (R_1ρ), relaxation dispersion, chemical exchange saturation transfer, and dark-state exchange saturation transfer data, we estimated the affinity of the protein to negatively charged NDs, the dissociation rate, and apparent association rate. We further showed that the protein in the ND-bound state adopts a conformation different from the native structure and the second helix is very likely involved in interactions with NDs. We also found a membrane-induced FABP conformational state that exists only in the presence of NDs. This state is native-like, different from other conformational states in structure, unbound to NDs, and in dynamic equilibrium with the ND-bound state.     

Energetics of inclusion-induced bilayer deformations.

The material properties of lipid bilayers can affect membrane protein function whenever conformational changes in the membrane-spanning proteins perturb the structure of the surrounding bilayer. This coupling between the protein and the bilayer arises from hydrophobic interactions between the protein and the bilayer. We analyze the free energy cost associated with a hydrophobic mismatch, i.e., a difference between the length of the protein's hydrophobic exterior surface and the average thickness of the bilayer's hydrophobic core, using a (liquid-crystal) elastic model of bilayer deformations. The free energy of the deformation is described as the sum of three contributions: compression-expansion, splay-distortion, and surface tension. When evaluating the interdependence among the energy components, one modulus renormalizes the other: e.g., a change in the compression-expansion modulus affects not only the compression-expansion energy but also the splay-distortion energy. The surface tension contribution always is negligible in thin solvent-free bilayers. When evaluating the energy per unit distance (away from the inclusion), the splay-distortion component dominates close to the bilayer/inclusion boundary, whereas the compression-expansion component is more prominent further away from the boundary. Despite this complexity, the bilayer deformation energy in many cases can be described by a linear spring formalism. The results show that, for a protein embedded in a membrane with an initial hydrophobic mismatch of only 1 A, an increase in hydrophobic mismatch to 1.3 A can increase the Boltzmann factor (the equilibrium distribution for protein conformation) 10-fold due to the elastic properties of the bilayer.     

In vitro synthesis and oligomerization of the mechanosensitive channel of large conductance, MscL, into a functional ion channel

Elucidation of high-resolution structures of integral membrane proteins is drastically lagging behind that of cytoplasmic proteins. In vitro synthesis and insertion of membrane proteins into synthetic membranes could circumvent bottlenecks associated with the overexpression of membrane proteins, producing sufficient membrane-inserted, correctly folded protein for structural studies. Using the mechanosensitive channel of large conductance, MscL, as a model protein we show that in vitro synthesized MscL inserts into YidC-containing proteoliposomes and oligomerizes to form a homopentamer. Using planar membrane bilayers, we show that MscL forms functional ion channels capable of ion transport. These data demonstrate that membrane insertion of MscL is YidC mediated, whereas subsequent oligomerization into a functional homopentamer is a spontaneous event.     

Unwinding of the Substrate Transmembrane Helix in Intramembrane Proteolysis

Intramembrane-cleaving proteases (I-CLiPs) activate pools of single-pass helical membrane protein signaling precursors that are key in the physiology of prokaryotic and eukaryotic cells. Proteases typically cleave peptide bonds within extended or flexible regions of their substrates, and thus the mechanism underlying the ability of I-CLiPs to hydrolyze the presumably α-helical transmembrane domain (TMD) of these membrane proteins is unclear. Using deep-ultraviolet resonance Raman spectroscopy in combination with isotopic labeling, we show that although predominantly in canonical α-helical conformation, the TMD of the established I-CLiP substrate Gurken displays 3₁₀-helical geometry. As measured by microscale thermophoresis, this substrate binds with high affinity to the I-CLiPs GlpG rhomboid and MCMJR1 presenilin homolog in detergent micelles. Binding results in deep-ultraviolet resonance Raman spectra, indicating conformational changes consistent with unwinding of the 3₁₀-helical region of the substrate’s TMD. This 3₁₀-helical conformation is key for intramembrane proteolysis, as the substitution of a single proline residue in the TMD of Gurken by alanine suppresses 3₁₀-helical content in favor of α-helical geometry and abolishes cleavage without affecting binding to the I-CLiP. Complemented by molecular dynamics simulations of the TMD of Gurken, our vibrational spectroscopy data provide biophysical evidence in support of a model in which the transmembrane region of cleavable I-CLiP substrates displays local deviations in canonical α-helical conformation characterized by chain flexibility, and binding to the enzyme results in conformational changes that facilitate local unwinding of the transmembrane helix for cleavage.     

Force Transduction and Lipid Binding in MscL: A Continuum-Molecular Approach

The bacterial mechanosensitive channel MscL, a small protein mainly activated by membrane tension, is a central model system to study the transduction of mechanical stimuli into chemical signals. Mutagenic studies suggest that MscL gating strongly depends on both intra-protein and interfacial lipid-protein interactions. However, there is a gap between this detailed chemical information and current mechanical models of MscL gating. Here, we investigate the MscL bilayer-protein interface through molecular dynamics simulations, and take a combined continuum-molecular approach to connect chemistry and mechanics. We quantify the effect of membrane tension on the forces acting on the surface of the channel, and identify interactions that may be critical in the force transduction between the membrane and MscL. We find that the local stress distribution on the protein surface is largely asymmetric, particularly under tension, with the cytoplasmic side showing significantly larger and more localized forces, which pull the protein radially outward. The molecular interactions that mediate this behavior arise from hydrogen bonds between the electronegative oxygens in the lipid headgroup and a cluster of positively charged lysine residues on the amphipathic S1 domain and the C-terminal end of the second trans-membrane helix. We take advantage of this strong interaction (estimated to be 10–13 kT per lipid) to actuate the channel (by applying forces on protein-bound lipids) and explore its sensitivity to the pulling magnitude and direction. We conclude by highlighting the simple motif that confers MscL with strong anchoring to the bilayer, and its presence in various integral membrane proteins including the human mechanosensitive channel K2P1 and bovine rhodopsin.     

A finite element framework for studying the mechanical response of macromolecules: application to the gating of the mechanosensitive channel MscL

The gating pathways of mechanosensitive channels of large conductance (MscL) in two bacteria (Mycobacterium tuberculosis and Escherichia coli) are studied using the finite element method. The phenomenological model treats transmembrane helices as elastic rods and the lipid membrane as an elastic sheet of finite thickness; the model is inspired by the crystal structure of MscL. The interactions between various continuum components are derived from molecular-mechanics energy calculations using the CHARMM all-atom force field. Both bacterial MscLs open fully upon in-plane tension in the membrane and the variation of pore diameter with membrane tension is found to be essentially linear. The estimated gating tension is close to the experimental value. The structural variations along the gating pathway are consistent with previous analyses based on structural models with experimental constraints and biased atomistic molecular-dynamics simulations. Upon membrane bending, neither MscL opens substantially, although there is notable and nonmonotonic variation in the pore radius. This emphasizes that the gating behavior of MscL depends critically on the form of the mechanical perturbation and reinforces the idea that the crucial gating parameter is lateral tension in the membrane rather than the curvature of the membrane. Compared to popular all-atom-based techniques such as targeted or steered molecular-dynamics simulations, the finite element method-based continuum-mechanics framework offers a unique alternative to bridge detailed intermolecular interactions and biological processes occurring at large spatial scales and long timescales. It is envisioned that such a hierarchical multiscale framework will find great value in the study of a variety of biological processes involving complex mechanical deformations such as muscle contraction and mechanotransduction.     

Membrane-Mediated Aggregation of Curvature-Inducing Nematogens and?Membrane Tubulation

The shapes of cell membranes are largely regulated by membrane-associated, curvature-active proteins. Herein, we use a numerical model of the membrane, recently developed by us, with elongated membrane inclusions possessing spontaneous directional curvatures that could be different along, and perpendicular to, the membrane’s long axis. We show that, due to membrane-mediated interactions, these curvature-inducing membrane-nematogens can aggregate spontaneously, even at low concentrations, and change the local shape of the membrane. We demonstrate that for a large group of such inclusions, where the two spontaneous curvatures have equal sign, the tubular conformation and sometimes the sheet conformation of the membrane are the common equilibrium shapes. We elucidate the factors necessary for the formation of these protein lattices. Furthermore, the elastic properties of the tubes, such as their compressional stiffness and persistence length, are calculated. Finally, we discuss the possible role of nematic disclination in capping and branching of the tubular membranes.     

Solid-state NMR spectroscopy to study protein–lipid interactions

The appropriate lipid environment is crucial for the proper function of membrane proteins. There is a tremendous variety of lipid molecules in the membrane and so far it is often unclear which component of the lipid matrix is essential for the function of a respective protein. Lipid molecules and proteins mutually influence each other; parameters such as acyl chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation are strongly modulated by proteins. More recent data also indicates that the influence of proteins goes beyond a single annulus of next-neighbor boundary lipids. Therefore, a mesoscopic approach to membrane lipid–protein interactions in terms of elastic membrane deformations has been developed. Solid-state NMR has greatly contributed to the understanding of lipid–protein interactions and the modern view of biological membranes. Methods that detect the influence of proteins on the membrane as well as direct lipid–protein interactions have been developed and are reviewed here. Examples for solid-state NMR studies on the interaction of Ras proteins, the antimicrobial peptide protegrin-1, the G protein-coupled receptor rhodopsin, and the K⁺ channel KcsA are discussed. This article is part of a Special Issue entitled Tools to study lipid functions.     

Neuroligin-2 dependent conformational activation of collybistin reconstituted in supported hybrid membranes

The assembly of the postsynaptic transmitter sensing machinery at inhibitory nerve cell synapses requires the intimate interplay between cell adhesion proteins, scaffold and adaptor proteins, and γ-aminobutyric acid (GABA) or glycine receptors. We developed an in vitro membrane system to reconstitute this process, to identify the essential protein components, and to define their mechanism of action, with a specific focus on the mechanism by which the cytosolic C terminus of the synaptic cell adhesion protein Neuroligin-2 alters the conformation of the adaptor protein Collybistin-2 and thereby controls Collybistin-2-interactions with phosphoinositides (PtdInsPs) in the plasma membrane. Supported hybrid membranes doped with different PtdInsPs and 1,2-dioleoyl-sn-glycero-3-{[N-(5-amino-1-carboxypentyl)iminodiacetic acid]succinyl} nickel salt (DGS-NTA(Ni)) to allow for the specific adsorption of the His₆-tagged intracellular domain of Neuroligin-2 (His-cytNL2) were prepared on hydrophobically functionalized silicon dioxide substrates via vesicle spreading. Two different collybistin variants, the WT protein (CB2_SH3) and a mutant that adopts an intrinsically ‘open’ and activated conformation (CB2_{SH3/W24A-E262A}), were bound to supported membranes in the absence or presence of His-cytNL2. The corresponding binding data, obtained by reflectometric interference spectroscopy, show that the interaction of the C terminus of Neuroligin-2 with Collybistin-2 induces a conformational change in Collybistin-2 that promotes its interaction with distinct membrane PtdInsPs.     

Phospholipid membranes affect tertiary structure of the soluble cytochrome b5 heme-binding domain

The influence of charged phospholipid membranes on the conformational state of the water-soluble fragment of cytochrome b₅ has been investigated by a variety of techniques at neutral pH. The results of this work provide the first evidence that aqueous solutions with high phospholipid/protein molar ratios (pH 7.2) induce the cytochrome to undergo a structural transition from the native conformation to an intermediate state with molten-globule like properties that occur in the presence of an artificial membrane surface and that leads to binding of the protein to the membrane. At other phospholipid/protein ratios, equilibrium was observed between cytochrome free in solution and cytochrome bound to the surface of vesicles. Inhibition of protein binding to the vesicles with increasing ionic strength indicated for the most part an electrostatic contribution to the stability of cytochrome b₅vesicle interactions at pH 7.2. The possible physiological role of membrane-induced conformational change in the structure of cytochrome b₅ upon the interaction with its redox partners is discussed.     

Contributions of Gaussian curvature and nonconstant lipid volume to protein deformation of lipid bilayers

An elastic model for membrane deformations induced by integral membrane proteins is presented. An earlier theory is extended to account for nonvanishing saddle splay modulus within lipid monolayers and perturbations to lipid volume proximal to the protein. Analytical results are derived for the deformation profile surrounding a single cylindrical protein inclusion, which compare favorably to coarse-grained simulations over a range of protein sizes. Numerical results for multi-protein systems indicate that membrane-mediated interactions between inclusions are strongly affected by Gaussian curvature and display nonpairwise additivity. Implications for the aggregation of proteins are discussed.     

Taming membranes: functional immobilization of biological membranes in hydrogels

Single molecule studies on membrane proteins embedded in their native environment are hampered by the intrinsic difficulty of immobilizing elastic and sensitive biological membranes without interfering with protein activity. Here, we present hydrogels composed of nano-scaled fibers as a generally applicable tool to immobilize biological membrane vesicles of various size and lipid composition. Importantly, membrane proteins immobilized in the hydrogel as well as soluble proteins are fully active. The triggered opening of the mechanosensitive channel of large conductance (MscL) reconstituted in giant unilamellar vesicles (GUVs) was followed in time on single GUVs. Thus, kinetic studies of vectorial transport processes across biological membranes can be assessed on single, hydrogel immobilized, GUVs. Furthermore, protein translocation activity by the membrane embedded protein conducting channel of bacteria, SecYEG, in association with the soluble motor protein SecA was quantitatively assessed in bulk and at the single vesicle level in the hydrogel. This technique provides a new way to investigate membrane proteins in their native environment at the single molecule level by means of fluorescence microscopy.     

Probing the Transition State of the Allosteric Pathway of the Shaker Kv Channel Pore by Linear Free-Energy Relations

Long-range coupling between distant functional elements of proteins may rely on allosteric communication trajectories lying along the protein structure, as described in the case of the Shaker voltage-activated potassium (Kv) channel model allosteric system. Communication between the distant Kv channel activation and slow inactivation pore gates was suggested to be mediated by a network of local pairwise and higher-order interactions among the functionally unique residues that constitute the allosteric trajectory. The mechanism by which conformational changes propagate along the Kv channel allosteric trajectory to achieve pore opening, however, remains unclear. Such conformational changes may propagate in either a concerted or a sequential manner during the reaction coordinate of channel opening. Residue-level structural information on the transition state of channel gating is required to discriminate between these possibilities. Here, we combine patch-clamp electrophysiology recordings of Kv channel gating and analysis using linear free-energy relations, focusing on a select set of residues spanning the allosteric trajectory of the Kv channel pore. We show that all allosteric trajectory residues tested exhibit an open-like conformation in the transition state of channel opening, implying that coupling interactions occur along the trajectory break in a concerted manner upon moving from the closed to the open state. Energetic coupling between the Kv channel gates thus occurs in a concerted fashion in both the spatial and the temporal dimensions, strengthening the notion that such trajectories correspond to pathways of mechanical deformation along which conformational changes propagate.     

Binding and interactions of L-BABP to lipid membranes studied by molecular dynamic simulations

Chicken liver bile acid-binding protein (L-BABP) is a member of the fatty acid-binding proteins super family. The common fold is a β-barrel of ten strands capped with a short helix-loop-helix motif called portal region, which is involved in the uptake and release of non-polar ligands. Using multiple-run molecular dynamics simulations we studied the interactions of L-BABP with lipid membranes of anionic and zwitterionic phospholipids. The simulations were in agreement with our experimental observations regarding the electrostatic nature of the binding and the conformational changes of the protein in the membrane. We observed that L-BABP migrated from the initial position in the aqueous bulk phase to the interface of anionic lipid membranes and established contacts with the head groups of phospholipids through the side of the barrel that is opposite to the portal region. The conformational changes in the protein occurred simultaneously with the binding to the membrane. Remarkably, these conformational changes were observed in the portal region which is opposite to the zone where the protein binds directly to the lipids. The protein was oriented with its macrodipole aligned in the configuration of lowest energy within the electric field of the anionic membrane, which indicates the importance of the electrostatic interactions to determine the preferred orientation of the protein. We also identified this electric field as the driving force for the conformational change. For all the members of the fatty acid-binding protein family, the interactions with lipid membranes is a relevant process closely related to the uptake, release and transfer of the ligand. The observations presented here suggest that the ligand transfer might not necessarily occur through the domain that directly interacts with the lipid membrane. The interactions with the membrane electric field that determine orientation and conformational changes described here can also be relevant for other peripheral proteins.     

Effects of preparation method on reversible conformational changes induced by neutral salts in spin-labeled erythrocyte membranes

Addition of “chaotropic” neutral salts to spin-labeled erythrocyte membranes produced changes in membrane conformation proportional to the known activities of the salts: KSCN > CaCl₂ ≥ KNO₃ ? NaCl. The effects had a threshold of about 0.10 m and increased through 0.60 m for all salts. Above 0.6 m other changes due to dielectric effects and/or protein loss began.The chaotropic effects were reversible on removal of the chaotropic ion, for human (or bovine) membranes prepared with a buffer ionic strength of 0.02. For membranes prepared with a buffer ionic strength of 0.005 the changes were irreversible. Exposure to salt concentration greater than 1 m also produced irreversible changes, probably due to solubilization. The results appear to show that the water-soluble membrane proteins are essential to maintain the “native” membrane conformation.