Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Challenges of Protein-Protein Docking of the Membrane Proteins

Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.

Machine learning in computational modelling of membrane protein sequences and structures: From methodologies to applications

Membrane proteins mediate a wide spectrum of biological processes, such as signal transduction and cell communication. Due to the arduous and costly nature inherent to the experimental process, membrane proteins have long been devoid of well-resolved atomic-level tertiary structures and, consequently, the understanding of their functional roles underlying a multitude of life activities has been hampered. Currently, computational tools dedicated to furthering the structure-function understanding are primarily focused on utilizing intelligent algorithms to address a variety of site-wise prediction problems (e.g., topology and interaction sites), but are scattered across different computing sources. Moreover, the recent advent of deep learning techniques has immensely expedited the development of computational tools for membrane protein-related prediction problems. Given the growing number of applications optimized particularly by manifold deep neural networks, we herein provide a review on the current status of computational strategies mainly in membrane protein type classification, topology identification, interaction site detection, and pathogenic effect prediction. Meanwhile, we provide an overview of how the entire prediction process proceeds, including database collection, data pre-processing, feature extraction, and method selection. This review is expected to be useful for developing more extendable computational tools specific to membrane proteins.

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.

Research progress of reduced amino acid alphabets in protein analysis and prediction

•

A comprehensive summary of the literature on the reduced amino acid alphabets.

•

A systematic review of the development history of reduced amino acid alphabets.

•

Rich application cases of amino acid reduction alphabets are described in the article.

•

A detailed analysis of the properties and uses of the reduced amino acid alphabets.

Proteins are the executors of cellular physiological activities, and accurate structural and function elucidation are crucial for the refined mapping of proteins. As a feature engineering method, the reduction of amino acid composition is not only an important method for protein structure and function analysis, but also opens a broad horizon for the complex field of machine learning. Representing sequences with fewer amino acid types greatly reduces the complexity and noise of traditional feature engineering in dimension, and provides more interpretable predictive models for machine learning to capture key features. In this paper, we systematically reviewed the strategy and method studies of the reduced amino acid (RAA) alphabets, and summarized its main research in protein sequence alignment, functional classification, and prediction of structural properties, respectively. In the end, we gave a comprehensive analysis of 672 RAA alphabets from 74 reduction methods.

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.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

Improved methanol tolerance of Rhizomucor miehei lipase based on N‑glycosylation within the α-helix region and its application in biodiesel production

BACKGROUND

Liquid lipases are widely used to convert oil into biodiesel. Methanol-resistant lipases with high catalytic activity are the first choice for practical production. Rhizomucor miehei lipase (RML) is a single-chain α/β-type protein that is widely used in biodiesel preparation. Improving the catalytic activity and methanol tolerance of RML is necessary to realise the industrial production of biodiesel.

RESULTS

In this study, a semi-rational design method was used to optimise the catalytic activity and methanol tolerance of ProRML. After N-glycosylation modification of the α-helix of the mature peptide in ProRML, the resulting mutants N218, N93, N115, N260, and N183 increased enzyme activity by 66.81, 13.54, 10.33, 3.69, and 2.39 times than that of WT, respectively. The residual activities of N218 and N260 were 88.78% and 86.08% after incubation in 50% methanol for 2.5 h, respectively. In addition, the biodiesel yield of all mutants was improved when methanol was added once and reacted for 24 h with colza oil as the raw material. N260 and N218 increased the biodiesel yield from 9.49% to 88.75% and 90.46%, respectively.

CONCLUSIONS

These results indicate that optimising N-glycosylation modification in the α-helix structure is an effective strategy for improving the performance of ProRML. This study provides an effective approach to improve the design of the enzyme and the properties of lipase mutants, thereby rendering them suitable for industrial biomass conversion.

Structural cavities are critical to balancing stability and activity of a membrane-integral enzyme

Packing interaction is a critical driving force in the folding of helical membrane proteins. Despite the importance, packing defects (i.e., cavities including voids, pockets, and pores) are prevalent in membrane-integral enzymes, channels, transporters, and receptors, playing essential roles in function. Then, a question arises regarding how the two competing requirements, packing for stability vs. cavities for function, are reconciled in membrane protein structures. Here, using the intramembrane protease GlpG of Escherichia coli as a model and cavity-filling mutation as a probe, we tested the impacts of native cavities on the thermodynamic stability and function of a membrane protein. We find several stabilizing mutations which induce substantial activity reduction without distorting the active site. Notably, these mutations are all mapped onto the regions of conformational flexibility and functional importance, indicating that the cavities facilitate functional movement of GlpG while compromising the stability. Experiment and molecular dynamics simulation suggest that the stabilization is induced by the coupling between enhanced protein packing and weakly unfavorable lipid desolvation, or solely by favorable lipid solvation on the cavities. Our result suggests that, stabilized by the relatively weak interactions with lipids, cavities are accommodated in membrane proteins without severe energetic cost, which, in turn, serve as a platform to fine-tune the balance between stability and flexibility for optimal activity.

Bacterial lyso-form lipoproteins are synthesized via an intramolecular acyl chain migration

All bacterial lipoproteins share a variably acylated N-terminal cysteine residue. Gram-negative bacterial lipoproteins are triacylated with a thioether-linked diacylglycerol moiety and an N-acyl chain. The latter is transferred from a membrane phospholipid donor to the α-amino terminus by the enzyme lipoprotein N-acyltransferase (Lnt), using an active-site cysteine thioester covalent intermediate. Many Gram-positive Firmicutes also have N-acylated lipoproteins, but the enzymes catalyzing N-acylation remain uncharacterized. The integral membrane protein Lit (lipoprotein intramolecular transacylase) from the opportunistic nosocomial pathogen Enterococcus faecalis synthesizes a specific lysoform lipoprotein (N-acyl S-monoacylglycerol) chemotype by an unknown mechanism that helps this bacterium evade immune recognition by the Toll-like receptor 2 family complex. Here, we used a deuterium-labeled lipoprotein substrate with reconstituted Lit to investigate intramolecular acyl chain transfer. We observed that Lit transfers the sn-2 ester-linked lipid from the diacylglycerol moiety to the α-amino terminus without forming a covalent thioester intermediate. Utilizing Mut-Seq to analyze an alanine scan library of Lit alleles, we identified two stretches of functionally important amino acid residues containing two conserved histidines. Topology maps based on reporter fusion assays and cysteine accessibility placed both histidines in the extracellular half of the cytoplasmic membrane. We propose a general acid base-promoted catalytic mechanism, invoking direct nucleophilic attack by the substrate α-amino group on the sn-2 ester to form a cyclic tetrahedral intermediate that then collapses to produce lyso-lipoprotein. Lit is a unique example of an intramolecular transacylase differentiated from that catalyzed by Lnt, and provides insight into the heterogeneity of bacterial lipoprotein biosynthetic systems.

De novo protein design, a retrospective

Proteins are molecular machines whose function depends on their ability to achieve complex folds with precisely defined structural and dynamic properties. The rational design of proteins from first-principles, or de novo, was once considered to be impossible, but today proteins with a variety of folds and functions have been realized. We review the evolution of the field from its earliest days, placing particular emphasis on how this endeavor has illuminated our understanding of the principles underlying the folding and function of natural proteins, and is informing the design of macromolecules with unprecedented structures and properties. An initial set of milestones in de novo protein design focused on the construction of sequences that folded in water and membranes to adopt folded conformations. The first proteins were designed from first-principles using very simple physical models. As computers became more powerful, the use of the rotamer approximation allowed one to discover amino acid sequences that stabilize the desired fold. As the crystallographic database of protein structures expanded in subsequent years, it became possible to construct proteins by assembling short backbone fragments that frequently recur in Nature. The second set of milestones in de novo design involves the discovery of complex functions. Proteins have been designed to bind a variety of metals, porphyrins, and other cofactors. The design of proteins that catalyze hydrolysis and oxygen-dependent reactions has progressed significantly. However, de novo design of catalysts for energetically demanding reactions, or even proteins that bind with high affinity and specificity to highly functionalized complex polar molecules remains an importnant challenge that is now being achieved. Finally, the protein design contributed significantly to our understanding of membrane protein folding and transport of ions across membranes. The area of membrane protein design, or more generally of biomimetic polymers that function in mixed or non-aqueous environments, is now becoming increasingly possible.

Packing of apolar side chains enables accurate design of highly stable membrane proteins

The features that stabilize the structures of membrane proteins remain poorly understood. Polar interactions contribute modestly, and the hydrophobic effect contributes little to the energetics of apolar side-chain packing in membranes. Disruption of steric packing can destabilize the native folds of membrane proteins, but is packing alone sufficient to drive folding in lipids? If so, then membrane proteins stabilized by this feature should be readily designed and structurally characterized-yet this has not been achieved. Through simulation of the natural protein phospholamban and redesign of variants, we define a steric packing code underlying its assembly. Synthetic membrane proteins designed using this code and stabilized entirely by apolar side chains conform to the intended fold. Although highly stable, the steric complementarity required for their folding is surprisingly stringent. Structural informatics shows that the designed packing motif recurs across the proteome, emphasizing a prominent role for precise apolar packing in membrane protein folding, stabilization, and evolution.

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

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

The HTLV-1 gp21 fusion peptide inhibits antigen specific T-cell activation in-vitro and in mice

The ability of the Lentivirus HIV-1 to inhibit T-cell activation by its gp41 fusion protein is well documented, yet limited data exists regarding other viral fusion proteins. HIV-1 utilizes membrane binding region of gp41 to inhibit T-cell receptor (TCR) complex activation. Here we examined whether this T-cell suppression strategy is unique to the HIV-1 gp41. We focused on T-cell modulation by the gp21 fusion peptide (FP) of the Human T-lymphotropic Virus 1 (HTLV-1), a Deltaretrovirus that like HIV infects CD4+ T-cells. Using mouse and human in-vitro T-cell models together with in-vivo T-cell hyper activation mouse model, we reveal that HTLV-1's FP inhibits T-cell activation and unlike the HIV FP, bypasses the TCR complex. HTLV FP inhibition induces a decrease in Th1 and an elevation in Th2 responses observed in mRNA, cytokine and transcription factor profiles. Administration of the HTLV FP in a T-cell hyper activation mouse model of multiple sclerosis alleviated symptoms and delayed disease onset. We further pinpointed the modulatory region within HTLV-1's FP to the same region previously identified as the HIV-1 FP active region, suggesting that through convergent evolution both viruses have obtained the ability to modulate T-cells using the same region of their fusion protein. Overall, our findings suggest that fusion protein based T-cell modulation may be a common viral trait.

Comparing side chain packing in soluble proteins, protein-protein interfaces, and transmembrane proteins

We compare side chain prediction and packing of core and non-core regions of soluble proteins, protein-protein interfaces, and transmembrane proteins. We first identified or created comparable databases of high-resolution crystal structures of these 3 protein classes. We show that the solvent-inaccessible cores of the 3 classes of proteins are equally densely packed. As a result, the side chains of core residues at protein-protein interfaces and in the membrane-exposed regions of transmembrane proteins can be predicted by the hard-sphere plus stereochemical constraint model with the same high prediction accuracies (>90%) as core residues in soluble proteins. We also find that for all 3 classes of proteins, as one moves away from the solvent-inaccessible core, the packing fraction decreases as the solvent accessibility increases. However, the side chain predictability remains high (80% within 30°) up to a relative solvent accessibility, rSASA≲0.3, for all 3 protein classes. Our results show that ≈40% of the interface regions in protein complexes are "core", that is, densely packed with side chain conformations that can be accurately predicted using the hard-sphere model. We propose packing fraction as a metric that can be used to distinguish real protein-protein interactions from designed, non-binding, decoys. Our results also show that cores of membrane proteins are the same as cores of soluble proteins. Thus, the computational methods we are developing for the analysis of the effect of hydrophobic core mutations in soluble proteins will be equally applicable to analyses of mutations in membrane proteins.

Measuring the Conformational Distance of GPCR-related Proteins Using a Joint-based Descriptor

Joint-based descriptor is a new level of macroscopic descriptor for protein structure using joints of secondary structures as a basic element. Here, we propose how the joint-based descriptor can be applied to examine the conformational distances or differences of transmembrane (TM) proteins. Specifically, we performed three independent studies that measured the global and conformational distances between GPCR A family and its related structures. First, the conformational distances of GPCR A family and other 7TM proteins were evaluated. This provided the information on the distant and close families or superfamilies to GPCR A family and permitted the identification of conserved local conformations. Second, computational models of GPCR A family proteins were validated, which enabled us to estimate how much they reproduce the native conformation of GPCR A proteins at global and local conformational level. Finally, the conformational distances between active and inactive states of GPCR proteins were estimated, which identified the difference of local conformation. The proposed macroscopic joint-based approach is expected to allow us to investigate structural features, evolutionary relationships, computational models and conformational changes of TM proteins in a more simplistic manner.

Hendra virus fusion protein transmembrane domain contributes to pre-fusion protein stability

Enveloped viruses utilize fusion (F) proteins studding the surface of the virus to facilitate membrane fusion with a target cell membrane. Fusion of the viral envelope with a cellular membrane is required for release of viral genomic material, so the virus can ultimately reproduce and spread. To drive fusion, the F protein undergoes an irreversible conformational change, transitioning from a metastable pre-fusion conformation to a more thermodynamically stable post-fusion structure. Understanding the elements that control stability of the pre-fusion state and triggering to the post-fusion conformation is important for understanding F protein function. Mutations in F protein transmembrane (TM) domains implicated the TM domain in the fusion process, but the structural and molecular details in fusion remain unclear. Previously, analytical ultracentrifugation was utilized to demonstrate that isolated TM domains of Hendra virus F protein associate in a monomer-trimer equilibrium (Smith, E. C., Smith, S. E., Carter, J. R., Webb, S. R., Gibson, K. M., Hellman, L. M., Fried, M. G., and Dutch, R. E. (2013) J. Biol. Chem. 288, 35726-35735). To determine factors driving this association, 140 paramyxovirus F protein TM domain sequences were analyzed. A heptad repeat of β-branched residues was found, and analysis of the Hendra virus F TM domain revealed a heptad repeat leucine-isoleucine zipper motif (LIZ). Replacement of the LIZ with alanine resulted in dramatically reduced TM-TM association. Mutation of the LIZ in the whole protein resulted in decreased protein stability, including pre-fusion conformation stability. Together, our data suggest that the heptad repeat LIZ contributed to TM-TM association and is important for F protein function and pre-fusion stability.

Computational Approaches for Revealing the Structure of Membrane Transporters: Case Study on Bilitranslocase

The structural and functional details of transmembrane proteins are vastly underexplored, mostly due to experimental difficulties regarding their solubility and stability. Currently, the majority of transmembrane protein structures are still unknown and this present a huge experimental and computational challenge. Nowadays, thanks to X-ray crystallography or NMR spectroscopy over 3000 structures of membrane proteins have been solved, among them only a few hundred unique ones. Due to the vast biological and pharmaceutical interest in the elucidation of the structure and the functional mechanisms of transmembrane proteins, several computational methods have been developed to overcome the experimental gap. If combined with experimental data the computational information enables rapid, low cost and successful predictions of the molecular structure of unsolved proteins. The reliability of the predictions depends on the availability and accuracy of experimental data associated with structural information. In this review, the following methods are proposed for in silico structure elucidation: sequence-dependent predictions of transmembrane regions, predictions of transmembrane helix–helix interactions, helix arrangements in membrane models, and testing their stability with molecular dynamics simulations. We also demonstrate the usage of the computational methods listed above by proposing a model for the molecular structure of the transmembrane protein bilitranslocase. Bilitranslocase is bilirubin membrane transporter, which shares similar tissue distribution and functional properties with some of the members of the Organic Anion Transporter family and is the only member classified in the Bilirubin Transporter Family. Regarding its unique properties, bilitranslocase is a potentially interesting drug target.

In Vivo Analysis of Infectivity, Fusogenicity, and Incorporation of a Mutagenic Viral Glycoprotein Library Reveals Determinants for Virus Incorporation

Enveloped viruses utilize transmembrane surface glycoproteins to gain entry into target cells. Glycoproteins from diverse viral families can be incorporated into nonnative viral particles in a process termed pseudotyping; however, the molecular mechanisms governing acquisition of these glycoproteins are poorly understood. For murine leukemia virus envelope (MLV Env) glycoprotein, incorporation into foreign viral particles has been shown to be an active process, but it does not appear to be caused by direct interactions among viral proteins. In this study, we coupled in vivo selection systems with Illumina next-generation sequencing (NGS) to test hundreds of thousands of MLV Env mutants for the ability to be enriched in viral particles and to perform other glycoprotein functions. NGS analyses on a subset of these mutants predicted that the residues important for incorporation are in the membrane-proximal external region (MPER), particularly W127 and W137, and the residues in the membrane-spanning domain (MSD) and also immediately flanking it (T140 to L163). These predictions were validated by directly measuring the impact of mutations in these regions on fusogenicity, infectivity, and incorporation. We suggest that these two regions dictate pseudotyping through interactions with specific lipid environments formed during viral assembly.

IMPORTANCE

Researchers from numerous fields routinely exploit the ability to manipulate viral tropism by swapping viral surface proteins. However, this process, termed pseudotyping, is poorly understood at the molecular level. For murine leukemia virus envelope (MLV Env) glycoprotein, incorporation into foreign viral particles is an active process, but it does not appear to occur through direct viral protein-protein interactions. In this study, we tested hundreds of thousands of MLV Env mutants for the ability to be enriched in viral particles as well as perform other glycoprotein functions. Our analyses on a subset of these mutants predict that the glycoprotein regions embedded in and immediately flanking the viral membrane dictate active incorporation into viral particles. We suggest that pseudotyping occurs through specific lipid-protein interactions at the viral assembly site.

Accurate Prediction of Contact Numbers for Multi-Spanning Helical Membrane Proteins

Prediction of the three-dimensional (3D) structures of proteins by computational methods is acknowledged as an unsolved problem. Accurate prediction of important structural characteristics such as contact number is expected to accelerate the otherwise slow progress being made in the prediction of 3D structure of proteins. Here, we present a dropout neural network-based method, TMH-Expo, for predicting the contact number of transmembrane helix (TMH) residues from sequence. Neuronal dropout is a strategy where certain neurons of the network are excluded from back-propagation to prevent co-adaptation of hidden-layer neurons. By using neuronal dropout, overfitting was significantly reduced and performance was noticeably improved. For multi-spanning helical membrane proteins, TMH-Expo achieved a remarkable Pearson correlation coefficient of 0.69 between predicted and experimental values and a mean absolute error of only 1.68. In addition, among those membrane protein-membrane protein interface residues, 76.8% were correctly predicted. Mapping of predicted contact numbers onto structures indicates that contact numbers predicted by TMH-Expo reflect the exposure patterns of TMHs and reveal membrane protein-membrane protein interfaces, reinforcing the potential of predicted contact numbers to be used as restraints for 3D structure prediction and protein-protein docking. TMH-Expo can be accessed via a Web server at www.meilerlab.org .

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/.

Characterizing the Murine Leukemia Virus Envelope Glycoprotein Membrane-Spanning Domain for Its Roles in Interface Alignment and Fusogenicity

The membrane-proximal region of murine leukemia virus envelope (Env) is a critical modulator of its functionality. We have previously shown that the insertion of one amino acid (+1 leucine) within the membrane-spanning domain (MSD) abolished protein functionality in infectivity assays. However, functionality could be restored to this +1 leucine mutant by either inserting two additional amino acids (+3 leucine) or by deleting the cytoplasmic tail domain (CTD) in the +1 leucine background. We inferred that the ectodomain and CTD have protein interfaces that have to be in alignment for Env to be functional. Here, we made single residue deletions to the Env mutant with the +1 leucine insertion to restore the interface alignment (gain of functionality) and therefore define the boundaries of the two interfaces. We identified the glycine-proline pairs near the N terminus (positions 147 and 148) and the C terminus (positions 159 and 160) of the MSD as being the boundaries of the two interfaces. Deletions between these pairs restored function, but deletions outside of them did not. In addition, the vast majority of the single residue deletions regained function if the CTD was deleted. The exceptions were four hydroxyl-containing amino acid residues (T139, T140, S143, and T144) that reside in the ectodomain interface and the proline at position 148, which were all indispensable for functionality. We hypothesize that the hydroxyl-containing residues at positions T139 and S143 could be a driving force for stabilizing the ectodomain interface through formation of a hydrogen-bonding network.

IMPORTANCE

The membrane-proximal external region (MPER) and membrane-spanning domains (MSDs) of viral glycoproteins have been shown to be critical for regulating glycoprotein fusogenicity. However, the roles of these two domains are poorly understood. We report here that point deletions and insertions within the MPER or MSD result in functionally inactive proteins. However, when the C-terminal tail domain (CTD) is deleted, the majority of the proteins remain functional. The only residues that were found to be critical for function regardless of the CTD were four hydroxyl-containing amino acids located at the C terminus of the MPER (T139 and T140) and at the N terminus of the MSD (S143 and T144) and a proline near the beginning of the MSD (P148). We demonstrate that hydrogen-bonding at positions T139 and S143 is critical for protein function. Our findings provide novel insights into the role of the MPER in regulating fusogenic activity of viral glycoproteins.

Structural Model of the Bilitranslocase Transmembrane Domain Supported by NMR and FRET Data

We present a 3D model of the four transmembrane (TM) helical regions of bilitranslocase (BTL), a structurally uncharacterized protein that transports organic anions across the cell membrane. The model was computed by considering helix-helix interactions as primary constraints, using Monte Carlo simulations. The interactions between the TM2 and TM3 segments have been confirmed by Förster resonance energy transfer (FRET) spectroscopy and nuclear magnetic resonance (NMR) spectroscopy, increasing our confidence in the model. Several insights into the BTL transport mechanism were obtained by analyzing the model. For example, the observed cis-trans Leu-Pro peptide bond isomerization in the TM3 fragment may indicate a key conformational change during anion transport by BTL. Our structural model of BTL may facilitate further studies, including drug discovery.

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.

Characteristics of protein residue-residue contacts and their application in contact prediction

Contact sites between amino acids characterize important structural features of a protein. We investigated characteristics of contact sites in a representative set of proteins and their relations between protein class or topology. For this purpose, we used a non-redundant set of 5872 protein domains, identically categorized by CATH and SCOP databases. The proteins represented alpha, beta, and alpha+beta classes. Contact maps of protein structures were obtained for a selected set of physical distances in the main backbone and separations in protein sequences. For each set a dependency between contact degree and distance parameters was quantified. We indicated residues forming contact sites most frequently and unique amino acid pairs which created contact sites most often within each structural class. Contact characteristics of specific topologies were compared to the characteristics of their protein classes showing protein groups with a distinguished contact characteristic. We showed that our results could be used to improve the performance of recent top contact predictor — direct coupling analysis. Our work provides values of contact site propensities that can be involved in bioinformatic databases.

Differential roles of tryptophan residues in the functional expression of human anion exchanger 1 (AE1, Band 3, SLC4A1)

Anion exchanger 1 (AE1) is a 95 kDa glycoprotein that facilitates Cl(-)=HCO(-)(3) exchange across the erythrocyte plasma membrane. This transport activity resides in the 52 kDa C-terminal membrane domain (Gly(361)-Val(911)) predicted to span the membrane 14 times. To explore the role of tryptophan (Trp) residues in AE1 function, the seven endogenous Trp residues in the membrane domain were mutated individually to alanine (Ala) and phenylalanine (Phe). Expression levels, cell surface abundance, inhibitor binding and transport activities of the mutants were measured upon expression in HEK-293 cells. The seven Trp residues divided into three classes according the impact of mutations on the functional expression of AE1: Class 1, dramatically decreased expression (Trp(492) and Trp(496)); Class 2, decreased expression by Ala substitution but not Phe (Trp(648), Trp(662) and Trp(723)); and Class 3, normal expression (Trp(831) and Trp(848)). The results indicate that Trp residues play differential roles in AE1 expression and function depending on their location in the protein and that Trp mutants with low expression are misfolded and retained in the endoplasmic reticulum.

General rules for the arrangements and gating motions of pore-lining helices in homomeric ion channels

The pore-lining helix (PLH) bundles are central to the function of all ion channels, as their conformational rearrangements dictate channel gating. Here we explore all plausible oligomeric arrangements of the PLH bundles within homomeric ion channels by building models using generic restraints. In particular, the distance between two neighbouring PLHs was bounded both below and above in order to avoid steric clash and allow proper packing. The resulting models provide a theoretical representation of the accessible space for oligomeric arrangements. While the represented space is confined, it encompasses nearly all the ion channel PLH bundles for which the structures are currently known. For a multitude of channels, gating models suggested by paths within the confined accessible space are in qualitative agreement with those established in previous structural and computational studies.

Rearrangements of the pore-lining helix (PLH) bundles of ion channels are central to their gating mechanisms. Here, Dai et al. use a modelling approach to define the general rules that govern the arrangements and gating motions of the PLHs in homomeric ion channels.

An analysis of oligomerization interfaces in transmembrane proteins

BACKGROUND

The amount of transmembrane protein (TM) structures solved to date is now large enough to attempt large scale analyses. In particular, extensive studies of oligomeric interfaces in the transmembrane region are now possible.

RESULTS

We have compiled the first fully comprehensive set of validated transmembrane protein interfaces in order to study their features and assess what differentiates them from their soluble counterparts.

CONCLUSIONS

The general features of TM interfaces do not differ much from those of soluble proteins: they are large, tightly packed and possess many interface core residues. In our set, membrane lipids were not found to significantly mediate protein-protein interfaces. Although no G protein-coupled receptor (GPCR) was included in the validated set, we analyzed the crystallographic dimerization interfaces proposed in the literature. We found that the putative dimer interfaces proposed for class A GPCRs do not show the usual patterns of stable biological interfaces, neither in terms of evolution nor of packing, thus they likely correspond to crystal interfaces. We cannot however rule out the possibility that they constitute transient or weak interfaces. In contrast we do observe a clear signature of biological interface for the proposed dimer of the class F human Smoothened receptor.

Quantitative residue-level structure-evolution relationships in the yeast membrane proteome

Membrane proteins exist in distinctly different environments than do soluble proteins, resulting in differences between their respective biophysical and evolutionary properties. In comparison with soluble proteins, relatively little is known about how the unique biophysical properties of membrane proteins affect their evolutionary properties at the residue level. In particular, transmembrane (TM) regions of membrane proteins tend to be more conserved than regions outside of the membrane (extramembrane [EM] regions), but the mechanisms underlying this phenomenon are not well understood. Here, we combine homology-based high-resolution three-dimensional protein models with rigorous evolutionary rate calculations to quantitatively assess residue-level structure-evolution relationships in the yeast membrane proteome. We find that residue evolutionary rate increases linearly with decreasing residue burial, regardless of the hydrophobic or hydrophilic nature of the solvent environment. This finding supports a direct relationship between a residue's selective constraint and the extent of its packing interactions with neighboring residues, independent of hydrophobic effects. Most importantly, for a fixed degree of burial, residues from TM regions tend to evolve more slowly than residues from EM regions. We attribute this difference to the increased importance of packing constraints and the decreased importance of hydrophobic effects in TM regions. This additional selective constraint on TM residues plays a dominant role in explaining why TM regions evolve more slowly than EM regions. In addition to revealing the universality of the linear relationship between residue burial and selective constraint across solvent environments, our work highlights the distinct residue-level evolutionary consequences imposed by the unique biophysical properties of the membrane environment.

Membrane protein native state discrimination by implicit membrane models

Four implicit membrane models [IMM1, generalized Born (GB)-surface area-implicit membrane (GBSAIM), GB with a simple switching (GBSW), and heterogeneous dielectric GB (HDGB)] were tested for their ability to discriminate the native conformation of five membrane proteins from 450 decoys generated by the Rosetta-Membrane program. The energy ranking of the native state and Z-scores were used to assess the performance of the models. The effect of membrane thickness was examined and was found to be substantial. Quite satisfactory discrimination was achieved with the all-atom IMM1 and GBSW models at 25.4 Å thickness and with the HDGB model at 28.5 Å thickness. The energy components by themselves were not discriminative. Both van der Waals and electrostatic interactions contributed to native state discrimination, to a different extent in each model. Computational efficiency of the models decreased in the order: extended-atom IMM1 > all-atom IMM1 > GBSAIM > GBSW > HDGB. These results encourage the further development and use of implicit membrane models for membrane protein structure prediction.

Comparison of fragments comprising the first two helices of the human Y4 and the yeast Ste2p G-protein-coupled receptors

Solution NMR techniques are used to determine the structure and the topology of micelle integration of a large fragment of the Y4 receptor, a human G-protein-coupled receptor, that contains the entire N-terminal domain plus the first two transmembrane (TM) segments. The structure calculations reveal that the putative TM helices are indeed helical to a large extent, but that interruptions of secondary structure occur close to internal polar or charged residues. This view is supported by (15)N relaxation data, amide-water exchange rates, and attenuations from micelle-integrating spin labels. No contacts between different helices are observed. This is in contrast to a similar TM1-TM2 fragment from the yeast Ste2p receptor for which locations of the secondary and the tertiary structure agreed well with the predictions from a homology model. The difference in structure is discussed in terms of principal biophysical properties of residues within central regions of the putative TM helices. Overall, using the biophysical scale of Wimley and White the TM regions of Ste2p display much more favorable free energies for membrane integration. Accordingly, the full secondary structure and the tertiary structure in TM1-TM2 of the Y4 receptor is likely to be formed only when tertiary contacts with other TM segments are created during folding of the receptor.

Replica exchange Monte-Carlo simulations of helix bundle membrane proteins: rotational parameters of helices

We propose a united-residue model of membrane proteins to investigate the structures of helix bundle membrane proteins (HBMPs) using coarse-grained (CG) replica exchange Monte-Carlo (REMC) simulations. To demonstrate the method, it is used to identify the ground state of HBMPs in a CG model, including bacteriorhodopsin (BR), halorhodopsin (HR), and their subdomains. The rotational parameters of transmembrane helices (TMHs) are extracted directly from the simulations, which can be compared with their experimental measurements from site-directed dichroism. In particular, the effects of amphiphilic interaction among the surfaces of TMHs on the rotational angles of helices are discussed. The proposed CG model gives a reasonably good structure prediction of HBMPs, as well as a clear physical picture for the packing, tilting, orientation, and rotation of TMHs. The root mean square deviation (RMSD) in coordinates of C(α) atoms of the ground state CG structure from the X-ray structure is 5.03 Å for BR and 6.70 Å for HR. The final structure of HBMPs is obtained from the all-atom molecular dynamics simulations by refining the predicted CG structure, whose RMSD is 4.38 Å for BR and 5.70 Å for HR.

Developing a high-quality scoring function for membrane protein structures based on specific inter-residue interactions

Membrane proteins are of particular biological and pharmaceutical importance, and computational modeling and structure prediction approaches play an important role in studies of membrane proteins. Developing an accurate model quality assessment program is of significance to the structure prediction of membrane proteins. Few such programs are proposed that can be applied to a broad range of membrane protein classes and perform with high accuracy. We developed a new model scoring function Interaction-based Quality assessment (IQ), based on the analysis of four types of inter-residue interactions within the transmembrane domains of helical membrane proteins. This function was tested using three high-quality model sets: all 206 models of GPCR Dock 2008, all 284 models of GPCR Dock 2010, and all 92 helical membrane protein models of the HOMEP set. For all three sets, the scoring function can select the native structures among all of the models with the success rates of 93, 85, and 100% respectively. For comparison, these three model sets were also adopted for a recently published model assessment program for membrane protein structures, ProQM, which gave the success rates of 85, 79, and 92% separately. These results suggested that IQ outperforms ProQM when only the transmembrane regions of the models are considered. This scoring function should be useful for the computational modeling of membrane proteins.

Mutational analysis of the West Nile virus NS4B protein

West Nile virus NS4B is a small hydrophobic nonstructural protein approximately 27 kDa in size whose function is poorly understood. Amino acid substitutions were introduced into the NS4B protein primarily targeting two distinct regions; the N-terminal domain (residues 35 through 60) and the central hydrophobic domain (residues 95 through 120). Only the NS4B P38G substitution was associated with both temperature-sensitive and small-plaque phenotypes. Importantly, this mutation was found to attenuate neuroinvasiveness greater than 10,000,000-fold and lower viremia titers compared to the wild-type NY99 virus in a mouse model. Full genome sequencing of the NS4B P38G mutant virus revealed two unexpected mutations at NS4B T116I and NS3 N480H (P38G/T116I/N480H), however, neither mutation alone was temperature sensitive or attenuated in mice. Following incubation of P38G/T116I/N480H at 41°C, five mutants encoding compensatory substitutions in the NS4B protein exhibited a reduction in the temperature-sensitive phenotype and reversion to a virulent phenotype in the mouse model.

Efficiency of detergents at maintaining membrane protein structures in their biologically relevant forms

High-resolution structural analysis of membrane proteins by X-ray crystallography or solution NMR spectroscopy often requires their solubilization in the membrane-mimetic environments of detergents. Yet the choice of a detergent suitable for a given study remains largely empirical. In the present work, we considered the micelle-crystallized structures of lactose permease (LacY), the sodium/galactose symporter (vSGLT), the vitamin B(12) transporter (BtuCD), and the arginine/agmatine antiporter (AdiC). Representative transmembrane (TM) segments were selected from these proteins based on their relative contact(s) with water, lipid, and/or within the protein, and were synthesized as Lys-tagged peptides. Each peptide was studied by circular dichroism and fluorescence spectroscopy in water, and in the presence of the detergents sodium dodecylsulfate (SDS, anionic); n-dodecyl phosphatidylcholine (DPC, zwitterionic); n-dodecyl-β-d-maltoside (DDM, neutral); and n-octyl-β-d-glucoside (OG, neutral, varying acyl tail length). We found that (i) the secondary structures of the TM segments were statistically indistinguishable in the four detergents studied; and (ii) a strong correlation exists between the extent of helical structure of each individual TM segment in detergents with its helicity level as it exists in the full-length protein, indicating that helix adoption is fundamentally the same in both environments. The denaturing properties of so-called 'harsh' detergents may thus largely be due to their interactions with non-membranous regions of proteins. Given the consistency of structural features observed for each TM segment in a variety of micellar media, the overall results suggest that the structure likely corresponds to its relevant biological form in the intact protein in its native lipid bilayer environment.

Identifying and measuring transmembrane helix-helix interactions by FRET

Specific interactions between helical transmembrane domains (TMs) play essential roles in the mechanisms governing the folding, stability and assembly of integral membrane proteins. Thus, it is appealing to identify helix-helix contacts and to seek the structural determinants of such interactions at the molecular level. Here, we provide a protocol for detecting and measuring specific helix-helix interactions in liposomes by Förster resonance energy transfer (FRET), using peptides corresponding to the TM domains of an integral membrane protein. We give a detailed procedure and practical guidelines on how to design, prepare, handle, and characterize fluorescently labeled TM peptides reconstituted in large unilamellar lipid vesicles. We also discuss some critical aspects of FRET measurements to ensure the correct analysis and interpretation of spectral data. Our method uses tryptophan/pyrene as the donor-acceptor FRET pair, but it can be easily adapted to other fluorescence pairs and to other membrane mimetic environments. The ability to identify crucial interhelical contacts is a valuable tool for the study of the stability, assembly, and function of the important and experimentally challenging helical membrane proteins.

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.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

Amino acid interaction preferences in helical membrane proteins

Membrane proteins are involved in a number of important biological functions. Yet, they are poorly understood from the structure and folding point of view. The external environment being drastically different from that of globular proteins, the intra-protein interactions in membrane proteins are also expected to be different. Hence, statistical potentials representing the features of inter-residue interactions based exclusively on the structures of membrane proteins are much needed. Currently, a reasonable number of structures are available, making it possible to undertake such an analysis on membrane proteins. In this study we have examined the inter-residue interaction propensities of amino acids in the membrane spanning regions of the alpha-helical membrane (HM) proteins. Recently we have shown that valuable information can be obtained on globular proteins by the evaluation of the pair-wise interactions of amino acids by classifying them into different structural environments, based on factors such as the secondary structure or the number of contacts that a residue can make. Here we have explored the possible ways of classifying the intra-protein environment of HM proteins and have developed scoring functions based on different classification schemes. On evaluation of different schemes, we find that the scheme which classifies amino acids to different intra-contact environment is the most promising one. Based on this classification scheme, we also redefine the hydrophobicity scale of amino acids in HM proteins.

Conserved GXXXG- and S/T-like motifs in the transmembrane domains of NS4B protein are required for hepatitis C virus replication

Hepatitis C virus (HCV) nonstructural protein 4B (NS4B) is an integral membrane protein, which plays an important role in the organization and function of the HCV replication complex (RC). Although much is understood about its amphipathic N-terminal and C-terminal domains, we know very little about the role of the transmembrane domains (TMDs) in NS4B function. We hypothesized that in addition to anchoring NS4B into host membranes, the TMDs are engaged in intra- and intermolecular interactions required for NS4B structure/function. To test this hypothesis, we have engineered a chimeric JFH1 genome containing the Con1 NS4B TMD region. The resulting virus titers were greatly reduced from those of JFH1, and further analysis indicated a defect in genome replication. We have mapped this incompatibility to NS4B TMD1 and TMD2 sequences, and we have defined putative TMD dimerization motifs (GXXXG in TMD2 and TMD3; the S/T cluster in TMD1) as key structural/functional determinants. Mutations in each of the putative motifs led to significant decreases in JFH1 replication. Like most of the NS4B chimeras, mutant proteins had no negative impact on NS4B membrane association. However, some mutations led to disruption of NS4B foci, implying that the TMDs play a role in HCV RC formation. Further examination indicated that the loss of NS4B foci correlates with the destabilization of NS4B protein. Finally, we have identified an adaptive mutation in the NS4B TMD2 sequence that has compensatory effects on JFH1 chimera replication. Taken together, these data underscore the functional importance of NS4B TMDs in the HCV life cycle.

Peptide model helices in lipid membranes: insertion, positioning, and lipid response on aggregation studied by X-ray scattering

Studying membrane active peptides or protein fragments within the lipid bilayer environment is particularly challenging in the case of synthetically modified, labeled, artificial, or recently discovered native structures. For such samples the localization and orientation of the molecular species or probe within the lipid bilayer environment is the focus of research prior to an evaluation of their dynamic or mechanistic behavior. X-ray scattering is a powerful method to study peptide/lipid interactions in the fluid, fully hydrated state of a lipid bilayer. For one, the lipid response can be revealed by observing membrane thickening and thinning as well as packing in the membrane plane; at the same time, the distinct positions of peptide moieties within lipid membranes can be elucidated at resolutions of up to several angstroms by applying heavy-atom labeling techniques. In this study, we describe a generally applicable X-ray scattering approach that provides robust and quantitative information about peptide insertion and localization as well as peptide/lipid interaction within highly oriented, hydrated multilamellar membrane stacks. To this end, we have studied an artificial, designed β-helical peptide motif in its homodimeric and hairpin variants adopting different states of oligomerization. These peptide lipid complexes were analyzed by grazing incidence diffraction (GID) to monitor changes in the lateral lipid packing and ordering. In addition, we have applied anomalous reflectivity using synchrotron radiation as well as in-house X-ray reflectivity in combination with iodine-labeling in order to determine the electron density distribution ρ(z) along the membrane normal (z axis), and thereby reveal the hydrophobic mismatch situation as well as the position of certain amino acid side chains within the lipid bilayer. In the case of multiple labeling, the latter technique is not only applicable to demonstrate the peptide's reconstitution but also to generate evidence about the relative peptide orientation with respect to the lipid bilayer.

Membrane anchoring and interaction between transmembrane domains are crucial for K+ channel function

The small viral channel Kcv is a Kir-like K(+) channel of only 94 amino acids. With this simple structure, the tetramer of Kcv represents the pore module of all complex K(+) channels. To examine the structural contribution of the transmembrane domains (TMDs) to channel function, we performed Ala scanning mutagenesis of the two domains and tested the functionality of the mutants in a yeast complementation assay. The data reveal, in combination with computational models, that the upper halves of both TMDs, which face toward the external medium, are rather rigid, whereas the inner parts are more flexible. The rigidity of the outer TMD is conferred by a number of essential aromatic amino acids that face the membrane and probably anchor this domain in the bilayer. The inner TMD is intimately connected with the rigid part of the outer TMD via π···π interactions between a pair of aromatic amino acids. This structural principle is conserved within the viral K(+) channels and also present in Kir2.2, implying a general importance of this architecture for K(+) channel function.

Non-specific binding and general cross-reactivity of Y receptor agonists are correlated and should importantly depend on their acidic sectors

Non-specific binding of Y receptor agonists to intact CHO cells, and to CHO cell or rat brain particulates, is much greater for human neuropeptide Y (hNPY) compared to porcine peptide Y (pPYY), and especially relative to human pancreatic polypeptide (hPP). This binding of hNPY is reduced by alkali cations in preference to non-ionic chaotrope urea, while the much lower non-specific binding of pPYY is more sensitive to urea. The difference could mainly be due to the 10-16 stretch in 36-residue Y agonists (residues 8-14 in N-terminally clipped 34-peptides), located in the sector that contains all acidic residues of physiological Y agonists. Anionic pairs containing aspartate in the 10-16 zone could be principally responsible for non-specific attachments, but may also aid the receptor site binding. Two such pairs are found in hNPY, one in pPYY, and none in hPP. The hydroxyl amino acid residue at position 13 in mammalian PYY and PP molecules could lower conformational plasticity and the non-selective binding via intrachain hydrogen bonding. The acidity of this tract could also be important in agonist selectivity of the Y receptor subtypes. The differences point to an evolutionary reduction of promiscuous protein binding from NPY to PP, and should also be important for Y agonist selectivity within NPY receptor group, and correlate with partial agonism and out-of group cross-reactivity with other receptors.