The transition state for integral membrane protein folding

Helical reorganization in the context of membrane protein folding: Insights from simulations with bacteriorhodopsin (BR) fragments

Membrane protein folding is distinct from folding of soluble proteins. Conformational acquisition in major membrane protein subclasses can be delineated into insertion and folding processes. An exception to the "two stage" folding, later developed to "three stage" folding, is observed within the last two helices in bacteriorhodopsin (BR), a system that serves as a model membrane protein. We employ a reductionist approach to understand interplay of molecular factors underlying the apparent defiance. Leveraging available solution NMR structures, we construct, sample in silico, and analyze partially (PIn) and fully inserted (FIn) BR membrane states. The membrane lateral C-terminal helix (CH) in PIn is markedly prone to transient structural distortions over microsecond timescales; a disorder prone region (DPR) is thereby identified. While clear transmembrane propensities are not acquired, the distortions induce alterations in local membrane curvature and area per lipid. Importantly, energetic decompositions reveal that overall, the N-terminal helix (NH) is thermodynamically more stable in the PIn. Higher overall stability of the FIn arises from favorable interactions between the NH and the CH. Our results establish lack of spontaneous transition of the PIn to the FIn, and attributes their partitioning to barriers that exceed those accessible with thermal fluctuations. This work paves the way for further detailed studies aimed at determining the thermo-kinetic roles of the initial five helices, or complementary external factors, in complete helical folding and insertion in BR. We comment that complementing such efforts with the growing field of machine learning assisted energy landscape searches may offer unprecedented insights.

Probing the energy barriers and stages of membrane protein unfolding using solid-state NMR spectroscopy

Understanding how the amino acid sequence dictates protein structure and defines its stability is a fundamental problem in molecular biology. It is especially challenging for membrane proteins that reside in the complex environment of a lipid bilayer. Here, we obtain an atomic-level picture of the thermally induced unfolding of a membrane-embedded α-helical protein, human aquaporin 1, using solid-state nuclear magnetic resonance spectroscopy. Our data reveal the hierarchical two-step pathway that begins with unfolding of a structured extracellular loop and proceeds to an intermediate state with a native-like helical packing. In the second step, the transmembrane domain unravels as a single unit, resulting in a heterogeneous misfolded state with high helical content but with nonnative helical packing. Our results show the importance of loops for the kinetic stabilization of the whole membrane protein structure and support the three-stage membrane protein folding model.

The full spectrum of SLC22 OCT1 mutations illuminates the bridge between drug transporter biophysics and pharmacogenomics

Mutations in transporters can impact an individual's response to drugs and cause many diseases. Few variants in transporters have been evaluated for their functional impact. Here, we combine saturation mutagenesis and multi-phenotypic screening to dissect the impact of 11,213 missense single-amino-acid deletions, and synonymous variants across the 554 residues of OCT1, a key liver xenobiotic transporter. By quantifying in parallel expression and substrate uptake, we find that most variants exert their primary effect on protein abundance, a phenotype not commonly measured alongside function. Using our mutagenesis results combined with structure prediction and molecular dynamic simulations, we develop accurate structure-function models of the entire transport cycle, providing biophysical characterization of all known and possible human OCT1 polymorphisms. This work provides a complete functional map of OCT1 variants along with a framework for integrating functional genomics, biophysical modeling, and human genetics to predict variant effects on disease and drug efficacy.

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.

Factors That Control the Force Needed to Unfold a Membrane Protein in Silico Depend on the Mode of Denaturation

Single-molecule force spectroscopy methods, such as AFM and magnetic tweezers, have proved extremely beneficial in elucidating folding pathways for soluble and membrane proteins. To identify factors that determine the force rupture levels in force-induced membrane protein unfolding, we applied our near-atomic-level Upside molecular dynamics package to study the vertical and lateral pulling of bacteriorhodopsin (bR) and GlpG, respectively. With our algorithm, we were able to selectively alter the magnitudes of individual interaction terms and identify that, for vertical pulling, hydrogen bond strength had the strongest effect, whereas other non-bonded protein and membrane-protein interactions had only moderate influences, except for the extraction of the last helix where the membrane-protein interactions had a stronger influence. The up-down topology of the transmembrane helices caused helices to be pulled out as pairs. The rate-limiting rupture event often was the loss of H-bonds and the ejection of the first helix, which then propagated tension to the second helix, which rapidly exited the bilayer. The pulling of the charged linkers across the membrane had minimal influence, as did changing the bilayer thickness. For the lateral pulling of GlpG, the rate-limiting rupture corresponded to the separation of the helices within the membrane, with the H-bonds generally being broken only afterward. Beyond providing a detailed picture of the rupture events, our study emphasizes that the pulling mode greatly affects the factors that determine the forces needed to unfold a membrane protein.

Methods to study folding of alpha-helical membrane proteins in lipids

How alpha-helical membrane proteins fold correctly in the highly hydrophobic membrane interior is not well understood. Their folding is known to be highly influenced by the lipids within the surrounding bilayer, but the majority of folding studies have focused on detergent-solubilized protein rather than protein in a lipid environment. There are different ways to study folding in lipid bilayers, and each method has its own advantages and disadvantages. This review will discuss folding methods which can be used to study alpha-helical membrane proteins in bicelles, liposomes, nanodiscs or native membranes. These folding methods include in vitro folding methods in liposomes such as denaturant unfolding studies, and single-molecule force spectroscopy studies in bicelles, liposomes and native membranes. This review will also discuss recent advances in co-translational folding studies, which use cell-free expression with liposomes or nanodiscs or are performed in vivo with native membranes.

How lipids affect the energetics of co-translational alpha helical membrane protein folding

Membrane proteins need to fold with precision in order to function correctly, with misfolding potentially leading to disease. The proteins reside within a hydrophobic lipid membrane and must insert into the membrane and fold correctly, generally whilst they are being translated by the ribosome. Favourable and unfavourable free energy contributions are present throughout each stage of insertion and folding. The unfavourable energy cost of transferring peptide bonds into the hydrophobic membrane interior is compensated for by the favourable hydrophobic effect of partitioning a hydrophobic transmembrane alpha-helix into the membrane. Native membranes are composed of many different types of lipids, but how these different lipids influence folding and the associated free energies is not well understood. Altering the lipids in the bilayer is known to affect the probability of transmembrane helix insertion into the membrane, and lipids also affect protein stability and can promote successful folding. This review will summarise the free energy contributions associated with insertion and folding of alpha helical membrane proteins, as well as how lipids can make these processes more or less favourable. We will also discuss the implications of this work for the free energy landscape during the co-translational folding of alpha helical membrane proteins.

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.

The high-energy transition state of the glutamate transporter homologue GltPh

Membrane transporters mediate cellular uptake of nutrients, signaling molecules, and drugs. Their overall mechanisms are often well understood, but the structural features setting their rates are mostly unknown. Earlier single-molecule fluorescence imaging of the archaeal model glutamate transporter homologue GltPh from Pyrococcus horikoshii suggested that the slow conformational transition from the outward- to the inward-facing state, when the bound substrate is translocated from the extracellular to the cytoplasmic side of the membrane, is rate limiting to transport. Here, we provide insight into the structure of the high-energy transition state of GltPh that limits the rate of the substrate translocation process. Using bioinformatics, we identified GltPh gain-of-function mutations in the flexible helical hairpin domain HP2 and applied linear free energy relationship analysis to infer that the transition state structurally resembles the inward-facing conformation. Based on these analyses, we propose an approach to search for allosteric modulators for transporters.

Watching helical membrane proteins fold reveals a common N-to-C-terminal folding pathway

To understand membrane protein biogenesis, we need to explore folding within a bilayer context. Here, we describe a single-molecule force microscopy technique that monitors the folding of helical membrane proteins in vesicle and bicelle environments. After completely unfolding the protein at high force, we lower the force to initiate folding while transmembrane helices are aligned in a zigzag manner within the bilayer, thereby imposing minimal constraints on folding. We used the approach to characterize the folding pathways of the Escherichia coli rhomboid protease GlpG and the human β₂-adrenergic receptor. Despite their evolutionary distance, both proteins fold in a strict N-to-C-terminal fashion, accruing structures in units of helical hairpins. These common features suggest that integral helical membrane proteins have evolved to maximize their fitness with cotranslational folding.

Revisiting the dynamic and thermodynamic processes driving the record-breaking January 2014 precipitation in the southern UK

Many attribution studies of precipitation extreme events have attempted to estimate the thermodynamic contribution (linked to temperature changes) and the dynamic contribution (linked to the atmospheric circulation). Those studies are based on statistical decompositions of atmospheric fields, and essentially focus on the horizontal motion of the atmosphere. This paper proposes a framework that decomposes those terms from first physical principles, which include the vertical atmospheric motion that has often been overlooked. The goal is to take into account the driving processes of the extreme event. We revisit a recent example of extreme precipitation that was extensively investigated through its relation with the atmospheric circulation. We find that although the horizontal motion plays a minor (but important) role, the vertical motion yields a dominating contribution to the event that is larger than the thermodynamic contribution. This analysis quantifies the processes leading to high winter precipitation rates, and can be extended for further attribution studies.

Characterization of Denatured States and Reversible Unfolding of Sensory Rhodopsin II

Our understanding on the folding of membrane proteins lags behind that of soluble proteins due to challenges posed by the exposure of hydrophobic regions during in vitro chemical denaturation and refolding experiments. While different folding models are accepted for soluble proteins, only the two-stage model and the long-range interactions model have been proposed so far for helical membrane proteins. To address our knowledge gap on how different membrane proteins traverse their folding pathways, we have systematically investigated the structural features of SDS-denatured states and the kinetics for reversible unfolding of sensory rhodopsin II (pSRII), a retinal-binding photophobic receptor from Natronomonas pharaonis. pSRII is difficult to denature, and only SDS can dislodge the retinal chromophore without rapid aggregation. Even in 30% SDS (0.998 Χ_SDS), pSRII retains the equivalent of six out of seven transmembrane helices, while the retinal-binding pocket is disrupted, with transmembrane residues becoming more solvent exposed. Folding of pSRII from an SDS-denatured state harboring a covalently bound retinal chromophore shows deviations from an apparent two-state behavior. SDS denaturation to form the sensory opsin apo-protein is reversible. We report pSRII as a new model protein which is suitable for membrane protein folding studies and has a unique folding mechanism that differs from those of bacteriorhodopsin and bovine rhodopsin.

Lipid bilayer composition modulates the unfolding free energy of a knotted α-helical membrane protein

α-Helical membrane proteins have eluded investigation of their thermodynamic stability in lipid bilayers. Reversible denaturation curves have enabled some headway in determining unfolding free energies. However, these parameters have been limited to detergent micelles or lipid bicelles, which do not possess the same mechanical properties as lipid bilayers that comprise the basis of natural membranes. We establish reversible unfolding of the membrane transporter LeuT in lipid bilayers, enabling the comparison of apparent unfolding free energies in different lipid compositions. LeuT is a bacterial ortholog of neurotransmitter transporters and contains a knot within its 12-transmembrane helical structure. Urea is used as a denaturant for LeuT in proteoliposomes, resulting in the loss of up to 30% helical structure depending upon the lipid bilayer composition. Urea unfolding of LeuT in liposomes is reversible, with refolding in the bilayer recovering the original helical structure and transport activity. A linear dependence of the unfolding free energy on urea concentration enables the free energy to be extrapolated to zero denaturant. Increasing lipid headgroup charge or chain lateral pressure increases the thermodynamic stability of LeuT. The mechanical and charge properties of the bilayer also affect the ability of urea to denature the protein. Thus, we not only gain insight to the long-sought-after thermodynamic stability of an α-helical protein in a lipid bilayer but also provide a basis for studies of the folding of knotted proteins in a membrane environment.

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.

Kinetics of peptide folding in lipid membranes

Despite our extensive understanding of water-soluble protein folding kinetics, much less is known about the folding dynamics and mechanisms of membrane proteins. However, recent studies have shown that for relatively simple systems, such as peptides that form a transmembrane α-helix, helical dimer, or helix-turn-helix, it is possible to assess the kinetics of several important steps, including peptide binding to the membrane from aqueous solution, peptide folding on the membrane surface, helix insertion into the membrane, and helix-helix association inside the membrane. Herein, we provide a brief review of these studies and also suggest new initiation and probing methods that could lead to improved temporal and structural resolution in future experiments.

A growing toolbox of techniques for studying β-barrel outer membrane protein folding and biogenesis

Great strides into understanding protein folding have been made since the seminal work of Anfinsen over 40 years ago, but progress in the study of membrane protein folding has lagged behind that of their water soluble counterparts. Researchers in these fields continue to turn to more advanced techniques such as NMR, mass spectrometry, molecular dynamics (MD) and single molecule methods to interrogate how proteins fold. Our understanding of β-barrel outer membrane protein (OMP) folding has benefited from these advances in the last decade. This class of proteins must traverse the periplasm and then insert into an asymmetric lipid membrane in the absence of a chemical energy source. In this review we discuss old, new and emerging techniques used to examine the process of OMP folding and biogenesis in vitro and describe some of the insights and new questions these techniques have revealed.

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.

Mechanistic studies of the biogenesis and folding of outer membrane proteins in vitro and in vivo: what have we learned to date?

Research into the mechanisms by which proteins fold into their native structures has been on-going since the work of Anfinsen in the 1960s. Since that time, the folding mechanisms of small, water-soluble proteins have been well characterised. By contrast, progress in understanding the biogenesis and folding mechanisms of integral membrane proteins has lagged significantly because of the need to create a membrane mimetic environment for folding studies in vitro and the difficulties in finding suitable conditions in which reversible folding can be achieved. Improved knowledge of the factors that promote membrane protein folding and disfavour aggregation now allows studies of folding into lipid bilayers in vitro to be performed. Consequently, mechanistic details and structural information about membrane protein folding are now emerging at an ever increasing pace. Using the panoply of methods developed for studies of the folding of water-soluble proteins. This review summarises current knowledge of the mechanisms of outer membrane protein biogenesis and folding into lipid bilayers in vivo and in vitro and discusses the experimental techniques utilised to gain this information. The emerging knowledge is beginning to allow comparisons to be made between the folding of membrane proteins with current understanding of the mechanisms of folding of water-soluble proteins.

Bacteriorhodopsin folds through a poorly organized transition state

The folding mechanisms of helical membrane proteins remain largely uncharted. Here we characterize the kinetics of bacteriorhodopsin folding and employ φ-value analysis to explore the folding transition state. First, we developed and confirmed a kinetic model that allowed us to assess the rate of folding from SDS-denatured bacteriorhodopsin (bRU) and provides accurate thermodynamic information even under influence of retinal hydrolysis. Next, we obtained reliable φ-values for 16 mutants of bacteriorhodopsin with good coverage across the protein. Every φ-value was less than 0.4, indicating the transition state is not uniquely structured. We suggest that the transition state is a loosely organized ensemble of conformations.

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.

Predictive energy landscapes for folding α-helical transmembrane proteins

We explore the hypothesis that the folding landscapes of membrane proteins are funneled once the proteins' topology within the membrane is established. We extend a protein folding model, the associative memory, water-mediated, structure, and energy model (AWSEM) by adding an implicit membrane potential and reoptimizing the force field to account for the differing nature of the interactions that stabilize proteins within lipid membranes, yielding a model that we call AWSEM-membrane. Once the protein topology is set in the membrane, hydrophobic attractions play a lesser role in finding the native structure, whereas polar-polar attractions are more important than for globular proteins. We examine both the quality of predictions made with AWSEM-membrane when accurate knowledge of the topology and secondary structure is available and the quality of predictions made without such knowledge, instead using bioinformatically inferred topology and secondary structure based on sequence alone. When no major errors are made by the bioinformatic methods used to assign the topology of the transmembrane helices, these two types of structure predictions yield roughly equivalent quality structures. Although the predictive energy landscape is transferable and not structure based, within the correct topological sector we find the landscape is indeed very funneled: Thermodynamic landscape analysis indicates that both the total potential energy and the contact energy decrease as native contacts are formed. Nevertheless the near symmetry of different helical packings with respect to native contact formation can result in multiple packings with nearly equal thermodynamic occupancy, especially at temperatures just below collapse.

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.

Relative domain folding and stability of a membrane transport protein

There is a limited understanding of the folding of multidomain membrane proteins. Lactose permease (LacY) of Escherichia coli is an archetypal member of the major facilitator superfamily of membrane transport proteins, which contain two domains of six transmembrane helices each. We exploit chemical denaturation to determine the unfolding free energy of LacY and employ Trp residues as site-specific thermodynamic probes. Single Trp LacY mutants are created with the individual Trps situated at mirror image positions on the two LacY domains. The changes in Trp fluorescence induced by urea denaturation are used to construct denaturation curves from which unfolding free energies can be determined. The majority of the single Trp tracers report the same stability and an unfolding free energy of approximately +2 kcal mol(-1). There is one exception; the fluorescence of W33 at the cytoplasmic end of helix I on the N domain is unaffected by urea. In contrast, the equivalent position on the first helix, VII, of the C-terminal domain exhibits wild-type stability, with the single Trp tracer at position 243 on helix VII reporting an unfolding free energy of +2 kcal mol(-1). This indicates that the region of the N domain of LacY at position 33 on helix I has enhanced stability to urea, when compared the corresponding location at the start of the C domain. We also find evidence for a potential network of stabilising interactions across the domain interface, which reduces accessibility to the hydrophilic substrate binding pocket between the two domains.

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.

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.

Malleability of the folding mechanism of the outer membrane protein PagP: parallel pathways and the effect of membrane elasticity

Understanding the interactions between membrane proteins and the lipid bilayer is key to increasing our ability to predict and tailor the folding mechanism, structure and stability of membrane proteins. Here, we have investigated the effects of changing the membrane composition and the relative concentrations of protein and lipid on the folding mechanism of the bacterial outer membrane protein PagP. The folding pathway, monitored by tryptophan fluorescence, was found to be characterized by a burst phase, representing PagP adsorption to the liposome surface, followed by a time course that reflects the folding and insertion of the protein into the membrane. In 1,2-dilauroyl-sn-glycero-3-phosphocholine (diC(12:0)PC) liposomes, the post-adsorption time course fits well to a single exponential at high lipid-to-protein ratios (LPRs), but at low LPRs, a second exponential phase with a slower folding rate constant is observed. Interrupted refolding assays demonstrated that the two exponential phases reflect the presence of parallel folding pathways. Partitioning between these pathways was found to be modulated by the elastic properties of the membrane. Folding into mixed 1,2-dilauroyl-sn-glycero-3-phosphoethanolamine:diC(12:0)PC liposomes resulted in a decrease in PagP adsorption to the liposomes and a switch to the slower folding pathway. By contrast, inclusion of 1,2-dilauroyl-sn-glycero-3-phosphoserine into diC(12:0)PC liposomes resulted in a decrease in the folding rate of the fast pathway. The results highlight the effect of lipid composition in tailoring the folding mechanism of a membrane protein, revealing that membrane proteins have access to multiple, competing folding routes to a unique native structure.

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.

Hydrogen exchange mass spectrometry of bacteriorhodopsin reveals light-induced changes in the structural dynamics of a biomolecular machine

Many proteins act as molecular machines that are fuelled by a nonthermal energy source. Examples include transmembrane pumps and stator-rotor complexes. These systems undergo cyclic motions (CMs) that are being driven along a well-defined conformational trajectory. Superimposed on these CMs are thermal fluctuations (TFs) that are coupled to stochastic motions of the solvent. Here we explore whether the TFs of a molecular machine are affected by the occurrence of CMs. Bacteriorhodopsin (BR) is a light-driven proton pump that serves as a model system in this study. The function of BR is based on a photocycle that involves trans/cis isomerization of a retinal chromophore, as well as motions of transmembrane helices. Hydrogen/deuterium exchange (HDX) mass spectrometry was used to monitor the TFs of BR, focusing on the monomeric form of the protein. Comparative HDX studies were conducted under illumination and in the dark. The HDX kinetics of BR are dramatically accelerated in the presence of light. The isotope exchange rates and the number of backbone amides involved in EX2 opening transitions increase roughly 2-fold upon illumination. In contrast, light/dark control experiments on retinal-free protein produced no discernible differences. It can be concluded that the extent of TFs in BR strongly depends on photon-driven CMs. The light-induced differences in HDX behavior are ascribed to protein destabilization. Specifically, the thermodynamic stability of the dark-adapted protein is estimated to be 5.5 kJ mol(-1) under the conditions of our work. This value represents the free energy difference between the folded state F and a significantly unfolded conformer U. Illumination reduces the stability of F by 2.2 kJ mol(-1). Mechanical agitation caused by isomerization of the chromophore is transferred to the surrounding protein scaffold, and subsequently, the energy dissipates into the solvent. Light-induced retinal motions therefore act analogously to an internal heat source that promotes the occurrence of TFs. Overall, our data highlight the potential of HDX methods for probing the structural dynamics of molecular machines under "engine on" and "engine off" conditions.

Folding and stability of membrane transport proteins in vitro

Transmembrane transporters are responsible for maintaining a correct internal cellular environment. The inherent flexibility of transporters together with their hydrophobic environment means that they are challenging to study in vitro, but recently significant progress been made. This review will focus on in vitro stability and folding studies of transmembrane alpha helical transporters, including reversible folding systems and thermal denaturation. The successful re-assembly of a small number of ATP binding cassette transporters is also described as this is a significant step forward in terms of understanding the folding and assembly of these more complex, multi-subunit proteins. The studies on transporters discussed here represent substantial advances for membrane protein studies as well as for research into protein folding. The work demonstrates that large flexible hydrophobic proteins are within reach of in vitro folding studies, thus holding promise for furthering knowledge on the structure, function and biogenesis of ubiquitous membrane transporter families. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Thermodynamic stability of bacteriorhodopsin mutants measured relative to the bacterioopsin unfolded state

The stability of bacteriorhodopsin (bR) has often been assessed using SDS unfolding assays that monitor the transition of folded bR (bR(f)) to unfolded (bR(u)). While many criteria suggest that the unfolding curves reflect thermodynamic stability, slow retinal (RET) hydrolysis during refolding makes it impossible to perform the most rigorous test for equilibrium, i.e., superimposable unfolding and refolding curves. Here we made a new equilibrium test by asking whether the refolding rate in the transition zone is faster than RET hydrolysis. We find that under conditions we have used previously, refolding is in fact slower than hydrolysis, strongly suggesting that equilibrium is not achieved. Instead, the apparent free energy values reported previously are dominated by unfolding rates. To assess how different the true equilibrium values are, we employed an alternative method by measuring the transition of bR(f) to unfolded bacterioopsin (bO(u)), the RET-free form of unfolded protein. The bR(f)-to-bO(u) transition is fully reversible, particular when we add excess RET. We compared the difference in unfolding free energies for 13 bR mutants measured by both assays. For 12 of the 13 mutants with a wide range of stabilities, the results are essentially the same within experimental error. The congruence of the results is fortuitous and suggests the energetic effects of most mutations may be focused on the folded state. The bR(f)-to-bO(u) reaction is inconvenient because many days are required to reach equilibrium, but it is the preferable measure of thermodynamic stability. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Linear rate-equilibrium relations arising from ion channel-bilayer energetic coupling

Linear rate-equilibrium (RE) relations, also known as linear free energy relations, are widely observed in chemical reactions, including protein folding, enzymatic catalysis, and channel gating. Despite the widespread occurrence of linear RE relations, the principles underlying the linear relation between changes in activation and equilibrium energy in macromolecular reactions remain enigmatic. When examining amphiphile regulation of gramicidin channel gating in lipid bilayers, we noted that the gating process could be described by a linear RE relation with a simple geometric interpretation. This description is possible because the gating process provides a well-understood reaction, in which structural changes in a bilayer-embedded model protein can be studied at the single-molecule level. It is thus possible to obtain quantitative information about the energetics of the reaction transition state and its position on a spatial coordinate. It turns out that the linear RE relation for the gramicidin monomer-dimer reaction can be understood, and the quantitative relation between changes in activation energy and equilibrium energy can be interpreted, by considering the effects of amphiphiles on the changes in bilayer elastic energy associated with channel gating. We are not aware that a similar simple mechanistic explanation of a linear RE relation has been provided for a chemical reaction in a macromolecule. RE relations generally should be useful for examining how amphiphile-induced changes in bilayer properties modulate membrane protein folding and function, and for distinguishing between direct (e.g., due to binding) and indirect (bilayer-mediated) effects.

Kinetic folding mechanism of an integral membrane protein examined by pulsed oxidative labeling and mass spectrometry

We report the application of pulsed oxidative labeling for deciphering the folding mechanism of a membrane protein. SDS-denatured bacteriorhodopsin (BR) was refolded by mixing with bicelles in the presence of free retinal. At various time points (20 ms to 1 day), the protein was exposed to a microsecond ·OH pulse that induces oxidative modifications at solvent-accessible methionine side chains. The extent of labeling was determined by mass spectrometry. These measurements were complemented by stopped-flow spectroscopy. Major time-dependent changes in solvent accessibility were detected for M20 (helix A) and M118 (helix D). Our kinetic data indicate a sequential folding mechanism, consistent with models previously suggested by others on the basis of optical data. Yet, ·OH labeling provides additional structural insights. An initial folding intermediate I(1) gets populated within 20 ms, concomitantly with formation of helix A. Subsequent structural consolidation leads to a transient species I(2). Noncovalent retinal binding to I(2) induces folding of helix D, thereby generating an intermediate I(R). In the absence of retinal, the latter transition does not take place. Hence, formation of helix D depends on retinal binding, whereas this is not the case for helix A. As the cofactor settles deeper into its binding pocket, a final transient species I(R) is generated. This intermediate converts into native BR within minutes by formation of the retinal-K216 Schiff base linkage. The combination of pulsed covalent labeling and optical spectroscopy employed here should also be suitable for exploring the folding mechanisms of other membrane proteins.

A folding pathway-dependent score to recognize membrane proteins

While various approaches exist to study protein localization, it is still a challenge to predict where proteins localize. Here, we consider a mechanistic viewpoint for membrane localization. Taking into account the steps for the folding pathway of α-helical membrane proteins and relating biophysical parameters to each of these steps, we create a score capable of predicting the propensity for membrane localization and call it FP(3)mem. This score is driven from the principal component analysis (PCA) of the biophysical parameters related to membrane localization. FP(3)mem allows us to rationalize the colocalization of a number of channel proteins with the Cav1.2 channel by their fewer propensities for membrane localization.

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.

Unfolding free energy of a two-domain transmembrane sugar transport protein

Understanding how an amino acid sequence folds into a functional, three-dimensional structure has proved to be a formidable challenge in biological research, especially for transmembrane proteins with multiple alpha helical domains. Mechanistic folding studies on helical membrane proteins have been limited to unusually stable, single domain proteins such as bacteriorhodopsin. Here, we extend such work to flexible, multidomain proteins and one of the most widespread membrane transporter families, the major facilitator superfamily, thus showing that more complex membrane proteins can be successfully refolded to recover native substrate binding. We determine the unfolding free energy of the two-domain, Escherichia coli galactose transporter, GalP; a bacterial homologue of human glucose transporters. GalP is reversibly unfolded by urea. Urea causes loss of substrate binding and a significant reduction in alpha helical content. Full recovery of helical structure and substrate binding occurs in dodecylmaltoside micelles, and the unfolding free energy can be determined. A linear dependence of this free energy on urea concentration allows the free energy of unfolding in the absence of urea to be determined as +2.5 kcal·mol(-1). Urea has often been found to be a poor denaturant for transmembrane helical structures. We attribute the denaturation of GalP helices by urea to the dynamic nature of the transporter structure allowing denaturant access via the substrate binding pocket, as well as to helical structure that extends beyond the membrane. This study gives insight into the final, critical folding step involving recovery of ligand binding for a multidomain membrane transporter.

The contribution of a covalently bound cofactor to the folding and thermodynamic stability of an integral membrane protein

The factors controlling the stability, folding, and dynamics of integral membrane proteins are not fully understood. The high stability of the membrane protein bacteriorhodopsin (bR), an archetypal member of the rhodopsin photoreceptor family, has been ascribed to its covalently bound retinal cofactor. We investigate here the role of this cofactor in the thermodynamic stability and folding kinetics of bR. Multiple spectroscopic probes were used to determine the kinetics and energetics of protein folding in mixed lipid/detergent micelles in the presence and absence of retinal. The presence of retinal increases extrapolated values for the overall unfolding free energy from 6.3 ± 0.4 kcal mol(-1) to 23.4 ± 1.5 kcal mol(-1) at zero denaturant, suggesting that the cofactor contributes 17.1 kcal mol(-1) towards the overall stability of bR. In addition, the cooperativity of equilibrium unfolding curves is markedly reduced in the absence of retinal with overall m-values decreasing from 31.0 ± 2.0 kcal mol(-1) to 10.9 ± 1.0 kcal mol(-1), indicating that the folded state of the apoprotein is less compact than the equivalent for the holoprotein. This change in the denaturant response means that the difference in the unfolding free energy at a denaturant concentration midway between the two unfolding curves is only ca 3-6 kcal mol(-1). Kinetic data show that the decrease in stability upon removal of retinal is associated with an increase in the apparent intrinsic rate constant of unfolding, k(u)(H2O), from ~1 × 10(-16) s(-1) to ~1 × 10(-4) s(-1) at 25 °C. This correlates with a decrease in the unfolding activation energy by 16.3 kcal mol(-1) in the apoprotein, extrapolated to zero SDS. These results suggest that changes in bR stability induced by retinal binding are mediated solely by changes in the activation barrier for unfolding. The results are consistent with a model in which bR is kinetically stabilized via a very slow rate of unfolding arising from protein-retinal interactions that increase the rigidity and compactness of the polypeptide chain.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.