Chapter 8 Methods for Measuring the Thermodynamic Stability of Membrane Proteins

Utility of FRET in studies of membrane protein oligomerization: The concept of the effective dissociation constant

The activity of many membrane receptors is controlled through their lateral association into dimers or higher-order oligomers. Although Förster resonance energy transfer (FRET) measurements have been used extensively to characterize the stability of receptor dimers, the utility of FRET in studies of larger oligomers has been limited. Here we introduce an effective equilibrium dissociation constant that can be extracted from FRET measurements for EphA2, a receptor tyrosine kinase (RTK) known to form active oligomers of heterogeneous distributions in response to its ligand ephrinA1-Fc. The newly introduced effective equilibrium dissociation constant has a well-defined physical meaning and biological significance. It denotes the receptor concentration for which half of the receptors are monomeric and inactive, and the other half are associated into oligomers and are active, irrespective of the exact oligomer size. This work introduces a new dimension to the utility of FRET in studies of membrane receptor association and signaling in the plasma membrane.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Free-energy changes of bacteriorhodopsin point mutants measured by single-molecule force spectroscopy

Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔG). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG, we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR's G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule-derived ΔΔG for mutant L223A (-2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔG for mutant V217A was 2.2-fold larger (-2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val²¹⁷ and a natively bound squalene lipid, highlighting the contribution of membrane protein-lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔG for a fully folded membrane protein embedded in its native bilayer.

Probing Membrane Protein Association Using Concentration-Dependent Number and Brightness

We introduce concentration-dependent number and brightness (cdN&B), a fluorescence fluctuation technique that can be implemented on a standard confocal microscope and can report on the thermodynamics of membrane protein association in the native plasma membrane. It uses transient transfection to enable measurements of oligomer size as a function of receptor concentration over a broad range, yielding the association constant. We discuss artifacts in cdN&B that are concentration-dependent and can distort the oligomerization curves, and we outline procedures that can correct for them. Using cdN&B, we characterize the association of neuropilin 1 (NRP1), a protein that plays a critical role in the development of the embryonic cardiovascular and nervous systems. We show that NRP1 associates into a tetramer in a concentration-dependent manner, and we quantify the strength of the association. This work demonstrates the utility of cdN&B as a powerful tool in biophysical chemistry.

Potential causes and consequences of rapid mitochondrial genome evolution in thermoacidophilic Galdieria (Rhodophyta)

BACKGROUND

The Cyanidiophyceae is an early-diverged red algal class that thrives in extreme conditions around acidic hot springs. Although this lineage has been highlighted as a model for understanding the biology of extremophilic eukaryotes, little is known about the molecular evolution of their mitochondrial genomes (mitogenomes).

RESULTS

To fill this knowledge gap, we sequenced five mitogenomes from representative clades of Cyanidiophyceae and identified two major groups, here referred to as Galdieria-type (G-type) and Cyanidium-type (C-type). G-type mitogenomes exhibit the following three features: (i) reduction in genome size and gene inventory, (ii) evolution of unique protein properties including charge, hydropathy, stability, amino acid composition, and protein size, and (iii) distinctive GC-content and skewness of nucleotides. Based on GC-skew-associated characteristics, we postulate that unidirectional DNA replication may have resulted in the rapid evolution of G-type mitogenomes.

CONCLUSIONS

The high divergence of G-type mitogenomes was likely driven by natural selection in the multiple extreme environments that Galdieria species inhabit combined with their highly flexible heterotrophic metabolism. We speculate that the interplay between mitogenome divergence and adaptation may help explain the dominance of Galdieria species in diverse extreme habitats.

Quantifying the Native Energetics Stabilizing Bacteriorhodopsin by Single-Molecule Force Spectroscopy

We quantified the equilibrium (un)folding free energy ΔG_{0} of an eight-amino-acid region starting from the fully folded state of the model membrane-protein bacteriorhodopsin using single-molecule force spectroscopy. Analysis of equilibrium and nonequilibrium data yielded consistent, high-precision determinations of ΔG_{0} via multiple techniques (force-dependent kinetics, Crooks fluctuation theorem, and inverse Boltzmann analysis). We also deduced the full 1D projection of the free-energy landscape in this region. Importantly, ΔG_{0} was determined in bacteriorhodopsin's native bilayer, an advance over traditional results obtained by chemical denaturation in nonphysiological detergent micelles.

Solid-state NMR spectroscopy based atomistic view of a membrane protein unfolding pathway

Membrane protein folding, structure, and function strongly depend on a cell membrane environment, yet detailed characterization of folding within a lipid bilayer is challenging. Studies of reversible unfolding yield valuable information on the energetics of folding and on the hierarchy of interactions contributing to protein stability. Here, we devise a methodology that combines hydrogen-deuterium (H/D) exchange and solid-state NMR (SSNMR) to follow membrane protein unfolding in lipid membranes at atomic resolution through detecting changes in the protein water-accessible surface, and concurrently monitoring the reversibility of unfolding. We obtain atomistic description of the reversible part of a thermally induced unfolding pathway of a seven-helical photoreceptor. The pathway is visualized through SSNMR-detected snapshots of H/D exchange patterns as a function of temperature, revealing the unfolding intermediate and its stabilizing factors. Our approach is transferable to other membrane proteins, and opens additional ways to characterize their unfolding and stabilizing interactions with atomic resolution.

Design of self-assembling transmembrane helical bundles to elucidate principles required for membrane protein folding and ion transport

Ion transporters and channels are able to identify and act on specific substrates among myriads of ions and molecules critical to cellular processes, such as homeostasis, cell signalling, nutrient influx and drug efflux. Recently, we designed Rocker, a minimalist model for Zn²⁺/H⁺ co-transport. The success of this effort suggests that de novo membrane protein design has now come of age so as to serve a key approach towards probing the determinants of membrane protein folding, assembly and function. Here, we review general principles that can be used to design membrane proteins, with particular reference to helical assemblies with transport function. We also provide new functional and NMR data that probe the dynamic mechanism of conduction through Rocker.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.

Quantifying the Interaction between EGFR Dimers and Grb2 in Live Cells

Adaptor proteins are a class of cytoplasmic proteins that bind to phosphorylated residues in receptor tyrosine kinases and trigger signaling cascades that control critically important cellular processes, such as cell survival, growth, differentiation, and motility. Here, we seek to characterize the interaction between epidermal growth factor receptor (EGFR) and the cytoplasmic adaptor protein growth factor receptor-bound protein 2 (Grb2) in a cellular context. To do so, we explore the utility of a highly biologically relevant model system, mammalian cells under reversible osmotic stress, and a recently introduced Förster resonance energy transfer microscopy method, fully quantified spectral imaging. We present a method that allows us to quantify the stoichiometry and the association constant of the EGFR-Grb2 binding interaction in the plasma membrane, in the presence and absence of activating ligand. The method that we introduce can have broad utility in membrane protein research, as it can be applied to different membrane protein-cytoplasmic protein pairs.

Engineering a Novel Porin OmpGF Via Strand Replacement from Computational Analysis of Sequence Motif

β-Barrelmembrane proteins (βMPs) form barrel-shaped pores in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. Because of the robustness of their barrel structures, βMPs have great potential as nanosensors for single-molecule detection. However, natural βMPs currently employed have inflexible biophysical properties and are limited in their pore geometry, hindering their applications in sensing molecules of different sizes and properties. Computational engineering has the promise to generate βMPs with desired properties. Here we report a method for engineering novel βMPs based on the discovery of sequence motifs that predominantly interact with the cell membrane and appear in more than 75% of transmembrane strands. By replacing β1-β6 strands of the protein OmpF that lack these motifs with β1-β6 strands of OmpG enriched with these motifs and computational verification of increased stability of its transmembrane section, we engineered a novel βMP called OmpGF. OmpGF is predicted to form a monomer with a stable transmembrane region. Experimental validations showed that OmpGF could refold in vitro with a predominant β-sheet structure, as confirmed by circular dichroism. Evidence of OmpGF membrane insertion was provided by intrinsic tryptophan fluorescence spectroscopy, and its pore-forming property was determined by a dye-leakage assay. Furthermore, single-channel conductance measurements confirmed that OmpGF function as a monomer and exhibits increased conductance than OmpG and OmpF. These results demonstrated that a novel and functional βMP can be successfully engineered through strand replacement based on sequence motif analysis and stability calculation.

A New Method to Study Heterodimerization of Membrane Proteins and Its Application to Fibroblast Growth Factor Receptors

The activity of receptor tyrosine kinases (RTKs) is controlled through their lateral association in the plasma membrane. RTKs are believed to form both homodimers and heterodimers, and the different dimers are believed to play unique roles in cell signaling. However, RTK heterodimers remain poorly characterized, as compared with homodimers, because of limitations in current experimental methods. Here, we develop a FRET-based methodology to assess the thermodynamics of hetero-interactions in the plasma membrane. To demonstrate the utility of the methodology, we use it to study the hetero-interactions between three fibroblast growth factor receptors-FGFR1, FGFR2, and FGFR3-in the absence of ligand. Our results show that all possible FGFR heterodimers form, suggesting that the biological roles of FGFR heterodimers may be as significant as the homodimer roles. We further investigate the effect of two pathogenic point mutations in FGFR3 (A391E and G380R) on heterodimerization. We show that each of these mutations stabilize most of the heterodimers, with the largest effects observed for FGFR3 wild-type/mutant heterodimers. We thus demonstrate that the methodology presented here can yield new knowledge about RTK interactions and can further our understanding of signal transduction across the plasma membrane.

Fully quantified spectral imaging reveals in vivo membrane protein interactions

Here we introduce the fully quantified spectral imaging (FSI) method as a new tool to probe the stoichiometry and stability of protein complexes in biological membranes. The FSI method yields two dimensional membrane concentrations and FRET efficiencies in native plasma membranes. It can be used to characterize the association of membrane proteins: to differentiate between monomers, dimers, or oligomers, to produce binding (association) curves, and to measure the free energies of association in the membrane. We use the FSI method to study the lateral interactions of Vascular Endothelial Growth Factor Receptor 2 (VEGFR2), a member of the receptor tyrosine kinase (RTK) superfamily, in plasma membranes, in vivo. The knowledge gained through the use of the new method challenges the current understanding of VEGFR2 signaling.

Characterization of membrane protein interactions in plasma membrane derived vesicles with quantitative imaging Förster resonance energy transfer

Here we describe an experimental tool, termed quantitative imaging Förster resonance energy transfer (QI-FRET), that enables the quantitative characterization of membrane protein interactions. The QI-FRET methodology allows us to acquire binding curves and calculate association constants for complex membrane proteins in the native plasma membrane environment. The method utilizes FRET detection, and thus requires that the proteins of interest are labeled with florescent proteins, either FRET donors or FRET acceptors. Since plasma membranes of cells have complex topologies precluding the acquisition of two-dimensional binding curves, the FRET measurements are performed in plasma membrane derived vesicles that bud off cells as a result of chemical or osmotic stress. The results overviewed here are acquired in vesicles produced with an osmotic vesiculation buffer developed in our laboratory, which does not utilize harsh chemicals. The concentrations of the donor-labeled and the acceptor-labeled proteins are determined, along with the FRET efficiencies, in each vesicle. The experiments utilize transient transfection, such that a wide variety of concentrations is sampled. Then, data from hundreds of vesicles are combined to yield dimerization curves. Here we discuss recent findings about the dimerization of receptor tyrosine kinases (RTKs), membrane proteins that control cell growth and differentiation via lateral dimerization in the plasma membrane. We focus on the dimerization of fibroblast growth factor receptor 3 (FGFR3), a RTK that plays a critically important role in skeletal development. We study the role of different FGFR3 domains in FGFR3 dimerization in the absence of ligand, and we show that FGFR3 extracellular domains inhibit unliganded dimerization, while contacts between the juxtamembrane domains, which connect the transmembrane domains to the kinase domains, stabilize the unliganded FGFR3 dimers. Since FGFR3 has been documented to harbor many pathogenic single amino acid mutations that cause skeletal and cranial dysplasias, as well as cancer, we also study the effects of these mutations on dimerization. First, we show that the A391E mutation, linked to Crouzon syndrome with acanthosis nigricans and to bladder cancer, significantly enhances FGFR3 dimerization in the absence of ligand and thus induces aberrant receptor interactions. Second, we present results about the effect of three cysteine mutations that cause thanatophoric dysplasia, a lethal phenotype. Such cysteine mutations have been hypothesized previously to cause constitutive dimerization, but we find instead that they have a surprisingly modest effect on dimerization. Most of the studied pathogenic mutations also altered FGFR3 dimer structure, suggesting that both increases in dimerization propensities and changes in dimer structure contribute to the pathological phenotypes. The results acquired with the QI-FRET method further our understanding of the interactions between FGFR3 molecules and RTK molecules in general. Since RTK dimerization regulates RTK signaling, our findings advance our knowledge of RTK activity in health and disease. The utility of the QI-FRET method is not restricted to RTKs, and we thus hope that in the future the QI-FRET method will be applied to other classes of membrane proteins, such as channels and G protein-coupled receptors.

Cooperative folding of a polytopic α-helical membrane protein involves a compact N-terminal nucleus and nonnative loops

Despite the ubiquity of helical membrane proteins in nature and their pharmacological importance, the mechanisms guiding their folding remain unclear. We performed kinetic folding and unfolding experiments on 69 mutants (engineered every 2-3 residues throughout the 178-residue transmembrane domain) of GlpG, a membrane-embedded rhomboid protease from Escherichia coli. The only clustering of significantly positive ϕ-values occurs at the cytosolic termini of transmembrane helices 1 and 2, which we identify as a compact nucleus. The three loops flanking these helices show a preponderance of negative ϕ-values, which are sometimes taken to be indicative of nonnative interactions in the transition state. Mutations in transmembrane helices 3-6 yielded predominantly ϕ-values near zero, indicating that this part of the protein has denatured-state-level structure in the transition state. We propose that loops 1-3 undergo conformational rearrangements to position the folding nucleus correctly, which then drives folding of the rest of the domain. A compact N-terminal nucleus is consistent with the vectorial nature of cotranslational membrane insertion found in vivo. The origin of the interactions in the transition state that lead to a large number of negative ϕ-values remains to be elucidated.

Cysteine residue is not essential for CPM protein thermal-stability assay

A popular thermal-stability assay developed especially for the study of membrane proteins uses a thiol-specific probe, 7-diethylamino-3-(4-maleimidophenyl)-4-methylcoumarin (CPM). The fluorescence emission of CPM surges when it forms a covalent bond with the side chain of a free Cys, which becomes more readily accessible upon protein thermal denaturation. Interestingly, the melting temperatures of membrane proteins determined using the CPM assay in literature are closely clustered in the temperature range 45-55 °C. A thorough understanding of the mechanism behind the observed signal change is critical for the accurate interpretation of the protein unfolding. Here we used two α-helical membrane proteins, AqpZ and AcrB, as model systems to investigate the nature of the fluorescence surge in the CPM assay. We found that the transition temperatures measured using circular-dichroism (CD) spectroscopy and the CPM assay were significantly different. To eliminate potential artifact that might arise from the presence of detergent, we monitored the unfolding of two soluble proteins. We found that, contrary to current understanding, the presence of a sulfhydryl group was not a prerequisite for the CPM thermal-stability assay. The observed fluorescence increase is probably caused by binding of the fluorophore to hydrophobic patches exposed upon protein unfolding.

The safety dance: biophysics of membrane protein folding and misfolding in a cellular context

Most biological processes require the production and degradation of proteins, a task that weighs heavily on the cell. Mutations that compromise the conformational stability of proteins place both specific and general burdens on cellular protein homeostasis (proteostasis) in ways that contribute to numerous diseases. Efforts to elucidate the chain of molecular events responsible for diseases of protein folding address one of the foremost challenges in biomedical science. However, relatively little is known about the processes by which mutations prompt the misfolding of α-helical membrane proteins, which rely on an intricate network of cellular machinery to acquire and maintain their functional structures within cellular membranes. In this review, we summarize the current understanding of the physical principles that guide membrane protein biogenesis and folding in the context of mammalian cells. Additionally, we explore how pathogenic mutations that influence biogenesis may differ from those that disrupt folding and assembly, as well as how this may relate to disease mechanisms and therapeutic intervention. These perspectives indicate an imperative for the use of information from structural, cellular, and biochemical studies of membrane proteins in the design of novel therapeutics and in personalized medicine.

Mechanisms of integral membrane protein insertion and folding

The biogenesis, folding, and structure of α-helical membrane proteins (MPs) are important to understand because they underlie virtually all physiological processes in cells including key metabolic pathways, such as the respiratory chain and the photosystems, as well as the transport of solutes and signals across membranes. Nearly all MPs require translocons--often referred to as protein-conducting channels--for proper insertion into their target membrane. Remarkable progress toward understanding the structure and functioning of translocons has been made during the past decade. Here, we review and assess this progress critically. All available evidence indicates that MPs are equilibrium structures that achieve their final structural states by folding along thermodynamically controlled pathways. The main challenge for cells is the targeting and membrane insertion of highly hydrophobic amino acid sequences. Targeting and insertion are managed in cells principally by interactions between ribosomes and membrane-embedded translocons. Our review examines the biophysical and biological boundaries of MP insertion and the folding of polytopic MPs in vivo. A theme of the review is the under-appreciated role of basic thermodynamic principles in MP folding and assembly. Thermodynamics not only dictates the final folded structure but also is the driving force for the evolution of the ribosome-translocon system of assembly. We conclude the review with a perspective suggesting a new view of translocon-guided MP insertion.

Membrane proteins bind lipids selectively to modulate their structure and function

Previous studies have established that the folding, structure and function of membrane proteins are influenced by their lipid environments and that lipids can bind to specific sites, for example, in potassium channels. Fundamental questions remain however regarding the extent of membrane protein selectivity towards lipids. Here we report a mass spectrometry approach designed to determine the selectivity of lipid binding to membrane protein complexes. We investigate the mechanosensitive channel of large conductance (MscL) from Mycobacterium tuberculosis and aquaporin Z (AqpZ) and the ammonia channel (AmtB) from Escherichia coli, using ion mobility mass spectrometry (IM-MS), which reports gas-phase collision cross-sections. We demonstrate that folded conformations of membrane protein complexes can exist in the gas phase. By resolving lipid-bound states, we then rank bound lipids on the basis of their ability to resist gas phase unfolding and thereby stabilize membrane protein structure. Lipids bind non-selectively and with high avidity to MscL, all imparting comparable stability; however, the highest-ranking lipid is phosphatidylinositol phosphate, in line with its proposed functional role in mechanosensation. AqpZ is also stabilized by many lipids, with cardiolipin imparting the most significant resistance to unfolding. Subsequently, through functional assays we show that cardiolipin modulates AqpZ function. Similar experiments identify AmtB as being highly selective for phosphatidylglycerol, prompting us to obtain an X-ray structure in this lipid membrane-like environment. The 2.3 Å resolution structure, when compared with others obtained without lipid bound, reveals distinct conformational changes that re-position AmtB residues to interact with the lipid bilayer. Our results demonstrate that resistance to unfolding correlates with specific lipid-binding events, enabling a distinction to be made between lipids that merely bind from those that modulate membrane protein structure and/or function. We anticipate that these findings will be important not only for defining the selectivity of membrane proteins towards lipids, but also for understanding the role of lipids in modulating protein function or drug binding.

Correlation between AcrB trimer association affinity and efflux activity

The majority of membrane proteins function as oligomers. However, it remains largely unclear how the oligomer stability of protein complexes correlates with their function. Understanding the relationship between oligomer stability and activity is essential to protein research and to virtually all cellular processes that depend on the function of protein complexes. Proteins make lasting or transient interactions as they perform their functions. Obligate oligomeric proteins exist and function exclusively at a specific oligomeric state. Although oligomerization is clearly critical for such proteins to function, a direct correlation between oligomer affinity and biological activity has not yet been reported. Here, we used an obligate trimeric membrane transporter protein, AcrB, as a model to investigate the correlation between its relative trimer affinity and efflux activity. AcrB is a component of the major multidrug efflux system in Escherichia coli. We created six AcrB constructs with mutations at the transmembrane intersubunit interface, and we determined their activities using both a drug susceptibility assay and an ethidium bromide accumulation assay. The relative trimer affinities of these mutants in detergent micelles were obtained using blue native polyacrylamide gel electrophoresis. A correlation between the relative trimer affinity and substrate efflux activity was observed, in which a threshold trimer stability was required to maintain efflux activity. The trimer affinity of the wild-type protein was approximately 3 kcal/mol more stable than the threshold value. Once the threshold was reached, an additional increase of stability in the range observed had no observable effect on protein activity.

Unfolding study of a trimeric membrane protein AcrB

The folding of a multi-domain trimeric α-helical membrane protein, Escherichia coli inner membrane protein AcrB, was investigated. AcrB contains both a transmembrane domain and a large periplasmic domain. Protein unfolding in sodium dodecyl sulfate (SDS) and urea was monitored using the intrinsic fluorescence and circular dichroism spectroscopy. The SDS denaturation curve displayed a sigmoidal profile, which could be fitted with a two-state unfolding model. To investigate the unfolding of separate domains, a triple mutant was created, in which all three Trp residues in the transmembrane domain were replaced with Phe. The SDS unfolding profile of the mutant was comparable to that of the wild type AcrB, suggesting that the observed signal change was largely originated from the unfolding of the soluble domain. Strengthening of trimer association through the introduction of an inter-subunit disulfide bond had little effect on the unfolding profile, suggesting that trimer dissociation was not the rate-limiting step in unfolding monitored by fluorescence emission. Under our experimental condition, AcrB unfolding was not reversible. Furthermore, we experimented with the refolding of a monomeric mutant, AcrBΔloop , from the SDS unfolded state. The CD spectrum of the refolded AcrBΔloop superimposed well onto the spectra of the original folded protein, while the fluorescence spectrum was not fully recovered. In summary, our results suggested that the unfolding of the trimeric AcrB started with a local structural rearrangement. While the refolding of secondary structure in individual monomers could be achieved, the re-association of the trimer might be the limiting factor to obtain folded wild-type AcrB.

Kinetics and thermodynamics of membrane protein folding

Understanding protein folding has been one of the great challenges in biochemistry and molecular biophysics. Over the past 50 years, many thermodynamic and kinetic studies have been performed addressing the stability of globular proteins. In comparison, advances in the membrane protein folding field lag far behind. Although membrane proteins constitute about a third of the proteins encoded in known genomes, stability studies on membrane proteins have been impaired due to experimental limitations. Furthermore, no systematic experimental strategies are available for folding these biomolecules in vitro. Common denaturing agents such as chaotropes usually do not work on helical membrane proteins, and ionic detergents have been successful denaturants only in few cases. Refolding a membrane protein seems to be a craftsman work, which is relatively straightforward for transmembrane β-barrel proteins but challenging for α-helical membrane proteins. Additional complexities emerge in multidomain membrane proteins, data interpretation being one of the most critical. In this review, we will describe some recent efforts in understanding the folding mechanism of membrane proteins that have been reversibly refolded allowing both thermodynamic and kinetic analysis. This information will be discussed in the context of current paradigms in the protein folding field.

Retinal proteins as model systems for membrane protein folding

Experimental folding studies of membrane proteins are more challenging than water-soluble proteins because of the higher hydrophobicity content of membrane embedded sequences and the need to provide a hydrophobic milieu for the transmembrane regions. The first challenge is their denaturation: due to the thermodynamic instability of polar groups in the membrane, secondary structures in membrane proteins are more difficult to disrupt than in soluble proteins. The second challenge is to refold from the denatured states. Successful refolding of membrane proteins has almost always been from very subtly denatured states. Therefore, it can be useful to analyze membrane protein folding using computational methods, and we will provide results obtained with simulated unfolding of membrane protein structures using the Floppy Inclusions and Rigid Substructure Topography (FIRST) method. Computational methods have the advantage that they allow a direct comparison between diverse membrane proteins. We will review here both, experimental and FIRST studies of the retinal binding proteins bacteriorhodopsin and mammalian rhodopsin, and discuss the extension of the findings to deriving hypotheses on the mechanisms of folding of membrane proteins in general. This article is part of a Special Issue entitled: Retinal Proteins-You can teach an old dog new tricks.

Role of a conserved residue R780 in Escherichia coli multidrug transporter AcrB

Multidrug efflux pumps play important roles in bacteria drug resistance. A major multidrug efflux system in Gram-negative bacteria is composed of the inner membrane transporter AcrB, outer membrane protein channel TolC, and membrane fusion protein AcrA. These three proteins form a large complex that spans both layers of cell membranes and the periplasmic space. AcrB exists and functions as a homotrimer. To identify residues at the trimer interface that play important roles in AcrB function, we conducted site directed mutagenesis and discovered a key residue, R780. Although R780K was partially functional, all other R780 mutants tested were completely nonfunctional. Replacement of R780 by other residues disrupted trimer association. However, a decrease of trimer stability was not the lone cause for the observed loss of activity, because the activity loss could not be restored by strengthening trimer interaction. Using both heat and chemical denaturation methods, we found that the mutation decreased protein stability. Finally, we identified a repressor mutation, M774K, through random mutagenesis. It restored the activity of AcrBR780A to a level close to that of the wild-type protein. To examine the mechanism of activity restoration, we monitored denaturation of AcrBR780A/M774K and found that the repressor mutation improved protein stability. These results suggest that R780 is critical for AcrB stability. When R780 was replaced by Ala, the protein retained the overall structure, still trimerized in the cell membrane, and interacted with AcrA. However, local structural rearrangement might have occurred and lead to the decrease of protein stability and loss of substrate efflux activity.

Reversible folding of human peripheral myelin protein 22, a tetraspan membrane protein

Misfolding of the α-helical membrane protein peripheral myelin protein 22 (PMP22) has been implicated in the pathogenesis of the common neurodegenerative disease known as Charcot-Marie-Tooth disease (CMTD) and also several other related peripheral neuropathies. Emerging evidence suggests that the propensity of PMP22 to misfold in the cell may be due to an intrinsic lack of conformational stability. Therefore, quantitative studies of the conformational equilibrium of PMP22 are needed to gain insight into the molecular basis of CMTD. In this work, we have investigated the folding and unfolding of wild type (WT) human PMP22 in mixed micelles. Both kinetic and thermodynamic measurements demonstrate that the denaturation of PMP22 by n-lauroyl sarcosine (LS) in dodecylphosphocholine (DPC) micelles is reversible. Assessment of the conformational equilibrium indicates that a significant fraction of unfolded PMP22 persists even in the absence of the denaturing detergent. However, we find the stability of PMP22 is increased by glycerol, which facilitates quantitation of thermodynamic parameters. To our knowledge, this work represents the first report of reversible unfolding of a eukaryotic multispan membrane protein. The results indicate that WT PMP22 possesses minimal conformational stability in micelles, which parallels its poor folding efficiency in the endoplasmic reticulum. Folding equilibrium measurements for PMP22 in micelles may provide an approach to assess the effects of cellular metabolites or potential therapeutic agents on its stability. Furthermore, these results pave the way for future investigation of the effects of pathogenic mutations on the conformational equilibrium of PMP22.

Measuring transmembrane helix interaction strengths in lipid bilayers using steric trapping

We have developed a method to measure strong transmembrane (TM) helix interaction affinities in lipid bilayers that are difficult to measure by traditional dilution methods. The method, called steric trapping, couples dissociation of biotinylated TM helices to a competitive binding by monovalent streptavidin (mSA), so that dissociation is driven by the affinity of mSA for biotin and mSA concentration. By adjusting the binding affinity of mSA through mutation, the method can obtain dissociation constants of TM helix dimers (K d,dimer) over a range of six orders of magnitudes. The K d,dimer limit of measurable target interaction is extended 3-4 orders of magnitude lower than possible by dilution methods. Thus, steric trapping opens up new opportunities to study the folding and assembly of α-helical membrane proteins in lipid bilayer environments. Here we provide detailed methods for applying steric trapping to a TM helix dimer.

Thermal unfolding of a mammalian pentameric ligand-gated ion channel proceeds at consecutive, distinct steps

Pentameric ligand-gated ion channels (LGICs) play an important role in fast synaptic signal transduction. Binding of agonists to the β-sheet-structured extracellular domain opens an ion channel in the transmembrane α-helical region of the LGIC. How the structurally distinct and distant domains are functionally coupled for such central transmembrane signaling processes remains an open question. To obtain detailed information about the stability of and the coupling between these different functional domains, we analyzed the thermal unfolding of a homopentameric LGIC, the 5-hydroxytryptamine receptor (ligand binding, secondary structure, accessibility of Trp and Cys residues, and aggregation), in plasma membranes as well as during detergent extraction, purification, and reconstitution into artificial lipid bilayers. We found a large loss in thermostability correlating with the loss of the lipid bilayer during membrane solubilization and purification. Thermal unfolding of the 5-hydroxytryptamine receptor occurred in consecutive steps at distinct protein locations. A loss of ligand binding was detected first, followed by formation of different transient low oligomeric states of receptor pentamers, followed by partial unfolding of helical parts of the protein, which finally lead to the formation receptor aggregates. Structural destabilization of the receptor in detergents could be partially reversed by reconstituting the receptor into lipid bilayers. Our results are important because they quantify the stability of LGICs during detergent extraction and purification and can be used to create stabilized receptor proteins for structural and functional studies.

Production of plasma membrane vesicles with chloride salts and their utility as a cell membrane mimetic for biophysical characterization of membrane protein interactions

Plasma membrane derived vesicles are used as a model system for the biochemical and biophysical investigations of membrane proteins and membrane organization. The most widely used vesiculation procedure relies on formaldehyde and dithiothreitol (DTT), but these active chemicals may introduce artifacts in the experimental results. Here we describe a procedure to vesiculate Chinese hamster ovary (CHO) cells, widely used for the expression of recombinant proteins, using a hypertonic vesiculation buffer containing chloride salts and no formaldehyde or DTT. We characterize the size distribution of the produced vesicles. We also show that these vesicles can be used for the biophysical characterization of interactions between membrane proteins.

Architectural and thermodynamic principles underlying intramembrane protease function

Intramembrane proteases hydrolyze peptide bonds within the membrane as a signaling paradigm universal to all life forms and with implications in disease. Deciphering the architectural strategies supporting intramembrane proteolysis is an essential but unattained goal. We integrated a new, quantitative and high-throughput thermal light-scattering technology, reversible equilibrium un/refolding, and quantitative protease assays to interrogate rhomboid architecture with 151 purified variants. Rhomboid proteases maintain low intrinsic thermodynamic stability (ΔG=2.1-4.5kcal/mol) resulting from a multitude of generally-weak transmembrane packing interactions, making them highly-responsive to their environment. Stability is consolidated by two buried glycines and several packing leucines, with a few multifaceted hydrogen bonds strategically-deployed to two peripheral regions. Opposite these regions lie transmembrane segment 5 and connected loops that are notably exempt of structural responsibility, suggesting intramembrane proteolysis involves considerable but localized protein dynamics. Our analyses provide a comprehensive ‘heat map’ of the physio-chemical anatomy underlying membrane-immersed enzyme function at unprecedented resolution.

Protein Structure in Membrane Domains

Of great interest to the academic and pharmaceutical research communities, helical transmembrane proteins are characterized by their ability to dissolve and fold in lipid bilayers—properties conferred by polypeptide spans termed transmembrane domains (TMDs). The apolar nature of TMDs necessitates the use of membrane-mimetic solvents for many structure and folding studies. This review examines the relationship between TMD structure and solvent environment, focusing on principles elucidated largely in membrane-mimetic environments with single-TMD protein and peptide models. Following a brief description of TMD sequence and conformational characteristics gleaned from the structural database, we present an overview of the conceptual models used to study folding in vitro. The impact of sequence and solvent context on the incorporation of TMDs into membranes, and its role in measurements of TMD self-assembly strengths, is then described. We conclude with a discussion of the nonspecific effects of membrane components on TMD stability.

Engineered oligomerization state of OmpF protein through computational design decouples oligomer dissociation from unfolding

Biogenesis of β-barrel membrane proteins is a complex, multistep, and as yet incompletely characterized process. The bacterial porin family is perhaps the best-studied protein family among β-barrel membrane proteins that allows diffusion of small solutes across the bacterial outer membrane. In this study, we have identified residues that contribute significantly to the protein-protein interaction (PPI) interface between the chains of outer membrane protein F (OmpF), a trimeric porin, using an empirical energy function in conjunction with an evolutionary analysis. By replacing these residues through site-directed mutagenesis either with energetically favorable residues or substitutions that do not occur in natural bacterial outer membrane proteins, we succeeded in engineering OmpF mutants with dimeric and monomeric oligomerization states instead of a trimeric oligomerization state. Moreover, our results suggest that the oligomerization of OmpF proceeds through a series of interactions involving two distinct regions of the extensive PPI interface: two monomers interact to form a dimer through the PPI interface near G19. This dimer then interacts with another monomer through the PPI interface near G135 to form a trimer. We have found that perturbing the PPI interface near G19 results in the formation of the monomeric OmpF only. Thermal denaturation of the designed dimeric OmpF mutant suggests that oligomer dissociation can be separated from the process of protein unfolding. Furthermore, the conserved site near G57 and G59 is important for the PPI interface and might provide the essential scaffold for PPIs.

Predicting three-dimensional structures of transmembrane domains of β-barrel membrane proteins

β-Barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. They are important for pore formation, membrane anchoring, and enzyme activity. These proteins are also often responsible for bacterial virulence. Due to difficulties in experimental structure determination, they are sparsely represented in the protein structure databank. We have developed a computational method for predicting structures of the transmembrane (TM) domains of β-barrel membrane proteins. Based on physical principles, our method can predict structures of the TM domain of β-barrel membrane proteins of novel topology, including those from eukaryotic mitochondria. Our method is based on a model of physical interactions, a discrete conformational state space, an empirical potential function, as well as a model to account for interstrand loop entropy. We are able to construct three-dimensional atomic structure of the TM domains from sequences for a set of 23 nonhomologous proteins (resolution 1.8-3.0 Å). The median rmsd of TM domains containing 75-222 residues between predicted and measured structures is 3.9 Å for main chain atoms. In addition, stability determinants and protein-protein interaction sites can be predicted. Such predictions on eukaryotic mitochondria outer membrane protein Tom40 and VDAC are confirmed by independent mutagenesis and chemical cross-linking studies. These results suggest that our model captures key components of the organization principles of β-barrel membrane protein assembly.

Revisiting the folding kinetics of bacteriorhodopsin

The elucidation of the physical principles that govern the folding and stability of membrane proteins is one of the greatest challenges in protein science. Several insights into the folding of α-helical membrane proteins have come from the investigation of the conformational equilibrium of H. halobium bacteriorhodopsin (bR) in mixed micelles using SDS as a denaturant. In an effort to confirm that folded bR and SDS-denatured bR reach the same conformational equilibrium, we found that bR folding is significantly slower than has been previously known. Interrogation of the effect of the experimental variables on folding kinetics reveals that the rate of folding is dependent not only on the mole fraction of SDS but also on the molar concentrations of mixed micelle components, a variable that was not controlled in the previous study of bR folding kinetics. Moreover, when the molar concentrations of mixed micelle components are fixed at the concentrations commonly employed for bR equilibrium studies, conformational relaxation in the transition zone is slower than hydrolysis of the retinal Schiff base. As a result, the conformational equilibrium between folded bR and SDS-denatured bR cannot be achieved under the conventional condition. Our finding suggests that the molar concentrations of mixed micelle components are important experimental variables in the investigation of the kinetics and thermodynamics of bR folding and should be accounted for to ensure the accurate assessment of the conformational equilibrium of bR without the interference of retinal hydrolysis.

Physical determinants of β-barrel membrane protein folding in lipid vesicles

The spontaneous folding of two Neisseria outer membrane proteins, opacity-associated (Opa)(60) and Opa(50) into lipid vesicles was investigated by systematically varying bulk and membrane properties. Centrifugal fractionation coupled with sodium dodecyl sulfate polyacrylamide gel electrophoresis mobility assays enabled the discrimination of aggregate, unfolded membrane-associated, and folded membrane-inserted protein states as well as the influence of pH, ionic strength, membrane surface potential, lipid saturation, and urea on each. Protein aggregation was reduced with increasing lipid chain length, basic pH, low salt, the incorporation of negatively charged guest lipids, or by the addition of urea to the folding reaction. Insertion from the membrane-associated form was improved in shorter chain lipids, with more basic pH and low ionic strength; it is hindered by unsaturated or ether-linked lipids. The isolation of the physical determinants of insertion suggests that the membrane surface and dipole potentials are driving forces for outer membrane protein insertion and folding into lipid bilayers.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

High-throughput selection of transmembrane sequences that enhance receptor tyrosine kinase activation

Dimerization is a critical requirement for the activation of the intracellular kinase domains of receptor tyrosine kinases (RTKs). The single transmembrane (TM) helices of RTKs contribute to dimerization, but the details are not well understood. Work with TM helices in various model systems has revealed a small number of specific dimerization sequence motifs, and it has been suggested that RTK dimerization is modulated by such motifs. Yet questions remain about the universality of these sequence motifs for RTK dimerization and about how TM domain dimerization in model systems relates to RTK activation in mammalian membranes. To investigate these questions, we designed a 3888-member combinatorial peptide library based on the TM domain of Neu (ErbB2) as a model RTK. The library contains many closely related, Neu-like sequences, including thousands of sequences with known dimerization motifs. We used an SDS-PAGE-based screen to select peptides that dimerize better than the native Neu sequence, and we assayed the activation of chimeric Neu receptors in mammalian cells with TM sequences selected in the screen. Despite the very high abundance of known dimerization motifs in the library, only a very few dimerizing sequences were identified by SDS-PAGE. About half of those sequences activated the Neu kinase significantly more than did the wild-type TM sequence. This work furthers our knowledge about the requirements for membrane protein interactions and the requirements for RTK activation in cells.

Dramatic destabilization of transmembrane helix interactions by features of natural membrane environments

Membrane proteins have evolved to fold and function in a lipid bilayer, so it is generally assumed that their stability should be optimized in a natural membrane environment. Yet optimal stability is not always in accord with optimization of function, so evolutionary pressure, occurring in a complex membrane environment, may favor marginal stability. Here, we find that the transmembrane helix dimer, glycophorin A (GpATM), is actually much less stable in the heterogeneous environment of a natural membrane than it is in model membranes and even common detergents. The primary destabilizing factors are electrostatic interactions between charged lipids and charged GpATM side chains, and nonspecific competition from other membrane proteins. These effects overwhelm stabilizing contributions from lateral packing pressure and excluded volume. Our work illustrates how evolution can employ membrane composition to modulate protein stability.

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.

Probing membrane protein unfolding with pulse proteolysis

Technical challenges have greatly impeded the investigation of membrane protein folding and unfolding. To develop a new tool that facilitates the study of membrane proteins, we tested pulse proteolysis as a probe for membrane protein unfolding. Pulse proteolysis is a method to monitor protein folding and unfolding, which exploits the significant difference in proteolytic susceptibility between folded and unfolded proteins. This method requires only a small amount of protein and, in many cases, may be used with unpurified proteins in cell lysates. To evaluate the effectiveness of pulse proteolysis as a probe for membrane protein unfolding, we chose Halobacterium halobium bacteriorhodopsin (bR) as a model system. The denaturation of bR in SDS has been investigated extensively by monitoring the change in the absorbance at 560 nm (A(560)). In this work, we demonstrate that denaturation of bR by SDS results in a significant increase in its susceptibility to proteolysis by subtilisin. When pulse proteolysis was applied to bR incubated in varying concentrations of SDS, the remaining intact protein determined by electrophoresis shows a cooperative transition. The midpoint of the cooperative transition (C(m)) shows excellent agreement with that determined by A(560). The C(m) values determined by pulse proteolysis for M56A and Y57A bRs are also consistent with the measurements made by A(560). Our results suggest that pulse proteolysis is a quantitative tool to probe membrane protein unfolding. Combining pulse proteolysis with Western blotting may allow the investigation of membrane protein unfolding in situ without overexpression or purification.