Membrane proteins can have high kinetic stability

The holocentric chromosome microevolution: From phylogeographic patterns to genomic associations with environmental gradients

Geographic isolation and chromosome evolution are two of the major drivers of diversification in eukaryotes in general, and specifically, in plants. On one hand, range shifts induced by Pleistocene glacial oscillations deeply shaped the evolutionary trajectories of species in the Northern Hemisphere. On the other hand, karyotype variability within species or species complexes may have adaptive potential as different karyotypes may represent different recombination rates and linkage groups that may be associated with locally adapted genes or supergenes. Organisms with holocentric chromosomes are ideal to study the link between local adaptation and chromosome evolution, due to their high cytogenetic variability, especially when it seems to be related to environmental variation. Here, we integrate the study of the phylogeography, chromosomal evolution and ecological requirements of a plant species complex distributed in the Western Euro-Mediterranean region (Carex gr. laevigata, Cyperaceae). We aim to clarify the relative influence of these factors on population differentiation and ultimately on speciation. We obtained a well-resolved RADseq phylogeny that sheds light on the phylogeographic patterns of molecular and chromosome number variation, which are compatible with south-to-north postglacial migration. In addition, landscape genomics analyses identified candidate loci for local adaptation, and also strong significant associations between the karyotype and the environment. We conclude that karyotype distribution in C. gr. laevigata has been constrained by both range shift dynamics and local adaptation. Our study demonstrates that chromosome evolution may be responsible, at least partially, for microevolutionary patterns of population differentiation and adaptation in Carex.

Effects of stand-alone polar residue on membrane protein stability and structure

Helical membrane proteins generally have a hydrophobic nature, with apolar side chains comprising the majority of the transmembrane (TM) helices. However, whenever polar side chains are present in the TM domain, they often exert a crucial role in structural interactions with other polar residues, such as TM helix associations and oligomerization. Moreover, polar residues in the TM region also often participate in protein functions, such as the Schiff base bonding between Lys residues and retinal in rhodopsin-like membrane proteins. Although many studies have focused on these functional polar residues, our understanding of stand-alone polar residues that are energetically unfavored in TM helixes is limited. Here, we adopted bacteriorhodopsin (bR) as a model system and systematically mutated 17 of its apolar Leu or Phe residues to polar Asn. Stability measurements of the resulting mutants revealed that all of these polar substitutions reduced bR stability to various extents, and the extent of destabilization of each mutant bR is also correlated to different structural factors, such as the relative accessible surface area and membrane depth of the mutation site. Structural analyses of these Asn residues revealed that they form sidechain-to-backbone hydrogen bonds that alleviate the unfavorable energetics in hydrophobic and apolar surroundings. Our results indicate that membrane proteins are able to accommodate certain stand-alone polar residues in the TM region without disrupting overall structures.

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.

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

A thermodynamic analysis of CLC transporter dimerization in lipid bilayers

The CLC-ec1 chloride/proton antiporter is a membrane-embedded homodimer with subunits that can dissociate and associate, but the thermodynamic driving forces favor the assembled dimer at biological densities. Yet, the physical reasons for this stability are confounding as dimerization occurs via the burial of hydrophobic interfaces away from the lipid solvent. For binding of nonpolar surfaces in aqueous solution, the driving force is often attributed to the hydrophobic effect, but this should not apply in the membrane since there is very little water. To investigate this further, we quantified the thermodynamic changes associated with CLC dimerization in membranes by carrying out a van 't Hoff analysis of the temperature dependency of the free energy of dimerization, ΔG°. To ensure that the reaction reached equilibrium at different temperatures, we utilized a Förster resonance energy transfer assay to report on relaxation kinetics of subunit exchange as a function of temperature. Equilibration times were then applied to measure CLC-ec1 dimerization isotherms at different temperatures using the single-molecule subunit-capture photobleaching analysis approach. The results demonstrate that the dimerization free energy of CLC in Escherichia coli-like membranes exhibits a nonlinear temperature dependency corresponding to a large, negative change in heat capacity, a signature of solvent ordering effects such as the hydrophobic effect. Consolidating this with our previous molecular analyses suggests that the nonbilayer defect required to solvate the monomeric state is one source of the observed change in heat capacity and indicates the existence of a generalizable driving force for protein association in membranes.

A thermodynamic analysis of CLC transporter dimerization in lipid bilayers

The CLC-ec1 chloride/proton antiporter is a membrane embedded homodimer where subunits can dissociate and associate, but the thermodynamic driving forces favor the assembled form at biological densities. Yet, the physical reasons for this stability are confounding since binding occurs via the burial of hydrophobic protein interfaces yet the hydrophobic effect should not apply since there is little water within the membrane. To investigate this further, we quantified the thermodynamic changes associated with CLC dimerization in membranes by carrying out a van 't Hoff analysis of the temperature dependency of the free energy of dimerization, ΔG° . To ensure that the reaction reached equilibrium under changing conditions, we utilized a Förster Resonance Energy Transfer based assay to report on the relaxation kinetics of subunit exchange as a function of temperature. These equilibration times were then applied to measure CLC-ec1 dimerization isotherms as a function of temperature using the single-molecule subunit-capture photobleaching analysis approach. The results demonstrate that the dimerization free energy of CLC in E. coli membranes exhibits a non-linear temperature dependency corresponding to a large, negative change in heat capacity, a signature of solvent ordering effects including the hydrophobic effect. Consolidating this with our previous molecular analyses suggests that the non-bilayer defect required to solvate the monomeric state is the molecular source of this large change in heat capacity and is a major and generalizable driving force for protein association in membranes.

Untangling the complexity of membrane protein folding

Delineating the folding steps of helical-bundle membrane proteins has been a challenging task. Many questions remain unanswered, including the conformation and stability of the states populated during folding, the shape of the energy barriers between the states, and the role of lipids as a solvent in mediating the folding. Recently, theoretical frames have matured to a point that permits detailed dissection of the folding steps, and advances in experimental techniques at both single-molecule and ensemble levels enable selective modulation of specific steps for quantitative determination of the folding energy landscapes. We also discuss how lipid molecules would play an active role in shaping the folding energy landscape of membrane proteins, and how folding of multi-domain membrane proteins can be understood based on our current knowledge. We conclude this review by offering an outlook for emerging questions in the study of membrane protein folding.

Revealing enzyme functional architecture via high-throughput microfluidic enzyme kinetics

Systematic and extensive investigation of enzymes is needed to understand their extraordinary efficiency and meet current challenges in medicine and engineering. We present HT-MEK (High-Throughput Microfluidic Enzyme Kinetics), a microfluidic platform for high-throughput expression, purification, and characterization of more than 1500 enzyme variants per experiment. For 1036 mutants of the alkaline phosphatase PafA (phosphate-irrepressible alkaline phosphatase of Flavobacterium), we performed more than 670,000 reactions and determined more than 5000 kinetic and physical constants for multiple substrates and inhibitors. We uncovered extensive kinetic partitioning to a misfolded state and isolated catalytic effects, revealing spatially contiguous regions of residues linked to particular aspects of function. Regions included active-site proximal residues but extended to the enzyme surface, providing a map of underlying architecture not possible to derive from existing approaches. HT-MEK has applications that range from understanding molecular mechanisms to medicine, engineering, and design.

Altering CLC stoichiometry by reducing non-polar side-chains at the dimerization interface

CLC-ec1 is a Cl^-/H⁺ antiporter that forms stable homodimers in lipid bilayers, with a free energy of -10.9 kcal/mol in 2:1 POPE/POPG lipid bilayers. The dimerization interface is formed by four transmembrane helices: H, I, P and Q, that are lined by non-polar side-chains that come in close contact, yet it is unclear as to whether their interactions drive dimerization. To investigate whether non-polar side-chains are required for dimer assembly, we designed a series of constructs where side-chain packing in the dimer state is significantly reduced by making 4-5 alanine substitutions along each helix (H-ala, I-ala, P-ala, Q-ala). All constructs are functional and three purify as stable dimers in detergent micelles despite the removal of significant side-chain interactions. On the other hand, H-ala shows the unique behavior of purifying as a mixture of monomers and dimers, followed by a rapid and complete conversion to monomers. In lipid bilayers, all four constructs are monomeric as examined by single-molecule photobleaching analysis. Further study of the H-helix shows that the single mutation L194A is sufficient to yield monomeric CLC-ec1 in detergent micelles and lipid bilayers. X-ray crystal structures of L194A reveal the protein re-assembles to form dimers, with a structure that is identical to wild-type. Altogether, these results demonstrate that non-polar membrane embedded side-chains play an important role in defining dimer stability, but the stoichiometry is highly contextual to the solvent environment. Furthermore, we discovered that L194 is a molecular hot-spot for defining dimerization of CLC-ec1.

Thermodynamic and kinetic stabilities of transmembrane helix bundles as revealed by single-pair FRET analysis: Effects of the number of membrane-spanning segments and cholesterol

The tertiary structures and conformational dynamics of transmembrane (TM) helical proteins are maintained by the interhelical interaction network in membranes, although it is complicated to analyze the underlying driving forces because the amino acid sequences can involve multiple and various types of interactions. To obtain insights into basal and common effects of the number of membrane-spanning segments and membrane cholesterol, we measured stabilities of helix bundles composed of simple TM helices (AALALAA)₃ (1TM) and (AALALAA)₃-G₅-(AALALAA)₃ (2TM). Association-dissociation dynamics for 1TM-1TM, 1TM-2TM, and 2TM-2TM pairs were monitored to compare stabilities of 2-, 3-, and 4-helical bundles, respectively, with single-pair fluorescence resonance energy transfer (sp-FRET) in liposome membranes. Both thermodynamic and kinetic stabilities of the helix bundles increased with a greater number of membrane-spanning segments in POPC. The presence of 30 mol% cholesterol strongly enhanced the formation of 1TM-1TM and 1TM-2TM bundles (~ - 9 kJ mol^-1), whereas it only weakly stabilized the 2TM-2TM bundle (~ - 3 kJ mol^-1). Fourier transform infrared-polarized attenuated total reflection (ATR-FTIR) spectroscopy revealed an ~30° tilt of the helix axis relative to bilayer normal for the 1TM-2TM pair in the presence of cholesterol, suggesting the formation of a tilted helix bundle to release high lateral pressure at the center of cholesterol-containing membranes. These results demonstrate that the number of membrane-spanning segments affects the stability and structure of the helix bundle, and their cholesterol-dependences. Such information is useful to understand the basics of folding and assembly of multispanning TM proteins.

Molecular Switch between Structural Compaction and Thermodynamic Stability by the Xxx-Pro Interface in Transmembrane β-Barrels

Transmembrane β-barrel scaffolds found in outer membrane proteins are formed and stabilized by a defined pattern of interstrand intraprotein H-bonds, in hydrophobic lipid bilayers. Introducing the conformationally constrained proline in β-barrels can cause significant destabilization of these structural regions that require H-bonding, with proline additionally acting as a secondary structure breaker. Membrane protein β-barrels are therefore expected to show poor tolerance to the presence of a transmembrane proline. Here, we assign a thermodynamic measure for the extent to which a single proline can be tolerated at the C-terminal interface of the model transmembrane β-barrel PagP. We find that proline drastically destabilizes PagP by 7.0 kcal mol^-1 with respect to wild-type PagP (F¹⁶¹ → P¹⁶¹). Interestingly, strategic modulation of the preceding residue can modify the measured energetics. Placing a hydrophobic or bulky side chain as the preceding residue increases the thermodynamic stability by ≤8.0 kcal mol^-1. While polar substituents at the preceding residue decrease the PagP stability, these residues demonstrate stronger tertiary packing interactions in the barrel and retain the catalytic activity of native PagP. This biophysical interplay between enhanced thermodynamic stability and attaining a structurally compact functional β-barrel scaffold highlights the detrimental effect caused by proline incorporation. Our findings also reveal alternative mechanisms that protein sequences can employ to salvage the structural integrity of transmembrane protein structures.

Reversible folding energetics of Yersinia Ail barrel reveals a hyperfluorescent intermediate

Deducing the molecular details of membrane protein folding has lately become an important area of research in biology. Using Ail, an outer membrane protein (OMP) from Yersina pestis as our model, we explore details of β-barrel folding, stability, and unfolding. Ail displays a simple transmembrane β-barrel topology. Here, we find that Ail follows a simple two-state mechanism in its folding and unfolding thermodynamics. Interestingly, Ail displays multi-step folding kinetics. The early kinetic intermediates in the folding pathway populate near the unfolded state (β_T ≈ 0.20), and do not display detectable changes in the local environment of the two interface indoles. Interestingly, tryptophans regulate the late events of barrel rearrangement, and Ail thermodynamic stability. We show that W¹⁴⁹ → Y/F/A substitution destabilizes Ail by ~0.13-1.7 kcal mol^-1, but retains path-independent thermodynamic equilibrium of Ail. In surprising contrast, substituting W⁴² and retaining W¹⁴⁹ shifts the thermodynamic equilibrium to an apparent kinetic retardation of only the unfolding process, which gives rise to an associated increase in scaffold stability by ~0.3-1.1 kcal mol^-1. This is accompanied by the formation of an unusual hyperfluorescent state in the unfolding pathway that is more structured, and represents a conformationally dynamic unfolding intermediate with the interface W¹⁴⁹ now lipid solvated. The defined role of each tryptophan and poorer folding efficiency of Trp mutants together presents compelling evidence for the importance of interface aromatics in the unique (un)folding pathway of Ail, and offers interesting insight on alternative pathways in generalized OMP assembly and unfolding mechanisms.

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.

Comparison of in vitro and in vivo oligomeric states of a wild type and mutant trimeric inner membrane multidrug transporter

Many membrane proteins exist and function as oligomers or protein complexes. Routine analytical methods involve extraction and solubilization of the proteins with detergents, which could disturb their actual oligomeric state. AcrB is a trimeric inner membrane multidrug transporter in E. coli. In previous studies, we created a mutant AcrBP223G, which behaves like a monomer when extracted from the cell membrane. However, the actual oligomeric state of AcrBP223G in cell membranes remained unclear, which complicated the interpretation of the mechanism by which the mutation affects function. Here we used several complementary methods to determine the oligomeric state of AcrBP223G in E. coli cell membranes. Two sets of quantitative fluorescent techniques were exploited. For these, we created fluorescent tagged AcrB, AcrB-CFP and AcrB-YPet. Fluorescence resonance energy transfer (FRET) and fluorescence recovery after photobleaching (FRAP) were employed to characterize independently the efficiency of energy transfer between co-expressed AcrB-CFP and AcrB-YPet, and the diffusion coefficient of AcrB-YPet and AcrBP223G-YPet in live E. coli cells. Second, we introduced Cys pairs at the inter-subunit interface and used controlled oxidation to probe inter-subunit distances. The results from all studies converge on the conclusion that AcrBP223G exists as a trimer in cell membranes, which dissociates during the purification steps. The small change in trimer affinity and structure leads to a significant loss of AcrB activity. In addition, throughout this study we developed protocols and established benchmark values, useful for further studies on membrane protein associations in cell membranes.

Unfolding of a ClC chloride transporter retains memory of its evolutionary history

ClC chloride channels and transporters are important for chloride homeostasis in species from bacteria to human. Mutations in ClC proteins cause genetically inherited diseases, some of which are likely to have folding defects. The ClC proteins present a challenging and unusual biological folding problem because they are large membrane proteins possessing a complex architecture with many re-entrant helices that go only part way through membrane and loop back out. Here we were able to examine the unfolding of the E. coli ClC transporter, ClC-ec1, using single-molecule forced unfolding methods. We find that the protein can be separated into two stable halves that unfold independently. The independence of the two domains is consistent with an evolutionary model in which the two halves arose from independent folding subunits that later fused together. Maintaining smaller folding domains of lesser complexity within large membrane proteins may be an advantageous strategy to avoid misfolding traps.

Folding-Degradation Relationship of a Membrane Protein Mediated by the Universally Conserved ATP-Dependent Protease FtsH

ATP-dependent protein degradation mediated by AAA+ proteases is one of the major cellular pathways for protein quality control and regulation of functional networks. While a majority of studies of protein degradation have focused on water-soluble proteins, it is not well understood how membrane proteins with abnormal conformation are selectively degraded. The knowledge gap stems from the lack of an in vitro system in which detailed molecular mechanisms can be studied as well as difficulties in studying membrane protein folding in lipid bilayers. To quantitatively define the folding-degradation relationship of membrane proteins, we reconstituted the degradation using the conserved membrane-integrated AAA+ protease FtsH as a model degradation machine and the stable helical-bundle membrane protein GlpG as a model substrate in the lipid bilayer environment. We demonstrate that FtsH possesses a substantial ability to actively unfold GlpG, and the degradation significantly depends on the stability and hydrophobicity near the degradation marker. We find that FtsH hydrolyzes 380-550 ATP molecules to degrade one copy of GlpG. Remarkably, FtsH overcomes the dual-energetic burden of substrate unfolding and membrane dislocation with the ATP cost comparable to that for water-soluble substrates by robust ClpAP/XP proteases. The physical principles elucidated in this study provide general insights into membrane protein degradation mediated by ATP-dependent proteolytic systems.

Biological Roles of Protein Kinetic Stability

A protein's stability may range from nonexistent, as in the case of intrinsically disordered proteins, to very high, as indicated by a protein's resistance to degradation, even under relatively harsh conditions. The stability of this latter group is usually under kinetic control because of a high activation energy for unfolding that virtually traps the protein in a specific conformation, thereby conferring resistance to proteolytic degradation and misfolding aggregation. The usual outcome of kinetic stability is a longer protein half-life. Thus, the protective role of protein kinetic stability is often appreciated, but relatively little is known about the extent of biological roles related to this property. In this Perspective, we will discuss several known or putative biological roles of protein kinetic stability, including protection from stressors to avoid aggregation or premature degradation, achieving long-term phenotypic change, and regulating cellular processes by controlling the trigger and timing of molecular motion. The picture that emerges from this analysis is that protein kinetic stability is involved in a myriad of known and yet to be discovered biological functions via its ability to confer degradation resistance and control the timing, extent, and permanency of molecular motion.

Interplay between Protein Thermal Flexibility and Kinetic Stability

Kinetic stability is a key parameter to comprehend protein behavior and it plays a central role to understand how evolution has reached the balance between function and stability in cell-relevant timescales. Using an approach that includes simulations, protein engineering, and calorimetry, we show that there is a clear correlation between kinetic stability determined by differential scanning calorimetry and protein thermal flexibility obtained from a novel method based on temperature-induced unfolding molecular dynamics simulations. Thermal flexibility quantitatively measures the increment of the conformational space available to the protein when energy in provided. The (β/α)₈ barrel fold of two closely related by evolution triosephosphate isomerases from two trypanosomes are used as model systems. The kinetic stability-thermal flexibility correlation has predictive power for the studied proteins, suggesting that the strategy and methodology discussed here might be applied to other proteins in biotechnological developments, evolutionary studies, and the design of protein based therapeutics.

Transmembrane β-barrels: Evolution, folding and energetics

The biogenesis of transmembrane β-barrels (outer membrane proteins, or OMPs) is an elaborate multistep orchestration of the nascent polypeptide with translocases, barrel assembly machinery, and helper chaperone proteins. Several theories exist that describe the mechanism of chaperone-assisted OMP assembly in vivo and unassisted (spontaneous) folding in vitro. Structurally, OMPs of bacterial origin possess even-numbered strands, while mitochondrial β-barrels are even- and odd-stranded. Several underlying similarities between prokaryotic and eukaryotic β-barrels and their folding machinery are known; yet, the link in their evolutionary origin is unclear. While OMPs exhibit diversity in sequence and function, they share similar biophysical attributes and structure. Similarly, it is important to understand the intricate OMP assembly mechanism, particularly in eukaryotic β-barrels that have evolved to perform more complex functions. Here, we deliberate known facets of β-barrel evolution, folding, and stability, and attempt to highlight outstanding questions in β-barrel biogenesis and proteostasis.

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.

Applications of Single-Molecule Methods to Membrane Protein Folding Studies

Protein folding is a fundamental life process with many implications throughout biology and medicine. Consequently, there have been enormous efforts to understand how proteins fold. Almost all of this effort has focused on water-soluble proteins, however, leaving membrane proteins largely wandering in the wilderness. The neglect has occurred not because membrane proteins are unimportant but rather because they present many theoretical and technical complications. Indeed, quantitative membrane protein folding studies are generally restricted to a handful of well-behaved proteins. Single-molecule methods may greatly alter this picture, however, because the ability to work at or near infinite dilution removes aggregation problems, one of the main technical challenges of membrane protein folding studies.

GXXXG-Mediated Parallel and Antiparallel Dimerization of Transmembrane Helices and Its Inhibition by Cholesterol: Single-Pair FRET and 2D IR Studies

Small-residue-mediated interhelical packings are ubiquitously found in helical membrane proteins, although their interaction dynamics and lipid dependence remain mostly uncharacterized. We used a single-pair FRET technique to examine the effect of a GXXXG motif on the association of de novo designed (AALALAA)₃ helices in liposomes. Dimerization occurred with sub-second lifetimes, which was abolished by cholesterol. Utilizing the nearly instantaneous time-resolution of 2D IR spectroscopy, parallel and antiparallel helix associations were identified by vibrational couplings across helices at their interface. Taken together, the data illustrate that the GXXXG motif controls helix packing but still allows for a dynamic and lipid-regulated oligomeric state.

The dimerization equilibrium of a ClC Cl(-)/H(+) antiporter in lipid bilayers

Interactions between membrane protein interfaces in lipid bilayers play an important role in membrane protein folding but quantification of the strength of these interactions has been challenging. Studying dimerization of ClC-type transporters offers a new approach to the problem, as individual subunits adopt a stable and functionally verifiable fold that constrains the system to two states - monomer or dimer. Here, we use single-molecule photobleaching analysis to measure the probability of ClC-ec1 subunit capture into liposomes during extrusion of large, multilamellar membranes. The capture statistics describe a monomer to dimer transition that is dependent on the subunit/lipid mole fraction density and follows an equilibrium dimerization isotherm. This allows for the measurement of the free energy of ClC-ec1 dimerization in lipid bilayers, revealing that it is one of the strongest membrane protein complexes measured so far, and introduces it as new type of dimerization model to investigate the physical forces that drive membrane protein association in membranes.

Selection for Protein Kinetic Stability Connects Denaturation Temperatures to Organismal Temperatures and Provides Clues to Archaean Life

The relationship between the denaturation temperatures of proteins (T_m values) and the living temperatures of their host organisms (environmental temperatures: T_ENV values) is poorly understood. Since different proteins in the same organism may show widely different T_m’s, no simple universal relationship between T_m and T_ENV should hold, other than T_m≥T_ENV. Yet, when analyzing a set of homologous proteins from different hosts, T_m’s are oftentimes found to correlate with T_ENV’s but this correlation is shifted upward on the T_m axis. Supporting this trend, we recently reported T_m’s for resurrected Precambrian thioredoxins that mirror a proposed environmental cooling over long geological time, while remaining a shocking ~50°C above the proposed ancestral ocean temperatures. Here, we show that natural selection for protein kinetic stability (denaturation rate) can produce a T_m↔T_ENV correlation with a large upward shift in T_m. A model for protein stability evolution suggests a link between the T_m shift and the in vivo lifetime of a protein and, more specifically, allows us to estimate ancestral environmental temperatures from experimental denaturation rates for resurrected Precambrian thioredoxins. The T_ENV values thus obtained match the proposed ancestral ocean cooling, support comparatively high Archaean temperatures, and are consistent with a recent proposal for the environmental temperature (above 75°C) that hosted the last universal common ancestor. More generally, this work provides a framework for understanding how features of protein stability reflect the environmental temperatures of the host organisms.

Steric trapping reveals a cooperativity network in the intramembrane protease GlpG

Membrane proteins are assembled through balanced interactions among protein, lipids and water. Studying their folding while maintaining the native lipid environment is necessary but challenging. Here we present methods for analyzing key elements in membrane protein folding including thermodynamic stability, compactness of the unfolded state and folding cooperativity under native conditions. The methods are based on steric trapping which couples unfolding of a doubly-biotinylated protein to binding of monovalent streptavidin (mSA). We further advanced this technology for general application by developing versatile biotin probes possessing spectroscopic reporters that are sensitized by mSA binding or protein unfolding. By applying these methods to an intramembrane protease GlpG of Escherichia coli, we elucidated a widely unraveled unfolded state, subglobal unfolding of the region encompassing the active site, and a network of cooperative and localized interactions to maintain the stability. These findings provide crucial insights into the folding energy landscape of membrane proteins.

The Effect of Detergent, Temperature, and Lipid on the Oligomeric State of MscL Constructs: Insights from Mass Spectrometry

The mechanosensitive channel of large conductance (MscL) acts as an emergency release valve for osmotic shock of bacteria preventing cell lysis. The large pore size, essential for function, requires the formation of oligomers with tetramers, pentamers, or hexamers observed depending on the species and experimental approach. We applied non-denaturing (native) mass spectrometry to five different homologs of MscL to determine the oligomeric state under more than 50 different experimental conditions elucidating lipid binding and subunit stoichiometry. We found equilibrium between pentameric and tetrameric species, which can be altered by detergent, disrupted by binding specific lipids, and perturbed by increasing temperature (37°C). We also established the presence of lipopolysaccharide bound to MscL and other membrane proteins expressed in Escherichia coli, revealing a potential source of heterogeneity. More generally, we highlight the use of mass spectrometry in probing membrane proteins under a variety of detergent-lipid environments relevant to structural biology.

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.

Personalized biochemistry and biophysics

Whole human genome sequencing of individuals is becoming rapid and inexpensive, enabling new strategies for using personal genome information to help diagnose, treat, and even prevent human disorders for which genetic variations are causative or are known to be risk factors. Many of the exploding number of newly discovered genetic variations alter the structure, function, dynamics, stability, and/or interactions of specific proteins and RNA molecules. Accordingly, there are a host of opportunities for biochemists and biophysicists to participate in (1) developing tools to allow accurate and sometimes medically actionable assessment of the potential pathogenicity of individual variations and (2) establishing the mechanistic linkage between pathogenic variations and their physiological consequences, providing a rational basis for treatment or preventive care. In this review, we provide an overview of these opportunities and their associated challenges in light of the current status of genomic science and personalized medicine, the latter often termed precision medicine.

Design and characterization of a membrane protein unfolding platform in lipid bilayers

Accurate measurement of membrane protein stability—and particularly how it may vary as a result of disease-phenotypic mutations—ideally requires a denaturant that can unfold a membrane-embedded structure while leaving the solubilizing environment unaffected. The steric trap method fulfills this requirement by using monovalent streptavidin (mSA) molecules to unfold membrane proteins engineered with two spatially close biotin tags. Here we adapted this method to an 87-residue helix-loop-helix (hairpin) construct derived from helices 3 and 4 in the transmembrane domain of the human cystic fibrosis transmembrane conductance regulator (CFTR), wherein helix-helix tertiary interactions are anticipated to confer a portion of construct stability. The wild type CFTR TM3/4 hairpin construct was modified with two accessible biotin tags for mSA-induced unfolding, along with two helix-terminal pyrene labels to monitor loss of inter-helical contacts by pyrene excimer fluorescence. A series of eight constructs with biotin tags at varying distances from the helix-terminal pyrene labels were expressed, purified and labeled appropriately; all constructs exhibited largely helical circular dichroism spectra. We found that addition of mSA to an optimized construct in lipid vesicles led to a complete and reversible loss in pyrene excimer fluorescence and mSA binding, and hence hairpin unfolding—results further supported by SDS-PAGE visualization of mSA bound and unbound species. While some dimeric/oligomeric populations persist that may affect quantitation of the unfolding step, our characterization of the design yields a promising prototype of a future platform for the systematic study of membrane protein folding in a lipid bilayer environment.

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.

Differential contribution of tryptophans to the folding and stability of the attachment invasion locus transmembrane β-barrel from Yersinia pestis

Attachment invasion locus (Ail) protein of Yersinia pestis is a crucial outer membrane protein for host invasion and determines bacterial survival within the host. Despite its importance in pathogenicity, surprisingly little is known on Ail biophysical properties. We investigate the contribution of micelle concentrations and interface tryptophans on the Ail β-barrel refolding and unfolding processes. Our results reveal that barrel folding is surprisingly independent of micelle amounts, but proceeds through an on-pathway intermediate that requires the interface W42 for cooperative barrel refolding. On the contrary, the unfolding event is strongly controlled by absolute micelle concentrations. We find that upon Trp → Phe substitution, protein stabilities follow the order W149F>WT>W42F for the refolding, and W42F>WT>W149F for unfolding. W42 confers cooperativity in barrel folding, and W149 clamps the post-folded barrel structure to its micelle environment. Our analyses reveal, for the first time, that interface tryptophan mutation can indeed render greater β-barrel stability. Furthermore, hysteresis in Ail stems from differential barrel-detergent interaction strengths in a micelle concentration-dependent manner, largely mediated by W149. The kinetically stabilized Ail β-barrel has strategically positioned tryptophans to balance efficient refolding and subsequent β-barrel stability, and may be evolutionarily chosen for optimal functioning of Ail during Yersinia pathogenesis.

Cysteine residues impact the stability and micelle interaction dynamics of the human mitochondrial β-barrel anion channel hVDAC-2

The anti-apoptotic 19-stranded transmembrane human voltage dependent anion channel isoform 2 (hVDAC-2) β-barrel stability is crucial for anion transport in mitochondria. The role of the unusually high number of cysteine residues in this isoform is poorly understood. Using a Cys-less construct of hVDAC-2, we haveinvestigated the contribution of cysteines to channel function, barrel stability and its influence on the strength of protein-micelle interactions. We observe that despite the overall preservation in barrel structure upon cysteine mutation, subtle local variations in the mode of interaction of the barrel with its refolded micellar environment arise, which may manifest itself in the channel activity of both the proteins.Fluorescence measurements of the Trp residues in hVDAC-2 point to possible differences in the association of the barrel with lauryldimethylamine oxide (LDAO) micelles. Upon replacement of cysteines in hVDAC-2, our data suggests greater barrel rigidity by way of intra-protein interactions. This, in turn, lowers the equilibrium barrel thermodynamic parameters in LDAOby perturbingthe stability of the protein-micelle complex. In addition to this, we also find a difference in the cooperativity of unfolding upon increasing the LDAO concentration, implying the importance of micelle concentration and micelle-protein ratios on the stability of this barrel. Our results indicate that the nine cysteine residues of hVDAC-2 are the key in establishing strong(er) barrel interactions with its environment and also impart additional malleability to the barrel scaffold.

Influence of protein-micelle ratios and cysteine residues on the kinetic stability and unfolding rates of human mitochondrial VDAC-2

Delineating the kinetic and thermodynamic factors which contribute to the stability of transmembrane β-barrels is critical to gain an in-depth understanding of membrane protein behavior. Human mitochondrial voltage-dependent anion channel isoform 2 (hVDAC-2), one of the key anti-apoptotic eukaryotic β-barrel proteins, is of paramount importance, owing to its indispensable role in cell survival. We demonstrate here that the stability of hVDAC-2 bears a strong kinetic contribution that is dependent on the absolute micellar concentration used for barrel folding. The refolding efficiency and ensuing stability is sensitive to the lipid-to-protein (LPR) ratio, and displays a non-linear relationship, with both low and high micellar amounts being detrimental to hVDAC-2 structure. Unfolding and aggregation process are sequential events and show strong temperature dependence. We demonstrate that an optimal lipid-to-protein ratio of 2600∶1 – 13000∶1 offers the highest protection against thermal denaturation. Activation energies derived only for lower LPRs are ∼17 kcal mol⁻¹ for full-length hVDAC-2 and ∼23 kcal mol⁻¹ for the Cys-less mutant, suggesting that the nine cysteine residues of hVDAC-2 impart additional malleability to the barrel scaffold. Our studies reveal that cysteine residues play a key role in the kinetic stability of the protein, determine barrel rigidity and thereby give rise to strong micellar association of hVDAC-2. Non-linearity of the Arrhenius plot at high LPRs coupled with observation of protein aggregation upon thermal denaturation indicates that contributions from both kinetic and thermodynamic components stabilize the 19-stranded β-barrel. Lipid-protein interaction and the linked kinetic contribution to free energy of the folded protein are together expected to play a key role in hVDAC-2 recycling and the functional switch at the onset of apoptosis.