| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Institute of Biochemistry Charité, University Medicine Berlin, 10117 Berlin, Germany
Correspondence: Address reprint requests to Peter Werner Hildebrand, Tel.: 49-30-450-528212; Fax: 49-30-450-528942; E-mail: peter.hildebrand{at}charite.de.
| ABSTRACT |
|---|
|
|
|---|
-helical membrane proteins. Using a modified Voronoi procedure, we calculated packing densities for helix-helix contacts in membrane spanning domains. Our results show that the transmembrane helices of protein channels and transporters are significantly more loosely packed compared with helices in globular proteins. The observed packing deficiencies of these membrane proteins are also reflected by a higher amount of cavities at functionally important sites. The cavities positioned along the gated pores of membrane channels and transporters are noticeably lined by polar amino acids that should be exposed to the aqueous medium when the protein is in the open state. In contrast, nonpolar amino acids surround the cavities in those protein regions where large rearrangements are supposed to take place, as near the hinge regions of transporters or at restriction sites of protein channels. We presume that the observed deficiencies of helix-helix packing are essential for the helical mobility that sustains the function of many membrane protein channels and transporters. | INTRODUCTION |
|---|
|
|
|---|
In membrane proteins, however, a somewhat different picture might be expected (MacKenzie and Engelman, 1998
; White and Wimley, 1999
). The distinctive surrounding milieu causes the weakening of the hydrophobic effect inside the protein within the lipid bilayer. In addition, the inclusion of prosthetic groups and the pervasion of ion channels and solute transporters with water-filled pores could result in deviations in molecular packing. Furthermore, it has been proposed that ion channels and solute transporters open throughout gating mechanisms that require broad molecular rearrangements of their transmembrane domains (Locher et al., 2003
; Perozo et al., 2002
; Swartz, 2004
). Thus the question arises whether this flexibility is facilitated by a reduced packing density.
Nevertheless, the validation of atomic packing densities in proteins is a difficult problem (Fleming and Richards, 2000
). Eilers and colleagues, for example, applied the occluded surface method (Pattabiraman et al., 1995
) to assess the packing densities of transmembrane helices by investigating 11 helical membrane protein structures (Eilers et al., 2002
). They presume, in conclusion, that helical membrane proteins are generally packed more densely than other proteins. This conclusion is surprising, because the hydrophobic effect as a driving force for helix-helix interaction is absent inside the lipid bilayer (MacKenzie and Engelman, 1998
). In addition, polar- or hydrogen-bonded interactions generally occur less frequently than in water-soluble globular proteins (DeGrado et al., 2003
). Leaving energetic considerations aside, the dense packing of transmembrane helices appears contradictory to the mobility that is thought to form the basis of the proper functioning of many membrane proteins (Jiang et al., 2002
).
To solve this conundrum, we used a modified Voronoi procedure to calculate the packing values and to evaluate the number of cavities found in the different groups of proteins (Goede et al., 1997
; Rother et al., 2003
). Hereby we concentrated in those atoms buried in the helix-helix interfaces leaving surface atoms aside. The advantage of our method is that it also works for atoms that reside in those protein regions, where large packing deficiencies occur, as in protein voids or pockets (cavities). As pointed out in a mathematical appendix, the occluded surface method fails to give exact packing values in those cases. The analysis of an updated set of 20 helical membrane proteins finally leads us to a point of view that is contrary to the proposed opinion that helical membrane proteins are generally packed more densely than other helical proteins. We show clearly that packing deficiencies often occur at functionally important protein regions. To comprehend the molecular base of the observed packing differences between helices in membrane and globular proteins, detailed statistics about the amino acid and atomic composition of helix-helix patches and cavities were collected.
| MATERIAL AND METHODS |
|---|
|
|
|---|
Data sets
We have recently shown that, according to their different functions and architecture, membrane proteins can be subdivided into two groups (Hildebrand et al., 2004
): membrane channels and transporters can be subsumed under the term "membrane gates", because of the regulated opening or closing of the pore underlying their function. On the other hand, proton pumps, receptors, and photosystems can be classified as "membrane coils", because the helix interaction motifs are very similar to those characteristic for soluble coiled-coils (Langosch and Heringa, 1998
). Transmembrane helices in membrane coils predominantly cross at left-handed angles, those in membrane gates predominantly at right-handed angles. We further demonstrated that the members of these two groups differ markedly in the torsion angles and intrahelical hydrogen bonds (Hildebrand et al., 2004
).
To establish whether the particular functionality is also reflected in different packing density values, the current data set of 20 high-resolution structures of nonrelated helical membrane proteins was specified into 10 membrane gates and 10 membrane coils (Table 1). All the calculations were carried out on monomers except for channels and transporters, where the functional multimers were used. We chose wild-type structures with highest resolution. Nonetheless, the structures of membrane coils are resolved at higher mean resolution than membrane gates. To qualify whether this could influence our analysis, the packing density values of 24 different structures of the photosynthetic reaction center was compared. As a result, we observed that the packing density (mean value = 0.81, standard deviation
= 0.006) does not depend on the crystallographic resolution (2.1 Å3.5 Å).
|
Packing densities
The local packing density of an atom is defined as the fraction PDloc = VVdW/Vse, where VVdW is the space inside the atoms' van der Waals sphere and Vse is the van der Waals sphere expanded by 1.4 Å, the radius of a solvent molecule. When two atoms are bonded or close to each other, parts of their spheres are cut off by separating panes, with impact on the local packing values. The amount of space remaining in each sphere is referred to as the atomic volume.
To calculate atomic volumes, the Voronoi cell procedure using hyperboloid interfaces between atoms was applied (Goede et al., 1997
). The interfaces between atoms are defined by all points that are equally distant from both van der Waals spheres of the two atoms. Atomic volumes were calculated numerically using a cubic lattice with a grid distance of 0.2 Å. We used the Stouten set of atom radii for volume calculations (Stouten et al., 1993
). Ligands and prosthetic groups were kept in the structure, but all lipid molecules found at the protein surfaces were removed.
For each protein structure, the average packing density PD was calculated as the mean value of the local packing densities of all contributing atoms having the standard deviation
and the variance
2. Atoms at the protein surface ("surface") are usually packed far less efficiently and were excluded from analysis. Surface atoms were defined as all atoms touched by a 1.4Å radius probe rolling along the VdW spheres at the whole protein surface. Thus, for the analysis of packing densities, solely "buried" atoms were considered.
Analyzing cavities
In protein structures, there are locations, where a 1.4 Å radius probe fits without intersecting any atoms' van der Waals sphere and without reaching the surface (Richmond, 1984
). These locations are defined as cavities. All buried atoms touched by the 1.4 Å probe were grouped as cavity neighbors. For this part of our analysis, the water molecules were removed. The relative amount of neighboring atoms of the cavities was estimated. Finally, the frequency of polar and nonpolar amino acids that are in contact with cavities was calculated, using the hydrophobicity scale of Eisenberg et al. (1984)
. Cavities enclosed by a fraction of nonpolar amino acids above average are termed "polar cavities"; those below the average "nonpolar cavities".
Molecular surface patches
Helical interfaces are defined as pairs of molecular surface patches between neighboring helices that are in direct contact with each other (Preissner et al., 1998
). Direct contact means that the atomic van der Waals surfaces are closer than a given cut-off distance. It was previously shown that molecular surface patches are generally flat atom assemblies with a length/width/depth ratio of
3:2:1. Nevertheless, the size and atomic composition of a molecular surface patch varies with the chosen cut-off (see Fig. 2). All packing values that will be mentioned were calculated at a medium cut-off value of 1.5 Å. At any rate, atoms being labeled as surface had been removed.
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
within the lipid bilayer (Chothia, 1975The comparison of the atomic packing density values reveals further details: The side-chain atoms buried between helices (PD = 0.77) are significantly (t-test, p-value of 0.005) less efficiently packed in the data set of membrane gates than the main-chain atoms (PD = 0.80). This is especially valid for the polar side-chain atoms (PD = 0.64). This difference is much smaller in the data sets of membrane coils and helices of globular proteins (Table 3). When a higher cut-off value was chosen, the helix patches expand in all three dimensions and the fraction of main-chain atoms at the interfaces increases (Fig. 2). Because the backbone atoms are packed more efficiently, the packing density of the investigated helix-helix interfaces accordingly rises (Fig. 3). This increase is plainest in the data set of membrane gates, because the side-chain atoms are packed most fluffily here.
|
|
It has been noted that the tight backbone-to-backbone packing is important for the stabilization of the tertiary structure of membrane proteins. This structural signature was successfully applied to the prediction of their tertiary structure (Fleishman and Ben-Tal, 2002
; Liu et al., 2004
). Accordingly, the C
-H
O hydrogen bond is a strong determinant of stability and specificity in transmembrane helix interactions (Senes et al., 2001
). Interestingly, multiple hydrogen bonds of this type are predominantly found between parallel transmembrane helices that cross at right-handed angles. In this regard, membrane gates again differ markedly from the remaining helical membrane proteins: Two-thirds of the helix crossings are right-handed, whereas only one-third of the transmembrane helices in membrane coils cross right-handed (Hildebrand et al., 2004
). Our approach thus supports the proposition that helices in membrane channels and transporters are stabilized throughout close approximations of their backbones, which again facilitates the formation of multiple interhelical hydrogen bonds, stabilizing the topology of the helices. The helical mobility that is necessary for the proper functioning of membrane channels and transporters (Swartz, 2004
) could again be realized by the loose packing of their side-chain atoms.
Finally, the question has to be addressed as to whether the observed packing deficiencies found in the three data sets are equally dispersed over the entire protein structure, or whether they cluster in certain protein regions (cavities). In the first case, the observed packing deficiencies could be traced back to intrinsic factors, such as the different values for the dielectric constant, the specific amino acid composition, or the different quality of the protein structures (Fleming and Richards, 2000
; Vaguine et al., 1999
). In the second case, the packing defects probably point toward functionally important sites of the protein. Subsequently, all atoms that are in direct contact with cavities were removed. As a result, the packing densities of the different data sets converge articulately (Fig. 1). This strongly favors the hypothesis that the packing deficiencies detected in membrane gates might be functionally important.
Cavities
Cavities are empty or solvent-filled spaces within a protein structure that are large enough to accommodate at least one water molecule. When the atoms exposed to cavities were excluded, the actual differences in packing densities are rather small (Fig. 1). Thus the cavities account for the reduced packing of membrane gates. The density distribution of membrane gates is therefore reflected by the localization and characterization of internal cavities.
Transmembrane domains of membrane coils resemble nearly the same content of internal cavities as helices of globular proteins (5.2% and 5.5% of the buried atoms are in direct contact with a cavity). By contrast, slightly more cavities were found in the data set of membrane gates (6.4%, Table 1). Most cavities are neighbored by 1525 atoms (
31%) or <5 (57%) atoms. The results therefore indicate that the size of cavities found in proteins depends neither on the surrounding milieu, nor on the fold of the proteins (Rother et al., 2003
).
Cavities in membrane gates are highly polar in contrast to cavities in membrane coils (Fig. 4). In globular proteins, most polar cavities are found close to the protein surface (Rother et al., 2003
). Unlike the helices in water-soluble globular proteins, transmembrane helices are encircled by hydrophobic lipid tails (see Material and Methods). Thus polar cavities within membrane proteins are positioned close to those protein regions that communicate with the extra- and intracellular environments, as the pores that pervade these proteins (Fig. 5). Therefore it seems feasible that the cavities of membrane gates are encircled by more polar amino acids than the cavities of membrane coils, and that many of them might be filled with water molecules.
|
|
However, it is assumed that deeply buried cavities contribute only little to this potential error, because they are hydrophobic, indicating that they are devoid of water (Rother et al., 2003
). This is supported by molecular dynamic simulations revealing that nonpolar cavities in the protein interior are even more compressible than other parts of the protein (Kocher et al., 1996
). Ongoing, we detected nonpolar cavities in protein regions that are supposed to be functionally important, too: around the restriction region and the so called "plug" of the translocon (1rh5, Fig. 6), and around the hinge regions of the structure of glycerol-3-phosphate transporter (1pw4, Fig. 5). We therefore conclude that focal packing defects allow for the structural flexibility that is required for the proper functioning of membrane channels and transporters.
|
| CONCLUSION |
|---|
|
|
|---|
The loose packing of membrane channels and transporters is in turn mainly caused by the frequent placement of polar side chains close to voids or pockets that are lined along the pores that pervade these proteins (Fig. 1). It is inferred by focal packing defects, i.e., cavities, rather than by steadily increased distances between midpoints of atoms. Nevertheless, if the proper functioning of membrane channels and transporters is based on a relatively loose packing of helical interfaces, how are these proteins stabilized compared to other membrane proteins? The comparatively close packing of the transmembrane backbones indicates that main-chain interactions probably compensate for the loose packing of side chains. The investigation of the packing density of transmembrane helices finally leads to the detection of structural details that promote a better understanding of the relation between stabilization and function of membrane proteins.
| APPENDIX: ESTIMATION OF THE MOLECULAR PACKING IN DEPENDENCE ON THE METHOD EMPLOYED AND THE DISTANCE TO THE CENTER OF A CAVITY |
|---|
|
|
|---|
from the center of the cavity (Fig. 7). If the atomic packing density (PD) is estimated by the Voronoi cell method, the situation will be the same regardless of whether the atom lies in the middle or more on the periphery of a cavity. If the occluded surface method (PV) is used, different values will be obtained. When calculating the packing density the following integral has to be solved:
![]() |
|
x1(
)) is the integral around the circle with radius x1, orthogonal to the x axis, and d(
) is the distance between the two spheres, measured from the center of the atom. The following equations must be combined, when the function is integrated:
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
Inserting Eqs. 3 and 4 into Eq. 5 results in
![]() |
![]() |
![]() |
This last quadratic equation can be solved by:
![]() |
![]() |
Because
should not be negative, we get
![]() |
The integral becomes
![]() |
For the packing value, we have to divide through the surface area of the atom 4
r2 and the maximal length 2.8 Å, getting
![]() |
Although the packing value PV (occluded surface method) estimates the mean distance to the next atoms, the packing density PK (Voronoi cell) rates the nonoccupied space around the atoms. This space increases with the power of 3 in relation to the distance to the center of the atom, whereas this distance has only linear influence on the packing value PV. Therefore, especially near cavities, the packing value PV will be larger compared to other cases with the same nonoccupied space around the atom and therefore the same packing density (Table 4). Because cavities are particularly frequent in protein channels and transporters, the molecular packing density of these proteins will be overestimated when calculated with the occluded surface method (PV).
|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported with grants from the Federal Ministry of Education and Research of Germany (BMBF), Berlin Center for Genome Based Bioinformatics (BCB), and the Charité.
Submitted on July 13, 2004; accepted for publication November 1, 2004.
| REFERENCES |
|---|
|
|
|---|
Chothia, C. 1975. Structural invariants in protein folding. Nature. 254:304308.[CrossRef][Medline]
Chothia, C., M. Levitt, and D. Richardson. 1981. Helix to helix packing in proteins. J. Mol. Biol. 145:215250.[CrossRef][Medline]
DeGrado, W. F., H. Gratkowski, and J. D. Lear. 2003. How do helix-helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo-oligomeric helical bundles. Protein Sci. 12:647665.
Eilers, M., A. B. Patel, W. Liu, and S. O. Smith. 2002. Comparison of helix interactions in membrane and soluble
-bundle proteins. Biophys. J. 82:27202736.
Eisenberg, D., R. M. Weiss, and T. C. Terwilliger. 1984. The hydrophobic moment detects periodicity in protein hydrophobicity. Proc. Natl. Acad. Sci. USA. 81:140144.
Fleishman, S. J., and N. Ben-Tal. 2002. A novel scoring function for predicting the conformations of tightly packed pairs of transmembrane
-helices. J. Mol. Biol. 321:363378.[CrossRef][Medline]
Fleming, P. J., and F. M. Richards. 2000. Protein packing: dependence on protein size, secondary structure and amino acid composition. J. Mol. Biol. 299:487498.[CrossRef][Medline]
Gille, C., A. Goede, R. Preissner, K. Rother, and C. Frömmel. 2000. Conservation of substructures in proteins: interfaces of secondary structural elements in proteasomal subunits. J. Mol. Biol. 299:11471154.[CrossRef][Medline]
Goede, A., R. Preissner, and C. Frömmel. 1997. Voronoi cell: new method for allocation of space among atoms: elimination of avoidable errors in calculation of atomic volume and density. J. Comput. Chem. 18:11131123.[CrossRef]
Hildebrand, P. W., R. Preissner, and C. Frömmel. 2004. Structural features of transmembrane helices. FEBS Lett. 559:145151.[CrossRef][Medline]
Huang, Y., M. J. Lemieux, J. Song, M. Auer, and D. N. Wang. 2003. Structure and mechanism of the glycerol-3-phosphate transporter from Escherichia coli. Science. 301:616620.
Jiang, Y., A. Lee, J. Chen, M. Cadene, B. T. Chait, and R. MacKinnon. 2002. The open pore conformation of potassium channels. Nature. 417:523526.[CrossRef][Medline]
Kim, S., and T. A. Cross. 2002. Uniformity, ideality, and hydrogen bonds in transmembrane
-helices. Biophys. J. 83:20842095.
Kocher, J. P., M. Prevost, S. J. Wodak, and B. Lee. 1996. Properties of the protein matrix revealed by the free energy of cavity formation. Structure. 4:15171529.[Medline]
Langosch, D., and J. Heringa. 1998. Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils. Proteins. 31:150159.[CrossRef][Medline]
Lemmon, M. A., and D. M. Engelman. 1994. Specificity and promiscuity in membrane helix interactions. FEBS Lett. 346:1720.[CrossRef][Medline]
Liu, W., M. Eilers, A. B. Patel, and S. O. Smith. 2004. Helix packing moments reveal diversity and conservation in membrane protein structure. J. Mol. Biol. 337:713729.[CrossRef][Medline]
Locher, K. P., R. B. Bass, and D. C. Rees. 2003. Structural biology. Breaching the barrier. Science. 301:603604.
MacKenzie, K. R., and D. M. Engelman. 1998. Structure-based prediction of the stability of transmembrane helix-helix interactions: the sequence dependence of glycophorin A dimerization. Proc. Natl. Acad. Sci. USA. 95:35833590.
Olivella, M., X. Deupi, C. Govaerts, and L. Pardo. 2002. Influence of the environment in the conformation of
-helices studied by protein database search and molecular dynamics simulations. Biophys. J. 82:32073213.
Pattabiraman, N., K. B. Ward, and P. J. Fleming. 1995. Occluded molecular surface: analysis of protein packing. J. Mol. Recognit. 8:334344.[CrossRef][Medline]
Perozo, E., D. M. Cortes, P. Sompornpisut, A. Kloda, and B. Martinac. 2002. Open channel structure of MscL and the gating mechanism of mechanosensitive channels. Nature. 418:942948.[CrossRef][Medline]
Popot, J. L., and D. M. Engelman. 2000. Helical membrane protein folding, stability, and evolution. Annu. Rev. Biochem. 69:881922.[CrossRef][Medline]
Preissner, R., A. Goede, and C. Frommel. 1999. Spare parts for helix-helix interaction. Protein Eng. 12:825832.
Preissner, R., A. Goede, and C. Frömmel. 1998. Dictionary of interfaces in proteins (DIP). Data bank of complementary molecular surface patches. J. Mol. Biol. 280:535550.[CrossRef][Medline]
Richmond, T. J. 1984. Solvent accessible surface area and excluded volume in proteins. Analytical equations for overlapping spheres and implications for the hydrophobic effect. J. Mol. Biol. 178:6389.[CrossRef][Medline]
Rother, K., R. Preissner, A. Goede, and C. Frömmel. 2003. Inhomogeneous molecular density: reference packing densities and distribution of cavities within proteins. Bioinformatics. 19:21122121.
Senes, A., I. Ubarretxena-Belandia, and D. M. Engelman. 2001. The C
-H
O hydrogen bond: a determinant of stability and specificity in transmembrane helix interactions. Proc. Natl. Acad. Sci. USA. 98:90569061.
Stouten, P. F. W., C. Frömmel, H. Nakamura, and C. Sander. 1993. An effective solvation term based on atomic occupancies for use in protein simulations. Mol. Simulat. 10:97120.[CrossRef]
Swartz, K. J. 2004. Opening the gate in potassium channels. Nat. Struct. Mol. Biol. 11:499501.[CrossRef][Medline]
Tsai, J., R. Taylor, C. Chothia, and M. Gerstein. 1999. The packing density in proteins: standard radii and volumes. J. Mol. Biol. 290:253266.[CrossRef][Medline]
Vaguine, A. A., J. Richelle, and S. J. Wodak. 1999. SFCHECK: a unified set of procedures for evaluating the quality of macromolecular structure-factor data and their agreement with the atomic model. Acta Crystallogr. D Biol. Crystallogr. 55:191205.[CrossRef][Medline]
Van den Berg, B., W. M. Clemons Jr., I. Collinson, Y. Modis, E. Hartmann, S. C. Harrison, and T. A. Rapoport. 2004. X-ray structure of a protein-conducting channel. Nature. 427:3644.[CrossRef][Medline]
White, S. H., and W. C. Wimley. 1999. Membrane protein folding and stability: physical principles. Annu. Rev. Biophys. Biomol. Struct. 28:319365.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
P. W. Hildebrand, S. Gunther, A. Goede, L. Forrest, C. Frommel, and R. Preissner Hydrogen-Bonding and Packing Features of Membrane Proteins: Functional Implications Biophys. J., March 15, 2008; 94(6): 1945 - 1953. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |