Higher-order interhelical spatial interactions in membrane proteins

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.

COMSAT: Residue contact prediction of transmembrane proteins based on support vector machines and mixed integer linear programming

In this article, we present COMSAT, a hybrid framework for residue contact prediction of transmembrane (TM) proteins, integrating a support vector machine (SVM) method and a mixed integer linear programming (MILP) method. COMSAT consists of two modules: COMSAT_SVM which is trained mainly on position-specific scoring matrix features, and COMSAT_MILP which is an ab initio method based on optimization models. Contacts predicted by the SVM model are ranked by SVM confidence scores, and a threshold is trained to improve the reliability of the predicted contacts. For TM proteins with no contacts above the threshold, COMSAT_MILP is used. The proposed hybrid contact prediction scheme was tested on two independent TM protein sets based on the contact definition of 14 Å between Cα-Cα atoms. First, using a rigorous leave-one-protein-out cross validation on the training set of 90 TM proteins, an accuracy of 66.8%, a coverage of 12.3%, a specificity of 99.3% and a Matthews' correlation coefficient (MCC) of 0.184 were obtained for residue pairs that are at least six amino acids apart. Second, when tested on a test set of 87 TM proteins, the proposed method showed a prediction accuracy of 64.5%, a coverage of 5.3%, a specificity of 99.4% and a MCC of 0.106. COMSAT shows satisfactory results when compared with 12 other state-of-the-art predictors, and is more robust in terms of prediction accuracy as the length and complexity of TM protein increase. COMSAT is freely accessible at http://hpcc.siat.ac.cn/COMSAT/.

Dimers of G-protein coupled receptors as versatile storage and response units

The status and use of transmembrane, extracellular and intracellular domains in oligomerization of heptahelical G-protein coupled receptors (GPCRs) are reviewed and for transmembrane assemblies also supplemented by new experimental evidence. The transmembrane-linked GPCR oligomers typically have as the minimal unit an asymmetric ~180 kDa pentamer consisting of receptor homodimer or heterodimer and a G-protein αβγ subunit heterotrimer. With neuropeptide Y (NPY) receptors, this assembly is converted to ~90 kDa receptor monomer-Gα complex by receptor and Gα agonists, and dimers/heteropentamers are depleted by neutralization of Gαi subunits by pertussis toxin. Employing gradient centrifugation, quantification and other characterization of GPCR dimers at the level of physically isolated and identified heteropentamers is feasible with labeled agonists that do not dissociate upon solubilization. This is demonstrated with three neuropeptide Y (NPY) receptors and could apply to many receptors that use large peptidic agonists.

Trimeric transmembrane domain interactions in paramyxovirus fusion proteins: roles in protein folding, stability, and function

Paramyxovirus fusion (F) proteins promote membrane fusion between the viral envelope and host cell membranes, a critical early step in viral infection. Although mutational analyses have indicated that transmembrane (TM) domain residues can affect folding or function of viral fusion proteins, direct analysis of TM-TM interactions has proved challenging. To directly assess TM interactions, the oligomeric state of purified chimeric proteins containing the Staphylococcal nuclease (SN) protein linked to the TM segments from three paramyxovirus F proteins was analyzed by sedimentation equilibrium analysis in detergent and buffer conditions that allowed density matching. A monomer-trimer equilibrium best fit was found for all three SN-TM constructs tested, and similar fits were obtained with peptides corresponding to just the TM region of two different paramyxovirus F proteins. These findings demonstrate for the first time that class I viral fusion protein TM domains can self-associate as trimeric complexes in the absence of the rest of the protein. Glycine residues have been implicated in TM helix interactions, so the effect of mutations at Hendra F Gly-508 was assessed in the context of the whole F protein. Mutations G508I or G508L resulted in decreased cell surface expression of the fusogenic form, consistent with decreased stability of the prefusion form of the protein. Sedimentation equilibrium analysis of TM domains containing these mutations gave higher relative association constants, suggesting altered TM-TM interactions. Overall, these results suggest that trimeric TM interactions are important driving forces for protein folding, stability and membrane fusion promotion.

Re-introduction of transmembrane serine residues reduce the minimum pore diameter of channelrhodopsin-2

Channelrhodopsin-2 (ChR2) is a microbial-type rhodopsin found in the green algae Chlamydomonas reinhardtii. Under physiological conditions, ChR2 is an inwardly rectifying cation channel that permeates a wide range of mono- and divalent cations. Although this protein shares a high sequence homology with other microbial-type rhodopsins, which are ion pumps, ChR2 is an ion channel. A sequence alignment of ChR2 with bacteriorhodopsin, a proton pump, reveals that ChR2 lacks specific motifs and residues, such as serine and threonine, known to contribute to non-covalent interactions within transmembrane domains. We hypothesized that reintroduction of the eight transmembrane serine residues present in bacteriorhodopsin, but not in ChR2, will restrict the conformational flexibility and reduce the pore diameter of ChR2. In this work, eight single serine mutations were created at homologous positions in ChR2. Additionally, an endogenous transmembrane serine was replaced with alanine. We measured kinetics, changes in reversal potential, and permeability ratios in different alkali metal solutions using two-electrode voltage clamp. Applying excluded volume theory, we calculated the minimum pore diameter of ChR2 constructs. An analysis of the results from our experiments show that reintroducing serine residues into the transmembrane domain of ChR2 can restrict the minimum pore diameter through inter- and intrahelical hydrogen bonds while the removal of a transmembrane serine results in a larger pore diameter. Therefore, multiple positions along the intracellular side of the transmembrane domains contribute to the cation permeability of ChR2.

An amino acid packing code for α-helical structure and protein design

This work demonstrates that all packing in α-helices can be simplified to repetitive patterns of a single motif: the knob-socket. Using the precision of Voronoi Polyhedra/Delauney Tessellations to identify contacts, the knob-socket is a four-residue tetrahedral motif: a knob residue on one α-helix packs into the three-residue socket on another α-helix. The principle of the knob-socket model relates the packing between levels of protein structure: the intra-helical packing arrangements within secondary structure that permit inter-helix tertiary packing interactions. Within an α-helix, the three-residue sockets arrange residues into a uniform packing lattice. Inter-helix packing results from a definable pattern of interdigitated knob-socket motifs between two α-helices. Furthermore, the knob-socket model classifies three types of sockets: (1) free, favoring only intra-helical packing; (2) filled, favoring inter-helical interactions; and (3) non, disfavoring α-helical structure. The amino acid propensities in these three socket classes essentially represent an amino acid code for structure in α-helical packing. Using this code, we used a novel yet straightforward approach for the design of α-helical structure to validate the knob-socket model. Unique sequences for three peptides were created to produce a predicted amount of α-helical structure: mostly helical, some helical, and no helix. These three peptides were synthesized, and helical content was assessed using CD spectroscopy. The measured α-helicity of each peptide was consistent with the expected predictions. These results and analysis demonstrate that the knob-socket motif functions as the basic unit of packing and presents an intuitive tool to decipher the rules governing packing in protein structure.

Bacterial mechanosensitive channels--MscS: evolution's solution to creating sensitivity in function

The discovery of mechanosensing channels has changed our understanding of bacterial physiology. The mechanosensitive channel of small conductance (MscS) is perhaps the most intensively studied of these channels. MscS has at least two states: closed, which does not allow solutes to exit the cytoplasm, and open, which allows rapid efflux of solvent and solutes. The ability to appropriately open or close the channel (gating) is critical to bacterial survival. We briefly review the science that led to the isolation and identification of MscS. We concentrate on the structure-function relationship of the channel, in particular the structural and biochemical approaches to understanding channel gating. We highlight the troubling discrepancies between the various models developed to understand MscS gating.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Enhanced Inter-helical Residue Contact Prediction in Transmembrane Proteins

In this paper, based on a recent work by McAllister and Floudas who developed a mathematical optimization model to predict the contacts in transmembrane alpha-helical proteins from a limited protein data set [1], we have enhanced this method by 1) building a more comprehensive data set for transmembrane alpha-helical proteins and this enhanced data set is then used to construct the probability sets, MIN-1N and MIN-2N, for residue contact prediction, 2) enhancing the mathematical model via modifications of several important physical constraints and 3) applying a new blind contact prediction scheme on different protein sets proposed from analyzing the contact prediction on 65 proteins from Fuchs et al. [2]. The blind contact prediction scheme has been tested on two different membrane protein sets. Firstly it is applied to five carefully selected proteins from the training set. The contact prediction of these five proteins uses probability sets built by excluding the target protein from the training set, and an average accuracy of 56% was obtained. Secondly, it is applied to six independent membrane proteins with complicated topologies, and the prediction accuracies are 73% for 2ZY9A, 21% for 3KCUA, 46% for 2W1PA, 64% for 3CN5A, 77% for 3IXZA and 83% for 3K3FA. The average prediction accuracy for the six proteins is 60.7%. The proposed approach is also compared with a support vector machine method (TMhit [3]) and it is shown that it exhibits better prediction accuracy.

Sensing bilayer tension: bacterial mechanosensitive channels and their gating mechanisms

Mechanosensitive channels sense and respond to changes in bilayer tension. In many respects, this is a unique property: the changes in membrane tension gate the channel, leading to the transient formation of open non-selective pores. Pore diameter is also high for the bacterial channels studied, MscS and MscL. Consequently, in cells, gating has severe consequences for energetics and homoeostasis, since membrane depolarization and modification of cytoplasmic ionic composition is an immediate consequence. Protection against disruption of cellular integrity, which is the function of the major channels, provides a strong evolutionary rationale for possession of such disruptive channels. The elegant crystal structures for these channels has opened the way to detailed investigations that combine molecular genetics with electrophysiology and studies of cellular behaviour. In the present article, the focus is primarily on the structure of MscS, the small mechanosensitive channel. The description of the structure is accompanied by discussion of the major sites of channel–lipid interaction and reasoned, but limited, speculation on the potential mechanisms of tension sensing leading to gating.

Driving forces for transmembrane alpha-helix oligomerization

We present what we believe to be a novel statistical contact potential based on solved structures of transmembrane (TM) alpha-helical bundles, and we use this contact potential to investigate the amino acid likelihood of stabilizing helix-helix interfaces. To increase statistical significance, we have reduced the full contact energy matrix to a four-flavor alphabet of amino acids, automatically determined by our methodology, in which we find that polarity is a more dominant factor of group identity than is size, with charged or polar groups most often occupying the same face, whereas polar/apolar residue pairs tend to occupy opposite faces. We found that the most polar residues strongly influence interhelical contact formation, although they occur rarely in TM helical bundles. Two-body contact energies in the reduced letter code are capable of determining native structure from a large decoy set for a majority of test TM proteins, at the same time illustrating that certain higher-order sequence correlations are necessary for more accurate structure predictions.

Structure elucidation of dimeric transmembrane domains of bitopic proteins

The interaction between transmembrane helices is of great interest because it directly determines biological activity of a membrane protein. Either destroying or enhancing such interactions can result in many diseases related to dysfunction of different tissues in human body. One much studied form of membrane proteins known as bitopic protein is a dimer containing two membrane-spanning helices associating laterally. Establishing structure-function relationship as well as rational design of new types of drugs targeting membrane proteins requires precise structural information about this class of objects. At present time, to investigate spatial structure and internal dynamics of such transmembrane helical dimers, several strategies were developed based mainly on a combination of NMR spectroscopy, optical spectroscopy, protein engineering and molecular modeling. These approaches were successfully applied to homo- and heterodimeric transmembrane fragments of several bitopic proteins, which play important roles in normal and in pathological conditions of human organism.

Helix-helix interaction patterns in membrane proteins

Membrane-spanning α-helices represent major sites of protein-protein interaction in membrane protein oligomerization and folding. As such, these interactions may be of exquisite specificity. Specificity often rests on a complex interplay of different types of residues forming the helix-helix interfaces via dense packing and different non-covalent forces, including van der Waal’s forces, hydrogen bonding, charge-charge interactions, and aromatic interactions. These interfaces often contain complex residue motifs where the contribution of constituent amino acids depends on the context of the surrounding sequence. Moreover, transmembrane helix-helix interactions are increasingly recognized as being dynamic and dependent on the functional state of a given protein.

End-to-end and end-to-middle interhelical interactions: new classes of interacting helix pairs in protein structures

Helix-helix interactions are important for the structure, stability and function of alpha-helical proteins. Helices that either cross in the middle or show extensive contacts between each other, such as coiled coils, have been investigated in previous studies. Interactions between two helices can also occur only at the terminal regions or between the terminal region of one helix and the middle region of another helix. Examples of such helix pairs are found in aquaporin, H(+)/Cl(-) transporter and Bcl-2 proteins. The frequency of the occurrence of such ;end-to-end' (EE) and ;end-to-middle' (EM) helix pairs in protein structures is not known. Questions regarding the residue preferences in the interface and the mode of interhelical interactions in such helix pairs also remain unanswered. In this study, high-resolution structures of all-alpha proteins from the PDB have been systematically analyzed and the helix pairs that interact only in EE or EM fashion have been extracted. EE and EM helix pairs have been categorized into five classes (N-N, N-C, C-C, N-MID and C-MID) depending on the region of interaction. Nearly 13% of 5725 helix pairs belonged to one of the five classes. Analysis of single-residue propensities indicated that hydrophobic and polar residues prefer to occur in the C-terminal and N-terminal regions, respectively. Hydrophobic C-terminal interacting residues and polar N-terminal interacting residues are also highly conserved. A strong correlation exists between some of the residue properties (surface area/volume and length of side chains) and their preferences for occurring in the interface of EE and EM helix pairs. In contrast to interacting non-EE/EM helix pairs, helices in EE and EM pairs are farther apart. In these helix pairs, residues with large surface area/volume and longer side chains are preferred in the interfacial region.

Complex patterns of histidine, hydroxylated amino acids and the GxxxG motif mediate high-affinity transmembrane domain interactions

Specific interactions of transmembrane helices play a pivotal role in the folding and oligomerization of integral membrane proteins. The helix-helix interfaces frequently depend on specific amino acid patterns. In this study, a heptad repeat pattern was randomized with all naturally occurring amino acids to uncover novel sequence motifs promoting transmembrane domain interactions. Self-interacting transmembrane domains were selected from the resulting combinatorial library by means of the ToxR/POSSYCCAT system. A comparison of the amino acid composition of high-and low-affinity sequences revealed that high-affinity transmembrane domains exhibit position-specific enrichment of histidine. Further, sequences containing His preferentially display Gly, Ser, and/or Thr residues at flanking positions and frequently contain a C-terminal GxxxG motif. Mutational analysis of selected sequences confirmed the importance of these residues in homotypic interaction. Probing heterotypic interaction indicated that His interacts in trans with hydroxylated residues. Reconstruction of minimal interaction motifs within the context of an oligo-Leu sequence confirmed that His is part of a hydrogen bonded cluster that is brought into register by the GxxxG motif. Notably, a similar motif contributes to self-interaction of the BNIP3 transmembrane domain.

Alpha-helical topology prediction and generation of distance restraints in membrane proteins

The field of protein structure prediction has seen significant advances in recent years. Researchers have followed a multitude of approaches, including methods based on comparative modeling, fold recognition and threading, and first-principles techniques. It is noteworthy that the structure prediction of membrane proteins is comparatively less studied by researchers in the field. A membrane protein is characterized by a protein structure that extends into or through the lipid-lipid bilayer of a cell. The structure is influenced by the combination of the hydrophobic bilayer region, the direct interaction with the bilayer, and the aqueous external environment. Due to the difficulty in obtaining reliable experimental structures, accurate computational prediction of membrane proteins is of paramount importance. An optimization model has been developed to predict the interhelical interactions in alpha-helical membrane proteins. A database of alpha-helical membrane proteins of known structure and limited sequence identity can be constructed to develop interaction probabilities. By then maximizing the occurrence of highly probable pairwise or three-residue interactions, realistic contacts can be predicted by imposing a number of geometrical constraints. The development of these low distance contacts can provide additional distance restraints for first principles-based approaches to the tertiary structure prediction problem. The proposed approach is shown to successfully predict interhelical contacts in several membrane protein systems, including bovine rhodopsin and the recently released human beta2 adrenergic receptor protein structure.

Surface recognition elements of membrane protein oligomerization

Although certain membrane proteins are functional as monomeric polypeptides, others must assemble into oligomers to carry out their biological roles. High-resolution membrane protein structures provide a valuable resource for examining the sequence features that facilitate-or preclude-assembly of membrane protein monomers into multimeric structures. Here we have utilized a data set of 28 high-resolution alpha-helical membrane protein structures comprising 32 nonredundant polypeptides to address this issue. The lipid-exposed surfaces of membrane proteins that have reached their fully assembled and functional biological units have been compared with those of the individual subunits that build quaternary structures. Though the overall amino acid composition of each set of surfaces is similar, a key distinction-the distribution of small-xxx-small motifs-delineates subunits from membrane proteins that have reached a functioning oligomeric state. Quaternary structure formation may therefore be dictated by small-xxx-small motifs that are not satisfied by intrachain contacts.

Phenylalanine promotes interaction of transmembrane domains via GxxxG motifs

Interactions of transmembrane helices play a crucial role in the folding and oligomerisation of integral membrane proteins. In order to uncover novel sequence motifs mediating these interactions, we randomised one face of a transmembrane helix with a set of non-polar or moderately polar amino acids. Those sequences capable of self-interaction upon integration into bacterial inner membranes were selected by means of the ToxR/POSSYCCAT system. A comparison between low/medium-affinity and high-affinity sequences reveals that high-affinity sequences are strongly enriched in phenylalanine residues that are frequently observed at the -3 position of GxxxG motifs, thus yielding FxxGxxxG motifs. Mutation of Phe or GxxxG in selected sequences significantly reduces self-interaction of the transmembrane domains without affecting their efficiency of membrane integration. Conversely, grafting FxxGxxxG onto unrelated transmembrane domains strongly enhances their interaction. Further, we find that FxxGxxxG is significantly over-represented in transmembrane domains of bitopic membrane proteins. The same motif contributes to self-interaction of the vesicular stomatitis virus G protein transmembrane domain. We conclude that Phe stabilises membrane-spanning GxxxG motifs. This is one example of how the role of certain side-chains in helix-helix interfaces is modulated by sequence context.

Empirical potential function for simplified protein models: combining contact and local sequence-structure descriptors

An effective potential function is critical for protein structure prediction and folding simulation. Simplified protein models such as those requiring only Calpha or backbone atoms are attractive because they enable efficient search of the conformational space. We show residue-specific reduced discrete-state models can represent the backbone conformations of proteins with small RMSD values. However, no potential functions exist that are designed for such simplified protein models. In this study, we develop optimal potential functions by combining contact interaction descriptors and local sequence-structure descriptors. The form of the potential function is a weighted linear sum of all descriptors, and the optimal weight coefficients are obtained through optimization using both native and decoy structures. The performance of the potential function in a test of discriminating native protein structures from decoys is evaluated using several benchmark decoy sets. Our potential function requiring only backbone atoms or Calpha atoms have comparable or better performance than several residue-based potential functions that require additional coordinates of side-chain centers or coordinates of all side-chain atoms. By reducing the residue alphabets down to size 10 for contact descriptors, the performance of the potential function can be further improved. Our results also suggest that local sequence-structure correlation may play important role in reducing the entropic cost of protein folding.

Prediction of transmembrane helix orientation in polytopic membrane proteins

Background

Membrane proteins compose up to 30% of coding sequences within genomes. However, their structure determination is lagging behind compared with soluble proteins due to the experimental difficulties. Therefore, it is important to develop reliable computational methods to predict structures of membrane proteins.

Results

We present a method for prediction of the TM helix orientation, which is an essential step in ab initio modeling of membrane proteins. Our method is based on a canonical model of the heptad repeat originally developed for coiled coils. We identify the helical surface patches that interface with lipid molecules at an accuracy of about 88% from the sequence information alone, using an empirical scoring function LIPS (LIPid-facing Surface), which combines lipophilicity and conservation of residues in the helix. We test and discuss results of prediction of helix-lipid interfaces on 162 transmembrane helices from 18 polytopic membrane proteins and present predicted orientations of TM helices in TRPV1 channel. We also apply our method to two structures of homologous cytochrome b₆f complexes and find discrepancy in the assignment of TM helices from subunits PetG, PetN and PetL. The results of LIPS calculations and analysis of packing and H-bonding interactions support the helix assignment found in the cytochrome b₆f structure from green alga but not the assignment of TM helices in the cyanobacterium b₆f structure.

Conclusion

LIPS calculations can be used for the prediction of helix orientation in ab initio modeling of polytopic membrane proteins. We also show with the example of two cytochrome b₆f structures that our method can identify questionable helix assignments in membrane proteins. The LIPS server is available online at .

Structural Role of (Bacterio)chlorophyll Ligated in the Energetically Unfavorable β-Position

Chlorophyll is attached to apoprotein in diastereotopically distinct ways, by beta- and alpha-ligation. Both the beta- and alpha-ligated chlorophylls of photosystem I are shown to have ample contacts to apoprotein within their proteinaceous binding sites, in particular, at C-13 of the isocyclic ring. The H-bonding patterns for the C-13(1) oxo groups, however, are clearly distinct for the beta-ligated and alpha-ligated chlorophylls. The beta-ligated chlorophylls frequently employ their C-13(1) oxo in H-bonds to neighboring helices and subunits. In contrast, the C-13(1) oxo of alpha-ligated chlorophylls are significantly less involved in H-bonding interactions, particularly to neighboring helices. Remarkably, in the peripheral antenna, light harvesting complex (LH2) from Rhodobacter sphaeroides, a single mutation in the alpha-subunit, introduced to eliminate H-bonding to the beta-bacteriochlorophyll-B850, which is ligated in the "beta-position," results in significant thermal destabilization of the LH2 in the membrane. In addition, in comparison with wild type LH2, the expression level of the LH2 lacking this H-bond is significantly reduced. These findings show that H-bonding to the C-13(1) keto group ofbeta-ligated (bacterio)-chlorophyll is a key structural motif and significantly contributes to the stability of bacteriochlorophyll proteins in the native membrane. Our analysis of photosystem I and II suggests that this hitherto unrecognized motif involving H-bonding to beta-ligated chlorophylls may be equally critical for the stable assembly of the inner core antenna of these multicomponent chlorophyll proteins.

Interstrand Pairing Patterns in β-Barrel Membrane Proteins: The Positive-outside Rule, Aromatic Rescue, and Strand Registration Prediction

beta-Barrel membrane proteins are found in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. Little is known about how residues in membrane beta-barrels interact preferentially with other residues on adjacent strands. We have developed probabilistic models to quantify propensities of residues for different spatial locations and for interstrand pairwise contact interactions involving strong H-bonds, side-chain interactions, and weak H-bonds. Using the reference state of exhaustive permutation of residues within the same beta-strand, the propensity values and p-values measuring statistical significance are calculated exactly by analytical formulae we have developed. Our findings show that there are characteristic preferences of residues for different membrane locations. Contrary to the "positive-inside" rule for helical membrane proteins, beta-barrel membrane proteins follow a significant albeit weaker "positive-outside" rule, in that the basic residues Arg and Lys are disproportionately favored in the extracellular cap region and disfavored in the periplasmic cap region. We find that different residue pairs prefer strong backbone H-bonded interstrand pairings (e.g. Gly-aromatic) or non-H-bonded pairings (e.g. aromatic-aromatic). In addition, we find that Tyr and Phe participate in aromatic rescue by shielding Gly from polar environments. We also show that these propensities can be used to predict the registration of strand pairs, an important task for the structure prediction of beta-barrel membrane proteins. Our accuracy of 44% is considerably better than random (7%). It also significantly outperforms a comparable registration prediction for soluble beta-sheets under similar conditions. Our results imply several experiments that can help to elucidate the mechanisms of in vitro and in vivo folding of beta-barrel membrane proteins. The propensity scales developed in this study will also be useful for computational structure prediction and for folding simulations.

Voronoi and Voronoi-related tessellations in studies of protein structure and interaction

The three-dimensional structure of a protein can be modeled by a set of polyhedra drawn around its atoms or residues. The tessellation invented by Voronoi in 1908, and other tessellations of space derived from it, provide versatile representations of three-dimensional structures. In recent years, they have been used to investigate a series of issues relating to proteins: atom and residue volumes, packing, folding, interactions and binding.

Geometric cooperativity and anticooperativity of three-body interactions in native proteins

Characterizing multibody interactions of hydrophobic, polar, and ionizable residues in protein is important for understanding the stability of protein structures. We introduce a geometric model for quantifying 3-body interactions in native proteins. With this model, empirical propensity values for many types of 3-body interactions can be reliably estimated from a database of native protein structures, despite the overwhelming presence of pairwise contacts. In addition, we define a nonadditive coefficient that characterizes cooperativity and anticooperativity of residue interactions in native proteins by measuring the deviation of 3-body interactions from 3 independent pairwise interactions. It compares the 3-body propensity value from what would be expected if only pairwise interactions were considered, and highlights the distinction of propensity and cooperativity of 3-body interaction. Based on the geometric model, and what can be inferred from statistical analysis of such a model, we find that hydrophobic interactions and hydrogen-bonding interactions make nonadditive contributions to protein stability, but the nonadditive nature depends on whether such interactions are located in the protein interior or on the protein surface. When located in the interior, many hydrophobic interactions such as those involving alkyl residues are anticooperative. Salt-bridge and regular hydrogen-bonding interactions, such as those involving ionizable residues and polar residues, are cooperative. When located on the protein surface, these salt-bridge and regular hydrogen-bonding interactions are anticooperative, and hydrophobic interactions involving alkyl residues become cooperative. We show with examples that incorporating 3-body interactions improves discrimination of protein native structures against decoy conformations. In addition, analysis of cooperative 3-body interaction may reveal spatial motifs that can suggest specific protein functions.

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.

Optimal bundling of transmembrane helices using sparse distance constraints

We present a two-step approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only sparse distance constraints, such as those derived from chemical cross-linking, dipolar EPR and FRET experiments. In Step 1, using an algorithm, we developed, the conformational space of membrane protein folds matching a set of distance constraints is explored to provide initial structures for local conformational searches. In Step 2, these structures refined against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. We begin by describing the statistical analysis of the solved membrane protein structures from which the theoretical portion of the penalty function was derived. We then describe the penalty function, and, using a set of six test cases, demonstrate that it is capable of distinguishing helical bundles that are close to the native bundle from those that are far from the native bundle. Finally, using a set of only 27 distance constraints extracted from the literature, we show that our method successfully recovers the structure of dark-adapted rhodopsin to within 3.2 A of the crystal structure.

Hydrogen bonding in a model bacteriochlorophyll-binding site drives assembly of light harvesting complex

In this study, the contribution of intramembrane hydrogen bonding at the interface between polypeptide and cofactor is explored in the native lipid environment by use of model bacteriochlorophyll proteins. In the peripheral antenna complex, LH2, large portions of the transmembrane helices, which make up the dimeric bacteriochlorophyll-binding site, are replaced by simplified, alternating alanine-leucine stretches. Replacement of either one of the two helices with the helices containing the model sequence at a time results in the assembly of complexes with nearly native light harvesting properties. In contrast, replacement of both helices results in the loss of antenna complexes from the membrane. The assembly of such doubly modified complexes is restored by a single intramembrane serine residue at position -4 relative to the liganding histidine of the alpha-subunit. In situ analysis of the spectral properties in a series of site-directed mutants reveals a critical dependence of the model complex assembly on the side chain of the residue at this position in the helix. A hydrogen bond between the hydroxy group of the serine and the 13(1) keto group of one of the central bacteriochlorophylls of the complexes is identified by Raman spectroscopy in the model antenna complex containing one of the alanine-leucine helices. The additional OH group of the serine residue, which participates in hydrogen bonding, increases the thermal stability of the model complexes in the native membrane. Intramembrane hydrogen bonding is thus shown to be a key factor for the binding of bacteriochlorophyll and assembly of this model cofactor-polypeptide site.

Sequence motifs, polar interactions and conformational changes in helical membrane proteins

The alpha helices of transmembrane proteins interact to form higher order structures. These interactions are frequently mediated by packing motifs (such as GxxxG) and polar residues. Recent structural data have revealed that small sidechains are able to both stabilize helical membrane proteins and allow conformational changes in the structure. The strong interactions involving polar sidechains often contribute to protein misfolding or malfunction.