Membrane hydrocarbon thickness modulates the dynamics of a membrane transport protein

Steric trapping strategy for studying the folding of helical membrane proteins

Elucidating the folding energy landscape of membrane proteins is essential to the understanding of the proteins' stabilizing forces, folding mechanisms, biogenesis, and quality control. This is not a trivial task because the reversible control of folding is inherently difficult in a lipid bilayer environment. Recently, novel methods have been developed, each of which has a unique strength in investigating specific aspects of membrane protein folding. Among such methods, steric trapping is a versatile strategy allowing a reversible control of membrane protein folding with minimal perturbation of native protein-water and protein-lipid interactions. In a nutshell, steric trapping exploits the coupling of spontaneous denaturation of a doubly biotinylated protein to the simultaneous binding of bulky monovalent streptavidin molecules. This strategy has been evolved to investigate key elements of membrane protein folding such as thermodynamic stability, spontaneous denaturation rates, conformational features of the denatured states, and cooperativity of stabilizing interactions. In this review, we describe the critical methodological advancement, limitation, and outlook of the steric trapping strategy.

Cardiolipin-Dependent Properties of Model Mitochondrial Membranes from Molecular Simulations

Cardiolipin is an anionic lipid found in the mitochondrial membranes of eukaryotes ranging from unicellular microorganisms to metazoans. This unique lipid contributes to various mitochondrial functions, including metabolism, mitochondrial membrane fusion and/or fission dynamics, and apoptosis. However, differences in cardiolipin content between the two mitochondrial membranes, as well as dynamic fluctuations in cardiolipin content in response to stimuli and cellular signaling events, raise questions about how cardiolipin concentration affects mitochondrial membrane structure and dynamics. Although cardiolipin’s structural and dynamic roles have been extensively studied in binary mixtures with other phospholipids, the biophysical properties of cardiolipin in higher number lipid mixtures are still not well resolved. Here, we used molecular dynamics simulations to investigate the cardiolipin-dependent properties of ternary lipid bilayer systems that mimic the major components of mitochondrial membranes. We found that changes to cardiolipin concentration only resulted in minor changes to bilayer structural features but that the lipid diffusion was significantly affected by those alterations. We also found that cardiolipin position along the bilayer surfaces correlated to negative curvature deflections, consistent with the induction of negative curvature stress in the membrane monolayers. This work contributes to a foundational understanding of the role of cardiolipin in altering the properties in ternary lipid mixtures composed of the major mitochondrial phospholipids, providing much-needed insights to help understand how cardiolipin concentration modulates the biophysical properties of mitochondrial membranes.

Phospholipid effects on SGLT1-mediated glucose transport in rabbit ileum brush border membrane vesicles

In small intestine, sodium-glucose cotransporter SGLT1 provides the main mechanism for sugar uptake. We investigated the effect of membrane phospholipids (PL) on this transport in rabbit ileal brush border membrane vesicles (BBMV). For this, PL of different charge, length, and saturation were incorporated into BBMV. Transport was measured related to (i) membrane surface charge (membrane-bound MC540 fluorescence), (ii) membrane thickness (PL incorporation of different acyl chain length), and (iii) membrane fluidity (r_12AS, fluorescence anisotropy of 12-AS). Compared to phosphatidylcholine (PC) carrying a neutral head group, inhibition of SGLT1 increased considerably with the acidic phosphatidic acid (PA) and phosphatidylinositol (PI) that increase membrane negative surface charge. The order of PL potency was PI>PA > PE = PS > PC. Inhibition by acidic PA-oleate was 5-times more effective than with neutral PE (phosphatidylethanolamine)-oleate. Lineweaver-Burk plot indicated uncompetitive inhibition of SGLT1 by PA. When membrane thickness was increased by neutral PC of varying acyl chain length, transport was increasingly inhibited by 16:1 PC to 22:1 PC. Even more pronounced inhibition was observed with mono-unsaturated instead of saturated acyl chains which increased membrane fluidity (indicated by decreased r_12AS). In conclusion, sodium-dependent glucose transport of rabbit ileal BBMV is modulated by (i) altered membrane surface charge, (ii) length of acyl chains via membrane thickness, and (iii) saturation of PL acyl chains altering membrane fluidity. Transport was attenuated by charged PL with longer and unsaturated acyl residues. Alterations of PL may provide a principle for attenuating dietary glucose uptake.

Role of the Native Outer-Membrane Environment on the Transporter BtuB

BtuB is a TonB-dependent transporter that permits the high-affinity binding and transport of cobalamin (CBL), or vitamin B₁₂, across the asymmetric outer membrane (OM) of Gram-negative bacteria. It has been shown that Ca²⁺ binding is necessary for high-affinity binding of CBL to BtuB, and earlier simulations suggested that calcium ions serve to stabilize key substrate-binding extracellular loops. However, those simulations did not account for the lipopolysaccharides in the OM. To illuminate the roles of both Ca²⁺ and lipopolysaccharides in protein functionality, we performed simulations of apo and Ca²⁺-loaded BtuB in symmetric and asymmetric bilayers. The simulations reveal that the oligosaccharides of LPS stabilize the extracellular loops to some degree, apparently obviating the need for Ca²⁺. However, it is shown that Ca²⁺ ions stabilize a key substrate-binding loop to an even greater degree, as well as reposition specific CBL-binding residues, bringing them closer to the organization found in the CBL-bound structure. These results indicate the importance of including realistic membrane models when simulating outer-membrane proteins.

Sensing membrane thickness: Lessons learned from cold stress

The lipid bilayer component of biological membranes is important for the distribution, organization, and function of bilayer spanning proteins. These physical barriers are subjected to bilayer perturbations. As a consequence, nature has evolved proteins that are able to sense changes in the bilayer properties and transform these lipid-mediated stimuli into intracellular signals. A structural feature that most signal-transducing membrane-embedded proteins have in common is one or more α-helices that traverse the lipid bilayer. Because of the interaction with the surrounding lipids, the organization of these transmembrane helices will be sensitive to membrane properties, like hydrophobic thickness. The helices may adapt to the lipids in different ways, which in turn can influence the structure and function of the intact membrane proteins. We review recent insights into the molecular basis of thermosensing via changes in membrane thickness and consider examples in which the hydrophobic matching can be demonstrated using reconstituted membrane systems. This article is part of a Special Issue entitled: The cellular lipid landscape edited by Tim P. Levine and Anant K. Menon.

Probing Structural Dynamics and Topology of the KCNE1 Membrane Protein in Lipid Bilayers via Site-Directed Spin Labeling and Electron Paramagnetic Resonance Spectroscopy

KCNE1 is a single transmembrane protein that modulates the function of voltage-gated potassium channels, including KCNQ1. Hereditary mutations in the genes encoding either protein can result in diseases such as congenital deafness, long QT syndrome, ventricular tachyarrhythmia, syncope, and sudden cardiac death. Despite the biological significance of KCNE1, the structure and dynamic properties of its physiologically relevant native membrane-bound state are not fully understood. In this study, the structural dynamics and topology of KCNE1 in bilayered lipid vesicles was investigated using site-directed spin labeling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy. A 53-residue nitroxide EPR scan of the KCNE1 protein sequence including all 27 residues of the transmembrane domain (45-71) and 26 residues of the N- and C-termini of KCNE1 in lipid bilayered vesicles was analyzed in terms of nitroxide side-chain motion. Continuous wave-EPR spectral line shape analysis indicated the nitroxide spin label side-chains located in the KCNE1 TMD are less mobile when compared to the extracellular region of KCNE1. The EPR data also revealed that the C-terminus of KCNE1 is more mobile when compared to the N-terminus. EPR power saturation experiments were performed on 41 sites including 18 residues previously proposed to reside in the transmembrane domain (TMD) and 23 residues of the N- and C-termini to determine the topology of KCNE1 with respect to the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) lipid bilayers. The results indicated that the transmembrane domain is indeed buried within the membrane, spanning the width of the lipid bilayer. Power saturation data also revealed that the extracellular region of KCNE1 is solvent-exposed with some of the portions partially or weakly interacting with the membrane surface. These results are consistent with the previously published solution NMR structure of KCNE1 in micelles.

A coiled coil switch mediates cold sensing by the thermosensory protein DesK

The thermosensor histidine kinase DesK from Bacillus subtilis senses changes in membrane fluidity initiating an adaptive response. Structural changes in DesK have been implicated in transmembrane signaling, but direct evidence is still lacking. On the basis of structure-guided mutagenesis, we now propose a mechanism of DesK-mediated signal sensing and transduction. The data indicate that stabilization/destabilization of a 2-helix coiled coil, which connects the transmembrane sensory domain of DesK to its cytosolic catalytic region, is crucial to control its signaling state. Computational modeling and simulations reveal couplings between protein, water and membrane mechanics. We propose that membrane thickening is the main driving force for signal sensing and that it acts by inducing helix stretching and rotation prompting an asymmetric kinase-competent state. Overall, the known structural changes of the sensor kinase, as well as further dynamic rearrangements that we now predict, consistently link structure determinants to activity modulation.

Characterizing solution surface loop conformational flexibility of the GM2 activator protein

GM2AP has a β-cup topology with numerous X-ray structures showing multiple conformations for some of the surface loops, revealing conformational flexibility that may be related to function, where function is defined as either membrane binding associated with ligand binding and extraction or interaction with other proteins. Here, site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy and molecular dynamic (MD) simulations are used to characterize the mobility and conformational flexibility of various structural regions of GM2AP. A series of 10 single cysteine amino acid substitutions were generated, and the constructs were chemically modified with the methanethiosulfonate spin label. Continuous wave (CW) EPR line shapes were obtained and subsequently simulated using the microscopic order macroscopic disorder (MOMD) program. Line shapes for sites that have multiple conformations in the X-ray structures required two spectral components, whereas spectra of the remaining sites were adequately fit with single-component parameters. For spin labeled sites L126C and I66C, spectra were acquired as a function of temperature, and simulations provided for the determination of thermodynamic parameters associated with conformational change. Binding to GM2 ligand did not alter the conformational flexibility of the loops, as evaluated by EPR and NMR spectroscopies. These results confirm that the conformational flexibility observed in the surface loops of GM2AP crystals is present in solution and that the exchange is slow on the EPR time scale (>ns). Furthermore, MD simulation results are presented and agree well with the conformational heterogeneity revealed by SDSL.

Membrane physical properties influence transmembrane helix formation

The pHLIP peptide has three states: (I) soluble in aqueous buffer, (II) bound to the bilayer surface at neutral pH, and (III) inserted as a transmembrane (TM) helix at acidic pH. The membrane insertion of pHLIP at low pH can be used to target the acidic tissues characteristic of different diseases, such as cancer. We find that the α-helix content of state II depends on lipid acyl chain length but not cholesterol, suggesting the helicity of the bound state may be controlled by the bilayer elastic bending modulus. Experiments with the P20G variant show the proline residue in pHLIP reduces the α-helix content of both states II and III. We also observe that the membrane insertion pKa is influenced by membrane physical properties, following a biphasic pattern similar to the membrane thickness optima observed for the function of eukaryotic membrane proteins. Because tumor cells exhibit altered membrane fluidity, we suggest this might influence pHLIP tumor targeting. We used a cell insertion assay to determine the pKa in live cells, observing that the properties in liposomes held in the more complex plasma membrane. Our results show that the formation of a TM helix is modulated by both the conformational propensities of the peptide and the physical properties of the bilayer. These results suggest a physical role for helix-membrane interactions in optimizing the function of more complex TM proteins.

Hydrophobic matching controls the tilt and stability of the dimeric platelet-derived growth factor receptor (PDGFR) β transmembrane segment

The platelet-derived growth factor receptor β is a member of the cell surface receptor tyrosine kinase family and dimerizes upon activation. We determined the structure of the transmembrane segment in dodecylphosphocholine micelles by liquid-state NMR and found that it forms a stable left-handed helical dimer. Solid-state NMR and oriented circular dichroism were used to measure the tilt angle of the helical segments in macroscopically aligned model membranes with different acyl chain lengths. Both methods showed that decreasing bilayer thickness (DEPC-POPC-DMPC) led to an increase in the helix tilt angle from 10° to 30° with respect to the bilayer normal. At the same time, reconstitution of the comparatively long hydrophobic segment became less effective, eventually resulting in complete protein aggregation in the short-chain lipid DLPC. Unrestrained molecular dynamics simulations of the dimer were carried out in explicit lipid bilayers (DEPC, POPC, DMPC, sphingomyelin), confirming the observed dependence of the helix tilt angle on bilayer thickness. Notably, molecular dynamics revealed that the left-handed dimer gets tilted en bloc, whereas conformational transitions to alternative (e.g. right-handed dimeric) states were not supported. The experimental data along with the simulation results demonstrate a pronounced interplay between the platelet-directed growth factor receptor β transmembrane segment and the bilayer thickness. The effect of hydrophobic mismatch might play a key role in the redistribution and activation of the receptor within different lipid microdomains of the plasma membrane in vivo.

Structural and dynamical analysis of an engineered FhuA channel protein embedded into a lipid bilayer or a detergent belt

Engineered channel proteins are promising nano-components with applications in nanodelivery and nanoreactors technology. Because few of the engineered channel proteins have been crystallized, solution studies based on Neutron Scattering, Circular Dichroism and NMR play a major role. Consequently, the understanding of membrane proteins dynamics in water/detergent solutions or when embedded in a lipid membrane, can clarify how the environment affects protein behavior. In this study, molecular dynamics simulations of the FhuA Escherichia coli outer membrane channel protein and its engineered FhuA Δ1-159 variant have been performed in two different environments: a DNPC (1,2-dinervonyl-sn-glycero-3-phosphocholine) lipid bilayer and a water/OES (N-octyl-2-hydroxyethyl sulfoxide) detergent solution. Furthermore the FhuA Δ1-159 variant has been simulated in the open and closed states, the last induced by the presence of six 3-(2-pyridyldithio)-propionic-acid in the channel inner core. Differences in protein structural and dynamical behavior between the two environments have been found. Considering the FhuA protein characterized by an elliptical-cylindrical symmetry: (a) neither variations on the secondary structure nor axial deformation have been observed in any of the systems; (b) the ellipticity of the channel section (open state) and its fluctuations are enhanced in presence of water/OES, while diminished or suppressed in the DNPC bilayer; (c) the insertion of hydrophobic pyridyl groups into the FhuA Δ1-159 channel (closed state) induces a higher ellipticity in water/OES solution, while shifting to a circular section in the DNPC membrane; (d) the cork domain represented by the first 159 amino acids does not play a major role for protein stability.

Molecular origin of electron paramagnetic resonance line shapes on β-barrel membrane proteins: the local solvation environment modulates spin-label configuration

In this work, electron paramagnetic resonance (EPR) spectroscopy and X-ray crystallography were used to examine the origins of EPR line shapes from spin-labels at the protein-lipid interface on the β-barrel membrane protein BtuB. Two atomic-resolution structures were obtained for the methanethiosulfonate spin-label derivatized to cysteines on the membrane-facing surface of BtuB. At one of these sites, position 156, the label side chain resides in a pocket formed by neighboring residues; however, it extends from the protein surface and yields a single-component EPR spectrum in the crystal that results primarily from fast rotation about the fourth and fifth bonds linking the spin-label to the protein backbone. In lipid bilayers, site 156 yields a multicomponent spectrum resulting from different rotameric states of the labeled side chain. Moreover, changes in the lipid environment, such as variations in bilayer thickness, modulate the EPR spectrum by modulating label rotamer populations. At a second site, position 371, the labeled side chain interacts with a pocket on the protein surface, leading to a highly immobilized single-component EPR spectrum that is not sensitive to hydrocarbon thickness. This spectrum is similar to that seen at other sites that are deep in the hydrocarbon, such as position 170. This work indicates that the rotameric states of spin-labels on exposed hydrocarbon sites are sensitive to the environment at the protein-hydrocarbon interface, and that this environment may modulate weak interactions between the labeled side chain and the protein surface. In the case of BtuB, lipid acyl chain packing is not symmetric around the β-barrel, and EPR spectra from labeled hydrocarbon-facing sites in BtuB may reflect this asymmetry. In addition to facilitating the interpretation of EPR spectra of membrane proteins, these results have important implications for the use of long-range distance restraints in protein structure refinement that are obtained from spin-labels.

Membrane thickness varies around the circumference of the transmembrane protein BtuB

BtuB is a large outer-membrane β-barrel protein that belongs to a class of active transport proteins that are TonB-dependent. These TonB-dependent transporters are based upon a 22-stranded antiparallel β-barrel, which is notably asymmetric in its length. Here, site-directed spin labeling and simulated annealing were used to locate the membrane lipid interface surrounding BtuB when reconstituted into phosphatidylcholine bilayers. Positions on the outer facing surface of the β-barrel and the periplasmic turns were spin-labeled and distances from the label to the membrane interface estimated by progressive power saturation of the electron paramagnetic resonance spectra. These distances were then used as atom-to-plane distance restraints in a simulated annealing routine, to dock the protein to two independent planes and produce a model representing the average position of the lipid phosphorus atoms at each interface. The model is in good agreement with the experimental data; however, BtuB is mismatched to the bilayer thickness and the resulting planes representing the bilayer interface are not parallel. In the model, the membrane thickness varies by 11 Å around the circumference of the protein, indicating that BtuB distorts the bilayer interface so that it is thinnest on the short side of the protein β-barrel.

Increased sensitivity and extended range of distance measurements in spin-labeled membrane proteins: Q-band double electron-electron resonance and nanoscale bilayers

We report a significant methodological advance in the application of double electron-electron resonance (DEER) spectroscopy to measure long-range distances in spin-labeled membrane proteins. In the pseudo two-dimensional environment of proteoliposomes, a steep intermolecular background shapes DEER signals leading to long accumulation times, complicating data analysis and reducing the maximal measurable distances from 70 A down to approximately 40-50 A. To eliminate these limitations, we took advantage of the homogeneity and monodispersity of a class of discoidal nanoscale phospholipid bilayers in conjunction with the micromolar DEER sensitivity at Q-band (34 GHz) microwave frequency. Spin-labeled mutants of the ABC transporter MsbA were functionally reconstituted at a ratio of one functional dimer per nanoscale apolipoprotein-bound bilayer (NABB). DEER echo intensities from NABB-reconstituted MsbA have linear baselines reflecting a three-dimensional spatial distribution. This results in an order-of-magnitude higher sensitivity at Q-band relative to proteoliposomes and restores the maximal observable distance effectively increasing experimental throughput. The advances described here set the stage for the use of DEER spectroscopy to analyze conformational dynamics of sample-limited eukaryotic membrane proteins.

Screening for the drug-phospholipid interaction: correlation to phospholipidosis

Phospholipid bilayers represent a complex, anisotropic environment fundamentally different from bulk oil or octanol, for instance. Even "simple" drug association to phospholipid bilayers can only be fully understood if the slab-of-hydrocarbon approach is abandoned and the complex, anisotropic properties of lipid bilayers reflecting the chemical structures and organization of the constituent phospholipids are considered. The interactions of drugs with phospholipids are important in various processes, such as drug absorption, tissue distribution, and subcellular distribution. In addition, drug-lipid interactions may lead to changes in lipid-dependent protein activities, and further, to functional and morphological changes in cells, a prominent example being the phospholipidosis (PLD) induced by cationic amphiphilic drugs. Herein we briefly review drug-lipid interactions in general and the significance of these interactions in PLD in particular. We also focus on a potential causal connection between drug-induced PLD and steatohepatitis, which is induced by some cationic amphiphilic drugs.

Coupling of calcium and substrate binding through loop alignment in the outer-membrane transporter BtuB

In Gram-negative bacteria, TonB-dependent outer-membrane transporters bind large, scarce organometallic substrates with high affinity preceding active transport. The cobalamin transporter BtuB requires the additional binding of two Ca(2+) ions before substrate binding can occur, but the underlying molecular mechanism is unknown. Using the crystallographic structures available for different bound states of BtuB, we have carried out extended molecular dynamics simulations of multiple functional states of BtuB to address the role of Ca(2+) in substrate recruitment. We find that Ca(2+) binding both stabilizes and repositions key extracellular loops of BtuB, optimizing interactions with the substrate. Interestingly, replacement by Mg(2+) abolishes this effect, in accordance with experiments. Using a set of new force-field parameters developed for cyanocobalamin, we also simulated the substrate-bound form of BtuB, where we observed interactions not seen in the crystal structure between the substrate and loops previously found to be important for binding and transport. Based on our results, we suggest that the large size of cobalamin compared to other TonB-dependent transporter substrates explains the requirement of Ca(2+) binding for high-affinity substrate recruitment in BtuB.