Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α‐helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α‐helices in both soluble and membrane proteins, rather than just transmembrane α‐helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length‐matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent. Proteins 2014; 82:1960–1970. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Helix capping.

Helix-capping motifs are specific patterns of hydrogen bonding and hydrophobic interactions found at or near the ends of helices in both proteins and peptides. In an alpha-helix, the first four >N-H groups and last four >C=O groups necessarily lack intrahelical hydrogen bonds. Instead, such groups are often capped by alternative hydrogen bond partners. This review enlarges our earlier hypothesis (Presta LG, Rose GD. 1988. Helix signals in proteins. Science 240:1632-1641) to include hydrophobic capping. A hydrophobic interaction that straddles the helix terminus is always associated with hydrogen-bonded capping. From a global survey among proteins of known structure, seven distinct capping motifs are identified-three at the helix N-terminus and four at the C-terminus. The consensus sequence patterns of these seven motifs, together with results from simple molecular modeling, are used to formulate useful rules of thumb for helix termination. Finally, we examine the role of helix capping as a bridge linking the conformation of secondary structure to supersecondary structure.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Phosphorylation stabilizes the N‐termini of α‐helices

The role of phosphorylation in stabilizing the N‐termini of α‐helices is examined using computer simulations of model peptides. The models comprise either a phosphorylated or unphosphorylated serine at the helix N‐terminus, followed by nine alanines. Monte Carlo/stochastic Dynamics simulations were performed on the model helices. The simulations revealed a distinct stabilization of the helical conformation at the N‐terminus after phosphorylation. The stabilization was attributable to favorable electrostatic interactions between the phosphate and the helix backbone. However, direct helix capping by the phosphorylated sidechain was not observed. The results of the calculations are consistent with experimental evidence on the stabilization of helices by phosphates and other anions. © 1999 John Wiley & Sons, Inc. Biopoly 49: 225–233, 1999     

Length‐dependent stability of α‐helical peptide upon adsorption to single‐walled carbon nanotube

Earlier studies have shown that the helical content of α‐helical peptide decreases upon its interaction with carbon nanotube (CNT). Further, the length of the α‐helix varies from few residues in the small globular protein to several number of residues in structural and membrane proteins. In structural and membrane proteins, helices are widely present as the supercoil i.e., helical bundles. Thus, in this study, the length‐dependent interaction pattern of α‐helical peptides with CNT and the stability of isolated α‐helical fragment versus supercoiled helical bundle upon interaction with CNT have been investigated using classical molecular dynamics (MD) simulation. Results reveal that the disruption in the helical motif on interaction with CNT is directly proportional to the length of the helix. Also it is found that the shorter helix does not undergo noticeable changes in the helicity upon adsorption with CNT. On the other hand, helicity of longer peptides is considerably affected by its interaction with CNT. In contrast to the known fact that the stability of the helix increases with its length, the disruption in the helical peptide increases with its length upon its interaction with CNT. Comparison of results shows that structural changes in the isolated helical fragment are higher than that in supercoiled helix. In fact, helical chain in supercoiled bundle does not undergo significant changes in the helicity upon interaction with CNT. Both the length of the helical peptide and the inherent stability of the helical unit in the supercoiled helix influence the interaction pattern with the CNT. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 357–369, 2013.     

Structural comparison and classification of alpha‐helical transmembrane domains based on helix interaction patterns

Structural classification of membrane proteins is still in its infancy due to the relative paucity of available three‐dimensional structures compared with soluble proteins. However, recent technological advances in protein structure determination have led to a significant increase in experimentally known membrane protein folds, warranting exploration of the structural universe of membrane proteins. Here, a new and completely membrane protein specific structural classification system is introduced that classifies α‐helical membrane proteins according to common helix architectures. Each membrane protein is represented by a helix interaction graph depicting transmembrane helices with their pairwise interactions resulting from individual residue contacts. Subsequently, proteins are clustered according to similarities among these helix interaction graphs using a newly developed structural similarity score called HISS. As HISS scores explicitly disregard structural properties of loop regions, they are more suitable to capture conserved transmembrane helix bundle architectures than other structural similarity scores. Importantly, we are able to show that a classification approach based on helix interaction similarity closely resembles conventional structural classification databases such as SCOP and CATH implying that helix interactions are one of the major determinants of α‐helical membrane protein folds. Furthermore, the classification of all currently available membrane protein structures into 20 recurrent helix architectures and 15 singleton proteins demonstrates not only an impressive variability of membrane helix bundles but also the conservation of common helix interaction patterns among proteins with distinctly different sequences. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Stabilization of conformationally dynamic helices by covalently attached acyl chains

Acylation of proteins is known to mediate membrane attachment and to influence subcellular sorting. Here, we report that acylation can stabilize secondary structure. Circular dichroism spectroscopy showed that N‐terminal attachment of acyl chains decreases the ability of an intrinsically flexible hydrophobic model peptide to refold from an α‐helical state to β‐sheet in response to changing solvent conditions. Acylation also stabilized the membrane‐embedded α‐helix. This increase of global helix stability did not result from decreased local conformational dynamics of the helix backbone as assessed by deuterium/hydrogen‐exchange experiments. We concluded that acylation can stabilize the structure of intrinsically dynamic helices and may thus prevent misfolding.     

Conformation of membrane‐bound proteins revealed by vacuum‐ultraviolet circular‐dichroism and linear‐dichroism spectroscopy

Knowledge of the conformations of a water‐soluble protein bound to a membrane is important for understanding the membrane‐interaction mechanisms and the membrane‐mediated functions of the protein. In this study we applied vacuum‐ultraviolet circular‐dichroism (VUVCD) and linear‐dichroism (LD) spectroscopy to analyze the conformations of α‐lactalbumin (LA), thioredoxin (Trx), and β‐lactoglobulin (LG) bound to phosphatidylglycerol liposomes. The VUVCD analysis coupled with a neural‐network analysis showed that these three proteins have characteristic helix‐rich conformations involving several helical segments, of which two amphiphilic or hydrophobic segments take part in interactions with the liposome. The LD analysis predicted the average orientations of these helix segments on the liposome: two amphiphilic helices parallel to the liposome surface for LA, two hydrophobic helices perpendicular to the liposome surface for Trx, and a hydrophobic helix perpendicular to and an amphiphilic helix parallel to the liposome surface for LG. This sequence‐level information about the secondary structures and orientations was used to formulate interaction models of the three proteins at the membrane surface. This study demonstrates the validity of a combination of VUVCD and LD spectroscopy in conformational analyses of membrane‐binding proteins, which are difficult targets for X‐ray crystallography and nuclear magnetic resonance spectroscopy. Proteins 2016; 84:349–359. © 2016 Wiley Periodicals, Inc.     

Current status of membrane protein structure classification

For over 2 decades, continuous efforts to organize the jungle of available protein structures have been underway. Although a number of discrepancies between different classification approaches for soluble proteins have been reported, the classification of membrane proteins has so far not been comparatively studied because of the limited amount of available structural data. Here, we present an analysis of α‐helical membrane protein classification in the SCOP and CATH databases. In the current set of 63 α‐helical membrane protein chains having between 1 and 13 transmembrane helices, we observed a number of differently classified proteins both regarding their domain and fold assignment. The majority of all discrepancies affect single transmembrane helix, two helix hairpin, and four helix bundle domains, while domains with more than five helices are mostly classified consistently between SCOP and CATH. It thus appears that the structural constraints imposed by the lipid bilayer complicate the classification of membrane proteins with only few membrane‐spanning regions. This problem seems to be specific for membrane proteins as soluble four helix bundles, not restrained by the membrane, are more consistently classified by SCOP and CATH. Our findings indicate that the structural space of small membrane helix bundles is highly continuous such that even minor differences in individual classification procedures may lead to a significantly different classification. Membrane proteins with few helices and limited structural diversity only seem to be reasonably classifiable if the definition of a fold is adapted to include more fine‐grained structural features such as helix–helix interactions and reentrant regions. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Structure of the hypothetical protein Ton 1535 from Thermococcus onnurineus NA1 reveals unique structural properties by a left‐handed helical turn in normal α‐solenoid protein

The crystal structure of Ton1535, a hypothetical protein from Thermococcus onnurineus NA1, was determined at 2.3 Å resolution. With two antiparallel α‐helices in a helix‐turn‐helix motif as a repeating unit, Ton1535 consists of right‐handed coiled N‐ and C‐terminal regions that are stacked together using helix bundles containing a left‐handed helical turn. One left‐handed helical turn in the right‐handed coiled structure produces two unique structural properties. One is the presence of separated concave grooves rather than one continuous concave groove, and the other is the contribution of α‐helices on the convex surfaces of the N‐terminal region to the extended surface of the concave groove of the C‐terminal region and vice versa. Proteins 2014; 82:1072–1078. © 2013 Wiley Periodicals, Inc.     

Structural basis of intramitochondrial phosphatidic acid transport mediated by Ups1‐Mdm35 complex

Ups1 forms a complex with Mdm35 and is critical for the transport of phosphatidic acid (PA) from the mitochondrial outer membrane to the inner membrane. We report the crystal structure of the Ups1‐Mdm35‐PA complex and the functional characterization of Ups1‐Mdm35 in PA binding and transfer. Ups1 features a barrel‐like structure consisting of an antiparallel β‐sheet and three α‐helices. Mdm35 adopts a three‐helical clamp‐like structure to wrap around Ups1 to form a stable complex. The β‐sheet and α‐helices of Ups1 form a long tunnel‐like pocket to accommodate the substrate PA, and a short helix α2 acts as a lid to cover the pocket. The hydrophobic residues lining the pocket and helix α2 are critical for PA binding and transfer. In addition, a hydrophilic patch on the surface of Ups1 near the PA phosphate‐binding site also plays an important role in the function of Ups1‐Mdm35. Our study reveals the molecular basis of the function of Ups1‐Mdm35 and sheds new light on the mechanism of intramitochondrial phospholipid transport by the MSF1/PRELI family proteins.     

Design of stable α‐helices using global sequence optimization

The rational design of peptide and protein helices is not only of practical importance for protein engineering but also is a useful approach in attempts to improve our understanding of protein folding. Recent modifications of theoretical models of helix‐coil transitions allow accurate predictions of the helix stability of monomeric peptides in water and provide new possibilities for protein design. We report here a new method for the design of α‐helices in peptides and proteins using AGADIR, the statistical mechanical theory for helix‐coil transitions in monomeric peptides and the tunneling algorithm of global optimization of multidimensional functions for optimization of amino acid sequences. CD measurements of helical content of peptides with optimized sequences indicate that the helical potential of protein amino acids is high enough to allow formation of stable α‐helices in peptides as short as of 10 residues in length. The results show the maximum achievable helix content (HC) of short peptides with fully optimized sequences at 5 °C is expected to be ～70–75%. Under certain conditions the method can be a powerful practical tool for protein engineering. Unlike traditional approaches that are often used to increase protein stability by adding a few favorable interactions to the protein structure, this method deals with all possible sequences of protein helices and selects the best one from them. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

An automated method for consistent helix assignment using turn information

A new automated helix assignment method is presented that leads to a more consistent definition of the helix termini, especially of the helix C-terminus. The method assigns a helix to segments of protein chain where adjacent helical turn structures are observed, capped by specific distorted turn types (e.g., open helical turns without a hydrogen bond) or capping motifs (e.g., the Schellman motif). Helix termini are detected by observing the behavior of the NH group in N-termini and the CO group in C-termini; in each case, the respective group must be free to interact with hydrogen bonding partners outside of the putative helix for a helix terminus to be assigned. The presented assignment method and SHAFT-assigned helices are part of Secbase and are made available with Relibase+ 3.0 and the free web version of Relibase 3.0. The method can also be used for the helix assignments of additional protein structures.     

Development of electron spin echo envelope modulation spectroscopy to probe the secondary structure of recombinant membrane proteins in a lipid bilayer

Membrane proteins conduct many important biological functions essential to the survival of organisms. However, due to their inherent hydrophobic nature, it is very difficult to obtain structural information on membrane‐bound proteins using traditional biophysical techniques. We are developing a new approach to probe the secondary structure of membrane proteins using the pulsed EPR technique of Electron Spin Echo Envelope Modulation (ESEEM) Spectroscopy. This method has been successfully applied to model peptides made synthetically. However, in order for this ESEEM technique to be widely applicable to larger membrane protein systems with no size limitations, protein samples with deuterated residues need to be prepared via protein expression methods. For the first time, this study shows that the ESEEM approach can be used to probe the local secondary structure of a ²H‐labeled d₈‐Val overexpressed membrane protein in a membrane mimetic environment. The membrane‐bound human KCNE1 protein was used with a known solution NMR structure to demonstrate the applicability of this methodology. Three different α‐helical regions of KCNE1 were probed: the extracellular domain (Val21), transmembrane domain (Val50), and cytoplasmic domain (Val95). These results indicated α‐helical structures in all three segments, consistent with the micelle structure of KCNE1. Furthermore, KCNE1 was incorporated into a lipid bilayer and the secondary structure of the transmembrane domain (Val50) was shown to be α‐helical in a more native‐like environment. This study extends the application of this ESEEM approach to much larger membrane protein systems that are difficult to study with X‐ray crystallography and/or NMR spectroscopy.     

PackHelix: a tool for helix-sheet packing during protein structure prediction

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Predicting the topology and constructing the geometry of structural motifs involving α‐helices and/or β‐strands are therefore key steps for accurate prediction of protein structure. While many efforts have focused on how to pack helices and on how to sample exhaustively the topologies and geometries of multiple strands forming a β‐sheet in a protein, there has been little progress on generating native‐like packings of helices on sheets. We describe a method that can generate the packing of multiple helices on a given β‐sheet for αβα sandwich type protein folds. This method mines the results of a statistical analysis of the conformations of αβ₂ motifs in protein structures to provide input values for the geometric attributes of the packing of a helix on a sheet. It then proceeds with a geometric builder that generates multiple arrangements of the helices on the sheet of interest by sampling through these values and performing consistency checks that guarantee proper loop geometry between the helices and the strands, minimal number of collisions between the helices, and proper formation of a hydrophobic core. The method is implemented as a module of ProteinShop. Our results show that it produces structures that are within 4–6 Å RMSD of the native one, regardless of the number of helices that need to be packed, though this number may increase if the protein has several helices between two consecutive strands in the sequence that pack on the sheet formed by these two strands. Proteins 2011; Published 2011 Wiley‐Liss, Inc.     

Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Folding of Aquaporin 1: Multiple evidence that helix 3 can shift out of the membrane core

The folding of most integral membrane proteins follows a two‐step process: initially, individual transmembrane helices are inserted into the membrane by the Sec translocon. Thereafter, these helices fold to shape the final conformation of the protein. However, for some proteins, including Aquaporin 1 (AQP1), the folding appears to follow a more complicated path. AQP1 has been reported to first insert as a four‐helical intermediate, where helix 2 and 4 are not inserted into the membrane. In a second step, this intermediate is folded into a six‐helical topology. During this process, the orientation of the third helix is inverted. Here, we propose a mechanism for how this reorientation could be initiated: first, helix 3 slides out from the membrane core resulting in that the preceding loop enters the membrane. The final conformation could then be formed as helix 2, 3, and 4 are inserted into the membrane and the reentrant regions come together. We find support for the first step in this process by showing that the loop preceding helix 3 can insert into the membrane. Further, hydrophobicity curves, experimentally measured insertion efficiencies and MD‐simulations suggest that the barrier between these two hydrophobic regions is relatively low, supporting the idea that helix 3 can slide out of the membrane core, initiating the rearrangement process.     

Structural features of transmembrane helices

A total of 160 transmembrane helices of 15 non-homologous high-resolution X-ray protein structures have been analyzed in respect of their structural features. The dihedral angles and hydrogen bonds of the helical sections that span the hydrophobic interior of the lipid bilayer have been investigated. The Ramachandran plot of protein channels and solute transporters exhibit a significant shift Delta (phi- and psi-angles) of Delta mean (+4.5 degrees and -5.4 degrees ), compared to a reference group of 151 alpha-helices of the same average length derived from water-soluble globular proteins. At the C-termini of transmembrane helices structural motifs equivalent to the Gly-caps of helices in globular proteins have been found, with two third of the transmembrane Gly-caps taking up a primary structure that is typically not found at helix termini exposed to a polar solvent. The structural particularities reported here are relevant for the three-dimensional modelling of membrane protein structures.     

Bihelix: Towards de novo structure prediction of an ensemble of G‐protein coupled receptor conformations

G‐Protein Coupled Receptors (GPCRs) play a critical role in cellular signal transduction pathways and are prominent therapeutic targets. Recently there has been major progress in obtaining experimental structures for a few GPCRs. Each GPCR, however, exhibits multiple conformations that play a role in their function and we have been developing methods aimed at predicting structures for all these conformations. Analysis of available structures shows that these conformations differ in relative helix tilts and rotations. The essential issue is, determining how to orient each of the seven helices about its axis since this determines how it interacts with the other six helices. Considering all possible helix rotations to ensure that no important packings are overlooked, and using rotation angle increments of 30° about the helical axis would still lead to 12⁷ or 35 million possible conformations each with optimal residue positions. We show in this paper how to accomplish this. The fundamental idea is to optimize the interactions between each pair of contacting helices while ignoring the other 5 and then to estimate the energies of all 35 million combinations using these pair‐wise interactions. This BiHelix approach dramatically reduces the effort to examine the complete set of conformations and correctly identifies the crystal packing for the experimental structures plus other near‐native packings we believe may play an important role in activation. This approach also enables a detailed structural analysis of functionally distinct conformations using helix‐helix interaction energy landscapes and should be useful for other helical transmembrane proteins as well. Proteins 2012. © 2011 Wiley Periodicals, Inc.