From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Dark nanodiscs for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern the structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently developed nonfluorescent membrane scaffold protein to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs loaded with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG) lipids show a complex biphasic thermal unfolding pattern with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5°C and 77.5°C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a DMPC and DMPG nanodisc. To extend the utility of this method, we evaluate and compare the thermostability of DsbB in different lipid environments. We introduce this method as a new tool that can be used to understand how compositionally and biophysically complex lipid environments drive membrane protein stability.

Dark nanodiscs as a model membrane for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently-developed non-fluorescent membrane scaffold protein (MSP) to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs show a complex biphasic thermal unfolding pattern in the presence of lipids with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5 °C and 77.5 °C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a nanodisc. We introduce this method as a new tool that can be used to understand how compositionally, and biophysically complex lipid environments drive membrane protein stability.

Observation of Strong Synergy in the Interfacial Water Response of Binary Ionic and Nonionic Surfactant Mixtures

Interfacial vibrational footprints of the binary mixture of sodium dodecyl sulfate (SDS) and hexaethylene glycol monododecyl ether (C₁₂E₆) were probed using heterodyne detected vibrational sum frequency generation (HDVSFG). Our results show that in the presence of C₁₂E₆ at CMC (70 μM) the effect of SDS on the orientation of interfacial water molecules is enhanced 10 times compared to just pure surfactants. The experimental results contest the traditional Langmuir adsorption model predictions. This is also evidenced by our molecular dynamics simulations that show a remarkable restructuring and enhanced orientation of the interfacial water molecules upon DS^- adsorption to the C₁₂E₆ surface. The simulations show that the adsorption free energy of DS^- ions to a water surface covered with C₁₂E₆ is an enthalpy-driven process and more attractive by ∼10 k_BT compared to the adsorption energy of DS^- to the surface of pure water.

How physical forces drive the process of helical membrane protein folding

Protein folding is a fundamental process of life with important implications throughout biology. Indeed, tens of thousands of mutations have been associated with diseases, and most of these mutations are believed to affect protein folding rather than function. Correct folding is also a key element of design. These factors have motivated decades of research on protein folding. Unfortunately, knowledge of membrane protein folding lags that of soluble proteins. This gap is partly caused by the greater technical challenges associated with membrane protein studies, but also because of additional complexities. While soluble proteins fold in a homogenous water environment, membrane proteins fold in a setting that ranges from bulk water to highly charged to apolar. Thus, the forces that drive folding vary in different regions of the protein, and this complexity needs to be incorporated into our understanding of the folding process. Here, we review our understanding of membrane protein folding biophysics. Despite the greater challenge, better model systems and new experimental techniques are starting to unravel the forces and pathways in membrane protein folding.

Folding and misfolding of potassium channel monomers during assembly and tetramerization

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τ_fold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization.

Cys-labeling kinetics of membrane protein GlpG: a role for specific SDS binding and micelle changes?

Empirically, α-helical membrane protein folding stability in surfactant micelles can be tuned by varying the mole fraction MF^SDS of anionic (sodium dodecyl sulfate (SDS)) relative to nonionic (e.g., dodecyl maltoside (DDM)) surfactant, but we lack a satisfying physical explanation of this phenomenon. Cysteine labeling (CL) has thus far only been used to study the topology of membrane proteins, not their stability or folding behavior. Here, we use CL to investigate membrane protein folding in mixed DDM-SDS micelles. Labeling kinetics of the intramembrane protease GlpG are consistent with simple two-state unfolding-and-exchange rates for seven single-Cys GlpG variants over most of the explored MF^SDS range, along with exchange from the native state at low MF^SDS (which inconveniently precludes measurement of unfolding kinetics under native conditions). However, for two mutants, labeling rates decline with MF^SDS at 0-0.2 MF^SDS (i.e., native conditions). Thus, an increase in MF^SDS seems to be a protective factor for these two positions, but not for the five others. We propose different scenarios to explain this and find the most plausible ones to involve preferential binding of SDS monomers to the site of CL (based on computational simulations) along with changes in size and shape of the mixed micelle with changing MF^SDS (based on SAXS studies). These nonlinear impacts on protein stability highlights a multifaceted role for SDS in membrane protein denaturation, involving both direct interactions of monomeric SDS and changes in micelle size and shape along with the general effects on protein stability of changes in micelle composition.

A complete picture of protein unfolding and refolding in surfactants

Interactions between proteins and surfactants are of relevance in many applications including food, washing powder formulations, and drug formulation. The anionic surfactant sodium dodecyl sulfate (SDS) is known to unfold globular proteins, while the non-ionic surfactant octaethyleneglycol monododecyl ether (C₁₂E₈) can be used to refold proteins from their SDS-denatured state. While unfolding have been studied in detail at the protein level, a complete picture of the interplay between protein and surfactant in these processes is lacking. This gap in our knowledge is addressed in the current work, using the β-sheet-rich globular protein β-lactoglobulin (bLG). We combined stopped-flow time-resolved SAXS, fluorescence, and circular dichroism, respectively, to provide an unprecedented in-depth picture of the different steps involved in both protein unfolding and refolding in the presence of SDS and C₁₂E₈. During unfolding, core-shell bLG-SDS complexes were formed within ∼10 ms. This involved an initial rapid process where protein and SDS formed aggregates, followed by two slower processes, where the complexes first disaggregated into single protein structures situated asymmetrically on the SDS micelles, followed by isotropic redistribution of the protein. Refolding kinetics (>100 s) were slower than unfolding (<30 s), and involved rearrangements within the mixing deadtime (∼5 ms) and transient accumulation of unfolded monomeric protein, differing in structure from the original bLG-SDS structure. Refolding of bLG involved two steps: extraction of most of the SDS from the complexes followed by protein refolding. These results reveal that surfactant-mediated unfolding and refolding of proteins are complex processes with rearrangements occurring on time scales from sub-milliseconds to minutes.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

Structure formation during translocon-unassisted co-translational membrane protein folding

Correctly folded membrane proteins underlie a plethora of cellular processes, but little is known about how they fold. Knowledge of folding mechanisms centres on reversible folding of chemically denatured membrane proteins. However, this cannot replicate the unidirectional elongation of the protein chain during co-translational folding in the cell, where insertion is assisted by translocase apparatus. We show that a lipid membrane (devoid of translocase components) is sufficient for successful co-translational folding of two bacterial α-helical membrane proteins, DsbB and GlpG. Folding is spontaneous, thermodynamically driven, and the yield depends on lipid composition. Time-resolving structure formation during co-translational folding revealed different secondary and tertiary structure folding pathways for GlpG and DsbB that correlated with membrane interfacial and biological transmembrane amino acid hydrophobicity scales. Attempts to refold DsbB and GlpG from chemically denatured states into lipid membranes resulted in extensive aggregation. Co-translational insertion and folding is thus spontaneous and minimises aggregation whilst maximising correct folding.

Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability

There is a compelling and growing need to accurately predict the impact of amino acid mutations on protein stability for problems in personalized medicine and other applications. Here the ability of 10 computational tools to accurately predict mutation-induced perturbation of folding stability (ΔΔG) for membrane proteins of known structure was assessed. All methods for predicting ΔΔG values performed significantly worse when applied to membrane proteins than when applied to soluble proteins, yielding estimated concordance, Pearson, and Spearman correlation coefficients of <0.4 for membrane proteins. Rosetta and PROVEAN showed a modest ability to classify mutations as destabilizing (ΔΔG < −0.5 kcal/mol), with a 7 in 10 chance of correctly discriminating a randomly chosen destabilizing variant from a randomly chosen stabilizing variant. However, even this performance is significantly worse than for soluble proteins. This study highlights the need for further development of reliable and reproducible methods for predicting thermodynamic folding stability in membrane proteins.

Topological constraints and modular structure in the folding and functional motions of GlpG, an intramembrane protease

We investigate the folding of GlpG, an intramembrane protease, using perfectly funneled structure-based models that implicitly account for the absence or presence of the membrane. These two models are used to describe, respectively, folding in detergent micelles and folding within a bilayer, which effectively constrains GlpG's topology in unfolded and partially folded states. Structural free-energy landscape analysis shows that although the presence of multiple folding pathways is an intrinsic property of GlpG's modular functional architecture, the large entropic cost of organizing helical bundles in the absence of the constraining bilayer leads to pathways that backtrack (i.e., local unfolding of previously folded substructures is required when moving from the unfolded to the folded state along the minimum free-energy pathway). This backtracking explains the experimental observation of thermodynamically destabilizing mutations that accelerate GlpG's folding in detergent micelles. In contrast, backtracking is absent from the model when folding is constrained within a bilayer, the environment in which GlpG has evolved to fold. We also characterize a near-native state with a highly mobile transmembrane helix 5 (TM5) that is significantly populated under folding conditions when GlpG is embedded in a bilayer. Unbinding of TM5 from the rest of the structure exposes GlpG's active site, consistent with studies of the catalytic mechanism of GlpG that suggest that TM5 serves as a substrate gate to the active site.

Mapping the energy landscape for second-stage folding of a single membrane protein

Membrane proteins are designed to fold and function in a lipid membrane, yet folding experiments within a native membrane environment are challenging to design. Here we show that single-molecule forced unfolding experiments can be adapted to study helical membrane protein folding under native-like bicelle conditions. Applying force using magnetic tweezers, we find that a transmembrane helix protein, Escherichia coli rhomboid protease GlpG, unfolds in a highly cooperative manner, largely unraveling as one physical unit in response to mechanical tension above 25 pN. Considerable hysteresis is observed, with refolding occurring only at forces below 5 pN. Characterizing the energy landscape reveals only modest thermodynamic stability (ΔG = 6.5 kBT) but a large unfolding barrier (21.3 kBT) that can maintain the protein in a folded state for long periods of time (t1/2 ∼3.5 h). The observed energy landscape may have evolved to limit the existence of troublesome partially unfolded states and impart rigidity to the structure.

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.

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.

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.

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.

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.