|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Department of Biochemistry, Weill Medical College of Cornell University, New York, New York 10021, USA
2 Bioinformatics Laboratory, Department of Applied Mathematics and Statistics, State University of New York at Stony Brook, Stony Brook, New York 11794-3600, USA
Reprint requests to: Ilya A. Vakser, Bioinformatics Laboratory, Department of Applied Mathematics and Statistics, State University of New York at Stony Brook, Stony Brook, NY 11794-3600, USA; e-mail: vakser{at}ams.sunysb.edu; fax: (631) 632-8490.
(RECEIVED November 3, 2003; FINAL REVISION January 16, 2004; ACCEPTED January 19, 2004)
| Abstract |
|---|
|
|
|---|
Keywords: protein modeling; structure prediction; protein folding; secondary structure packing; docking
Article published ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.03505804.
| Introduction |
|---|
|
|
|---|
-helices,
-strands, and loops). Low-resolution docking studies (Vakser 1996b; Vakser et al. 1999) of helixhelix interactions (Vakser 1996a; Vakser and Jiang 2002) indicated the existence of geometrically preferred sides/faces in helices for the interaction with other helices. Our earlier study of helix bundles in integral membrane proteins (Jiang and Vakser 2000) revealed that side chains inside helixhelix interfaces, on average, are shorter than those in the noninterface parts of the helices. Thus, transmembrane helices, if viewed along the main axis, on average, are not circles but rather ellipses. This effect brings main axes of helices closer to each other and maximizes the surface contact between the helices (thus enhancing the packing/stability of the protein). The membrane-spanning helices were chosen because they are well packed, aligned, and almost parallel to each other in the transmembrane bundles, thus being an ideal object for such a study. These results obtained on a few available crystal structures of transmembrane proteins raised a natural question: whether the same phenomenon applies to helixhelix pairs in all protein structures. In the present study, a comprehensive nonredundant set of interacting helices from the entire PDB was used for the analysis of side-chain length distribution. The results show that similar to the earlier analyzed transmembrane subset, on average, the side chains in helices are shorter inside the helixhelix interfaces than those outside the interface areas. This small but statistically significant difference places helices closer to each other and enhances the packing of the structure. This phenomenon provides a new insight into the packing and stability of protein structures.
| Data set of helices |
|---|
|
|
|---|
Selection of helices
Helical fragments in the protein structures were identified by DSSP (Kabsch and Sander 1983). To guarantee at least two turns in all faces of a helix, only those with
8 residues were selected for the next step. Each selected helix was checked against the rest of the helices in the protein for interaction. Two helices were considered interacting if each had
4 residues in contact with the other helix. Following our earlier study (Jiang and Vakser 2000), two residues were considered to be in contact if they had atoms within the sum of their van der Waals radii plus 0.2 Å.
Unlike most transmembrane helices that are roughly parallel to each other, helices in globular proteins are packed at the variety of angles (Reddy and Blundell 1993; Walther et al. 1996; Hespenheide and Kuhn 2003). The helix packing angle was used to determine the interacting segments in the helices. The outstanding parts of the helices were deleted. The definition of the interacting segments is shown in Figure 1A
. All pairs of helical segments satisfying the above criteria were considered, including more than one segmentsegment interaction per helix (where the helix interacted with more than one helix).
|
atom was in the interface. The criterion was applied regardless of the number of helixhelix interfaces per helix. If a segment of a helix was involved in more than one helixhelix interaction, a residue in that segment could be in the interface in one interaction and outside the interface in the other interaction. In such cases the residue was considered to be the interface one.
Results and Discussion
The sizes of the 20 residue types (Table 1
) taken from Levitt (1976) were defined as average distances from the C
atom to the side-chain center of mass in a database of protein structures. Although the number of available structures in that database was limited, the number of occurrences of each amino acid type was sufficient for a simple length quantification scheme. Sizes of interface side chains and noninterface side chains were averaged over all helices in the data set. The average size of 6749 interface residues was 1.92 Å, and the average size of 12,806 noninterface residues was 2.07 Å (Fig. 2
). The 0.15-Å difference is statistically significant, with a P-value = 1029 in the Student t-test, assuming equal variance.
|
|
The existence of a small but statistically significant difference in length between helixhelix interface and noninterface side chains in soluble proteins determined in this study is consistent with a similar effect discovered in the transmembrane helices (Jiang and Vakser 2000). The example of this effect is shown in Figure 3
. It is important to emphasize though that this difference is determined statistically, on a large number of helixhelix pairs; thus, the existence and the extent of such a difference in individual helixhelix pairs vary. There might be even entire classes of protein structures and structural motifs where this general rule does not necessarily apply (e.g., in leucine zippers the helixhelix interfaces are known to be formed mostly by residues with long side chainsLeu, Lys, Arg, Glu, Gln). However, on average, for the entire data set of protein structures the rule still holds.
|
Such interface/noninterface difference enhances the packing of the helices by flattening the interface and bringing the helices closer to each other (Fig. 4
), and thus, may have resulted from evolutionary pressure to create tighter packing of protein structures. The effect provides an interesting insight into aspects of protein structure and folding.
|
| Acknowledgments |
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| References |
|---|
|
|
|---|
Hobohm, U. and Sander, C. 1994. Enlarged representative set of protein structures. Protein Sci. 3: 522524.[Abstract]
Jiang, S. and Vakser, I.A. 2000. Side chains in transmembrane helices are shorter at helixhelix interfaces. Proteins 40: 429435.[CrossRef][Medline]
Jiang, S., Tovchigrechko, A., and Vakser, I.A. 2003. The role of geometric complementarity in secondary structure packing: A systematic docking study. Protein Sci. 12: 16461651.
Kabsch, W. and Sander, C. 1983. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 22: 25772637.[CrossRef][Medline]
Levitt, M. 1976. A simplified representation of protein conformations for rapid simulation of protein folding. J. Mol. Biol. 104: 59107.[CrossRef][Medline]
Li, R. and Woodward, C. 1999. The hydrogen exchange core and protein folding. Protein Sci. 8: 15711591.[Abstract]
Reddy, B.V.B. and Blundell, T.L. 1993. Packing of secondary structure elements in proteins. Analysis and prediction of inter-helix distances. J. Mol. Biol. 233: 464479.[CrossRef][Medline]
Vakser, I.A. 1996a. Long-distance potentials: An approach to the multiple-minima problem in ligandreceptor interaction. Protein Eng. 9: 3741.
. 1996b. Low-resolution docking: Prediction of complexes for underdetermined structures. Biopolymers 39: 455464.[CrossRef][Medline]
Vakser, I.A. and Jiang, S. 2002. Strategies for modeling the interactions of the transmembrane helices of G-protein coupled receptors by geometric complementarity using the GRAMM computer algorithm. Methods Enzymol. 343: 313328.[Medline]
Vakser, I.A., Matar, O.G., and Lam, C.F. 1999. A systematic study of low-resolution recognition in proteinprotein complexes. Proc. Natl. Acad. Sci. 96: 84778482.
Walther, D., Eisenhaber, F., and Argos, P. 1996. Principles of helixhelix packing in proteins: The helical lattice superposition model. J. Mol. Biol. 255: 536553.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |