Charge asymmetry in the proteins of the outer membrane

Protein-lipid charge interactions control the folding of outer membrane proteins into asymmetric membranes

Biological membranes consist of two leaflets of phospholipid molecules that form a bilayer, each leaflet comprising a distinct lipid composition. This asymmetry is created and maintained in vivo by dedicated biochemical pathways, but difficulties in creating stable asymmetric membranes in vitro have restricted our understanding of how bilayer asymmetry modulates the folding, stability and function of membrane proteins. In this study, we used cyclodextrin-mediated lipid exchange to generate liposomes with asymmetric bilayers and characterize the stability and folding kinetics of two bacterial outer membrane proteins (OMPs), OmpA and BamA. We found that excess negative charge in the outer leaflet of a liposome impedes their insertion and folding, while excess negative charge in the inner leaflet accelerates their folding relative to symmetric liposomes with the same membrane composition. Using molecular dynamics, mutational analysis and bioinformatics, we identified a positively charged patch critical for folding and stability. These results rationalize the well-known 'positive-outside' rule of OMPs and suggest insights into the mechanisms that drive OMP folding and assembly in vitro and in vivo.

General features of transmembrane beta barrels from a large database

Large datasets contribute new insights to subjects formerly investigated by exemplars. We used coevolution data to create a large, high-quality database of transmembrane β-barrels (TMBB). By applying simple feature detection on generated evolutionary contact maps, our method (IsItABarrel) achieves 95.88% balanced accuracy when discriminating among protein classes. Moreover, comparison with IsItABarrel revealed a high rate of false positives in previous TMBB algorithms. In addition to being more accurate than previous datasets, our database (available online) contains 1,938,936 bacterial TMBB proteins from 38 phyla, respectively, 17 and 2.2 times larger than the previous sets TMBB-DB and OMPdb. We anticipate that due to its quality and size, the database will serve as a useful resource where high-quality TMBB sequence data are required. We found that TMBBs can be divided into 11 types, three of which have not been previously reported. We find tremendous variance in proteome percentage among TMBB-containing organisms with some using 6.79% of their proteome for TMBBs and others using as little as 0.27% of their proteome. The distribution of the lengths of the TMBBs is suggestive of previously hypothesized duplication events. In addition, we find that the C-terminal β-signal varies among different classes of bacteria though its consensus sequence is LGLGYRF. However, this β-signal is only characteristic of prototypical TMBBs. The ten non-prototypical barrel types have other C-terminal motifs, and it remains to be determined if these alternative motifs facilitate TMBB insertion or perform any other signaling function.

Asymmetric Lipid Vesicles: Techniques, Applications, and Future Perspectives as an Innovative Drug Delivery System

Novel lipid-based nanosystems have been of interest in improving conventional drug release methods. Liposomes are the most studied nanostructures, consisting of lipid bilayers ideal for drug delivery, thanks to their resemblance to the cell plasma membrane. Asymmetric liposomes are vesicles with different lipids in their inner and outer layers; because of this, they can be configured to be compatible with the therapeutic drug while achieving biocompatibility and stability. Throughout this review, topics such as the applications, advantages, and synthesis techniques of asymmetric liposomes will be discussed. Further, an in silico analysis by computational tools will be examined as a helpful tool for designing and understanding asymmetric liposome mechanisms in pharmaceutical applications. The dual-engineered design of asymmetric liposomes makes them an ideal alternative for transdermal drug delivery because of the improved protection of pharmaceuticals without lowering adsorption rates and system biocompatibility.

Drug‑resistant Acinetobacter baumannii: From molecular mechanisms to potential therapeutics (Review)

Bacterial drug resistance is increasingly becoming an important problem that needs to be solved urgently in modern clinical practices. Infection caused by Acinetobacter baumannii is a serious threat to the life and health of patients. The drug resistance rate of Acinetobacter baumannii strains is increasing, thus research on the drug resistance of Acinetobacter baumannii has also seen an increase. When patients are infected with drug-resistant Acinetobacter baumannii, the availability of suitable antibiotics commonly used in clinical practices is becoming increasingly limited and the prognosis of patients is worsening. Studying the molecular mechanism of the drug resistance of Acinetobacter baumannii is fundamental to solving the problem of drug-resistant Acinetobacter baumannii and potentially other 'super bacteria'. Drug resistance mechanisms primarily include enzymes, membrane proteins, efflux pumps and beneficial mutations. Research on the underlying mechanisms provides a theoretical basis for the use and development of antibiotics and the development of novel treatment methods.

Putative Riemerella anatipestifer Outer Membrane Protein H Affects Virulence

Riemerella anatipestifer causes serious contagious disease in ducks, geese, and other fowl. However, as a harmful pathogen causing significant economic losses in the poultry industry, R. anatipestifer is still poorly understood for its pathogenesis mechanisms. In a previous study, we developed an indirect ELISA method for detecting R. anatipestifer infection using B739_0832 protein, a putative outer membrane protein H (OmpH) that is conserved among different serotypes of R. anatipestifer. Although OmpH in some pathogenic bacteria, such as Pasteurella, has been reported as a virulence factor, it is still not clear whether B739_0832 protein contributes to the virulence of R. anatipestifer. In this study, we confirmed that B739_0832 protein in R. anatipestifer localizes to the outer membrane. We constructed a B739_0832 deletion mutant strain (ΔB739_0832) and assayed various effects from the deletion of B739_0832. ΔB739_0832 strain had a similar growth rate to wild-type R. anatipestifer CH-1. However, the survival rate of ducklings in 10 days after infection from ΔB739_0832 strain was 50%, whereas no ducklings survived from wild-type R. anatipestifer infection. Furthermore, the median lethal dose (LD₅₀) of the ΔB739_0832 strain was approximately 150 times higher than that of the wild-type strain. Pathology examinations on infected ducklings found that, at 36 h after infection, bacterial loads in blood, liver, and brain tissues from ΔB739_0832-infected ducklings were considerably lower than those from wild-type infected ducklings. These results demonstrate that the B739_0832 protein contributes to the virulence of R. anatipestifer CH-1.

Principles and Methods in Computational Membrane Protein Design

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

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

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

Membrane Barrels Are Taller, Fatter, Inside-Out Soluble Barrels

Up-and-down β-barrel topology exists in both the membrane and soluble environment. By comparing features of these structurally similar proteins, we can determine what features are particular to the environment rather than the fold. Here we compare structures of membrane β-barrels to soluble β-barrels and evaluate their relative size, shape, amino acid composition, hydrophobicity, and periodicity. We find that membrane β-barrels are generally larger than soluble β-barrels, with more strands per barrel and more amino acids per strand, making them wider and taller. We also find that membrane β-barrels are inside-out soluble β-barrels. The inward region of membrane β-barrels has similar hydrophobicity to the outward region of soluble β-barrels, and the outward region of membrane β-barrels has similar hydrophobicity to the inward region of the soluble β-barrels. Moreover, even though both types of β-barrel have been assumed to have strands with amino acids that alternate in direction and hydrophobicity, we find that the membrane β-barrels have more regular alternation than soluble β-barrels. These features give insight into how membrane barrels maintain their fold and function in the membrane.

The Outer Membrane Proteins OmpA, CarO, and OprD of Acinetobacter baumannii Confer a Two-Pronged Defense in Facilitating Its Success as a Potent Human Pathogen

Of all the ESKAPE pathogens, carbapenem-resistant and multidrug-resistant Acinetobacter baumannii is the leading cause of hospital-acquired and ventilator-associated pneumonia. A. baumannii infections are notoriously hard to eradicate due to its propensity to rapidly acquire multitude of resistance determinants and the virulence factor cornucopia elucidated by the bacterium that help it fend off a wide range of adverse conditions imposed upon by host and environment. One such weapon in the arsenal of A. baumannii is the outer membrane protein (OMP) compendium. OMPs in A. baumannii play distinctive roles in facilitating the bacterial acclimatization to antibiotic- and host-induced stresses, albeit following entirely different mechanisms. OMPs are major immunogenic proteins in bacteria conferring bacteria host-fitness advantages including immune evasion, stress tolerance, and resistance to antibiotics and antibacterials. In this review, we summarize the current knowledge of major A. baumannii OMPs and discuss their versatile role in antibiotic resistance and virulence. Specifically, we explore how OmpA, CarO, and OprD-like porins mediate antibiotic and amino acid shuttle and host virulence.

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.

Peptoid drug discovery and optimization via surface X-ray scattering

Synthetic polymers mimicking antimicrobial peptides have drawn considerable interest as potential therapeutics. N-substituted glycines, or peptoids, are recognized by their in vivo stability and ease of synthesis. Peptoids are thought to act primarily on the negatively charged lipids that are abundant in bacterial cell membranes. A mechanistic understanding of lipid-peptoid interaction at the molecular level will provide insights for rational design and optimization of peptoids. Here, we highlight recent studies that utilize synchrotron liquid surface X-ray scattering to characterize the underlying peptoid interactions with bacterial and eukaryotic membranes. Cellular membranes are highly complex, and difficult to characterize at the molecular level. Model systems including Langmuir monolayers, are used in these studies to reduce system complexity. The general workflow of these systems and the corresponding data analysis techniques are presented alongside recent findings. These studies investigate the role of peptoid physicochemical characteristics on membrane activity. Specifically, the roles of cationic charge, conformational constraint via macrocyclization, and hydrophobicity are shown to correlate their membrane interactions to biological activities in vitro. These structure-activity relationships have led to new insights into the mechanism of action by peptoid antimicrobials, and suggest optimization strategies for future therapeutics based on peptoids.

Evolutionary pathways of repeat protein topology in bacterial outer membrane proteins

Outer membrane proteins (OMPs) are the proteins in the surface of Gram-negative bacteria. These proteins have diverse functions but a single topology: the β-barrel. Sequence analysis has suggested that this common fold is a β-hairpin repeat protein, and that amplification of the β-hairpin has resulted in 8-26-stranded barrels. Using an integrated approach that combines sequence and structural analyses, we find events in which non-amplification diversification also increases barrel strand number. Our network-based analysis reveals strand-number-based evolutionary pathways, including one that progresses from a primordial 8-stranded barrel to 16-strands and further, to 18-strands. Among these pathways are mechanisms of strand number accretion without domain duplication, like a loop-to-hairpin transition. These mechanisms illustrate perpetuation of repeat protein topology without genetic duplication, likely induced by the hydrophobic membrane. Finally, we find that the evolutionary trace is particularly prominent in the C-terminal half of OMPs, implicating this region in the nucleation of OMP folding.

Identification of a novel post-insertion step in the assembly of a bacterial outer membrane protein

Although the barrel assembly machinery (Bam) complex has been shown to facilitate the insertion of β barrel proteins into the bacterial outer membrane (OM), the stage at which β barrels fold is unknown. Here, we describe insights into β barrel assembly that emerged from an analysis of a member of the autotransporter family of OM proteins (EspP) in Escherichia coli. EspP contains an extracellular 'passenger' domain that is translocated across the OM and then released from the covalently linked β barrel domain in an intra-barrel cleavage reaction. We found that the mutation of an unusual lipid-exposed lysine residue impairs a previously unidentified late folding step that follows both the membrane insertion of the β barrel domain and the secretion of the passenger domain but that precedes proteolytic maturation. Our results demonstrate that β barrel assembly can be completed at a post-insertion stage and raise the possibility that interactions with membrane lipids can promote folding in vivo. Furthermore, by showing that the passenger domain is secreted before the β barrel domain is fully assembled, our results also provide evidence against the long-standing hypothesis that autotransporters are autonomous protein secretion systems.

Efflux Pumps Represent Possible Evolutionary Convergence onto the β-Barrel Fold

There are around 100 varieties of outer membrane proteins in each Gram-negative bacteria. All of these proteins have the same fold-an up-down β-barrel. It has been suggested that all membrane β-barrels excluding lysins are homologous. Here we suggest that β-barrels of efflux pumps have converged on this fold as well. By grouping structurally solved outer membrane β-barrels (OMBBs) by sequence we find that the membrane environment may have led to convergent evolution of the barrel fold. Specifically, the lack of sequence linkage to other barrels coupled with distinctive structural differences, such as differences in strand tilt and barrel radius, suggest that the outer membrane factor of efflux pumps evolutionarily converged on the barrel. Rather than being related to other OMBBs, sequence and structural similarity in the periplasmic region of the outer membrane factor of efflux pumps suggests an evolutionary link to the periplasmic subunit of the same pump complex.

Tight Turns of Outer Membrane Proteins: An Analysis of Sequence, Structure, and Hydrogen Bonding

As a structural class, tight turns can control molecular recognition, enzymatic activity, and nucleation of folding. They have been extensively characterized in soluble proteins but have not been characterized in outer membrane proteins (OMPs), where they also support critical functions. We clustered the 4 to 6 residue tight turns of 110 OMPs to characterize the phi/psi angles, sequence, and hydrogen bonding of these structures. We find significant differences between reports of soluble protein tight turns and OMP tight turns. Since OMP strands are less twisted than soluble strands, they favor different turn structures types. Moreover, the membrane localization of OMPs yields different sequence hallmarks for their tight turns relative to soluble protein turns. We also characterize the differences in phi/psi angles, sequence, and hydrogen bonding between OMP extracellular loops and OMP periplasmic turns. As previously noted, the extracellular loops tend to be much longer than the periplasmic turns. We find that this difference in length is due to the broader distribution of lengths of the extracellular loops not a large difference in the median length. Extracellular loops also tend to have more charged residues as predicted by the charge-out rule. Finally, in all OMP tight turns, hydrogen bonding between the side chain and backbone 2 to 4 residues away from that side chain plays an important role. These bonds preferentially use an Asp, Asn, Ser, or Thr residue in a beta or pro phi/psi conformation. We anticipate that this study will be applicable to future design and structure prediction of OMPs.

Preparation and Properties of Asymmetric Synthetic Membranes Based on Lipid and Polymer Self-Assembly

Cell membrane asymmetry is a common structural feature of all biological cells. Researchers have tried for decades to better study its formation and its function in membrane-regulated phenomena. In particular, there has been increasing interest in developing synthetic asymmetric membrane models in the laboratory, with the aim of studying basic physical chemistry properties that may be correlated to a relevant biological function. The present article aims to summarize the main presented approaches to prepare asymmetric membranes, which are most often made from lipids, polymers, or a combination of both.

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

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

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

Motivation

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

Results

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

Availability and Implementation

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

Contact

: jliang@uic.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

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

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

Antibacterial properties of TiO2 modified with reduced graphene oxide

In this paper, the antibacterial activity of titanium dioxide modified with reduced graphene oxide (rGO) was presented. TiO₂/rGO photocatalysts were prepared by the hydrothermal method under elevated pressure at 180°C and heated at 100°C in Ar flow. The obtained photocatalysts were characterized by means of XRD, FTIR/DRS, UV-vis/DR, Raman spectroscopy and scanning electron microscopy (SEM). The carbon content was also examined. FTIR/DRS and Raman analysis confirmed the presence of rGO in the TiO₂ structure, suggesting a successful modification. The antimicrobial photoactivity of photocatalysts was conducted against E. coli under an artificial solar light. The results show that all TiO₂/rGO photocatalysts exhibited an antibacterial activity higher than unmodified TiO₂. The best result was found for sample with 1.5wt% additive of reduced graphene oxide. In this case, total inactivation of E. coli was noticed after 75min of irradiation. It was found that the presence of rGO in sample improves the antimicrobial activity.

Selective pressure for rapid membrane integration constrains the sequence of bacterial outer membrane proteins

Almost all bacterial outer membrane proteins (OMPs) contain a β barrel domain that serves as a membrane anchor, but the assembly and quality control of these proteins are poorly understood. Here, we show that the introduction of a single lipid-facing arginine residue near the middle of the β barrel of the Escherichia coli OMPs OmpLA and EspP creates an energy barrier that impedes membrane insertion. Although several unintegrated OmpLA mutants remained insertion-competent, they were slowly degraded by the periplasmic protease DegP. Two EspP mutants were also gradually degraded by DegP but were toxic because they first bound to the Bam complex, an essential heteroligomer that catalyzes the membrane insertion of OMPs. Interestingly, another EspP mutant likewise formed a prolonged, deleterious interaction with the Bam complex but was protected from degradation and eventually inserted into the membrane in a native conformation. The different types of interactions between the EspP mutants and the Bam complex that we observed may correspond to distinct stages in OMP assembly. Our results show that sequences that significantly delay assembly are disfavored not only because unintegrated OMPs are subjected to degradation, but also because OMPs that assemble slowly can form dominant-negative interactions with the Bam complex.

Role of the Native Outer-Membrane Environment on the Transporter BtuB

BtuB is a TonB-dependent transporter that permits the high-affinity binding and transport of cobalamin (CBL), or vitamin B₁₂, across the asymmetric outer membrane (OM) of Gram-negative bacteria. It has been shown that Ca²⁺ binding is necessary for high-affinity binding of CBL to BtuB, and earlier simulations suggested that calcium ions serve to stabilize key substrate-binding extracellular loops. However, those simulations did not account for the lipopolysaccharides in the OM. To illuminate the roles of both Ca²⁺ and lipopolysaccharides in protein functionality, we performed simulations of apo and Ca²⁺-loaded BtuB in symmetric and asymmetric bilayers. The simulations reveal that the oligosaccharides of LPS stabilize the extracellular loops to some degree, apparently obviating the need for Ca²⁺. However, it is shown that Ca²⁺ ions stabilize a key substrate-binding loop to an even greater degree, as well as reposition specific CBL-binding residues, bringing them closer to the organization found in the CBL-bound structure. These results indicate the importance of including realistic membrane models when simulating outer-membrane proteins.

Structural transformation of vesicles formed by a polystyrene-b-poly(acrylic acid)/polystyrene-b-poly(4-vinyl pyridine) mixture: from symmetric to asymmetric membranes

Asymmetric vesicles with different inner and outer corona compositions are applicable in microreactors, drug delivery, and biomimics because of their unique functions in membrane permeability and protein localization. In this study, we develop a novel approach to construct asymmetric vesicles and demonstrate the first structural transformation of polymeric vesicles from symmetric to asymmetric membranes. Experimental results and Monte Carlo simulation results clearly reveal that increased intercorona repulsion and enhanced hydrophobic chain mobility are essential to realize this transformation. Moreover, similar transformation processes are observed where either HCl or NaOH is added to change the intercorona interaction. This finding indicates that the observed structural transformation is dominated by physical interactions rather than chemical environment. The constructed asymmetric vesicles can be selectively decorated with gold nanoparticles on the outer corona. This study introduces a novel approach to prepare asymmetric vesicles and provides insights into the mechanism underlying the structural transformation of polymeric vesicles from symmetric to asymmetric membranes.

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

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

Computing structure-based lipid accessibility of membrane proteins with mp_lipid_acc in RosettaMP

BACKGROUND

Membrane proteins are underrepresented in structural databases, which has led to a lack of computational tools and the corresponding inappropriate use of tools designed for soluble proteins. For membrane proteins, lipid accessibility is an essential property. Although programs are available for sequence-based prediction of lipid accessibility and structure-based identification of solvent-accessible surface area, the latter does not distinguish between water accessible and lipid accessible residues in membrane proteins.

RESULTS

Here we present mp_lipid_acc, the first method to identify lipid accessible residues from the protein structure, implemented in the RosettaMP framework and available as a webserver. Our method uses protein structures transformed in membrane coordinates, for instance from PDBTM or OPM databases, and a defined membrane thickness to classify lipid accessibility of residues. mp_lipid_acc is applicable to both α-helical and β-barrel membrane proteins of diverse architectures with or without water-filled pores and uses a concave hull algorithm for surface-residue classification. We further provide a manually curated benchmark dataset that can be used for further method development.

CONCLUSIONS

We present a novel tool to classify lipid accessibility from the protein structure, which is applicable to proteins of diverse architectures and achieves prediction accuracies of 90% on a manually curated database. mp_lipid_acc is part of the Rosetta software suite, available at www.rosettacommons.org . The webserver is available at http://rosie.graylab.jhu.edu/mp_lipid_acc/submit and the benchmark dataset is available at http://tinyurl.com/mp-lipid-acc-dataset .

Predicting Beta Barrel Transmembrane Proteins Using HMMs

Transmembrane beta-barrels (TMBBs) constitute an important structural class of membrane proteins located in the outer membrane of gram-negative bacteria, and in the outer membrane of chloroplasts and mitochondria. They are involved in a wide variety of cellular functions and the prediction of their transmembrane topology, as well as their discrimination in newly sequenced genomes is of great importance as they are promising targets for antimicrobial drugs and vaccines. Several methods have been applied for the prediction of the transmembrane segments and the topology of beta barrel transmembrane proteins utilizing different algorithmic techniques. Hidden Markov Models (HMMs) have been efficiently used in the development of several computational methods used for this task. In this chapter we give a brief review of different available prediction methods for beta barrel transmembrane proteins pointing out sequence and structural features that should be incorporated in a prediction method. We then describe the procedure of the design and development of a Hidden Markov Model capable of predicting the transmembrane beta strands of TMBBs and discriminating them from globular proteins.

Outer membrane protein design

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

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

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

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

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

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

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

Living on the edge: Simulations of bacterial outer-membrane proteins

Gram-negative bacteria are distinguished in part by a second, outer membrane surrounding them. This membrane is distinct from others, possessing an outer leaflet composed not of typical phospholipids but rather large, highly charged molecules known as lipopolysaccharides. Therefore, modeling the structure and dynamics of proteins embedded in the outer membrane requires careful consideration of their native environment. In this review, we examine how simulations of such outer-membrane proteins have evolved over the last two decades, culminating most recently in detailed, highly accurate atomistic models of the outer membrane. We also draw attention to how the simulations have coupled with experiments to produce novel insights unattainable through a single approach. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

An Integrated Framework Advancing Membrane Protein Modeling and Design

Membrane proteins are critical functional molecules in the human body, constituting more than 30% of open reading frames in the human genome. Unfortunately, a myriad of difficulties in overexpression and reconstitution into membrane mimetics severely limit our ability to determine their structures. Computational tools are therefore instrumental to membrane protein structure prediction, consequently increasing our understanding of membrane protein function and their role in disease. Here, we describe a general framework facilitating membrane protein modeling and design that combines the scientific principles for membrane protein modeling with the flexible software architecture of Rosetta3. This new framework, called RosettaMP, provides a general membrane representation that interfaces with scoring, conformational sampling, and mutation routines that can be easily combined to create new protocols. To demonstrate the capabilities of this implementation, we developed four proof-of-concept applications for (1) prediction of free energy changes upon mutation; (2) high-resolution structural refinement; (3) protein-protein docking; and (4) assembly of symmetric protein complexes, all in the membrane environment. Preliminary data show that these algorithms can produce meaningful scores and structures. The data also suggest needed improvements to both sampling routines and score functions. Importantly, the applications collectively demonstrate the potential of combining the flexible nature of RosettaMP with the power of Rosetta algorithms to facilitate membrane protein modeling and design.

Computational redesign of the lipid-facing surface of the outer membrane protein OmpA

Advances in computational design methods have made possible extensive engineering of soluble proteins, but designed β-barrel membrane proteins await improvements in our understanding of the sequence determinants of folding and stability. A subset of the amino acid residues of membrane proteins interact with the cell membrane, and the design rules that govern this lipid-facing surface are poorly understood. We applied a residue-level depth potential for β-barrel membrane proteins to the complete redesign of the lipid-facing surface of Escherichia coli OmpA. Initial designs failed to fold correctly, but reversion of a small number of mutations indicated by backcross experiments yielded designs with substitutions to up to 60% of the surface that did support folding and membrane insertion.

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1-2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug-specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans-membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α-helical MPs as well as β-barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge-based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade.

Phase behavior under a noncentrosymmetric interaction: shifted-charge colloids investigated by Monte Carlo simulation

Using Monte Carlo simulations, we investigate the structural characteristics of an interacting hard-sphere system with shifted charge to elucidate the effect of the noncentrosymmetric interaction on its phase behavior. Two different phase transitions are identified for this model system. With increasing volume fraction, an abrupt liquid-to-crystal transition first occurs at a significantly lower volume fraction than in centrosymmetrically charged systems. This is due to the stronger effective interparticle repulsion caused by the additional charge anisotropy. Moreover, within the crystal state at higher volume fraction, the system further undergoes a continuous disorder-to-order transition with respect to charge orientation. Detailed analyses in this work disclose the nature of these transitions, and orientation fluctuation can cause noncentrosymmetric unit cells. The dependence of crystal formation and orientational ordering on temperature was also examined. These findings indicate that the noncentrosymmetric interaction in this work results in additional freedoms to fine-tune the phase diagram and increase the functionalities of materials. Moreover, these model studies are essential to advance the future understanding regarding the fundamental physiochemical properties of novel Janus colloidal particles and protein crystallization conditions.

E. coli outer membrane and interactions with OmpLA

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.