Influence of Protein Scaffold on Side-Chain Transfer Free Energies

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.

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

Interpreting the molecular mechanisms of disease variants in human transmembrane proteins

Next-generation sequencing of human genomes reveals millions of missense variants, some of which may lead to loss of protein function and ultimately disease. Here, we investigate missense variants in membrane proteins-key drivers in cell signaling and recognition. We find enrichment of pathogenic variants in the transmembrane region across 19,000 functionally classified variants in human membrane proteins. To accurately predict variant consequences, one fundamentally needs to understand the underlying molecular processes. A key mechanism underlying pathogenicity in missense variants of soluble proteins has been shown to be loss of stability. Membrane proteins, however, are widely understudied. Here, we interpret variant effects on a larger scale by performing structure-based estimations of changes in thermodynamic stability using a membrane-specific energy function and analyses of sequence conservation during evolution of 15 transmembrane proteins. We find evidence for loss of stability being the cause of pathogenicity in more than half of the pathogenic variants, indicating that this is a driving factor also in membrane-protein-associated diseases. Our findings show how computational tools aid in gaining mechanistic insights into variant consequences for membrane proteins. To enable broader analyses of disease-related and population variants, we include variant mappings for the entire human proteome.

Diverse Scientific Benchmarks for Implicit Membrane Energy Functions

Energy functions are fundamental to biomolecular modeling. Their success depends on robust physical formalisms, efficient optimization, and high-resolution data for training and validation. Over the past 20 years, progress in each area has advanced soluble protein energy functions. Yet, energy functions for membrane proteins lag behind due to sparse and low-quality data, leading to overfit tools. To overcome this challenge, we assembled a suite of 12 tests on independent data sets varying in size, diversity, and resolution. The tests probe an energy function's ability to capture membrane protein orientation, stability, sequence, and structure. Here, we present the tests and use the franklin2019 energy function to demonstrate them. We then identify areas for energy function improvement and discuss potential future integration with machine-learning-based optimization methods. The tests are available through the Rosetta Benchmark Server (https://benchmark.graylab.jhu.edu/) and GitHub (https://github.com/rfalford12/Implicit-Membrane-Energy-Function-Benchmark).

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.

Characterization of the Features of Water Inside the SecY Translocon

Despite extended experimental and computational studies, the mechanism regulating membrane protein folding and stability in cell membranes is not fully understood. In this review, I will provide a personal and partial account of the scientific efforts undertaken by Dr. Stephen White to shed light on this topic. After briefly describing the role of water and the hydrophobic effect on cellular processes, I will discuss the physical chemistry of water confined inside the SecY translocon pore. I conclude with a review of recent literature that attempts to answer fundamental questions on the pathway and energetics of translocon-guided membrane protein insertion.

Membrane Protein Engineering with Rosetta

Protein engineering can yield new molecular tools for nanotechnology and therapeutic applications through modulating physiochemical and biological properties. Engineering membrane proteins is especially attractive because they perform key cellular processes including transport, nutrient uptake, removal of toxins, respiration, motility, and signaling. In this chapter, we describe two protocols for membrane protein engineering with the Rosetta software: (1) ΔΔG calculations for single point mutations and (2) sequence optimization in different membrane lipid compositions. These modular protocols are easily adaptable for more complex problems and serve as a foundation for efficient membrane protein engineering calculations.

SurA is a cryptically grooved chaperone that expands unfolded outer membrane proteins

The periplasmic chaperone network ensures the biogenesis of bacterial outer membrane proteins (OMPs) and has recently been identified as a promising target for antibiotics. SurA is the most important member of this network, both due to its genetic interaction with the β-barrel assembly machinery complex as well as its ability to prevent unfolded OMP (uOMP) aggregation. Using only binding energy, the mechanism by which SurA carries out these two functions is not well-understood. Here, we use a combination of photo-crosslinking, mass spectrometry, solution scattering, and molecular modeling techniques to elucidate the key structural features that define how SurA solubilizes uOMPs. Our experimental data support a model in which SurA binds uOMPs in a groove formed between the core and P1 domains. This binding event results in a drastic expansion of the rest of the uOMP, which has many biological implications. Using these experimental data as restraints, we adopted an integrative modeling approach to create a sparse ensemble of models of a SurA•uOMP complex. We validated key structural features of the SurA•uOMP ensemble using independent scattering and chemical crosslinking data. Our data suggest that SurA utilizes three distinct binding modes to interact with uOMPs and that more than one SurA can bind a uOMP at a time. This work demonstrates that SurA operates in a distinct fashion compared to other chaperones in the OMP biogenesis network.

Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane

Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.

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.

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.

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.

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.

A Membrane Burial Potential with H-Bonds and Applications to Curved Membranes and Fast Simulations

We use the statistics of a large and curated training set of transmembrane helical proteins to develop a knowledge-based potential that accounts for the dependence on both the depth of burial of the protein in the membrane and the degree of side-chain exposure. Additionally, the statistical potential includes depth-dependent energies for unsatisfied backbone hydrogen bond donors and acceptors, which are found to be relatively small, ∼2 RT. Our potential accurately places known proteins within the bilayer. The potential is applied to the mechanosensing MscL channel in membranes of varying thickness and curvature, as well as to the prediction of protein structure. The potential is incorporated into our new Upside molecular dynamics algorithm. Notably, we account for the exchange of protein-lipid interactions for protein-protein interactions as helices contact each other, thereby avoiding overestimating the energetics of helix association within the membrane. Simulations of most multimeric complexes find that isolated monomers and the oligomers retain the same orientation in the membrane, suggesting that the assembly of prepositioned monomers presents a viable mechanism of oligomerization.

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.

Building Blocks of the Outer Membrane: Calculating a General Elastic Energy Model for β-Barrel Membrane Proteins

The outer membranes of Gram negative bacteria are the first points of contact these organisms make with their environment. Understanding how composition determines the mechanical properties of this essential barrier is of paramount importance. Therefore, we developed a new computational method to measure the elasticity of transmembrane proteins found in the outer membrane. Using all-atom molecular dynamics simulations of these proteins, we apply a set of external forces to mechanically stress the transmembrane β-barrels. Our results from four representative β-barrels show that outer membrane proteins display elastic properties that are approximately 70 to 190 times stiffer than neat lipid membranes. These findings suggest that outer membrane β-barrels are a significant source of mechanical stability in bacteria. Our all-atom approach further reveals that resistance to radial stress is encoded by a general mechanism that includes stretching of backbone hydrogen bonds and tilting of β-strands with respect to the bilayer normal. This computational framework facilitates an increased theoretical understanding of how varying lipid and protein amounts affect the mechanical properties of the bacterial outer membrane.

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.

Position-Specific contribution of interface tryptophans on membrane protein energetics

Interface tryptophans are key residues that facilitate the folding and stability of membrane proteins. Escherichia coli OmpX possesses two unique interface tryptophans, namely Trp76, which is present at the interface and is solvent-exposed, and Trp140, which is relatively more lipid solvated than Trp76 in symmetric lipid membranes. Here, we address the requirement for tryptophan and the consequences of aromatic amino acid substitutions on the folding and stability of OmpX. Using spectroscopic measurements of OmpX-Trp/Tyr/Phe mutants, we show that the specific mutation W76→Y allows barrel assembly >1.5-fold faster than native OmpX, and increases stability by ~0.4kcalmol^-1. In contrast, mutating W140→F/Y lowers OmpX thermodynamic stability by ~0.4kcalmol^-1, without affecting the folding kinetics. We conclude that the stabilizing effect of tryptophan at the membrane interface can be position-and local environment-specific. We propose that the thermodynamic contributions for interface residues be interpreted with caution.