Role of Aromatic Side Chains in the Folding and Thermodynamic Stability of Integral Membrane Proteins

Solvent Tolerance Improvement of Lipases Enhanced Their Applications: State of the Art

Lipases, crucial catalysts in biochemical synthesis, find extensive applications across industries such as food, medicine, and cosmetics. The efficiency of lipase-catalyzed reactions is significantly influenced by the choice of solvents. Polar organic solvents often result in a decrease, or even loss, of lipase activity. Conversely, nonpolar organic solvents induce excessive rigidity in lipases, thereby affecting their activity. While the advent of new solvents like ionic liquids and deep eutectic solvents has somewhat improved the activity and stability of lipases, it fails to address the fundamental issue of lipases' poor solvent tolerance. Hence, the rational design of lipases for enhanced solvent tolerance can significantly boost their industrial performance. This review provides a comprehensive summary of the structural characteristics and properties of lipases in various solvent systems and emphasizes various strategies of protein engineering for non-aqueous media to improve lipases' solvent tolerance. This study provides a theoretical foundation for further enhancing the solvent tolerance and industrial properties of lipases.

Conformational fingerprinting with Raman spectroscopy reveals protein structure as a translational biomarker of muscle pathology

Neuromuscular disorders are a group of conditions that can result in weakness of skeletal muscles. Examples include fatal diseases such as amyotrophic lateral sclerosis and conditions associated with high morbidity such as myopathies (muscle diseases). Many of these disorders are known to have abnormal protein folding and protein aggregates. Thus, easy to apply methods for the detection of such changes may prove useful diagnostic biomarkers. Raman spectroscopy has shown early promise in the detection of muscle pathology in neuromuscular disorders and is well suited to characterising the conformational profiles relating to protein secondary structure. In this work, we assess if Raman spectroscopy can detect differences in protein structure in muscle in the setting of neuromuscular disease. We utilise in vivo Raman spectroscopy measurements from preclinical models of amyotrophic lateral sclerosis and the myopathy Duchenne muscular dystrophy, together with ex vivo measurements of human muscle samples from individuals with and without myopathy. Using quantitative conformation profiling and matrix factorisation we demonstrate that quantitative 'conformational fingerprinting' can be used to identify changes in protein folding in muscle. Notably, myopathic conditions in both preclinical models and human samples manifested a significant reduction in α-helix structures, with concomitant increases in β-sheet and, to a lesser extent, nonregular configurations. Spectral patterns derived through non-negative matrix factorisation were able to identify myopathy with a high accuracy (79% in mouse, 78% in human tissue). This work demonstrates the potential of conformational fingerprinting as an interpretable biomarker for neuromuscular disorders.

MPA-Pred: A machine learning approach for predicting the binding affinity of membrane protein-protein complexes

Membrane protein-protein interactions are essential for several functions including cell signaling, ion transport, and enzymatic activity. These interactions are mainly dictated by their binding affinities. Although several methods are available for predicting the binding affinity of protein-protein complexes, there exists no specific method for membrane protein-protein complexes. In this work, we collected the experimental binding affinity data for a set of 114 membrane protein-protein complexes and derived several structure and sequence-based features. Our analysis on the relationship between binding affinity and the features revealed that the important factors mainly depend on the type of membrane protein and the functional class of the protein. Specifically, aromatic and charged residues at the interface, and aromatic-aromatic and electrostatic interactions are found to be important to understand the binding affinity. Further, we developed a method, MPA-Pred, for predicting the binding affinity of membrane protein-protein complexes using a machine learning approach. It showed an average correlation and mean absolute error of 0.83 and 0.91 kcal/mol, respectively, using the jack-knife test on a set of 114 complexes. We have also developed a web server and it is available at https://web.iitm.ac.in/bioinfo2/MPA-Pred/. This method can be used for predicting the affinity of membrane protein-protein complexes at a large scale and aid to improve drug design strategies.

Aspergillus fumigatus ZrfC Zn(II) transporter scavengers zincophore-bound Zn(II)

Aspergillus fumigatus is an opportunistic pathogen that is able to invade and grow in the lungs of immunosuppressed patients and cause invasive pulmonary aspergillosis. The concentration of free Zn(II) in living tissues is much lower than that required for optimal fungal growth; thus, to obtain Zn(II) from the host, Aspergillus fumigatus uses highly specified Zn(II) transporters: ZrfA, ZrfB and ZrfC. The ZrfC transporter plays the main role in Zn(II) acquisition from the host in neutral and mildly alkaline environment via interacting with the secreted Aspf2 zincophore. Understanding the Aspf2-ZrfC interactions is therefore necessary for explaining the process of Zn(II) acquisition by Aspergillus fumigatus, and identifying Zn(II) binding sites in its transporter and describing the thermodynamics of such binding are the fundamental steps to achieve this goal. We focus on two probable ZrfC Zn(II) binding sites and show that the Ac-MNCHFHAGVEHCIGAGESESGSSQ-NH2 region binds Zn(II) with higher affinity than the Ac-TGCHSHGS-NH2 one and that this binding is much stronger than the binding of Zn(II) to the Aspf2 zincophore, allowing efficient Zn(II) transport from the Aspf2 zincophore to the ZrfC transporter. The same ZrfC fragments also able to bind Ni(II), another metal ion essential for fungi that could also compete with Zn(II) binding, with comparable affinity.

Sculpting conducting nanopore size and shape through de novo protein design

Transmembrane β-barrels (TMBs) are widely used for single molecule DNA and RNA sequencing and have considerable potential for a broad range of sensing and sequencing applications. Current engineering approaches for nanopore sensors are limited to naturally occurring channels such as CsgG, which have evolved to carry out functions very different from sensing, and hence provide sub-optimal starting points. In contrast, de novo protein design can in principle create an unlimited number of new nanopores with any desired properties. Here we describe a general approach to the design of transmembrane β-barrel pores with different diameter and pore geometry. NMR and crystallographic characterization shows that the designs are stably folded with structures close to the design models. We report the first examples of de novo designed TMBs with 10, 12 and 14 stranded β-barrels. The designs have distinct conductances that correlate with their pore diameter, ranging from 110 pS (~0.5 nm pore diameter) to 430 pS (~1.1 nm pore diameter), and can be converted into sensitive small-molecule sensors with high signal to noise ratio. The capability to generate on demand β-barrel pores of defined geometry opens up fundamentally new opportunities for custom engineering of sequencing and sensing technologies.

Metal coordination to solute binding proteins - exciting chemistry with potential biological meaning

Zn(II) is essential for bacterial survival and virulence. In host cells, its abundance is extremely limited, thus, bacteria have evolved transport mechanisms that enable them to take up this essential metal nutrient. Paracoccus denitrificans encodes two solute binding proteins (SBPs) - ZnuA and AztC, which are responsible for zinc acquisition from the host cells. We focus on understanding the interactions of Zn(II) and Ni(II) (zinc's potential competitor, which is a biologically relevant metal ion essential for various bacterial enzymes) with the extracellular ZnuA and AztC's loops from P. denitrificans that are expected to be possible Zn(II) binding sites. In the case of Zn(II) complexes with ZnuA outercellular loop regions, the numerous histidines act as anchoring donors, forming complexes with up to four coordinated His residues, while in the AztC region, three imidazole nitrogens and one water molecule are involved in Zn(II) binding. In Zn(II) complexes with ZnuA His-rich loop regions, so-called polymorphic binding sites are observed. The large number of available imidazoles and carboxylic side chains also strongly affects the structure of Ni(II) complexes; the more histidines in the studied peptide, the higher the affinity to bind Ni(II) and the higher the pH value at which amide nitrogens start to participate in Ni(II) binding. Additionally, for Ni(II)-ZnuA complexes, a more rare octahedral geometry is observed and such complexes are more stable than the corresponding Zn(II) ones, in contrast to what was observed in the AztC region, suggesting that the numerous histidyl and glutamic acid side chains are more tempting for Ni(II) than for Zn(II).The general strong affinity of Zn(II)-zincophore complexes is also discussed.

In-silico prediction of TGF-β1 non-synonymous variants and their impact on binding affinity to Fresolimumab

TGF-β1 is a potent immunoregulatory cytokine that plays diverse roles in development, bone healing, fibrosis, and cancer. However, characterizing TGF-β1 gene variants is challenging because the structural and functional consequences of these variants are still undetermined. In this study, we aimed to perform an in-silico analysis of TGF-β1 non-synonymous variants and their pathogenic effects on the TGF-β1 protein. A total of 10,252 TGF-β1 SNPs were collected from the NCBI dbSNP database and in-silico tools (SIFT, PROVEAN, Mutation Taster, ClinVar, PolyPhen-2, CScape, MutPred, and ConSurf) were used. The in-silico predicted potential variants were further investigated for their binding to the TGF-β1 targeting drug "Fresolimumab". Molecular docking was performed using HADDOCK and confirmed by PRODIGY and PDBsum. The in-silico analysis predicted four potential TGF-β1 nsSNPs: E47G in the LAP domain of the propeptide and I22T, L28F, and E35D in the mature TGF-β1 peptide. HADDOCK and molecular dynamics simulations revealed that the I22T and E35D variants have higher binding affinity for Fresolimumab as compared to the wild type and L28F variants. Molecular dynamics simulations (100 ns) and principal component analysis showed that TGF-β1 variants influenced the protein structure and caused variations in the internal dynamics of protein complexes with the antibody. Among them, the E35D variant significantly destabilized the TGF-β1 protein structure, resulting in rearrangement in the binding site and affecting the interactions with the Fresolimumab. This study identified four variants that can affect the TGF-β1 protein structure and result in functional consequences such as impaired response to Fresolimumab.Communicated by Ramaswamy H. Sarma.

Cell-Free Expression of De Novo Designed Peptides That Form β-Barrel Nanopores

Nanopore sensing has attracted much attention as a rapid, simple, and label-free single-molecule detection technology. To apply nanopore sensing to extensive targets including polypeptides, nanopores are required to have a size and structure suitable for the target. We recently designed a de novo β-barrel peptide nanopore (SVG28) that constructs a stable and monodispersely sized nanopore. To develop the sizes and functionality of peptide nanopores, systematic exploration is required. Here we attempt to use a cell-free synthesis system that can readily express peptides using transcription and translation. Hydrophilic variants of SVG28 were designed and expressed by the PURE system. The peptides form a monodispersely sized nanopore, with a diameter 1.1 or 1.5 nm smaller than that of SVG28. Such cell-free synthesizable peptide nanopores have the potential to enable the systematic custom design of nanopores and comprehensive sequence screening of nanopore-forming peptides.

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.

Thermodynamic analysis of the GASright transmembrane motif supports energetic model of dimerization

The GASright motif, best known as the fold of the glycophorin A transmembrane dimer, is one of the most common dimerization motifs in membrane proteins, characterized by its hallmark GxxxG-like sequence motifs (GxxxG, AxxxG, GxxxS, and similar). Structurally, GASright displays a right-handed crossing angle and short interhelical distance. Contact between the helical backbones favors the formation of networks of weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O=C). To understand the factors that modulate the stability of GASright, we previously presented a computational and experimental structure-based analysis of 26 predicted dimers. We found that the contributions of van der Waals packing and Cα-H hydrogen bonding to stability, as inferred from the structural models, correlated well with relative dimerization propensities estimated experimentally with the in vivo assay TOXCAT. Here we test this model with a quantitative thermodynamic analysis. We used Förster resonance energy transfer (FRET) to determine the free energy of dimerization of a representative subset of seven of the 26 original TOXCAT dimers using FRET. To overcome the technical issue arising from limited sampling of the dimerization isotherm, we introduced a globally fitting strategy across a set of constructs comprising a wide range of stabilities. This strategy yielded precise thermodynamic data that show strikingly good agreement between the original propensities and ΔG° of association in detergent, suggesting that TOXCAT is a thermodynamically driven process. From the correlation between TOXCAT and thermodynamic stability, the predicted free energy for all the 26 GASright dimers was calculated. These energies correlate with the in silico ΔE scores of dimerization that were computed on the basis of their predicted structure. These findings corroborate our original model with quantitative thermodynamic evidence, strengthening the hypothesis that van der Waals and Cα-H hydrogen bond interactions are the key modulators of GASright stability.

Investigating the functional role of a buried interchain aromatic cluster in Escherichia coli GrpE dimer

Aromatic clusters in the core of proteins are often involved in imparting structural stability to proteins. However, their functional importance is not always clear. In this study, we investigate the thermosensing role of a phenylalanine cluster present in the GrpE homodimer. GrpE, which acts as a nucleotide exchange factor for the molecular chaperone DnaK, is well known for its thermosensing activity resulting from temperature-dependent structural changes that allow control of chaperone function. Using mutational analysis, we show that an interchain phenylalanine cluster in a four-helix bundle of the GrpE homodimer assists in the thermosensing ability of the co-chaperone. Substitution of aromatic residues with hydrophobic ones in the core of the four-helix bundle reduces the thermal stability of the bundle and that of a connected coiled-coil domain, which impacts thermosensing. Cell growth assays and SEM images of the mutants show filamentous growth of Escherichia coli cells at 42°C, which corroborates with the defect in thermosensing. Our work suggests that the interchain edge-to-face aromatic cluster is important for the propagation of the structural signal from the coiled-coil domain to the four-helical bundle of GrpE, thus facilitating GrpE-mediated thermosensing in bacteria.

Multiplex Chemical Labeling of Amino Acids for Protein Footprinting Structure Assessment

Protein footprinting with mass spectrometry is an established structural biology technique for mapping solvent accessibility and assessing molecular-level interactions of proteins. In hydroxyl radical protein footprinting (HRPF), hydroxyl (OH) radicals generated by water radiolysis or other methods covalently label protein side chains. Because of the wide dynamic range of OH reactivity, not all side chains are easily detected in a single experiment. Novel reagent development and the use of radical chain reactions for labeling, including trifluoromethyl radicals, is a potential approach to normalize the labeling across a diverse set of residues. HRPF in the presence of a trifluoromethylation reagent under the right conditions could provide a "one-pot" reaction for multiplex labeling of protein side chains. Toward this goal, we have systematically evaluated amino acid labeling with the recently investigated Langlois' reagent (LR) activated by X-ray-mediated water radiolysis, followed by three different mass spectrometry methods. We compared the reactivity of CF₃ and OH radical labeling for all 20 protein side chains in a competition-free environment. We found that all 20 amino acids exhibited CF₃ or OH labeling in LR. Our investigations provide the evidence and knowledge set to perfect hydroxyl radical-activated trifluoromethyl chemistry as "one-pot" reaction for multiplex labeling of protein side chains to achieve higher resolution in HRPF.

How physical forces drive the process of helical membrane protein folding

Protein folding is a fundamental process of life with important implications throughout biology. Indeed, tens of thousands of mutations have been associated with diseases, and most of these mutations are believed to affect protein folding rather than function. Correct folding is also a key element of design. These factors have motivated decades of research on protein folding. Unfortunately, knowledge of membrane protein folding lags that of soluble proteins. This gap is partly caused by the greater technical challenges associated with membrane protein studies, but also because of additional complexities. While soluble proteins fold in a homogenous water environment, membrane proteins fold in a setting that ranges from bulk water to highly charged to apolar. Thus, the forces that drive folding vary in different regions of the protein, and this complexity needs to be incorporated into our understanding of the folding process. Here, we review our understanding of membrane protein folding biophysics. Despite the greater challenge, better model systems and new experimental techniques are starting to unravel the forces and pathways in membrane protein folding.

Transmembrane β-barrel proteins of bacteria: From structure to function

The outer membrane of Gram-negative bacteria is a specialized organelle conferring protection to the cell against various environmental stresses and resistance to many harmful compounds. The outer membrane has a number of unique features, including an asymmetric lipid bilayer, the presence of lipopolysaccharides and an individual proteome. The vast majority of the integral transmembrane proteins in the outer membrane belongs to the family of β-barrel proteins. These evolutionarily related proteins share a cylindrical, anti-parallel β-sheet core fold spanning the outer membrane. The loops and accessory domains attached to the β-barrel allow for a remarkable versatility in function for these proteins, ranging from diffusion pores and transporters to enzymes and adhesins. We summarize the current knowledge on β-barrel structure and folding and give an overview of their functions, evolution, and potential as drug targets.

A pH-responsive, endosomolytic liposome functionalized with membrane-anchoring, comb-like pseudopeptides for enhanced intracellular delivery and cancer treatment

Low intracellular delivery efficiency and multidrug resistance are among major barriers to effective cancer therapy.

De novo design of a nanopore for single-molecule detection that incorporates a β-hairpin peptide

The amino-acid sequence of a protein encodes information on its three-dimensional structure and specific functionality. De novo design has emerged as a method to manipulate the primary structure for the development of artificial proteins and peptides with desired functionality. This paper describes the de novo design of a pore-forming peptide, named SV28, that has a β-hairpin structure and assembles to form a stable nanopore in a bilayer lipid membrane. This large synthetic nanopore is an entirely artificial device for practical applications. The peptide forms multidispersely sized nanopore structures ranging from 1.7 to 6.3 nm in diameter and can detect DNAs. To form a monodispersely sized nanopore, we redesigned the SV28 by introducing a glycine-kink mutation. The resulting redesigned peptide forms a monodisperse pore with a diameter of 1.7 nm leading to detection of a single polypeptide chain. Such de novo design of a β-hairpin peptide has the potential to create artificial nanopores, which can be size adjusted to a target molecule.

Diethylpyrocarbonate Footprints a Membrane Protein in Micelles

Membrane proteins play crucial roles in cell signaling and transport and, thus, are the targets of many small molecule drugs. The characterization of membrane protein structures poses challenges for the high-resolution biophysical tools because the transmembrane (TM) domain is hydrophobic, opening an opportunity for mass spectrometry (MS)-based footprinting. The hydrophobic reagent diethylpyrocarbonate (DEPC), a heavily studied footprinter for water-soluble proteins, can label up to 30% of surface residues via a straightforward protocol, streamlining the MS-based footprinting workflow. To test its applicability to membrane proteins, we footprinted vitamin K epoxide reductase (VKOR) membrane protein with DEPC. The results demonstrate that besides labeling the hydrophilic extracellular (extramembrane (EM)) domain, DEPC can also diffuse into the hydrophobic TM domain and subsequently label that region. The labeling process was facilitated by tip sonication to enhance reagent diffusion into micelles. We then analyzed the correlation between the residue modification extent and the theoretical accessible surface area percentage (%ASA); the data generally show good correlation with the residue location. Compared with conventional hydrophilic footprinters, the relatively hydrophobic DEPC can map a membrane protein's TM domain, suggesting that the reagent's hydrophobicity can be exploited to obtain structural information on the membrane-spanning region. This encouraging result should assist in the development of more efficient footprinters for membrane protein TM domain footprinting, enabled by further understanding the relationship between a reagent's hydrophobicity and its preferred labeling sites.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ₁, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Bacterial Outer Membrane Proteins Are Targeted to the Bam Complex by Two Parallel Mechanisms

Membrane proteins that are integrated into the outer membrane of Gram-negative bacteria typically contain a unique "β barrel" structure that serves as a membrane spanning segment. A conserved "β signal" motif is located at the C terminus of the β barrel of many outer membrane proteins (OMPs), but the function of this sequence is unclear. We found that mutations in the β signal slightly delayed the assembly of three model Escherichia coli OMPs by reducing their affinity for the barrel assembly machinery (Bam) complex, a heterooligomer that catalyzes β barrel insertion, and led to the degradation of a fraction of the protein in the periplasm. Interestingly, the absence of the periplasmic chaperone SurA amplified the effect of the mutations and caused the complete degradation of the mutant proteins. In contrast, the absence of another periplasmic chaperone (Skp) suppressed the effect of the mutations and considerably enhanced the efficiency of assembly. Our results reveal the existence of two parallel OMP targeting mechanisms that rely on a cis-acting peptide (the β signal) and a trans-acting factor (SurA), respectively. Our results also challenge the long-standing view that periplasmic chaperones are redundant and provide evidence that they have specialized functions.IMPORTANCE Proteins that are embedded in the outer membrane of Gram-negative bacteria (OMPs) play an important role in protecting the cell from harmful chemicals. OMPs share a common architecture and often contain a conserved sequence motif (β motif) of unknown function. Although OMPs are escorted to the outer membrane by proteins called chaperones, the exact function of the chaperones is also unclear. Here, we show that the β motif and the chaperone SurA both target OMPs to the β barrel insertion machinery in the outer membrane. In contrast, the chaperone Skp delivers unintegrated OMPs to protein degradation complexes. Our results challenge the long-standing view that chaperones are functionally redundant and strongly suggest that they have specialized roles in OMP targeting and quality control.

The antibiotic darobactin mimics a β-strand to inhibit outer membrane insertase

Antibiotics that target Gram-negative bacteria in new ways are needed to resolve the antimicrobial resistance crisis^1-3. Gram-negative bacteria are protected by an additional outer membrane, rendering proteins on the cell surface attractive drug targets^4,5. The natural compound darobactin targets the bacterial insertase BamA⁶-the central unit of the essential BAM complex, which facilitates the folding and insertion of outer membrane proteins^7-13. BamA lacks a typical catalytic centre, and it is not obvious how a small molecule such as darobactin might inhibit its function. Here we resolve the mode of action of darobactin at the atomic level using a combination of cryo-electron microscopy, X-ray crystallography, native mass spectrometry, in vivo experiments and molecular dynamics simulations. Two cyclizations pre-organize the darobactin peptide in a rigid β-strand conformation. This creates a mimic of the recognition signal of native substrates with a superior ability to bind to the lateral gate of BamA. Upon binding, darobactin replaces a lipid molecule from the lateral gate to use the membrane environment as an extended binding pocket. Because the interaction between darobactin and BamA is largely mediated by backbone contacts, it is particularly robust against potential resistance mutations. Our results identify the lateral gate as a functional hotspot in BamA and will allow the rational design of antibiotics that target this bacterial Achilles heel.

Free-Radical Membrane Protein Footprinting by Photolysis of Perfluoroisopropyl Iodide Partitioned to Detergent Micelle by Sonication

A free-radical footprinting approach is described for integral membrane protein (IMP) that extends, significantly, the "fast photochemical oxidation of proteins" (FPOP) platform. This new approach exploits highly hydrophobic perfluoroisopropyl iodide (PFIPI) together with tip sonication to ensure efficient transport into the micelle interior, allowing laser dissociation and footprinting of the transmembrane domains. In contrast to water soluble footprinters, PFIPI footprints both the hydrophobic intramembrane and the hydrophilic extramembrane domains of the IMP vitamin K epoxide reductase (VKOR). The footprinting is fast, giving high coverage for Tyr (100 %) and Trp. The incorporation of the reagent with sonication does not significantly affect VKOR's enzymatic function, and tyrosine iodination does not compromise protease digestion and the subsequent analysis. The locations for the modifications are largely consistent with the corresponding solvent accessibilities, recommending this approach for future membrane protein footprinting.

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.

Outer membrane protein evolution

Outer membrane proteins have remarkably homogeneous structure. They are all up down β-barrels. Up down barrels themselves are composed of repeated sets of β-hairpins. The consistency of the usage of the β-hairpin throughout the outer membrane milieu allows for interrogation of the evolution of these repetitive structures. Here we describe recent investigations of outer membrane protein evolution and how evolutionary precepts have been used for novel outer membrane protein design.

The Uniqueness of Tryptophan in Biology: Properties, Metabolism, Interactions and Localization in Proteins

Tryptophan (Trp) holds a unique place in biology for a multitude of reasons. It is the largest of all twenty amino acids in the translational toolbox. Its side chain is indole, which is aromatic with a binuclear ring structure, whereas those of Phe, Tyr, and His are single-ring aromatics. In part due to these elaborate structural features, the biosynthetic pathway of Trp is the most complex and the most energy-consuming among all amino acids. Essential in the animal diet, Trp is also the least abundant amino acid in the cell, and one of the rarest in the proteome. In most eukaryotes, Trp is the only amino acid besides Met, which is coded for by a single codon, namely UGG. Due to the large and hydrophobic π-electron surface area, its aromatic side chain interacts with multiple other side chains in the protein, befitting its strategic locations in the protein structure. Finally, several Trp derivatives, namely tryptophylquinone, oxitriptan, serotonin, melatonin, and tryptophol, have specialized functions. Overall, Trp is a scarce and precious amino acid in the cell, such that nature uses it parsimoniously, for multiple but selective functions. Here, the various aspects of the uniqueness of Trp are presented in molecular terms.

Tyrosine in the hinge region of the pore-forming motif regulates oligomeric β-barrel pore formation by Vibrio cholerae cytolysin

β-barrel pore-forming toxins perforate cell membranes by forming oligomeric β-barrel pores. The most crucial step is the membrane-insertion of the pore-forming motifs that create the transmembrane β-barrel scaffold. Molecular mechanism that regulates structural reorganization of these pore-forming motifs during β-barrel pore-formation still remains elusive. Using Vibrio cholerae cytolysin as an archetypical example of the β-barrel pore-forming toxin, we show that a key tyrosine residue (Y321) in the hinge region of the pore-forming motif plays crucial role in this process. Mutation of Y321 abrogates oligomerization of the membrane-bound toxin protomers, and blocks subsequent steps of pore-formation. Our study suggests that the presence of Y321 in the hinge region of the pore-forming motif is crucial for the toxin molecule to sense membrane-binding, and to trigger essential structural rearrangements required for the subsequent oligomerization and pore-formation process. Such a regulatory mechanism of pore-formation by V. cholerae cytolysin has not been documented earlier in the structurally related β-barrel pore-forming toxins.

How bilayer properties influence membrane protein folding

The question of how proteins manage to organize into a unique three-dimensional structure has been a major field of study since the first protein structures were determined. For membrane proteins, the question is made more complex because, unlike water-soluble proteins, the solvent is not homogenous or even unique. Each cell and organelle has a distinct lipid composition that can change in response to environmental stimuli. Thus, the study of membrane protein folding requires not only understanding how the unfolded chain navigates its way to the folded state, but also how changes in bilayer properties can affect that search. Here we review what we know so far about the impact of lipid composition on bilayer physical properties and how those properties can affect folding. A better understanding of the lipid bilayer and its effects on membrane protein folding is not only important for a theoretical understanding of the folding process, but can also have a practical impact on our ability to work with and design membrane proteins.

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.

Interaction of aliphatic amino acids with zwitterionic and charged lipid membranes: hydration and dehydration phenomena

Aliphatic amino acids interact differently in order to induce gelation or fluidization in zwitterionic and charged lipid membranes as a result of hydration or dehydration of the membrane surface.

Single-residue physicochemical characteristics kinetically partition membrane protein self-assembly and aggregation

Ninety-five percent of all transmembrane proteins exist in kinetically trapped aggregation-prone states that have been directly linked to neurodegenerative diseases. Interestingly, the primary sequence almost invariably avoids off-pathway aggregate formation, by folding reliably into its native, thermodynamically stabilized structure. However, with the rising incidence of protein aggregation diseases, it is now important to understand the underlying mechanism(s) of membrane protein aggregation. Micromolecular physicochemical and biochemical alterations in the primary sequence that trigger the formation of macromolecular cross-β aggregates can be measured only through combinatorial spectroscopic experiments. Here, we developed spectroscopic thermal perturbation with 117 experimental variables to assess how subtle protein sequence variations drive the molecular transition of the folded protein to oligomeric aggregates. Using the Yersinia pestis outer transmembrane β-barrel Ail as a model, we delineated how a single-residue substitution that alters the membrane-anchoring ability of Ail significantly contributes to the kinetic component of Ail stability. We additionally observed a stabilizing role for interface aliphatics, and that interface aromatics physicochemically contribute to Ail self-assembly and aggregation. Moreover, our method identified the formation of structured oligomeric intermediates during Ail aggregation. We show that the self-aggregation tendency of Ail is offset by the evolution of a thermodynamically compromised primary sequence that balances folding, stability, and oligomerization. Our approach provides critical information on how subtle changes in protein primary sequence trigger cross-β fibril formation, with insights that have direct implications for deducing the molecular progression of neurodegeneration and amyloidogenesis in humans.

Planarity and out-of-plane vibrational modes of tryptophan and tyrosine in biomolecular modeling

Tryptophan and tyrosine are amino acids that play significant roles in the folding processes of proteins at water-membrane interfaces because of their amphipathic heteroaromatic rings. Employing appropriate heteroaromatic molecular structures is essential for obtaining accurate dynamics and predictive capabilities in molecular simulations of these amino acids. In this study, molecular dynamics simulations that applied the most recent version of the CHARMM36 force field were conducted on aqueous solutions of tryptophan and of tyrosine. Geometric analysis and dynamics quantified how aromatic rings deviated from planar structures and exhibited out-of-plane fluctuations. Radial distribution functions showed possible biological significance because the extent of ring planarity slightly affected local water concentrations near aromatic rings. Instantaneous all-atom normal mode analysis (NMA) and Fourier transformation of time autocorrelation functions of out-of-plane displacements were applied to study out-of-plane vibrations of atoms in these rings. The NMA started with minimum energy configurations and then averaged over fluctuations in aqueous solution. The frequencies and frequency patterns that were obtained for tryptophan and tyrosine with CHARMM36 differed from literature reports of Raman spectra, infrared spectra, and frequencies calculated using quantum mechanics, with some out-of-plane modes found at higher frequencies. Effects of imposing improper torsion potentials and changing torsion angle force constants were investigated for all atoms in the rings of tryptophan and tyrosine. Results show that these coarse force field variations only affect planarity and out-of-plane vibrations of atoms within the rings, and not other vibrations. Although increasing improper torsion force constants reduced deviations from aromatic ring planarity significantly, it increased out-of-plane mode frequencies. Reducing torsion angle force constants (with and without improper torsions) shifted modes to lower frequencies. A combination of decreasing most torsion angle force constants for ring atoms in both amino acids and including improper torsion forces attained frequencies and frequency patterns for out-of-plane normal modes that were more similar to the literature spectra. These force field variations decreased the extents of out-of-plane vibrations within the heteroaromatic rings of tryptophan, especially around the nitrogen atom in the ring, but not within the heteroaromatic ring of tyrosine. Conclusions were unaffected by the peptide endgroup, water, or simulation ensemble.

Reversible folding energetics of Yersinia Ail barrel reveals a hyperfluorescent intermediate

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

A comprehensive computational study of amino acid interactions in membrane proteins

Transmembrane proteins play a fundamental role in a wide series of biological processes but, despite their importance, they are less studied than globular proteins, essentially because their embedding in lipid membranes hampers their experimental characterization. In this paper, we improved our understanding of their structural stability through the development of new knowledge-based energy functions describing amino acid pair interactions that prevail in the transmembrane and extramembrane regions of membrane proteins. The comparison of these potentials and those derived from globular proteins yields an objective view of the relative strength of amino acid interactions in the different protein environments, and their role in protein stabilization. Separate potentials were also derived from α-helical and β-barrel transmembrane regions to investigate possible dissimilarities. We found that, in extramembrane regions, hydrophobic residues are less frequent but interactions between aromatic and aliphatic amino acids as well as aromatic-sulfur interactions contribute more to stability. In transmembrane regions, polar residues are less abundant but interactions between residues of equal or opposite charges or non-charged polar residues as well as anion-π interactions appear stronger. This shows indirectly the preference of the water and lipid molecules to interact with polar and hydrophobic residues, respectively. We applied these new energy functions to predict whether a residue is located in the trans- or extramembrane region, and obtained an AUC score of 83% in cross validation, which demonstrates their accuracy. As their application is, moreover, extremely fast, they are optimal instruments for membrane protein design and large-scale investigations of membrane protein stability.

Nano-viscosimetry analysis of the membrane disrupting action of the bee venom peptide melittin

Melittin is one of the most studied α-helical cationic membrane disrupting peptides. It is the main component of bee venom, however it is considered an antimicrobial peptide for its ability to kill bacteria. Melittin is believed to act by opening large toroidal pores in the plasma membrane of the targeted cells/bacteria, although this is questioned by some authors. Little is known, however, about the molecular mechanism leading to this activity. In this study the mechanism of action of melittin was studied by dye leakage and quartz crystal microbalance fingerprinting analysis in biomimetic model membranes. The results revealed the existence of multiple stages in the membrane disrupting action with characteristic differences between different membrane types. In bacterial-mimetic (charged) lipid mixtures the viscoelastic fingerprints suggest a surface-acting mechanism, whereas in mammalian-mimetic (neutral) membranes melittin appears to penetrate the bilayer already at low concentrations. In domain-forming mixed membranes melittin shows a preference for the domain containing predominantly zwitterionic lipids. The results confirm membrane poration but are inconsistent with the insertion-to-toroidal pore pathway. Therefore hypotheses of the two membrane disrupting pathways were developed, describing the membrane disruption as either surface tension modulation leading to toroidal pore formation, or linear aggregation leading to fissure formation in the membrane.

A conserved and buried edge-to-face aromatic interaction in small ubiquitin-like modifier (SUMO) has a role in SUMO stability and function

Aromatic amino acids buried at a protein's core are often involved in mutual paired interactions. Ab initio energy calculations have highlighted that the conformational orientations and the effects of substitutions are important for stable aromatic interactions among aromatic rings, but studies in the context of a protein's fold and function are elusive. Small ubiquitin-like modifier (SUMO) is a common post-translational modifier that affects diverse cellular processes. Here, we report that a highly conserved aromatic triad of three amino acids, Phe36-Tyr51-Phe64, is a unique SUMO signature that is absent in other ubiquitin-like homologous folds. We found that a specific edge-to-face conformation between the Tyr51-Phe64 pair of interacting aromatics is vital to the fold and stability of SUMO. Moreover, the noncovalent binding of SUMO-interacting motif (SIM) at the SUMO surface was critically dependent on the paired aromatic interactions buried at the core. NMR structural studies revealed that perturbation of the Tyr51-Phe64 conformation disrupts several long-range tertiary contacts in SUMO, leading to a heterogeneous and dynamic protein with attenuated SUMOylation both in vitro and in cells. A subtle perturbation of the edge-to-face conformation by a Tyr to Phe substitution significantly decreased stability, SUMO/SIM affinity, and the rate of SUMOylation. Our results highlight that absolute co-conservation of specific aromatic pairs inside the SUMO protein core has a role in its stability and function.

Covalently Linking Oligomerization-Impaired GlpF Protomers Does Not Completely Re-establish Wild-Type Channel Activity

Integral membrane proteins of the aquaporin family facilitate rapid water flux across cellular membranes in all domains of life. Although the water-conducting pore is clearly defined in an aquaporin monomer, all aquaporins assemble into stable tetramers. In order to investigate the role of protomer–protomer interactions, we analyzed the activity of heterotetramers containing increasing fractions of mutated monomers, which have an impaired oligomerization propensity and activity. In order to enforce interaction between the protomers, we designed and analyzed a genetically fused homotetramer of GlpF, the aquaglyceroporin of the bacterium Escherichia coli (E. coli). However, increasing fractions of the oligomerization-impaired mutant GlpF E43A affected the activity of the GlpF heterotetramer in a nearly linear manner, indicating that the reduced protein activity, caused by the introduced mutations, cannot be fully compensated by simply covalently linking the monomers. Taken together, the results underline the importance of exactly positioned monomer–monomer contacts in an assembled GlpF tetramer.

Filling the Void: Introducing Aromatic Interactions into Solvent Tunnels To Enhance Lipase Stability in Methanol

An enhanced stability of enzymes in organic solvents is desirable under industrial conditions. The potential of lipases as biocatalysts is mainly limited by their denaturation in polar alcohols. In this study, we focused on selected solvent tunnels in lipase from Geobacillus stearothermophilus T6 to improve its stability in methanol during biodiesel synthesis. Using rational mutagenesis, bulky aromatic residues were incorporated to occupy solvent channels and induce aromatic interactions leading to a better inner core packing. The chemical and structural characteristics of each solvent tunnel were systematically analyzed. Selected residues were replaced with Phe, Tyr, or Trp. Overall, 16 mutants were generated and screened in 60% methanol, from which 3 variants showed an enhanced stability up to 81-fold compared with that of the wild type. All stabilizing mutations were found in the longest tunnel detected in the "closed-lid" X-ray structure. The combination of Phe substitutions in an A187F/L360F double mutant resulted in an increase in unfolding temperature (T_m ) of 7°C in methanol and a 3-fold increase in biodiesel synthesis yield from waste chicken oil. A kinetic analysis with p-nitrophenyl laurate revealed that all mutants displayed lower hydrolysis rates (k _cat), though their stability properties mostly determined the transesterification capability. Seven crystal structures of different variants were solved, disclosing new π-π or CH/π intramolecular interactions and emphasizing the significance of aromatic interactions for improved solvent stability. This rational approach could be implemented for the stabilization of other enzymes in organic solvents.IMPORTANCE Enzymatic synthesis in organic solvents holds increasing industrial opportunities in many fields; however, one major obstacle is the limited stability of biocatalysts in such a denaturing environment. Aromatic interactions play a major role in protein folding and stability, and we were inspired by this to redesign enzyme voids. The rational protein engineering of solvent tunnels of lipase from Geobacillus stearothermophilus is presented here, offering a promising approach to introduce new aromatic interactions within the enzyme core. We discovered that longer tunnels leading from the surface to the enzyme active site were more beneficial targets for mutagenesis for improving lipase stability in methanol during biodiesel biosynthesis. A structural analysis of the variants confirmed the generation of new interactions involving aromatic residues. This work provides insights into stability-driven enzyme design by targeting the solvent channel void.

The rhomboid protease GlpG has weak interaction energies in its active site hydrogen bond network

Rhomboid proteases are membrane-integrated enzymes that hydrolyze peptide bonds in the transmembrane domains of protein substrates. Gaffney and Hong experimentally determine interaction energies between active site residues to reveal weak coupling, which may explain the slow proteolysis mediated by GlpG.

Intramembrane rhomboid proteases are of particular interest because of their function to hydrolyze a peptide bond of a substrate buried in the membrane. Crystal structures of the bacterial rhomboid protease GlpG have revealed a catalytic dyad (Ser201-His254) and oxyanion hole (His150/Asn154/the backbone amide of Ser201) surrounded by the protein matrix and contacting a narrow water channel. Although multiple crystal structures have been solved, the catalytic mechanism of GlpG is not completely understood. Because it is a serine protease, hydrogen bonding interactions between the active site residues are thought to play a critical role in the catalytic cycle. Here, we dissect the interaction energies among the active site residues His254, Ser201, and Asn154 of Escherichia coli GlpG, which form a hydrogen bonding network. We combine double mutant cycle analysis with stability measurements using steric trapping. In mild detergent, the active site residues are weakly coupled with interaction energies (ΔΔG_Inter) of ‒1.4 kcal/mol between His254 and Ser201 and ‒0.2 kcal/mol between Ser201 and Asn154. Further, by analyzing the propagation of single mutations of the active site residues, we find that these residues are important not only for function but also for the folding cooperativity of GlpG. The weak interaction between Ser and His in the catalytic dyad may partly explain the unusually slow proteolysis by GlpG compared with other canonical serine proteases. Our result suggests that the weak hydrogen bonds in the active site are sufficient to carry out the proteolytic function of rhomboid proteases.

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.

Lipid bilayer composition modulates the unfolding free energy of a knotted α-helical membrane protein

α-Helical membrane proteins have eluded investigation of their thermodynamic stability in lipid bilayers. Reversible denaturation curves have enabled some headway in determining unfolding free energies. However, these parameters have been limited to detergent micelles or lipid bicelles, which do not possess the same mechanical properties as lipid bilayers that comprise the basis of natural membranes. We establish reversible unfolding of the membrane transporter LeuT in lipid bilayers, enabling the comparison of apparent unfolding free energies in different lipid compositions. LeuT is a bacterial ortholog of neurotransmitter transporters and contains a knot within its 12-transmembrane helical structure. Urea is used as a denaturant for LeuT in proteoliposomes, resulting in the loss of up to 30% helical structure depending upon the lipid bilayer composition. Urea unfolding of LeuT in liposomes is reversible, with refolding in the bilayer recovering the original helical structure and transport activity. A linear dependence of the unfolding free energy on urea concentration enables the free energy to be extrapolated to zero denaturant. Increasing lipid headgroup charge or chain lateral pressure increases the thermodynamic stability of LeuT. The mechanical and charge properties of the bilayer also affect the ability of urea to denature the protein. Thus, we not only gain insight to the long-sought-after thermodynamic stability of an α-helical protein in a lipid bilayer but also provide a basis for studies of the folding of knotted proteins in a membrane environment.

Membrane Anchoring of α-Helical Proteins: Role of Tryptophan

The function of membrane proteins relies on a defined orientation of protein relative to lipid. In apparent correlation to protein anchoring, tryptophan residues are enriched in the lipid headgroup region. To characterize the thermodynamic and structural basis of this relationship in α-helical membrane proteins, we examined the role of three conserved tryptophans in the folding of the heterodimeric integrin αIIbβ3 transmembrane (TM) complex in phospholipid bicelles and mammalian membranes. In the homogenous lipid environment of bicelles, tryptophan was replaceable by residues of distinct polarities. The appropriate polarity was guided by the electrostatic potential of the tryptophan surrounding, suggesting that tryptophan can complement diverse environments by adjusting the orientation of its anisotropic side chain to achieve site-specific anchoring. As a sole membrane anchor, tryptophan made a contribution of 0.4 kcal/mol to TM complex stability in bicelles. In membranes, it proved more difficult to replace tryptophan even by tyrosine, indicating a superior capacity to interact with heterogeneous lipids of biological membranes. Interestingly, at intracellular TM helix ends, where integrin activation is initiated, sequence motifs that interact with lipids via opposing polarity patterns were found to restrict TM helix orientations beyond tryptophan anchoring. In contrast to bicelles, phenylalanine became the least accepted substitute in membranes, demonstrating an increased role of the hydrophobic effect. Altogether, our study implicates a wide amphiphilic range of tryptophan, membrane complexity, and the hydrophobic effect to be important factors in tryptophan membrane anchoring.

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.

The GlpF residue Trp219 is part of an amino-acid cluster crucial for aquaglyceroporin oligomerization and function

The vestibule loop regions of aquaglyceroporins are involved in accumulation of glycerol inside the channel pore. Even though most loop regions do not show high sequence similarity among aquaglyceroporins, loop E is highly conserved in aquaglyceroporins, but not in members of the homologous aquaporins. Specifically, a tryptophan residue is extremely conserved within this loop. We have investigated the role of this residue (Trp219) that deeply protrudes into the protein and potentially interacts with adjacent loops, using the E. coli aqualgyeroporin GlpF as a model. Replacement of Trp219 affects the activity of GlpF and impairs the stability of the tetrameric protein. Furthermore, we have identified an amino acid cluster involving Trp219 that stabilizes the GlpF tetramer. Based on our results we propose that Trp219 is key for formation of a defined vestibule structure, which is crucial for glycerol accumulation as well as for the stability of the active GlpF tetramer.

How Many Protein Sequences Fold to a Given Structure? A Coevolutionary Analysis

Quantifying the relationship between protein sequence and structure is key to understanding the protein universe. A fundamental measure of this relationship is the total number of amino acid sequences that can fold to a target protein structure, known as the "sequence capacity," which has been suggested as a proxy for how designable a given protein fold is. Although sequence capacity has been extensively studied using lattice models and theory, numerical estimates for real protein structures are currently lacking. In this work, we have quantitatively estimated the sequence capacity of 10 proteins with a variety of different structures using a statistical model based on residue-residue co-evolution to capture the variation of sequences from the same protein family. Remarkably, we find that even for the smallest protein folds, such as the WW domain, the number of foldable sequences is extremely large, exceeding the Avogadro constant. In agreement with earlier theoretical work, the calculated sequence capacity is positively correlated with the size of the protein, or better, the density of contacts. This allows the absolute sequence capacity of a given protein to be approximately predicted from its structure. On the other hand, the relative sequence capacity, i.e., normalized by the total number of possible sequences, is an extremely tiny number and is strongly anti-correlated with the protein length. Thus, although there may be more foldable sequences for larger proteins, it will be much harder to find them. Lastly, we have correlated the evolutionary age of proteins in the CATH database with their sequence capacity as predicted by our model. The results suggest a trade-off between the opposing requirements of high designability and the likelihood of a novel fold emerging by chance.