A team of chaperones play to win in the bacterial periplasm.

The survival and virulence of Gram-negative bacteria require proper biogenesis and maintenance of the outer membrane (OM), which is densely packed with ß-barrel OM proteins (OMPs). Before reaching the OM, precursor unfolded OMPs (uOMPs) must cross the whole cell envelope. A network of periplasmic chaperones and proteases maintains unfolded but folding-competent conformations of these membrane proteins in the aqueous periplasm while simultaneously preventing off-pathway aggregation. These periplasmic proteins utilize different strategies, including conformational heterogeneity, oligomerization, multivalency, and kinetic partitioning, to perform and regulate their functions. Redundant and unique characteristics of the individual periplasmic players synergize to create a protein quality control team capable responding to changing environmental stresses.

Generation of unfolded outer membrane protein ensembles defined by hydrodynamic properties.

Outer membrane proteins (OMPs) must exist as an unfolded ensemble while interacting with a chaperone network in the periplasm of Gram-negative bacteria. Here, we developed a method to model unfolded OMP (uOMP) conformational ensembles using the experimental properties of two well-studied OMPs. The overall sizes and shapes of the unfolded ensembles in the absence of a denaturant were experimentally defined by measuring the sedimentation coefficient as a function of urea concentration. We used these data to model a full range of unfolded conformations by parameterizing a targeted coarse-grained simulation protocol. The ensemble members were further refined by short molecular dynamics simulations to reflect proper torsion angles. The final conformational ensembles have polymer properties different from unfolded soluble and intrinsically disordered proteins and reveal inherent differences in the unfolded states that necessitate further investigation. Building these uOMP ensembles advances the understanding of OMP biogenesis and provides essential information for interpreting structures of uOMP-chaperone complexes.

FkpA enhances membrane protein folding using an extensive interaction surface.

Outer membrane protein (OMP) biogenesis in gram-negative bacteria is managed by a network of periplasmic chaperones that includes SurA, Skp, and FkpA. These chaperones bind unfolded OMPs (uOMPs) in dynamic conformational ensembles to suppress aggregation, facilitate diffusion across the periplasm, and enhance folding. FkpA primarily responds to heat-shock stress, but its mechanism is comparatively understudied. To determine FkpA chaperone function in the context of OMP folding, we monitored the folding of three OMPs and found that FkpA, unlike other periplasmic chaperones, increases the folded yield but decreases the folding rate of OMPs. The results indicate that FkpA behaves as a chaperone and not as a folding catalyst to influence the OMP folding trajectory. Consistent with the folding assay results, FkpA binds all three uOMPs as determined by sedimentation velocity (SV) and photo-crosslinking experiments. We determine the binding affinity between FkpA and uOmpA171 by globally fitting SV titrations and find it to be intermediate between the known affinities of Skp and SurA for uOMP clients. Notably, complex formation steeply depends on the urea concentration, suggesting an extensive binding interface. Initial characterizations of the complex using photo-crosslinking indicate that the binding interface spans the entire FkpA molecule. In contrast to prior findings, folding and binding experiments performed using subdomain constructs of FkpA demonstrate that the full-length chaperone is required for full activity. Together these results support that FkpA has a distinct and direct effect on OMP folding that it achieves by utilizing an extensive chaperone-client interface to tightly bind clients.

Revisiting macromolecular hydration with HullRadSAS.

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

A proteome-wide map of chaperone-assisted protein refolding in a cytosol-like milieu.

The journey by which proteins navigate their energy landscapes to their native structures is complex, involving (and sometimes requiring) many cellular factors and processes operating in partnership with a given polypeptide chain's intrinsic energy landscape. The cytosolic environment and its complement of chaperones play critical roles in granting many proteins safe passage to their native states; however, it is challenging to interrogate the folding process for large numbers of proteins in a complex background with most biophysical techniques. Hence, most chaperone-assisted protein refolding studies are conducted in defined buffers on single purified clients. Here, we develop a limited proteolysis-mass spectrometry approach paired with an isotope-labeling strategy to globally monitor the structures of refolding Escherichia coli proteins in the cytosolic medium and with the chaperones, GroEL/ES (Hsp60) and DnaK/DnaJ/GrpE (Hsp70/40). GroEL can refold the majority (85%) of the E. coli proteins for which we have data and is particularly important for restoring acidic proteins and proteins with high molecular weight, trends that come to light because our assay measures the structural outcome of the refolding process itself, rather than binding or aggregation. For the most part, DnaK and GroEL refold a similar set of proteins, supporting the view that despite their vastly different structures, these two chaperones unfold misfolded states, as one mechanism in common. Finally, we identify a cohort of proteins that are intransigent to being refolded with either chaperone. We suggest that these proteins may fold most efficiently cotranslationally, and then remain kinetically trapped in their native conformations.

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.

Local Bilayer Hydrophobicity Modulates Membrane Protein Stability.

Through the insertion of nonpolar side chains into the bilayer, the hydrophobic effect has long been accepted as a driving force for membrane protein folding. However, how the changing chemical composition of the bilayer affects the magnitude of the side-chain transfer free energies (ΔGsc°) has historically not been well understood. A particularly challenging region for experimental interrogation is the bilayer interfacial region that is characterized by a steep polarity gradient. In this study, we have determined the ΔGsc° for nonpolar side chains as a function of bilayer position using a combination of experiment and simulation. We discovered an empirical correlation between the surface area of the nonpolar side chain, the transfer free energies, and the local water concentration in the membrane that allows for ΔGsc° to be accurately estimated at any location in the bilayer. Using these water-to-bilayer ΔGsc° values, we calculated the interface-to-bilayer transfer free energy (ΔGi,b°). We find that the ΔGi,b° are similar to the "biological", translocon-based transfer free energies, indicating that the translocon energetically mimics the bilayer interface. Together these findings can be applied to increase the accuracy of computational workflows used to identify and design membrane proteins as well as bring greater insight into our understanding of how disease-causing mutations affect membrane protein folding and function.

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.

Domain interactions determine the conformational ensemble of the periplasmic chaperone SurA.

SurA is thought to be the most important periplasmic chaperone for outer membrane protein (OMP) biogenesis. Its structure is composed of a core region and two peptidylprolyl isomerase domains, termed P1 and P2, connected by flexible linkers. As such these three independent folding units are able to adopt a number of distinct spatial positions with respect to each other. The conformational dynamics of these domains are thought to be functionally important yet are largely unresolved. Here we address this question of the conformational ensemble using sedimentation equilibrium, small-angle neutron scattering, and folding titrations. This combination of orthogonal methods converges on a SurA population that is monomeric at physiological concentrations. The conformation that dominates this population has the P1 and core domains docked to one another, for example, "P1-closed" and the P2 domain extended in solution. We discovered that the distribution of domain orientations is defined by modest and favorable interactions between the core domain and either the P1 or the P2 domains. These two peptidylprolyl domains compete with each other for core-binding but are thermodynamically uncoupled. This arrangement implies two novel insights. Firstly, an open conformation must exist to facilitate P1 and P2 exchange on the core, indicating that the open client-binding conformation is populated at low levels even in the absence of client unfolded OMPs. Secondly, competition between P1 and P2 binding paradoxically occludes the client binding site on the core, which may serve to preserve the reservoir of binding-competent apo-SurA in the periplasm.

Membrane Proteins Have Distinct Fast Internal Motion and Residual Conformational Entropy.

The internal motions of integral membrane proteins have largely eluded comprehensive experimental characterization. Here the fast side-chain dynamics of the α-helical sensory rhodopsin II and the ß-barrel outer membrane protein W have been investigated in lipid bilayers and detergent micelles by solution NMR relaxation techniques. Despite their differing topologies, both proteins have a similar distribution of methyl-bearing side-chain motion that is largely independent of membrane mimetic. The methyl-bearing side chains of both proteins are, on average, more dynamic in the ps-ns timescale than any soluble protein characterized to date. Accordingly, both proteins retain an extraordinary residual conformational entropy in the folded state, which provides a counterbalance to the absence of the hydrophobic effect. Furthermore, the high conformational entropy could greatly influence the thermodynamics underlying membrane-protein functions, including ligand binding, allostery, and signaling.

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.

Backbone Hydrogen Bond Energies in Membrane Proteins Are Insensitive to Large Changes in Local Water Concentration.

A hallmark feature of biological lipid bilayer structure is a depth-dependent polarity gradient largely resulting from the change in water concentration over the angstrom length scale. This gradient is particularly steep as it crosses the membrane interfacial regions where the water concentration drops at least a million-fold along the direction of the bilayer normal. Although local water content is often assumed to be a major determinant of membrane protein stability, the effect of the water-induced polarity gradient upon backbone hydrogen bond strength has not been systematically investigated. We addressed this question by measuring the free energy change for a number of backbone hydrogen bonds in the transmembrane protein OmpW. These values were obtained at 33 backbone amides from hydrogen/deuterium fractionation factors by nuclear magnetic resonance spectroscopy. We surprisingly found that OmpW backbone hydrogen bond energies do not vary over a wide range of water concentrations that are characteristic of the solvation environment in the bilayer interfacial region. We validated the interpretation of our results by determining the hydrodynamic and solvation properties of our OmpW-micelle complex using analytical ultracentrifugation and molecular dynamics simulations. The magnitudes of the backbone hydrogen bond free energy changes in our study are comparable to those observed in water-soluble proteins, the H-segment of the leader peptidase helix used in the von Heijne and White biological scale experiments, and several interfacial peptides. Our results agree with those reported for the transmembrane α-helical portion of the amyloid precursor protein after the latter values were adjusted for kinetic isotope effects. Overall, our work suggests that backbone hydrogen bonds provide modest thermodynamic stability to membrane protein structures and that many amides are unaffected by dehydration within the bilayer.

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.

Plasticity and transient binding are key ingredients of the periplasmic chaperone network.

SurA, Skp, FkpA, and DegP constitute a chaperone network that ensures biogenesis of outer membrane proteins (OMPs) in Gram-negative bacteria. Both Skp and FkpA are holdases that prevent the self-aggregation of unfolded OMPs, whereas SurA accelerates folding and DegP is a protease. None of these chaperones is essential, and we address here how functional plasticity is manifested in nine known null strains. Using a comprehensive computational model of this network termed OMPBioM, our results suggest that a threshold level of steady state holdase occupancy by chaperones is required, but the cell is agnostic to the specific holdase molecule fulfilling this function. In addition to its foldase activity, SurA moonlights as a holdase when there is no expression of Skp and FkpA. We further interrogate the importance of chaperone-client complex lifetime by conducting simulations using lifetime values for Skp complexes that range in length by six orders of magnitude. This analysis suggests that transient occupancy of durations much shorter than the Escherichia coli doubling time is required. We suggest that fleeting chaperone occupancy facilitates rapid sampling of the periplasmic conditions, which ensures that the cell can be adept at responding to environmental changes. Finally, we calculated the network effects of adding multivalency by computing populations that include two Skp trimers per unfolded OMP. We observe only modest perturbations to the system. Overall, this quantitative framework of chaperone-protein interactions in the periplasm demonstrates robust plasticity due to its dynamic binding and unbinding behavior.

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.

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties.

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.