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

Folding of Circularly Permuted and Split Outer Membrane Protein F via Electrostatic Interactions with Terminal Residues

Membrane proteins are essential targets in drug design, biosensing, and catalysis. In this study, we explored the folding of engineered outer membrane protein F (OmpF), an abundant and functional β-barrel protein expressed in Gram-negative bacteria. We carried out circular permutation, splitting and self-complementation, and point mutation. The folding efficiency and kinetic analyses demonstrated that the N- and C-terminal residues of OmpF played critical roles in folding via electrostatic interactions with lipid headgroups. Our results indicate that native porins without charged terminal residues may be tightly downregulated to retain the integrity of the outer membrane, and this modification may facilitate the insertion and folding of modified membrane proteins under in vitro and in vivo conditions for various applications.

Comparing Interfacial Trp, Interfacial His and pH Dependence for the Anchoring of Tilted Transmembrane Helical Peptides

Charged and aromatic amino acid residues, being enriched toward the terminals of membrane-spanning helices in membrane proteins, help to stabilize particular transmembrane orientations. Among them, histidine is aromatic and can be positively charge at low pH. To enable investigations of the underlying protein-lipid interactions, we have examined the effects of single or pairs of interfacial histidine residues using the constructive low-dynamic GWALP23 (acetyl-GG²ALW⁵LALALALALALALW¹⁹LAG²²A-amide) peptide framework by incorporating individual or paired histidines at locations 2, 5, 19 or 22. Analysis of helix orientation by means of solid-state ²H NMR spectra of labeled alanine residues reveals marked differences with H^2,22 compared to W^2,22. Nevertheless, the properties of membrane-spanning H^2,22WALP23 helices show little pH dependence and are similar to those having Gly, Arg or Lys at positions 2 and 22. The presence of H5 or H19 influences the helix rotational preference but not the tilt magnitude. H5 affects the helical integrity, as residue 7 unwinds from the core helix; yet once again the helix orientation and dynamic properties show little sensitivity to pH. The overall results reveal that the detailed properties of transmembrane helices depend upon the precise locations of interfacial histidine residues.

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.

Hydrophobic Characteristic Is Energetically Preferred for Cysteine in a Model Membrane Protein

The naturally occurring amino acid cysteine has often been implicated with a crucial role in maintaining protein structure and stability. An intriguing duality in the intrinsic hydrophobicity of the cysteine side chain is that it exhibits both polar as well as hydrophobic characteristics. Here, we have utilized a cysteine-scanning mutational strategy on the transmembrane β-barrel PagP to examine the membrane depth-dependent energetic contribution of the free cysteine side chain (thiolate) versus the parent residue at an experimental pH of 9.5 in phosphatidylcholine vesicles. We find that introduction of cysteine causes destabilization at several of the 26 lipid-facing sites of PagP that we mutated in this study. The destabilization is minimal (0.5-1.5 kcal/mol) when the mutation is toward the bilayer midplane, whereas it is higher in magnitude (3.0-5.0 kcal/mol) near the bilayer interface. These observations suggest that cysteine forms more favorable interactions with the hydrophobic lipid core as compared to the amphiphilic water-lipid interface. The destabilizing effect is more pronounced when cysteine replaces the interfacial aromatics, which are known to participate in tertiary interaction networks in transmembrane β-barrels. Our observations from experiments involving the introduction of cysteine at the bilayer midplane further strengthen previous views that the free cysteine side chain does possess strongly apolar characteristics. Additionally, the free energy changes observed upon cysteine incorporation show a depth-dependent correlation with the estimated energetic cost of partitioning derived from reported hydrophobicity scales. Our results and observations from the thermodynamic analysis of the PagP barrel may explain why cysteine, despite possessing a polar sulfhydryl group, tends to behave as a hydrophobic (rather than polar) residue in folded protein structures.

Salvaging the Thermodynamic Destabilization of Interface Histidine in Transmembrane β-Barrels

The ability of histidine to participate in a wide range of stabilizing polar interactions preferentially populates this residue in functionally important sites of proteins. Histidine possesses an amphiphilic and electrostatic nature that is essential for amino acids residing at membrane interfaces. However, the frequency of occurrence of histidine at membrane interfaces, particularly transmembrane β-barrels, is lower than those of other aromatic residues. Here, we carry out comprehensive energetic measurements using equilibrium folding of the outer membrane enzyme PagP to address the contribution of a C-terminal interface histidine to barrel stability. We show that placing histidine at the C-terminus universally destabilizes PagP by 4.0-8.0 kcal mol^-1 irrespective of the neighboring residue. Spectroscopic and electrophoretic measurements indicate that the altered stability may arise from a loss of barrel compaction. Isoleucine, methionine, and valine salvage this destabilization marginally (in addition to tyrosine, which shows an exceptionally high folding free energy value), when placed at the penultimate position, at the expense of an altered folding pathway. Double-mutant cycle analysis indicates that the coupling energy between the terminal and penultimate residues in PagP-X¹⁶⁰H¹⁶¹ increases when the level of intrinsic destabilization by the terminal H¹⁶¹ is high. Our observations that neighboring residues cannot salvage the energetic destabilization of histidine may explain why histidine is less abundant at membrane interfaces.

Comparative Analysis of TM and Cytoplasmic β-barrel Conformations Using Joint Descriptor

Macroscopic descriptors have become valuable as coarse-grained features of complex proteins and are complementary to microscopic descriptors. Proteins macroscopic geometric features provide effective clues in the quantification of distant similarity and close dissimilarity searches for structural comparisons. In this study, we performed a systematic comparison of β-barrels, one of the important classes of protein folds in various transmembrane (TM) proteins against cytoplasmic barrels to estimate the conformational features using a joint-based descriptor. The approach uses joint coordinates and dihedral angles (β and γ) based on the β-strand joints and loops to determine the arrangements and propensities at the local and global levels. We then confirmed that there is a clear preference in the overall β and γ distribution, arrangements of β-strands and loops, signature patterns, and the number of strand effects between TM and cytoplasmic β-barrel geometries. As a robust and simple approach, we determine that the joint-based descriptor could provide a reliable static structural comparison aimed at macroscopic level between complex protein conformations.

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.

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.

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.

High-resolution structure prediction of β-barrel membrane proteins

[Formula: see text]-Barrel membrane proteins ([Formula: see text]MPs) play important roles, but knowledge of their structures is limited. We have developed a method to predict their 3D structures. We predict strand registers and construct transmembrane (TM) domains of [Formula: see text]MPs accurately, including proteins for which no prediction has been attempted before. Our method also accurately predicts structures from protein families with a limited number of sequences and proteins with novel folds. An average main-chain rmsd of 3.48 Å is achieved between predicted and experimentally resolved structures of TM domains, which is a significant improvement ([Formula: see text]3 Å) over a recent study. For [Formula: see text]MPs with NMR structures, the deviation between predictions and experimentally solved structures is similar to the difference among the NMR structures, indicating excellent prediction accuracy. Moreover, we can now accurately model the extended [Formula: see text]-barrels and loops in non-TM domains, increasing the overall coverage of structure prediction by [Formula: see text]%. Our method is general and can be applied to genome-wide structural prediction of [Formula: see text]MPs.

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.

Energetics of side-chain partitioning of β-signal residues in unassisted folding of a transmembrane β-barrel protein

The free energy of water-to-interface amino acid partitioning is a major contributing factor in membrane protein folding and stability. The interface residues at the C terminus of transmembrane β-barrels form the β-signal motif required for assisted β-barrel assembly in vivo but are believed to be less important for β-barrel assembly in vitro. Here, we experimentally measured the thermodynamic contribution of all 20 amino acids at the β-signal motif to the unassisted folding of the model β-barrel protein PagP. We obtained the partitioning free energy for all 20 amino acids at the lipid-facing interface (ΔΔG⁰_w,i(φ)) and the protein-facing interface (ΔΔG⁰_w,i(π)) residues and found that hydrophobic amino acids are most favorably transferred to the lipid-facing interface, whereas charged and polar groups display the highest partitioning energy. Furthermore, the change in non-polar surface area correlated directly with the partitioning free energy for the lipid-facing residue and inversely with the protein-facing residue. We also demonstrate that the interface residues of the β-signal motif are vital for in vitro barrel assembly, because they exhibit a side chain–specific energetic contribution determined by the change in nonpolar accessible surface. We further establish that folding cooperativity and hydrophobic collapse are balanced at the membrane interface for optimal stability of the PagP β-barrel scaffold. We conclude that the PagP C-terminal β-signal motif influences the folding cooperativity and stability of the folded β-barrel and that the thermodynamic contributions of the lipid- and protein-facing residues in the transmembrane protein β-signal motif depend on the nature of the amino acid side chain.

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

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

Outer membrane protein design

Membrane proteins are the gateway to the cell. These proteins are also a control center of the cell, as information from the outside is passed through membrane proteins as signals to the cellular machinery. The design of membrane proteins seeks to harness the power of these gateways and signal carriers. This review will focus on the design of the membrane proteins that are in the outer membrane, a membrane which only exists for gram negative bacteria, mitochondria, and chloroplasts. Unlike other membrane proteins, outer membrane proteins are uniquely shaped as β-barrels. Herein, I describe most known examples of membrane β-barrel design to date, focusing particularly on categorizing designs as: Firstly, structural deconstruction; secondly, structural changes; thirdly, chemical function design; and finally, the creation of new folds.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

Aromatic Side Chain Water-to-Lipid Transfer Free Energies Show a Depth Dependence across the Membrane Normal

Quantitating and understanding the physical forces responsible for the interactions of biomolecules are fundamental to the biological sciences. This is especially challenging for membrane proteins because they are embedded within cellular bilayers that provide a unique medium in which hydrophobic sequences must fold. Knowledge of the energetics of protein-lipid interactions is thus vital to understand cellular processes involving membrane proteins. Here we used a host-guest mutational strategy to calculate the Gibbs free energy changes of water-to-lipid transfer for the aromatic side chains Trp, Tyr, and Phe as a function of depth in the membrane. This work reveals an energetic gradient in the transfer free energies for Trp and Tyr, where transfer was most favorable to the membrane interfacial region and comparatively less favorable into the bilayer center. The transfer energetics follows the concentration gradient of polar atoms across the bilayer normal that naturally occurs in biological membranes. Additional measurements revealed nearest-neighbor coupling in the data set are influenced by a network of aromatic side chains in the host protein. Taken together, these results show that aromatic side chains contribute significantly to membrane protein stability through either aromatic-aromatic interactions or placement at the membrane interface.