Overcoming hysteresis to attain reversible equilibrium folding for outer membrane phospholipase A in phospholipid bilayers

Folding pathway of a discontinuous two-domain protein

It is estimated that two-thirds of all proteins in higher organisms are composed of multiple domains, many of them containing discontinuous folds. However, to date, most in vitro protein folding studies have focused on small, single-domain proteins. As a model system for a two-domain discontinuous protein, we study the unfolding/refolding of a slow-folding double mutant of the maltose binding protein (DM-MBP) using single-molecule two- and three-color Förster Resonance Energy Transfer experiments. We observe a dynamic folding intermediate population in the N-terminal domain (NTD), C-terminal domain (CTD), and at the domain interface. The dynamic intermediate fluctuates rapidly between unfolded states and compact states, which have a similar FRET efficiency to the folded conformation. Our data reveals that the delayed folding of the NTD in DM-MBP is imposed by an entropic barrier with subsequent folding of the highly dynamic CTD. Notably, accelerated DM-MBP folding is routed through the same dynamic intermediate within the cavity of the GroEL/ES chaperone system, suggesting that the chaperonin limits the conformational space to overcome the entropic folding barrier. Our study highlights the subtle tuning and co-dependency in the folding of a discontinuous multi-domain protein.

Hysteretic Pressure Dependence of Ca2+ Binding in LH1 Bacterial Membrane Chromoproteins

Much of the thermodynamic parameter values that support life are set by the properties of proteins. While the denaturing effects of pressure and temperature on proteins are well documented, their precise structural nature is rarely revealed. This work investigates the destabilization of multiple Ca²⁺ binding sites in the cyclic LH1 light-harvesting membrane chromoprotein complexes from two Ca-containing sulfur purple bacteria by hydrostatic high-pressure perturbation spectroscopy. The native (Ca-saturated) and denatured (Ca-depleted) phases of these complexes are well distinguishable by much-shifted bacteriochlorophyll a exciton absorption bands serving as innate optical probes in this study. The pressure-induced denaturation of the complexes related to the failure of the protein Ca-binding pockets and the concomitant breakage of hydrogen bonds between the pigment chromophores and protein environment were found cooperative, involving all or most of the Ca²⁺ binding sites, but irreversible. The strong hysteresis observed in the spectral and kinetic characteristics of phase transitions along the compression and decompression pathways implies asymmetry in the relevant free energy landscapes and activation free energy distributions. A phase transition pressure equal to about 1.9 kbar was evaluated for the complexes from Thiorhodovibrio strain 970 from the pressure dependence of biphasic kinetics observed in the minutes to 100 h time range.

The Thermodynamic Stability of Membrane Proteins in Micelles and Lipid Bilayers Investigated with the Ferrichrom Receptor FhuA

Extraction of integral membrane proteins into detergents for structural and functional studies often leads to a strong loss in protein stability. The impact of the lipid bilayer on the thermodynamic stability of an integral membrane protein in comparison to its solubilized form in detergent was examined and compared for FhuA from Escherichia coli and for a mutant, FhuAΔ5-160, lacking the N-terminal cork domain. Urea-induced unfolding was monitored by fluorescence spectroscopy to determine the effective free energies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta G{^\text{o}_{\rm u}} $$\end{document}ΔGuo of unfolding. To obtain enthalpic and entropic contributions of unfolding of FhuA, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta G{^\text{o}_{\rm u}} $$\end{document}ΔGuo were determined at various temperatures. When solubilized in LDAO detergent, wt-FhuA and FhuAΔ5-160 unfolded in a single step. The 155-residue cork domain stabilized wt-FhuA by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta\Delta G{^\text{o}_{\rm u}} $$\end{document}ΔΔGuo~ 40 kJ/mol. Reconstituted into lipid bilayers, wt-FhuA unfolded in two steps, while FhuAΔ5-160 unfolded in a single step, indicating an uncoupled unfolding of the cork domain. For FhuAΔ5-160 at 35 °C, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta G{^\text{o}_{\rm u}} $$\end{document}ΔGuo increased from ~ 5 kJ/mol in LDAO micelles to about ~ 20 kJ/mol in lipid bilayers, while the temperature of unfolding increased from T_M ~ 49 °C in LDAO micelles to T_M ~ 75 °C in lipid bilayers. Enthalpies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta H{_{\rm M}^\text{o}}$$\end{document}ΔHMowere much larger than free energies \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta G{^\text{o}_{\rm u}} $$\end{document}ΔGuo, for FhuAΔ5-160 and for wt-FhuA, and compensated by a large gain of entropy upon unfolding. The gain in conformational entropy is expected to be similar for unfolding of FhuA from micelles or bilayers. The strongly increased T_M and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta H{_{\rm M}^\text{o}}$$\end{document}ΔHMo observed for the lipid bilayer-reconstituted FhuA in comparison to the LDAO-solubilized forms, therefore, very likely arise from a much-increased solvation entropy of FhuA in bilayers.

Amphiphilic Membrane Environments Regulate Enzymatic Behaviors of Salmonella Outer Membrane Protease

The role of an amphiphilic environment in the functional regulation of integral membrane proteins is well appreciated but how specific amphiphilic surrounding influences the conformational plasticity and function of a protein is less obvious. We focus on the Salmonella phosphoglycerate transport system (pgt)-encoded outer membrane protease E (PgtE), which plays an important role in tissue infiltration and survival of Salmonella enterica. Despite our understanding of its physiological functions, elucidation of its enzymatic behavior in response to the immediate amphiphilic surrounding is lacking. We monitor the proteolytic activity of PgtE reconstituted in Zwittergent 3-12 detergent micelles or a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayer and examine factors that influence its activity. We find, to our surprise, that PgtE, which is thought to elicit a rapid response toward various substrates, showed hysteretic enzymatic behavior, characterized by a prominent lag phase prior to achieving the exponential steady state in its detergent-stabilized form as well as in the outer membrane embedded native state in live bacteria. The lag phase was abolished under three conditions: preformation of an inactive detergent-stabilized PgtE-substrate complex without lipopolysaccharide (LPS), LPS-bound detergent-stabilized PgtE that had reached steady state velocity, or PgtE reconstituted into a POPC bilayer environment. Interestingly, detergent- and bilayer-stabilized PgtE showed comparable steady-state activity. And strikingly, lipopolysaccharide (LPS) becomes nonessential for the activation of PgtE when the protein is reconstituted in the phospholipid bilayer, contrasting a long-standing notion that LPS is required for proteases belonging to the omptin family to be proteolytically active. These findings suggest intriguing biological nuances for the proteolytic function of PgtE that were not well appreciated previously and offer new perspectives that may generally be applicable for omptins.

Principles and Methods in Computational Membrane Protein Design

After decades of progress in computational protein design, the design of proteins folding and functioning in lipid membranes appears today as the next frontier. Some notable successes in the de novo design of simplified model membrane protein systems have helped articulate fundamental principles of protein folding, architecture and interaction in the hydrophobic lipid environment. These principles are reviewed here, together with the computational methods and approaches that were used to identify them. We provide an overview of the methodological innovations in the generation of new protein structures and functions and in the development of membrane-specific energy functions. We highlight the opportunities offered by new machine learning approaches applied to protein design, and by new experimental characterization techniques applied to membrane proteins. Although membrane protein design is in its infancy, it appears more reachable than previously thought.

Membrane proteins enter the fold

Membrane proteins have historically been recalcitrant to biophysical folding studies. However, recent adaptations of methods from the soluble protein folding field have found success in their applications to transmembrane proteins composed of both α-helical and β-barrel conformations. Avoiding aggregation is critical for the success of these experiments. Altogether these studies are leading to discoveries of folding trajectories, foundational stabilizing forces and better-defined endpoints that enable more accurate interpretation of thermodynamic data. Increased information on membrane protein folding in the cell shows that the emerging biophysical principles are largely recapitulated even in the complex biological environment.

De novo design of transmembrane β barrels

Transmembrane β-barrel proteins (TMBs) are of great interest for single-molecule analytical technologies because they can spontaneously fold and insert into membranes and form stable pores, but the range of pore properties that can be achieved by repurposing natural TMBs is limited. We leverage the power of de novo computational design coupled with a "hypothesis, design, and test" approach to determine TMB design principles, notably, the importance of negative design to slow β-sheet assembly. We design new eight-stranded TMBs, with no homology to known TMBs, that insert and fold reversibly into synthetic lipid membranes and have nuclear magnetic resonance and x-ray crystal structures very similar to the computational models. These advances should enable the custom design of pores for a wide range of applications.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Determining the Free Energies of Outer Membrane Proteins in Lipid Bilayers

The thermodynamic stabilities of membrane proteins are of fundamental interest to provide a biophysical description of their structure-function relationships because energy determines conformational populations. In addition, structure-energy relationships can be exploited in membrane protein design and in synthetic biology. To determine the thermodynamic stability of a membrane protein, it is not sufficient to be able to unfold and refold the molecule: establishing path independence of this reaction is essential. Here we describe the procedures required to measure and verify path independence for the folding of outer membrane proteins in large unilamellar vesicles.

GeTFEP: A general transfer free energy profile of transmembrane proteins

Free energy of transferring amino acid side-chains from aqueous environment into lipid bilayers, known as transfer free energy (TFE), provides important information on the thermodynamic stability of membrane proteins. In this study, we derived a TFE profile named General Transfer Free Energy Profile (GeTFEP) based on computation of the TFEs of 58 β-barrel membrane proteins (βMPs). The GeTFEP agrees well with experimentally measured and computationally derived TFEs. Analysis based on the GeTFEP shows that residues in different regions of the transmembrane (TM) segments of βMPs have different roles during the membrane insertion process. Results further reveal the importance of the sequence pattern of TM strands in stabilizing βMPs in the membrane environment. In addition, we show that GeTFEP can be used to predict the positioning and the orientation of βMPs in the membrane. We also show that GeTFEP can be used to identify structurally or functionally important amino acid residue sites of βMPs. Furthermore, the TM segments of α-helical membrane proteins can be accurately predicted with GeTFEP, suggesting that the GeTFEP is of general applicability in studying membrane protein.

Single-Molecule Mechanistic Study of Enzyme Hysteresis

Hysteresis is an important feature of enzyme-catalyzed reactions, as it reflects the influence of enzyme regulation in the presence of ligands such as substrates or allosteric molecules. In typical kinetic studies of enzyme activity, hysteretic behavior is observed as a "lag" or "burst" in the time course of the catalyzed reaction. These lags and bursts are due to the relatively slow transition from one state to another state of the enzyme molecule, with different states having different kinetic properties. However, it is difficult to understand the underlying mechanism of hysteresis by observing bulk reactions because the different enzyme molecules in the population behave stochastically. In this work, we studied the hysteretic behavior of mutant β-glucuronidase (GUS) using a high-throughput single-molecule array platform and investigated the effect of thermal treatment on the hysteresis.

Statistically derived asymmetric membrane potentials from α-helical and β-barrel membrane proteins

Modeling membrane protein (MP) folding, insertion, association and their interactions with other proteins, lipids, and drugs requires accurate transfer free energies (TFEs). Various TFE scales have been derived to quantify the energy required or released to insert an amino acid or protein into the membrane. Experimental measurement of TFEs is challenging, and only few scales were extended to depth-dependent energetic profiles. Statistical approaches can be used to derive such potentials; however, this requires a sufficient number of MP structures. Furthermore, MPs are tightly coupled to bilayers that are heterogeneous in terms of lipid composition, asymmetry, and protein content between organisms and organelles. Here we derived asymmetric implicit membrane potentials from β-barrel and α-helical MPs and use them to predict topology, depth and orientation of proteins in the membrane. Our data confirm the ‘charge-outside’ and ‘positive-inside’ rules for β-barrels and α-helical proteins, respectively. We find that the β-barrel profiles have greater asymmetry than the ones from α-helical proteins, as a result of the different membrane architecture of gram-negative bacterial outer membranes and the existence of lipopolysaccharide in the outer leaflet. Our data further suggest that pore-facing residues in β-barrels have a larger contribution to membrane insertion and stability than previously suggested.

Balancing Force Field Protein-Lipid Interactions To Capture Transmembrane Helix-Helix Association

Atomistic simulations have recently been shown to be sufficiently accurate to reversibly fold globular proteins and have provided insights into folding mechanisms. Gaining similar understanding from simulations of membrane protein folding and association would be of great medical interest. All-atom simulations of the folding and assembly of transmembrane protein domains are much more challenging, not least due to very slow diffusion within the lipid bilayer membrane. Here, we focus on a simple and well-characterized prototype of membrane protein folding and assembly, namely the dimerization of glycophorin A, a homodimer of single transmembrane helices. We have determined the free energy landscape for association of the dimer using the CHARMM36 force field. We find that the native structure is a metastable state, but not stable as expected from experimental estimates of the dissociation constant and numerous experimental structures obtained under a variety of conditions. We explore two straightforward approaches to address this problem and demonstrate that they result in stable dimers with dissociation constants consistent with experimental data.

Efficient computation of transfer free energies of amino acids in beta-barrel membrane proteins

Motivation

Transmembrane beta-barrel proteins (TMBs) serve a multitude of essential cellular functions in Gram-negative bacteria, mitochondria and chloroplasts. Transfer free energies (TFEs) of residues in the transmembrane (TM) region provides fundamental quantifications of thermodynamic stabilities of TMBs, which are important for the folding and the membrane insertion processes, and may help in understanding the structure-function relationship. However, experimental measurement of TFEs of TMBs is challenging. Although a recent computational method can be used to calculate TFEs, the results of which are in excellent agreement with experimentally measured values, this method does not scale up, and is limited to small TMBs.

Results

We have developed an approximation method that calculates TFEs of TM residues in TMBs accurately, with which depth-dependent transfer free energy profiles can be derived. Our results are in excellent agreement with experimental measurements. This method is efficient and applicable to all bacterial TMBs regardless of the size of the protein.

Availability and Implementation

An online webserver is available at http://tanto.bioe.uic.edu/tmb-tfe .

Contact

: jliang@uic.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Outer membrane protein folding from an energy landscape perspective

The cell envelope is essential for the survival of Gram-negative bacteria. This specialised membrane is densely packed with outer membrane proteins (OMPs), which perform a variety of functions. How OMPs fold into this crowded environment remains an open question. Here, we review current knowledge about OMP folding mechanisms in vitro and discuss how the need to fold to a stable native state has shaped their folding energy landscapes. We also highlight the role of chaperones and the β-barrel assembly machinery (BAM) in assisting OMP folding in vivo and discuss proposed mechanisms by which this fascinating machinery may catalyse OMP folding.

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.

Influence of Protein Scaffold on Side-Chain Transfer Free Energies

The process by which membrane proteins fold involves the burial of side chains into lipid bilayers. Both structure and function of membrane proteins depend on the magnitudes of side-chain transfer free energies (ΔΔG_sc^o). In the absence of other interactions, ΔΔG_sc^o is an independent property describing the energetics of an isolated side chain in the bilayer. However, in reality, side chains are attached to the peptide backbone and surrounded by other side chains in the protein scaffold in biology, which may alter the apparent ΔΔG_sc^o. Previously we reported a whole protein water-to-bilayer hydrophobicity scale using the transmembrane β-barrel Escherichia coli OmpLA as a scaffold protein. To investigate how a different protein scaffold can modulate these energies, we measured ΔΔG_sc^o for all 20 amino acids using the transmembrane β-barrel E. coli PagP as a scaffold protein. This study represents, to our knowledge, the first instance of ΔΔG_sc^o measured in the same experimental conditions in two structurally and sequentially distinct protein scaffolds. Although the two hydrophobicity scales are strongly linearly correlated, we find that there are apparent scaffold induced changes in ΔΔG_sc^o for more than half of the side chains, most of which are polar residues. We propose that the protein scaffold affects the ΔΔG_sc^o of side chains that are buried in unfavorable environments by dictating the mechanisms by which the side chain can reach a more favorable environment and thus modulating the magnitude of ΔΔG_sc^o.

Solubilization of Membrane Proteins into Functional Lipid-Bilayer Nanodiscs Using a Diisobutylene/Maleic Acid Copolymer

Once removed from their natural environment, membrane proteins depend on membrane-mimetic systems to retain their native structures and functions. To this end, lipid-bilayer nanodiscs that are bounded by scaffold proteins or amphiphilic polymers such as styrene/maleic acid (SMA) copolymers have been introduced as alternatives to detergent micelles and liposomes for in vitro membrane-protein research. Herein, we show that an alternating diisobutylene/maleic acid (DIBMA) copolymer shows equal performance to SMA in solubilizing phospholipids, stabilizes an integral membrane enzyme in functional bilayer nanodiscs, and extracts proteins of various sizes directly from cellular membranes. Unlike aromatic SMA, aliphatic DIBMA has only a mild effect on lipid acyl-chain order, does not interfere with optical spectroscopy in the far-UV range, and does not precipitate in the presence of low millimolar concentrations of divalent cations.

Concatemers of Outer Membrane Protein A Take Detours in the Folding Landscape

Outer membrane protein A (OmpA) is the most abundant protein in the outer membrane of Escherichia coli. The N-terminal domain forms an eight-stranded membrane-embedded β-barrel that is widely used as a model protein for in vitro folding into the membrane and into surfactant micelles. Under conditions that include a low surfactant concentration, OmpA can form stable higher-order structures by intermolecular association. Other β-barrel membrane proteins also associate to form noncovalently linked trimers in vivo. This inspired us to test how topological constraints imposed by intramolecular links between individual OmpA molecules affect this process. Here we report on the properties of concatemers consisting of two and three copies of the transmembrane part of OmpA. Both concatemers could be folded to a native state in surfactant micelles according to spectroscopy and electrophoretic band shifts. This native state had the same thermodynamic stability against chemical denaturation as the original OmpA. Above 1.5 M GdmCl, concatemerization increased both refolding and unfolding rates, which we attribute to entropic effects. However, below 1.5 M GdmCl, folding kinetics were 2-3 orders of magnitude slower and more complex, involving a greater degree of parallel folding steps and species that could be classified as off-pathway. Only OmpA2 could quantitatively be folded into vesicles (though to an extent lower than that of OmpA), while OmpA3 formed three species with different levels of folding. Thus, close spatial and sequential proximity of OmpA domains on the same polypeptide chain have a strong tendency to trap the protein in different misfolded states.

Negative Charge Neutralization in the Loops and Turns of Outer Membrane Phospholipase A Impacts Folding Hysteresis at Neutral pH

Hysteresis in equilibrium protein folding titrations is an experimental barrier that must be overcome to extract meaningful thermodynamic quantities. Traditional approaches to solving this problem involve testing a spectrum of solution conditions to find ones that achieve path independence. Through this procedure, a specific pH of 3.8 was required to achieve path independence for the water-to-bilayer equilibrium folding of outer membrane protein OmpLA. We hypothesized that the neutralization of negatively charged side chains (Asp and Glu) at pH 3.8 could be the physical basis for path-independent folding at this pH. To test this idea, we engineered variants of OmpLA with Asp → Asn and Glu → Gln mutations to neutralize the negative charges within various regions of the protein and tested for reversible folding at neutral pH. Although not fully resolved, our results show that these mutations in the periplasmic turns and extracellular loops are responsible for 60% of the hysteresis in wild-type folding. Overall, our study suggests that negative charges impact the folding hysteresis in outer membrane proteins and their neutralization may aid in protein engineering applications.

Outer Membrane Protein Folding and Topology from a Computational Transfer Free Energy Scale

Knowledge of the transfer free energy of amino acids from aqueous solution to a lipid bilayer is essential for understanding membrane protein folding and for predicting membrane protein structure. Here we report a computational approach that can calculate the folding free energy of the transmembrane region of outer membrane β-barrel proteins (OMPs) by combining an empirical energy function with a reduced discrete state space model. We quantitatively analyzed the transfer free energies of 20 amino acid residues at the center of the lipid bilayer of OmpLA. Our results are in excellent agreement with the experimentally derived hydrophobicity scales. We further exhaustively calculated the transfer free energies of 20 amino acids at all positions in the TM region of OmpLA. We found that the asymmetry of the Gram-negative bacterial outer membrane as well as the TM residues of an OMP determine its functional fold in vivo. Our results suggest that the folding process of an OMP is driven by the lipid-facing residues in its hydrophobic core, and its NC-IN topology is determined by the differential stabilities of OMPs in the asymmetrical outer membrane. The folding free energy is further reduced by lipid A and assisted by general depth-dependent cooperativities that exist between polar and ionizable residues. Moreover, context-dependency of transfer free energies at specific positions in OmpLA predict regions important for protein function as well as structural anomalies. Our computational approach is fast, efficient and applicable to any OMP.

N-helix and Cysteines Inter-regulate Human Mitochondrial VDAC-2 Function and Biochemistry

Human voltage-dependent anion channel-2 (hVDAC-2) functions primarily as the crucial anti-apoptotic protein in the outer mitochondrial membrane, and additionally as a gated bidirectional metabolite transporter. The N-terminal helix (NTH), involved in voltage sensing, bears an additional 11-residue extension (NTE) only in hVDAC-2. In this study, we assign a unique role for the NTE as influencing the chaperone-independent refolding kinetics and overall thermodynamic stability of hVDAC-2. Our electrophysiology data shows that the N-helix is crucial for channel activity, whereas NTE sensitizes this isoform to voltage gating. Additionally, hVDAC-2 possesses the highest cysteine content, possibly for regulating reactive oxygen species content. We identify interdependent contributions of the N-helix and cysteines to channel function, and the measured stability in micellar environments with differing physicochemical properties. The evolutionary demand for the NTE in the presence of cysteines clearly emerges from our biochemical and functional studies, providing insight into factors that functionally demarcate hVDAC-2 from the other VDACs.

The negative charge of the membrane has opposite effects on the membrane entry and exit of pH-low insertion peptide

The pH-low insertion peptide (pHLIP) targets acidic diseases such as cancer. The acidity of the environment causes key aspartic acids in pHLIP to become protonated, causing the peptide to insert into membranes. Here we investigate how the negative charge of the membrane influences how pHLIP enters and exits the lipid bilayer. We found that electrostatic repulsion affected differently the membrane entry and exit of pHLIP for negatively charged membranes. As a consequence, a large hysteresis was observed. We propose this is not a consequence of structural changes but results from local changes in the environment of aspartic acids, shifting their pK values.

Thermodynamic, structural and functional properties of membrane protein inclusion bodies are analogous to purified counterparts: case study from bacteria and humans

Purification-free transmembrane protein inclusion body preparations for rapid and cost-effective biophysical, functional and structural studies.

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.

Folded monomers and hexamers of the ectodomain of the HIV gp41 membrane fusion protein: potential roles in fusion and synergy between the fusion peptide, hairpin, and membrane-proximal external region

HIV is an enveloped virus and fusion between the HIV and host cell membranes is catalyzed by the ectodomain of the HIV gp41 membrane protein. Both the N-terminal fusion peptide (FP) and C-terminal membrane-proximal external region (MPER) are critical for fusion and are postulated to bind to the host cell and HIV membranes, respectively. Prior to fusion, the gp41 on the virion is a trimer in noncovalent complex with larger gp120 subunits. The gp120 bind host cell receptors and move away or dissociate from gp41 which subsequently catalyzes fusion. In the present work, large gp41 ectodomain constructs were produced and biophysically and structurally characterized. One significant finding is observation of synergy between the FP, hairpin, and MPER in vesicle fusion. The ectodomain-induced fusion can be very efficient with only ∼15 gp41 per vesicle, which is comparable to the number of gp41 on a virion. Conditions are found with predominant monomer or hexamer but not trimer and these may be oligomeric states during fusion. Monomer gp41 ectodomain is hyperthermostable and has helical hairpin structure. A new HIV fusion model is presented where (1) hemifusion is catalyzed by folding of gp41 ectodomain monomers into hairpins and (2) subsequent fusion steps are catalyzed by assembly into a hexamer with FPs in an antiparallel β sheet. There is also significant interest in the gp41 MPER because it is the epitope of several broadly neutralizing antibodies. Two of these antibodies bind our gp41 ectodomain constructs and support investigation of the gp41 ectodomain as an immunogen in HIV vaccine development.

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.

Model for coupled insertion and folding of membrane-spanning proteins

Current understanding of the forces directing the folding of integral membrane proteins is very limited compared to the detailed picture available for water-soluble proteins. While mechanistic studies of the folding process in vitro have been conducted for only a small number of membrane proteins, the available evidence indicates that their folding process is thermodynamically driven like that of soluble proteins. In vivo, however, the majority of integral membrane proteins are installed in membranes by dedicated machinery, suggesting that the cellular systems may act to facilitate and regulate the spontaneous physical process of folding. Both the in vitro folding process and the in vivo pathway must navigate an energy landscape dominated by the energetically favorable burial of hydrophobic segments in the membrane interior and the opposition to folding due to the need for passage of polar segments across the membrane. This manuscript describes a simple, exactly solvable model which incorporates these essential features of membrane protein folding. The model is used to compare the folding time under conditions which depict both the in vitro and in vivo pathways. It is proposed that the cellular complexes responsible for insertion of membrane proteins act by lowering the energy barrier for passage of polar regions through the membrane, thereby allowing the chain to more rapidly achieve the folded state.

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.

Methionine mutations of outer membrane protein X influence structural stability and beta-barrel unfolding

We report the biochemical and biophysical characterization of outer membrane protein X (OmpX), an eight-stranded transmembrane β-barrel from E. coli, and compare the barrel behavior with a mutant devoid of methionine residues. Transmembrane outer membrane proteins of bacterial origin are known to display high tolerance to sequence rearrangements and mutations. Our studies with the triple mutant of OmpX that is devoid of all internal methionine residues (M18L; M21L; M118L) indicate that Met replacement has no influence on the refolding efficiency and structural characteristics of the protein. Surprisingly, the conserved substitution of Met→Leu leads to barrel destabilization and causes a lowering of the unfolding free energy by a factor of ∼8.5 kJ/mol, despite the mutations occurring at the loop regions. We report that the barrel destabilization is accompanied by a loss in cooperativity of unfolding in the presence of chemical denaturants. Furthermore, we are able to detect an unfolding intermediate in the Met-less barrel, whereas the parent protein exhibits a classic two-state unfolding. Thermal denaturation measurements also suggest a greater susceptibility of the OmpX barrel to heat, in the Met-less construct. Our studies reveal that even subtle variations in the extra-membrane region of rigid barrel structures such as OmpX, may bear severe implications on barrel stability. We propose that methionines contribute to efficient barrel structuring and protein-lipid interactions, and are therefore important elements of OmpX stability.

Dissecting the effects of periplasmic chaperones on the in vitro folding of the outer membrane protein PagP

Although many periplasmic folding factors have been identified, the mechanisms by which they interact with unfolded outer membrane proteins (OMPs) to promote correct folding and membrane insertion remain poorly understood. Here, we have investigated the effect of two chaperones, Skp and SurA, on the folding kinetics of the OMP, PagP. Folding kinetics of PagP into both zwitterionic diC12:0PC (1,2-dilauroyl-sn-glycero-3-phosphocholine) liposomes and negatively charged 80:20 diC12:0PC:diC12:0PG [1,2-dilauroyl-sn-glycero-3-phospho-(1'-rac-glycerol)] liposomes were investigated using a combination of spectroscopic and SDS-PAGE assays. The results indicate that Skp modulates the observed rate of PagP folding in a manner that is dependent on the composition of the membrane and the ionic strength of the buffer used. These data suggest that electrostatic interactions play an important role in Skp-assisted substrate delivery to the membrane. In contrast, SurA showed no effect on the observed folding rates of PagP, consistent with the view that these chaperones act by distinct mechanisms in partially redundant parallel chaperone pathways that facilitate OMP assembly. In addition to delivery of the substrate protein to the membrane, the ability of Skp to prevent OMP aggregation was investigated. The results show that folding and membrane insertion of PagP can be restored, in part, by Skp in conditions that strongly favour PagP aggregation. These results illustrate the utility of in vitro systems for dissecting the complex folding environment encountered by OMPs in the periplasm and demonstrate the key role of Skp in holding aggregation-prone OMPs prior to their direct or indirect delivery to the membrane.

Membrane protein thermodynamic stability may serve as the energy sink for sorting in the periplasm

Thermodynamic stabilities are pivotal for understanding structure-function relationships of proteins, and yet such determinations are rare for membrane proteins. Moreover, the few measurements that are available have been conducted under very different experimental conditions, which compromises a straightforward extraction of physical principles underlying stability differences. Here, we have overcome this obstacle and provided structure-stability comparisons for multiple membrane proteins. This was enabled by measurements of the free energies of folding and the m values for the transmembrane proteins PhoP/PhoQ-activated gene product (PagP) and outer membrane protein W (OmpW) from Escherichia coli. Our data were collected in the same lipid bilayer and buffer system we previously used to determine those parameters for E. coli outer membrane phospholipase A (OmpLA). Biophysically, our results suggest that the stabilities of these proteins are strongly correlated to the water-to-bilayer transfer free energy of the lipid-facing residues in their transmembrane regions. We further discovered that the sensitivities of these membrane proteins to chemical denaturation, as judged by their m values, was consistent with that previously observed for water-soluble proteins having comparable differences in solvent exposure between their folded and unfolded states. From a biological perspective, our findings suggest that the folding free energies for these membrane proteins may be the thermodynamic sink that establishes an energy gradient across the periplasm, thus driving their sorting by chaperones to the outer membranes in living bacteria. Binding free energies of these outer membrane proteins with periplasmic chaperones support this energy sink hypothesis.

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.