Side-chain conformations in 4-{alpha}-helical bundles

Vasiliki E. Fadouloglou1,2, Nicholas M. Glykos2 and Michael Kokkinidis1,2,3

1 Department of Biology, University of Crete, P.O. Box 2208, GR-71409 Heraklion and 2 Foundation for Research and Technology–Hellas, Institute of Molecular Biology and Biotechnology (IMBB), P.O. Box 1527, GR-71110 Heraklion, Crete, Greece


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
The distribution of the {chi}1, {chi}2 dihedral angles in a dataset consisting of 12 unrelated 4-{alpha}-helical bundle proteins was determined and qualitatively compared with that observed in globular proteins. The analysis suggests that the 4-{alpha}-helical bundle motif could occasionally impose steric constraints on side chains: (i) the side-chain conformations are limited to only a subset of the conformations observed in globular proteins and for some amino acids they are sterically more constrained than those in helical regions of globular proteins; (ii) aspartic acid and asparagine occasionally adopt rotamers that have not been previously reported for globular or helical proteins; (iii) some rotamers of tyrosine and isoleucine are predominantly or exclusively associated with hydrophobic core positions (a, d); (iv) mutations in the hydrophobic core occur preferentially between residue types which among other physicochemical properties also share a predominant rotamer.

Keywords: 4-{alpha}-helical bundles/{chi}-dihedral angles/rotamer/side chain


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
The conformation of side chains is an essential feature of protein architecture. Consequently, knowledge of the factors that affect the side-chain conformations is significant, both for the understanding of protein folding and for the successful design of mutated proteins.

Side chains in proteins prefer certain conformations as shown by the non-uniform distribution of the {chi}-dihedral angles (Janin et al., 1978Go). The preferred conformations correspond to energy minima that are generally represented by three regions of {chi}1 (around 60°, 180° and –60°; Janin et al., 1978). Analyses of the distribution of {chi}-dihedral angles by two groups (James and Sielecki, 1983Go; Ponder and Richards, 1987Go) led to the definition of a rotamer as a dense cluster of points in the {chi}-angle space; furthermore, a rotamer library was developed based upon 19 well-determined protein structures. The relationship between secondary structure and side-chain conformations was subsequently investigated (McGregor et al., 1987Go; Summers et al., 1987Go). These studies revealed that the rotamer preferences of side chains are strongly affected by the secondary structure. A significant correlation between the backbone {phi},{Psi} values and the side-chain dihedral angles was also found (Dunbrack and Karplus, 1993Go). Based on a significantly enlarged data set, Schrauber et al. (Schrauber et al., 1993Go) further refined the original rotamer library for globular proteins introduced by Ponder and Richards (Ponder and Richards, 1987Go). The influence of backbone conformations was taken into account by grouping the rotamer distributions for each amino acid according to several secondary structure-based classes. Furthermore, the term `rotamericity' of an amino acid was introduced (Schrauber et al., 1993Go), defined as the ratio of the total number of occurrences of the specific amino acid in any of the possible rotamers to the total number of occurrences of this amino acid in the sample.

In the present study, the role of the geometric constraints posed by a specific topology to the side-chain dihedral angles was investigated. As a model for protein topology, the 4-{alpha}-helical bundle motif was used. This simple, recurrent tertiary motif consists of four {alpha}-helices packed against each other in an antiparallel manner at an angle of about 20° (Figure 1aGo) (Weber and Salemme, 1980Go; Cohen and Parry, 1986Go). The {alpha}-helices are usually connected together with loop regions; alternatively, the bundle is formed as an assembly of helices belonging to different polypeptide chains, as is the case with the ColE1 Rop protein (Banner et al., 1987Go; Presnell and Cohen, 1989Go; Harris et al., 1994Go). The amino acid sequences of the helices follow a specific pattern of hydrophilic and hydrophobic residues of the type (a,b,c,d,e,f,g)n (Crick, 1953Go). This pattern is repeated every seven residues (heptads). Positions a and d form the core of the bundle and are generally occupied by hydrophobic amino acids (Figure 1bGo).



View larger version (60K):
[in this window]
[in a new window]
 
Fig. 1. (a) Schematic representation of the 4-{alpha}-helical bundle motif. Each cylinder represents an {alpha}-helix. (b) Transversal section of a canonical 4-{alpha}-helical bundle. The positions of the heptad repeat are marked with the letters ag. The hydrophobic packing interactions among the residues of the core are illustrated with dotted lines.

 
The amino acid frequencies for the seven topologically distinct positions of the heptad repeat revealed highly specific relationships between topology and sequence preferences (Paliakasis and Kokkinidis, 1992Go). These preferences reflect the constraints imposed by topology to the amino acid sequence. A question of interest is whether steric hindrances posed by the 4-{alpha}-helical bundle topology are reflected not only on the pattern of amino acid sequence but also in the conformations of side chains as expressed in terms of {chi}-angles. In this study, the distribution of {chi}1, {chi}2 dihedral angles found in a sample of 12 4-{alpha}-helical bundles is compared with the conformational preferences of side chains observed in globular proteins or specifically in {alpha}-helices. However, this work can only provide qualitative information since the small size of the sample makes it difficult to draw definite conclusions.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
To analyze the conformations of the {chi}-dihedral angles, an initial sample of 443 4-{alpha}-helical bundle structures (326 of which were lysozyme molecules) was obtained from the Brookhaven Protein Data Bank (Bernstein et al., 1977; http://pdb-browsers.ebi.ac.uk//) and the FSSP database (Holm and Sander, 1996; http://www2.ebi.ac.uk/dali/fssp/fssp.html). To avoid bias due to sequence homologies, a subset of 12 structures with <25% sequence identity was used (Table IGo). The analysis was restricted to {alpha}-helical residues only (588 out of a total of 741 residues in the sample). Gly, Pro and Ala residues were excluded. Assignment of the heptad positions (ag) in the sequences of our sample was carried out manually by inspecting their structures with the program O (Jones et al., 1991Go). This assignment was generally unambiguous because of the high regularity of the motif. Side-chain dihedral angles were also determined with the program O. Owing to the high variability of the torsion angles beyond {chi}2, the analysis was restricted to {chi}1, {chi}2 only. For the classification of {chi}1, {chi}2 combinations, the conventions of Schrauber et al. (1993) were used, i.e. a side-chain conformation was assigned to a specific rotamer if the dihedral angles did not deviate by more than 20° from the values reported for this rotamer. To investigate whether the tightly packed hydrophobic cores impose special constraints to side chains, each amino acid type was analyzed according to its topological position in the bundle using a classification of residues into two groups. The first group consisted of internal residues (a and d positions) while the second group comprised more exposed residues (positions b, c, e, f and g).


View this table:
[in this window]
[in a new window]
 
Table I. The 4-{alpha}-helical bundles used for the analysis of the side-chain conformations
 
Furthermore, the constraints imposed by 4-{alpha}-helical bundle cores on the evolution of their sequences were examined from the perspective of rotamer conservation. More specifically, the internal positions (a and d) in a series of aligned sequences of homologous 4-{alpha}-helical bundles were compared and the pattern of amino acid substitutions was examined from the aspect of rotamer conservation. To perform this analysis, representative sets of homologous sequences for eight of the proteins of our sample were identified using the Swiss Prot data bank (Bairoch and Apweiler, 1997Go; see also Table IIGo) and the FASTA program (Pearson and Lipman, 1988Go). For each of the eight resulting groups, the sequence of the PDB structure was used as a reference. Within each group, the sequence alignment program ALIGN (http://www2.igh.cnrs.fr/bin/align-guess.cgi) was used to perform the pairwise alignments of all sequences with the reference sequence. Sequence homology combined with the characteristics of the heptad pattern allowed an unambiguous identification of the equivalent a, d positions within each group.


View this table:
[in this window]
[in a new window]
 
Table II. The eight groups of homologous sequences and their Swiss Prot entries
 

    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Owing to the small size of the sample (see Materials and methods section), only qualitative information can be obtained by this work. Nevertheless, as far as we can assert by estimated standard deviations, our results (rotamer preferences, etc.) are generally consistent with the statistics of the much larger sample of globular proteins performed by other authors (Schrauber et al., 1993Go).

Amino acid composition of the sample

The amino acid composition of the helical parts of the proteins used in this work is presented in Figure 2Go. Compared with an earlier analysis (Paliakasis and Kokkinidis, 1992Go) it shows only minor differences, i.e. (i) a higher occurrence of Leu and (ii) an increase in the Leu to Ala ratio. The composition of internal (a, d) positions (Figure 2Go) shows a clear predominance of Leu, Ala, Ile and Val residues, which agrees well with the preferences found by Paliakasis and Kokkinidis.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 2. Composition of the sample of 4-{alpha}-helical bundles. Gray bars, over all amino acids; black bars, composition of the hydrophobic cores (positions a and d).

 
Rotamer preferences

The clusters of {chi}1, {chi}2 dihedral angles found for the side chains of our sample are presented in Table IIIGo. The most striking feature of this distribution for the majority of amino acid types is that a large fraction of the side chains belong to a single rotamer. Ile, Val, Thr and Cys are extreme examples, where the dominant cluster contains at least five times more members than the next one. On the other hand, Leu and Arg have two densely populated rotamers whereas Ser and Glu do not show pronounced preferences for a particular rotamer.


View this table:
[in this window]
[in a new window]
 
Table III. Summary of the rotamers found in 4-{alpha}-helical bundles and rotamers of globular proteins reported by other authors

 
Table IIIGo also includes the rotamers reported for globular proteins by Schrauber et al. (Schrauber et al., 1993Go) and Ponder and Richards (Ponder and Richards, 1987Go), and also those which are most frequently found in {alpha}-helical regions of globular proteins (in bold in Table IIIGo). Comparison between our sample and globular proteins shows that, as a rule, the preferred side-chain conformations of 4-{alpha}-helical bundles are limited to a small subset of the rotamers found in globular proteins. Furthermore, the side-chain conformations of Tyr, Met, Thr and Cys appear to be more constrained in our sample compared with {alpha}-helical regions of globular proteins because there is a strong preference for one out of several possible `helical rotamers'. As shown in Table IIIGo, at least the 50% of the side-chain population of the above amino acid types adopt just one of the rotamers found in {alpha}-helical regions. A novel rotamer with {chi}1 = 178° and {chi}2 = 64° which to our knowledge has not been reported earlier, was found for Asp (see Table IIIGo and Figure 3aGo). Similarly, in the case of Asn, there are indications of a novel cluster with {chi}1 = –170° and {chi}2 = 59° (Table IIIGo and Figure 3bGo).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 3. Distribution of the {chi}1, {chi}2 dihedral angles for (a) Asp and (b) Asn residues. The circles represent the data points of our sample. Stars mark the positions of the rotamers listed in the Ponder and Richards rotamer library (Ponder and Richards, 1987Go) [no values were reported for Asp and Asn by Schrauber et al. (Schrauber et al., 1993Go)]. The symmetry around the {chi}2 angle is taken into account. Graphs were prepared using the program Xmgr (http://plasma-gate.weizmann.ac.il/Xmgr/).

 
A possible interpretation of these findings could be that the 4-{alpha}-helical bundle motif imposes special constraints on side-chain conformations. If this is indeed so, one would expect these constraints to be more pronounced for residues of the tightly packed hydrophobic core. For this reason, the residues occupying the hydrophobic core positions a and d of the heptad repeat were treated separately. For each amino acid type, the residues were distinguished in two groups according to their location in (a, d) and (b, c, e, f, g) positions and the {chi}1, {chi}2 plots were obtained. Inspection of these diagrams showed that for Ile and Tyr residues the two groups are segregated to different rotamers (Figure 4a and bGo). Two clusters which appear to be better compatible with internal residues (a and d positions) were found for these cases. These clusters. with {chi}1 = –62°, {chi}2 = –60° for Ile and {chi}1 = –70°, {chi}2 = 107° for Tyr, correspond to rotamers which are also present in the globular proteins, but which have not been reported yet in connection with protein cores.



View larger version (27K):
[in this window]
[in a new window]
 
Fig. 4. {chi}1, {chi}2 angle distribution for (a) Tyr and (b) Ile residues. Internal positions (a, d) are indicated by black circles, semi-exposed or exposed residues are presented as white circles. Crosses denote rotamers for globular proteins (Schrauber et al., 1993Go). For Tyr the symmetry around the {chi}2 angle is taken into account.

 
Another useful indicator for assessing the influence of position-specific effects to the side-chain conformations is the rotamericity (Schrauber et al., 1993Go). For each amino acid type, the rotamericities were calculated separately for the internal and external positions (data not shown). A tendency of the internal residues for higher rotamericity values is observed compared with the external positions.

Correlation between frequency of natural mutations and rotamer preferences in the hydrophobic core of the 4-{alpha}-helical bundles

Natural mutations in the hydrophobic core of 4-{alpha}-helical bundles were analyzed from the aspect of rotamer conservation [rotamer conservation in proteins has been reported (Summers et al., 1987Go) and forms the basis of molecular modeling by homology]. In this analysis, the sequence alignments described in the Materials and methods section were examined and the pattern of sequence variation at the a and d positions was studied. From a total of 226 pairs of aligned amino acids in the core, 81 substitutions were found. In Table IVGo the observed substitutions and the corresponding frequencies for the residues commonly found in a and d positions are presented. An interesting observation is the relatively high frequency of Val->Leu and Met->Leu substitutions. It is worth noting that only four types of substitution account for 40% of the total observed (Met->Leu, Val->Leu, Leu->Ile and Leu->Val). Comparison of Tables III and IVGoGo shows that substitutions occur overwhelmingly between residue types that have their major rotamers in common. For example, the main rotamer of Ile and Met coincides with one of the main rotamers of Leu. As an exception to the above observations, Phe->Ile substitution occurs between residues that do not share any common rotamer while some relatively rare substitutions (e.g. Met->Phe, Tyr->Asp and Tyr->Met) occur between residues which do not share a major rotamer but have at least one less populated {chi}1, {chi}2 cluster in common. For example, the dominant rotamer of Tyr and Phe ({chi}1 = 176°, {chi}2 = 78° and {chi}1 = 180°, {chi}2 = 75°, respectively) coincide with one of the `rare' rotamers of Met ({chi}1 = –175°, {chi}2 = 64°).


View this table:
[in this window]
[in a new window]
 
Table IV. Substitutions observed in a and d positions of 4-{alpha}-helical bundles
 

    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
The conclusions of this work can be summarized as follows. (i) In addition to the known influence of {alpha}-helical secondary structure on the side-chain conformations (Schrauber et al., 1993Go), the 4-{alpha}-helical bundle topology probably imposes on some residues additional constraints and restricts the permitted conformations overwhelmingly to a unique rotamer. (ii) The side chains of aspartic acid and probably asparagine occasionally adopt a conformation that is associated with a novel rotamer. (iii) For isoleucine and tyrosine, rotamers were found which appear to be more compatible with residues in the hydrophobic core. (iv) Natural mutations in a and d positions of the 4-{alpha}-helical bundles tend to occur between amino acids which among other properties have their main rotamers in common.

Previous studies have demonstrated that the secondary structure significantly affects the side-chain conformations by restricting them to a subset of those found in globular proteins (McGregor et al., 1987Go; Schrauber et al., 1993Go). This is confirmed and reinforced by our analysis, which in addition showed that for some amino acids (Tyr, Met, Thr, Cys) the preference of specific side-chain conformations is much more pronounced even when compared to {alpha}-helices. This could be an effect of the 4-{alpha}-helical bundle topology, which imposes additional constraints to side chains limiting the permitted conformations to a subset of those adopted by {alpha}-helices.

Generally, the majority of side chains in the 4-{alpha}-helical bundles do not adopt novel or unusual conformations; side chains are clustered in some of the already known regions (from the globular proteins) of the {chi}1, {chi}2 space. A notable exception to this statement is the Asp residue, several occurrences of which have been classified as a novel rotamer (Figure 5aGo). This result was unexpected given that Asp is not subject to the specific constraints of the motif, since it systematically occurs in external or semi-buried positions. Indications of a new rotamer are also present for Asn. An analysis of all Asp and Asn residues in our sample shows that their temperature factors are fairly low and very close to or even lower than the average temperature factor of the structure to which they belong. Thus, the observed behavior is probably not an artifact due to an increased side-chain flexibility. The novel rotamers are not associated with a particular position of Asp/Asn residues (e.g. capping residues) or a particular backbone conformation (all Asp/Asn residues studied have helical {phi},{Psi} angles). Detailed inspection and comparison with the known rotamers of the globular proteins, showed that the novel rotamer places the side chain in an optimal position relative to the protein backbone and the side chains of the neighboring {alpha}-helices, so that electrostatic repulsions between charged groups and steric hindrances between the side chains are minimized. However, the size of our sample is too small to go beyond a qualitative discussion of these novel rotamers.



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 5. Examples of side-chain conformations. (a) Asp residues (1nfn). A represents a motif-specific conformation ({chi}1 = 177°, {chi}2 = 69°) and B a non-specific conformation ({chi}1 = –62°, {chi}2 = 113°). (b) Tyr residues (1rhg). C represents a position-dependent rotamer ({chi}1 = 159°, {chi}2 = 62°) and D a position-independent rotamer ({chi}1 = –78°, {chi}2 = –76°). The figure was prepared using the programs O (Jones et al., 1991Go) and Rasmol (Collaborative Computational Project Number 4, 1994Go).

 
A related observation is that Ile and Tyr exhibit a rotamer that is exclusively occupied by internal residues (positions a and d). Taking into account our limited data we could only assert that these conformations appear to be less well accepted by external residues; however, they are adopted by a significant fraction of the total number of internal residues (~50% for Tyr and 20% for Ile). For Ile, in particular, this position-specific conformation does not belong to the preferred rotamer of {alpha}-helices (Schrauber et al., 1993Go) and is systematically adopted by residues located at the terminal caps of the bundles. However, when the individual amino acids that exhibit this special conformation were systematically examined in comparison with their local environments, no obvious reason that could provide an explanation for this positional dependence was found. On the other hand, for the case of Tyr it is clear that the position-dependent rotamer arranges the bulky, aromatic group out of the hydrophobic core of the bundle and unfavorable packing of the side chains is prevented (Figure 5bGo).

Side-chain shape, volume, polarity, packing density and cavity volume have been reported as factors that affect the pattern of substitutions in the protein interior (Bordo and Argos, 1990Go; Vlassi et al., 1999Go). Conservation of these structural parameters is important for hydrophobic core mutations. Our present study provides evidence that the pattern of side-chain substitutions in 4-{alpha}-helix bundles is also consistent with the conservation of highly populated rotamers. In a previous study, Summers et al. (1987) concluded that for structurally and functionally homologous proteins, there is a high probability that the side-chain conformations are conserved. For amino acid substitutions, the conservation of orientation of both C{gamma} and C{delta} atoms, was estimated to be in order of 35–75% (Summers et al., 1987Go). In the case of 4-{alpha}-helical bundles, the frequency of the amino acid substitutions that occur between residues with at least one rotamer in common was estimated to be ~75% and indicates an extensive conservation of rotamers in the core of homologous bundles. This is consistent with recent work by Vlassi et al. (1999), which found a pronounced tendency of 4-{alpha}-helical bundles to preserve hydrophobic core packing interactions upon mutations. Exceptions to this behavior (e.g. Phe->Ile substitutions), where rotamer conservation is not possible, occur at the end of the bundles. It is reasonable to assume that the constraints of the motif at those positions are weaker.

The aim of this work, as mentioned, was to examine qualitatively whether steric hindrances posed by the 4-{alpha}-helical bundle topology are reflected in the conformations of side chains. Within this framework, it has been shown that a tertiary motif can affect to some extent the conformations adopted by the side chains of some amino acids. This happens first through additional restrictions imposed to the side-chain conformations, second through the formation of novel {chi}1, {chi}2 clusters and third through some rotamers that appear to be more compatible with internal residues. A natural extension of this work would be a systematic analysis of additional known tertiary motifs. Such a study would provide valuable information for the understanding of protein folding and would find applications in the successful design of mutations and in the homology modeling of new structures.


    Notes
 
3 To whom correspondence should be addressed. E-mail: kokkinid{at}imbb.forth.gr Back


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Bairoch,A. and Apweiler,P. (1997) J. Mol. Med., 75, 312–316.[ISI][Medline]

Banner,D.W., Kokkinidis,M. and Tsernoglou,D. (1987) J. Mol. Biol., 196, 657–675.[ISI][Medline]

Bernstein,F.C., Koetzle,T.F., Williams,G.J.B., Meyer,E.F., Brice,M.D., Rogers,J.R., Kennard,O., Shimanouchi,T. and Tsumi,M. (1997) J. Mol. Biol., 112, 535–542.

Bordo,D. and Argos,P. (1990) J. Mol. Biol., 211, 975–988.[ISI][Medline]

Cohen,C. and Parry,D.A.D. (1986) Trends Biochem. Sci., 11, 245–248.[ISI]

Collaborative Computational Project Number 4 (1994) Acta Crystallogr., D50, 760–763.[ISI]

Crick,F.H.C. (1953) Acta Crystallogr., 6, 689–697.[ISI]

Dunbrack,L.R. and Karplus,M. (1993) J. Mol. Biol., 230, 543–574.[ISI][Medline]

Harris,N.L., Presnell,S.R. and Cohen,F.E. (1994) J. Mol. Biol., 236, 1356–1368.[ISI][Medline]

Holm,L. and Sander,C. (1996) Science, 273, 595–602.[Abstract/Free Full Text]

James,M.N.G. and Sielecki,A.R. (1983) J. Mol. Biol., 183, 299–361.

Janin,J., Wodak,S., Levitt,M. and Maigret,B. (1978) J. Mol. Biol., 125, 357–386.[ISI][Medline]

Jones,T.A., Zou,J-Y., Cowan,S.W. and Kjeldgaard,M. (1991) Acta Crystallogr., A47, 110–119.[ISI]

McGregor,M.J., Islam,S.A. and Sternberg,M.J.E. (1987) J. Mol. Biol., 198, 295–310.[ISI][Medline]

Paliakasis,C.D. and Kokkinidis,M. (1992) Protein Eng., 5, 739–748.[Abstract]

Pearson,W.R. and Lipman,D.J. (1988) Proc. Natl Acad. Sci. USA, 85, 2444–2448.[Abstract]

Ponder,J.W. and Richards,M. (1987) J. Mol. Biol., 193, 775–791.[ISI][Medline]

Presnell,S.R. and Cohen,F.E. (1989) Proc. Natl Acad. Sci. USA, 86, 6592–6596.[Abstract]

Schrauber,H., Eisenhaber,F. and Argos,P. (1993) J. Mol. Biol., 230, 592–612.[ISI][Medline]

Summers,N.L., Carlson,W.D. and Karplus,M. (1987) J. Mol. Biol., 196, 175–198.[ISI][Medline]

Vlassi,M., Cesareni,G. and Kokkinidis,M. (1999) J. Mol. Biol., 285, 817–827.[ISI][Medline]

Weber,P.C. and Salemme,F.R. (1980) Nature, 287, 82–84.[ISI][Medline]

Received September 14, 2000; accepted February 15, 2001.





This Article
Abstract
FREE Full Text (PDF)
Alert me when this article is cited
Alert me if a correction is posted
Services
Email this article to a friend
Similar articles in this journal
Similar articles in ISI Web of Science
Similar articles in PubMed
Alert me to new issues of the journal
Add to My Personal Archive
Download to citation manager
Search for citing articles in:
ISI Web of Science (1)
Request Permissions
Google Scholar
Articles by Fadouloglou, V. E.
Articles by Kokkinidis, M.
PubMed
PubMed Citation
Articles by Fadouloglou, V. E.
Articles by Kokkinidis, M.