Article Information

• PDF (343 kb)
• Full Text with Large Figures
• Cited by in Scopus (5)

PubMed

• Articles by Peter Werner Hildebrand
• Articles by Robert Preissner

Related Articles

• Direct Simulation of Transmembrane Helix Association: Role o...

  Direct Simulation of Transmembrane Helix Association: Role of Asparagines Biophysical Journal, Volume 87, Issue 3, 1 September 2004, Pages 1650-1656 Thomas Stockner, Walter L. Ash, Justin L. MacCallum and D. Peter Tieleman Abstract The forces contributing to the association of transmembrane helices in folded membrane proteins have received considerable attention recently. In this study we investigate the importance of hydrogen bonding by studying the effect of a single Asn residue in the center of an otherwise hydrophobic transmembrane peptide using computer simulations. We use the model peptide MS1 which has been derived from the leucine zipper coiled-coil dimer of the transcription factor peptide GCN4-P1. We follow the trajectory of 36 initially monomeric MS1 transmembrane helical peptides in a membrane-mimicking octane layer as they associate into larger structures. These peptides predominately form dimers. The interaction between the polar asparagine residues, capable of simultaneously being a hydrogen-bond donor and acceptor, contributes strongly to the stability of associated helices. Only dimers with interhelical hydrogen bonds form stable structures, whereas aggregates without any hydrogen-bonding interactions form very transient structures. We examine the hydrogen-bonding patterns and find that there are two forms of dimer, one with symmetric hydrogen bonds and one with asymmetric hydrogen bonds. Based on the structures in our simulation we propose a model with a monomer ↔ symmetric dimer ↔ asymmetric dimer ↔ trimer equilibrium. Abstract | Full Text | PDF (324 kb)

• Comparison of Helix Interactions in Membrane and Soluble α-B...

  Comparison of Helix Interactions in Membrane and Soluble α-Bundle Proteins Biophysical Journal, Volume 82, Issue 5, 1 May 2002, Pages 2720-2736 Markus Eilers, Ashish B. Patel, Wei Liu and Steven O. Smith Abstract Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble -bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability. Abstract | Full Text | PDF (523 kb)

• Complete Predicted Three-Dimensional Structure of the Facili...

  Complete Predicted Three-Dimensional Structure of the Facilitator Transmembrane Protein and Hepatitis C Virus Receptor CD81: Conserved and Variable Structural Domains in the Tetraspanin Superfamily Biophysical Journal, Volume 90, Issue 1, 1 January 2006, Pages 212-227 Michel Seigneuret Abstract Tetraspanins are a superfamily of transmembrane proteins implicated in cellular development, motility, and activation through their interactions with a large range of proteins and with specific membrane microdomains. The complete three-dimensional structure of the tetraspanin CD81 has been predicted by molecular modeling and from the crystallographic structure of the EC2 large extracellular domain. Periodicity of sequence conservation, homology modeling, secondary structure prediction, and protein docking were used. The transmembrane domain appears organized as a four-stranded left-handed coiled coil directly connecting to two helices of the EC2. A smaller extracellular loop EC1 contains a small largely hydrophobic -strand that packs in a conserved hydrophobic groove of the EC2. The palmitoylable intracellular N-terminal segment forms an amphipathic membrane-parallel helix. Structural variability occurs mainly in an hypervariable subdomain of the EC2 and in intracellular regions. Therefore, the variable interaction selectivity of tetraspanins originates both from sequence variability within structurally conserved domains and from the occurrence of small structurally variable domains. In CD81 and other tetraspanins, the numerous membrane-exposed aromatic residues are asymmetrically clustered and protrude on one side of the transmembrane domain. This may represent a functional specialization of these two sides for interactions with cholesterol, proteins, or membrane microdomains. Abstract | Full Text | PDF (1088 kb)

• …more

Copyright © 2008 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 94, Issue 6, 1945-1953, 15 March 2008

doi:10.1529/biophysj.107.110395

Biophysical Theory and Modeling

Hydrogen-Bonding and Packing Features of Membrane Proteins: Functional Implications

Peter Werner Hildebrand^*, ,  , Stefan Günther^*, Andrean Goede^*, Lucy Forrest^†, Cornelius Frömmel^‡ and Robert Preissner^*

^* Institute of Molecular Biology and Bioinformatics, Charité University Medicine, Berlin, Germany
^† Howard Hughes Medical Institute, Department of Biochemistry and Molecular Biology and Center for Computational Biology and Bioinformatics, Columbia University, New York, New York 10032 USA
^‡ Medical Faculty, Georg-August-Universität Göttingen, Göttingen, Germany

Address reprint requests to Dr. Peter Werner Hildebrand, University Medicine Berlin, Institute of Biochemistry Charité, Monbijoustr. 1a, Berlin D-10117, Germany.

Membrane proteins are embedded into the hydrophobic environment of lipid bilayers where they modulate the exchange of information and mass between the different partitions of cells and tissues. This class of proteins includes highly biomedically and pharmaceutically relevant proteins such as G-protein coupled receptors, channels, and transporters. Most human membrane proteins are assemblies of hydrophobic transmembrane helices that bind coenzymes or ligands or, alternatively, form protein channels. This type of membrane channel can be distinguished from the β-stranded membrane proteins that form the rather unselective pores in membranes that originate from bacterial or mitochondrial outer membranes 1. These β-barrels are constructed from β-sheets in which the polar amino acids lining the pores and the hydrophobic residues that face the membrane are arranged in an alternating manner. By contrast, the amino acid composition of helical membrane proteins is primarily hydrophobic, with ∼90%–95% of the membrane-spanning residues being nonpolar 2.

Another effect of the lipid bilayer is to cause a weakening of the hydrophobic effect for the protein, which is the main driving force for the folding of water-soluble globular proteins 3. Consequently, it has been proposed that optimized van der Waals interactions between transmembrane helices could compensate for the lack of the hydrophobic effect and provide the driving force for membrane protein folding 4. However, recent analysis has revealed that helical membrane proteins are not more tightly packed than water-soluble proteins 5,6. Therefore, it remains an important outstanding question how the assembly of transmembrane helices is actually accomplished 7.

Stable transmembrane helices are formed as a result of regular patterns of hydrogen bonds between polar main-chain atoms 2,3. However, since the helical backbone is only partially shielded by the side chains 8, electrostatic interactions between polar backbone atoms are likely also to play an important role in the interactions between transmembrane helices. For example, packing interfaces containing the small amino acids Gly or Ala allow for the formation of Cα–H…O hydrogen bonds in helical membrane proteins 9. Since there is no entropic cost in burying the backbone atoms of small amino acids, such positions may even serve as initial points for the folding of helical membrane proteins. Such Cα–H…O bonds probably have a stabilizing effect of up to −1 kcal/mol 10 and accordingly could be nearly as strong as classical amide hydrogen bonds of globular proteins 11, where the strength of hydrogen bonds is diminished by the high dielectric effect in water and the competition from water for hydrogen bonds 7. Nevertheless, the role of Cα–H…O hydrogen bonds in the stability of the tertiary structure of helical membrane proteins has been somewhat controversial 10,12,13.

For classical hydrogen bonds, the free energy of formation is much higher within the hydrophobic core of the lipid bilayer, because of the low dielectric constant: estimates of this value range from −2 to −5kcalmol⁻¹ per bond 3,14. Consistent with this observation, it has been shown that a single polar amino acid such as Asn or Glu can drive homomeric association of model transmembrane peptides 15,16. However, such charged and polar amino acids are rare in helical membrane proteins, probably because of the free energy cost of desolvating polar side chains. Nevertheless, about half of all helix pairs in membrane proteins are stabilized by classical hydrogen bonds 5,17, and most of these are found in motifs involving medium polar residues such as Ser or Thr. Indeed, it has been shown that such Ser/Thr motifs can also drive the association of model transmembrane helices 18. Thus, classical hydrogen bonds are clearly important in helix-helix interactions in membrane proteins.

The strength of both classical and Cα–H…O hydrogen bonds depends on the distance, the chemistry of the donor and acceptor atoms, the relative arrangement of donor and acceptor atoms, and the nature of the surrounding milieu 19,20. As a consequence, even small conformational changes may result in the breaking of hydrogen bonds 13. Conformational changes such as the shifting, rotating, or tilting of helices are believed to occur frequently in some membrane proteins and are suggested to be facilitated by specific structural characteristics such as helix kinks, the smoothness of helical surfaces, or local packing defects 6,21,22. The analysis of the packing and the geometrical features of helix pairs that are also involved in hydrogen bonds could therefore provide clues to understanding the structure-function relationship of membrane proteins.

We have recently shown that the prevalent motif of helix-helix interactions in channels and transporters (collectively referred to as channels) is the right-handed motif, whereas this type of contact is clearly underrepresented in other membrane proteins (collectively referred to as membrane coils) 6. It is also known that multiple Cα–H…O hydrogen bonds tend to cluster at such right-handed contacts 9, which is consistent with the fact that channels are rich in sequence patterns containing small and medium polar amino acids 23. Moreover, channels are packed significantly less densely than membrane coils 6. These packing defects are likely to be functionally important, because they cluster at the proposed hinge regions of transporters or in the pores of channels. Together, these results suggest that the distribution of hydrogen bonds in membrane proteins may correlate (via specific architectural features such as handedness of helix-helix interactions or packing density) with the function of the protein.

To test this hypothesis we report here a comprehensive analysis of the abundance, types, and location of hydrogen bonds, including Cα–H…O hydrogen bonds, in known membrane protein structures. We first correlate differences in hydrogen-bond distributions with the two classes of membrane proteins: channels and membrane coils. In an earlier study, we observed that many crystal structures of membrane proteins include cavities that presumably contain water molecules not resolved in the electron density 6. Thus, we also inspect the regions around hydrogen bonds for large cavities to identify whether the presence of hydrogen bonds correlates with nearby putative buried waters. Finally, we show how the differences in helix-helix interactions (i.e., crossing angles) in the two protein classes correlate with the hydrogen-bonding patterns. More broadly, because of the importance of hydrogen bonding in the membrane interior, the detailed description of hydrogen bonds obtained from this analysis can be considered an important step toward a solution to the membrane protein folding problem 7.

We used a nonredundant data set of 27 high resolution membrane protein structures for the analysis of hydrogen bonds; these can be divided into 13 channels and 14 membrane coils structures (Table 1). To differentiate between those amino acids located within or outside the lipophilic environment, two parallel planes were drawn to delineate the hydrophobic part of the lipid bilayer, applying the criteria described in our previous analysis 2. Polar interactions were considered for this analysis only if both hydrogen-bonding partners are located between these two planes. The Ca²⁺-ATPase is the only member of the channel group for which the transport mechanism has been revealed by the structural elucidation of four different states 24. The ATP-bound (Protein Data Bank (PDB) code: 1vfp) form of the channel was taken for the statistics of hydrogen bonds. To estimate the number of hydrogen bonds preserved during transport, the Ca²⁺-bound and -unbound forms (1su4, 1iwo) and the intermediate state with hydrolyzed ATP (1wpg) were also studied 24.

The Cα–H…O and classical (N, O, C_β−ζ)–H…O hydrogen bonds were calculated using the programs TMCHbond and HBexplore, respectively 9,20. The hydrogen coordinates were generated according to standard geometrical rules in a first step, taking into account the hybridization of the donor atoms and the atomic environment 25. For the calculation of classical hydrogen bonds, only those interactions were considered with distances between the acceptor O and the donor H equal to or shorter than 2.5Å and with bonding angles ζ larger than 90°. We used the standard parameters recommended by the authors to calculate the Cα–H…O bonds with TMCHbond 9. That is, the maximum distance between the acceptor O and the donor H was 3.5Å with bond angles ζ>120° or 3.0Å with ζ>90°.

For the calculation of packing densities, the Voronoi cell method was applied 26. This method uses curved instead of planar interfaces between atoms to create more reasonable assignments. Another advantage of this method is that it also works for atoms located in protein regions with large packing defects, which are frequently found in channels 6. Hydrogen bonds were referred to as “in close proximity to” such an internal cavity if at least one atom lying within a radius of 6Å around the center of the bond is in contact with the cavity. To evaluate the atomic packing densities, all buried atoms of two neighboring transmembrane helices with their van der Waals radii closer than 1.5Å were considered. Atoms are classified as buried when less than 40% of the atoms’ surface contacts the surrounding milieu 6.

For the evaluation of crossing angles, we first identified pairs of helices with at least two residues with atoms less than 1.5Å from one another. Transmembrane helices are regularly kinked or curved 2. The crossing angle was therefore determined using the local axes of the two helical sections that are actually involved in the contact after extending each helical section by an additional turn at both termini. Thus helices are only considered in their entirety for helix pairs with virtually parallel arrangements.

The Shannon entropy H is a well-established measure to estimate the conservation of a certain amino acid position during evolution 27. It was taken from the corresponding position of the multiple sequence alignment of the PDB-associated HSSP file 28. The absolute value of the Shannon entropy depends on the average similarity of all the sequences used for the alignment. Thus, only homologous proteins with a sequence similarity above 30% were considered. To estimate the variability of a certain amino acid, a reasonable amount of distantly related proteins is required for the alignment. Protein alignments with an average pair-wise sequence identity above 50% (1ppj and 2occ) were not taken into account. The entropy of residues exposed to lipid, H_exposed, is generally high in all helical membrane proteins 29,30,31. The average value of H_exposed over all residues in each protein was used as a reference to facilitate comparison across different proteins. Thus, for an amino acid, a, the conservation of the buried amino acids is defined as

Analysis of interhelical hydrogen bonding in membrane proteins reveals significant differences between channels and membrane coils. Specifically, in channel proteins, there are 0.29 hydrogen bonds per amino acid, whereas residues in membrane coils are stabilized on average by 0.19 hydrogen bonds (Table 1). Hence, channels are stabilized by 50% more hydrogen bonds than membrane coils. This statistically highly significant difference (t-test: 99.9% probability) is mainly caused by the higher number of Cα–H…O hydrogen bonds in channels. Almost two times as many Cα–H…O hydrogen bonds per residue are formed in channels (0.14) compared to membrane coils (0.08). In addition, the number of classical hydrogen bonds per residue is slightly increased in channels (0.15 vs. 0.12).

Another way to interpret these results is to say that approximately every third amino acid in channels and every fifth amino acid in membrane coils are involved in interhelical hydrogen bonding. Since approximately half of the residues of the investigated membrane proteins are exposed to lipid 23, this would mean that, on average, approximately every second buried amino acid is hydrogen bonded. It is important to note that we considered only interactions with both donor and acceptor atoms located inside the hydrophobic part of the lipid bilayer. In that region, the value of the dielectric constant is low and electrostatic interactions are of considerable strength 3. Thus, hydrogen bonding appears to play a key role in the stabilization of the tertiary structure of helical membrane proteins.

To estimate the strengths of the hydrogen bonds in the two data sets, the distances between acceptor and donor atoms were compared. At longer distances (2.6Å) the number of classical hydrogen bonds begins to decrease in both data sets (Fig. 1). This distribution is also characteristic of hydrogen bonding in water-soluble globular proteins 19,32, suggesting that the number of classical hydrogen bonds does not depend on the surrounding milieu 5.

Display large version of this figure

Figure 1

Distribution (%) of different lengths (Å) centered at each value of classical (solid) and Cα–H…O (hatched) hydrogen bonds in channels (black columns) and membrane coils (white columns).

At distances shorter than 2.0Å, we observe that the fraction of all hydrogen-bond types is higher in membrane coils than in channels. Conversely, the fraction of hydrogen bonds longer than 2.0Å is considerably higher in channels. This trend is enhanced by the higher occurrence of Cα–H…O hydrogen bonds in channels, also starting at bond lengths of ∼2.0Å. As a consequence, the fraction of helix-helix contacts stabilized by the shortest—and hence the strongest—hydrogen bonds is considerably lower in channels than in membrane coils (Fig. 1).

The analysis of interhelical hydrogen bonding of membrane proteins therefore suggests that channels are enriched in hydrogen bonds, and particularly in bonds with weak to medium strength.

The composition of residues forming interhelical hydrogen bonds varies significantly between channels and membrane coils. In channels, the smallest amino acids (Gly and Ala) form the majority of the interhelical hydrogen bonds (Fig. 2), whereas in membrane coils the medium polar amino acids Thr and Ser form most bonds. This difference is caused by the higher content of these residues in the respective protein types 2 as well as the higher likelihood that Gly will form Cα–H…O hydrogen bonds in channels (t-test: 95% probability). Specifically, 37% of Gly residues are involved in Cα–H…O hydrogen bonds in these proteins. By contrast, in membrane coils only 22% of the Gly residues contribute to Cα–H…O hydrogen bonds. Although the statistical significance is not very high, these results suggest that Gly may be of particular importance for stabilizing channels via hydrogen bonding.

Display large version of this figure

Figure 2

Proportion (%) of amino acid types involved in classical (solid) and Cα–H…O (hatched) hydrogen bonds in channels (black columns; total=342) and in membrane coils (white columns; total=358).

The variability of a certain amino acid position during evolution indicates whether it plays a vital role for the protein in terms of folding, stability, or function. The Shannon entropy, calculated from a sample of homologous protein sequences, is a well-established measure for this variability. Indeed, since exposed residues are generally less conserved than buried residues, this measure has been used to effectively predict buried versus exposed residues 30,31,33. Since the Shannon entropy largely depends on the quality of each alignment, the differential to an internal standard has to be used to compare the residual variability of different proteins (see Materials and Methods). Using this measure we found that buried residues in helical membrane proteins are conserved up to 25% more than exposed residues (Fig. 3). This is in agreement with previous analyses 29,30,31.

Display large version of this figure

Figure 3

Relative decrease in Shannon entropy of buried or hydrogen-bonded residues versus exposed residues in channels (black columns) and membrane coils (white columns). Error bars signify the 95% confidence levels.

When considering the Shannon entropy of residues that form hydrogen bonds, we find that in both protein types these residues are noticeably more conserved than buried residues as a whole (Fig. 3). In fact, hydrogen-bonding residues are ∼40% more conserved than exposed residues. We previously observed similar results in an analysis of classical hydrogen bonds in membrane proteins 5. The impact of Cα–H…O hydrogen bonds on the stability of helical membrane proteins, however, is somewhat more controversial 12,13. Therefore, it is important to note that, in both protein types, amino acids involved in this type of hydrogen bond are significantly more conserved (t-test: 99% probability) than buried residues in general. Indeed, in channels, residues taking part in Cα–H···O hydrogen bonding are slightly more conserved even than residues involved in classical hydrogen bonds (Fig. 3). The analysis of the per-residue variability, therefore, supports the notion that Cα–H…O hydrogen bonds are important for the stabilization and function of helical membrane proteins 9.

It has previously been stated that the presence of hydrogen bonds correlates with higher density helix-helix contacts in all types of membrane proteins 17. This finding is only partially supported by our analysis. In membrane coils (which account for 10 of the 13 membrane protein structures considered in the analysis of Adamian and Liang 17), this increase is indeed highly significant for those helix pairs that are stabilized by multiple hydrogen bonds (Figure 4AB). By contrast, in channels the packing density of hydrogen-bonded helix-helix contacts is only marginally increased compared to non-hydrogen-bonded pairs. We have pointed out recently that helix-helix contacts in membrane coils are packed more densely than in channels 6. This observation also holds for hydrogen-bonded helix pairs of membrane coils, which are packed significantly (t-test: 99.5% probability) more densely than those of channels.

Display large version of this figure

Figure 4

Comparison of packing densities of non-hydrogen-bonded and hydrogen-bonded helix pairs (A, classical; B, Cα–H…O; AB, both types) in channels (black columns) and membrane coils (white columns). Error bars signify the 99% confidence levels.

The lower packing density in channels may be due to the presence of cavities that contain internal waters 6. Since the stability of interhelical hydrogen bonds does not depend only on geometry but also on the chemical environment 3, we investigated the portion of hydrogen bonds that is located in close proximity to protein voids or pockets (see Materials and Methods). In membrane coils, only about one quarter (28%) of the hydrogen bonds are located closer than 6Å to such cavities (Table 1). By contrast, one half (51%) of the hydrogen bonds of channels are found in close proximity to a cavity. The breaking of interhelical hydrogen bonds in channels during their function could therefore be facilitated by interactions with water molecules in nearby internal cavities.

In the Ca²⁺-ATPase, we found that 59% of the hydrogen bonds are in close proximity to water-sized cavities (Table 1, ). We detected a total of 130 different hydrogen bonds that are formed during the reaction cycle (data not shown). The only two classical and six Cα–H…O hydrogen bonds that are preserved during the entire reaction cycle are located in the anchor domain, in which the helices do not change conformation 24. Furthermore, water-sized cavities in the anchor domain are practically absent (only 1 of 38 cavities of the Ca²⁺-bound state is located there). Accordingly, in the Ca²⁺-ATPase, there seems to be a correlation between the presence of water-sized cavities and the ability of interhelical hydrogen bonds to break.

An analysis of helix crossing angles reveals that 69% of the hydrogen-bonded helix pairs of channels cross at right-handed angles (Table 1). This is significantly more than in membrane coils, where only 27% of the hydrogen-bonded helix pairs cross at right-handed angles (χ²: 99.9% probability). Another striking difference is that the helix crossing angles of hydrogen-bonded helices cluster around |τ|≈40° for channels and |τ|≈20° for membrane coils (Fig. 5). The preference for certain packing modes in helical membrane proteins has also been observed by others 5,34,35. However, the subdivision of the membrane proteins into channels and membrane coils shows that these differences are correlated with protein function.

Display large version of this figure

Figure 5

Proportion (%) of different helix crossing angles (τ) stabilized by classical (solid) and Cα–H…O (hatched) hydrogen bonds in channels (black columns) and membrane coils (white columns).

The right-handed and left-handed packing motifs prevalent in channels and membrane coils match the class a and class c “knobs-into-holes” packing motifs, respectively, derived from a helical lattice superposition model 36. Walther and co-workers stated in that work that there should be greater flexibility in class a helix packing since small amino acids are frequently involved in the packing core. Consistent with that hypothesis, the distribution of helix-helix crossing angles for the antiparallel right-handed (class a) helix pairs in channels is very broad (Fig. 5). Moreover, hydrogen-bonded helix pairs with |τ|≈40° are preferred in both parallel and antiparallel helix-helix interactions of channels. By contrast, antiparallel (class c) helix-helix interactions are strongly preferred in membrane coils for which the distribution of helix crossing angles is extremely narrow.

The two different types of helix-helix interactions characteristic for channels and membrane coils are illustrated schematically in Fig. 6. The class c packing 36 of helices in membrane coils (|τ|≈20°) results in a parallel arrangement of the helix turns and the clear preference for antiparallel helix packing motifs, where the bulky side chains of one helix interlock between two bulky residues of the other helix (Fig. 7). In the class a packing 36 of helices in channels (|τ|≈40°), the turns of the contacting helices are positioned nearly orthogonal to each other, where the Cα-Cβ vectors of the buried residues point away from the packing core. The preferred burial of small residues at helix-helix interfaces of channels 23,37 together with the geometry of the packing of side chains in right-handed contacts (Fig. 7) allows the helices to arrange in both parallel and antiparallel motifs.

Display large version of this figure

Figure 6

Schematic representation of the left-handed knobs-into-hole packing of membrane coils (left) with parallel arrangement of the helix turns and the right-handed knobs-into-holes packing of channels (right) with orthogonal arrangement of the helix turns. The parallel arrangement of the helix turns and the resultant arrangements of side chains (Fig. 7) lead to a preference for antiparallel interactions in membrane coils. The dotted lines indicate the angle of the turn on the rear of the helices.

Display large version of this figure

Figure 7

Schematic representation of the knobs-into-holes packing of an antiparallel left-handed helix pair of membrane coils (PDB code: 1c3w) (left) and of a right-handed parallel helix pair of channels (1kpl) (right). The Cα-Cβ vectors of the core residues (in magenta) in right-handed interactions point away from the contact, where small residues are again preferred. The knobs are therefore articulately flattened, and the sterical restrictions imposed on the conformation of these helix pairs are limited. For the interdigitations of side chains, the conformation of left-handed helix pairs is sterically constrained and antiparallel arrangements are strongly preferred.

Finally, there is a clear preference for parallel class a packing (τ=−40°±10°) and antiparallel class c packing (τ=−160°±10°) in hydrogen-bonded helix pairs of channels and membrane coils, respectively, compared to non-hydrogen-bonded helix pairs. Exactly 50% of all hydrogen-bonded helix pairs of membrane coils are found in the antiparallel class c packing motif (Fig. 5). This is two times as many as in non-hydrogen-bonded helix pairs, where the same motif accounts for only 22% of the pairs (data not shown). In channels, three times as many (or 23%) of the hydrogen-bonded helix pairs are found in the parallel class a packing motif than in non-hydrogen-bonded helix pairs (8%). Our results therefore highlight that the different types of helix packing in membrane proteins correlate with the presence of hydrogen bonds.

Using a comprehensive analysis of geometrical features and hydrogen bonding, we found that there are significant differences between the two classes of membrane proteins studied here. In terms of helix packing, as shown also in our previous work 23, it is clear that helices in channels form right-handed contacts, whereas helices in membrane coils tend to form left-handed crossing angles. Interestingly, DeGrado and co-workers recently illustrated that if helix pairs of membrane proteins are categorized according to their structural similarity 35, clusters very similar to those obtained from the classifications in channels and membrane coils are achieved (Fig. 5). That is, the antiparallel left-handed helix pairs (typical for membrane coils) and the parallel and antiparallel right-handed helix pairs (typical for channels) form the biggest clusters of structurally similar helix pairs. The subdivisions of helical membrane proteins according to structural 35 or functional (Fig. 5) characteristics, therefore, refer to the same building blocks.

Here, we have expanded on this analysis to show how these different helix-packing arrangements may contribute to the function of helical membrane proteins. In particular, in right-handed packing arrangements the helices adopt a much broader range of crossing angles than in the left-handed packing motifs (Fig. 5). Moreover, due to the nature of the side-chain packing, these right-handed helix pairs—which predominate in channels—are not interlocked the same way as the left-handed helix pairs of membrane coils (Fig. 7). It therefore seems that there is a higher conformational freedom in the right-handed helix pairs that are clearly overrepresented in channels. Since the function of channels frequently involves large structural rearrangements, we propose that this higher conformational freedom is important for their function.

The low packing density of the helix pairs of channels (Fig. 4) may also facilitate greater conformational freedom. It has recently been proposed that the energetic barrier in separating the surfaces of two hydrophobic helices may be very high 38. This separation involves energetically highly unfavorable transition states, as void volumes form and steric dewetting occurs 39,40. This result suggests that a “tight” packing of transmembrane helices 4, i.e., via optimized van der Waals contacts, would therefore lead to rather rigid contacts. This energetic barrier, however, does not exist for helix pairs separated by individual water molecules 38. We have previously shown that solvated helix pairs are typically found in channels 6. It is therefore possible that the water that accumulates at the functionally important regions of channels is necessary to facilitate the required conformational switching.

One of the most significant findings of this analysis is that the different types of helix packing in membrane proteins correlate with the presence of hydrogen bonds of specific types (Fig. 5). This is notable since hydrogen bonds are particularly strong in the hydrophobic environment of the lipid bilayer. Specifically, in right-handed crossings the “knobs” of one helix that pack into the “holes” of the other helix are typically flattened, resulting in a close packing of backbones (Fig. 7). Channels are therefore clearly enriched in Cα–H…O hydrogen bonds (Table 1). Moreover, helix pairs stabilized by hydrogen bonds are also found to arrange in “knobs-onto-knobs” packing patterns 36 in which the side chains of adjacent helices are typically not interlocked. Such a motif where Gly and Ala residues of contacting helices are directly opposed is vital for the assembly and the opening mechanism of the mechanosensitive channel of small conductance 41. Indeed, the conservation of residues involved in hydrogen bonds (Fig. 3) argues for a central role of hydrogen bonds for the folding, stabilization, and function of helical membrane proteins in general 9,17,42.

Hydrogen bonds are frequently found close to water-sized cavities (Table 1). By providing alternative hydrogen-bonding partners, the nearby water may be important to stabilize transition states in which the hydrogen-bonding network is significantly altered, as shown here for the Ca²⁺-ATPase. The requirement for easy breaking and reformation of hydrogen bonds is also consistent with the fact that channels are enriched in hydrogen bonds of medium strengths, such as Cα–H…O hydrogen bonds, rather than the short, and thus presumably strong, hydrogen bonds found in membrane coils (Fig. 1). Accordingly, by the use of hydrogen bonds, channels may achieve a high flexibility while preserving energetic stability.

Finally, in addition to providing a functional explanation for the presence of specific helix packing and hydrogen-bonding types, the results provide new approaches for the prediction of the structure and function of helical membrane proteins. For example, the observation that certain helix-packing architectures are observed more frequently reduces the conformational search space required for low-resolution modeling. Furthermore, it will be interesting to see whether such helix-packing arrangements can be predicted via the detection of specific sequence patterns. Alternatively, the identification of water-containing cavities and particular structure motifs in proteins of known structure may enable identification of regions that are involved in conformational changes.

We have described several basic geometrical features characteristic of helical membrane proteins with different functions. Hydrogen bonds are clearly overrepresented in the structural motifs that are common to channels and tend to originate from highly conserved residues. The fact that channels are significantly enriched in Cα–H…O hydrogen bonds is consistent with the proposal that the enrichment of short interhelical distances between transmembrane helices (due to an enrichment of small amino acid residue types) may be important for energetic stabilization of membrane proteins 5,6,9. The presumed weakness of the hydrogen bonds in channels, along with the prevalence of nearby cavities that may contain water molecules, together support a hypothesis that breaking and reformation of different networks of hydrogen bonds may facilitate conformational changes in channels. As more structures of channels in different functional states become available, it should become clear whether this hypothesis holds as a general rule.