| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
/ß)8 Barrel Proteins
Computational Biology Research Center (CBRC), Institute of Advanced Industrial Science and Technology (AIST) 2-41-6 Aomi, Koto-ku, Tokyo 135-0064, Japan; and * Department of Physics, Bharathidasan University, Tiruchirapalli 620 024, Tamilnadu, India
Correspondence: Address reprint requests to M. Michael Gromiha, Tel.: +81-3-3599-8046; Fax: +81-3-3599-8081; E-mail: michael-gromiha{at}aist.go.jp.
| ABSTRACT |
|---|
|
|
|---|
/ß)8 barrel proteins provides ample light to understand the factors that are responsible for directing and maintaining their common fold. In this work, the hydrophobically enriched clusters are identified in 92% of the considered (
/ß)8 barrel proteins. The residue segments with hydrophobic clusters have high thermal stability. Further, these clusters are formed and stabilized through long-range interactions. Specifically, a network of long-range contacts connects adjacent ß-strands of the (
/ß)8 barrel domain and the hydrophobic clusters. The implications of hydrophobic clusters and long-range networks in providing a feasible common mechanism for the folding of (
/ß)8 barrel proteins are proposed. | INTRODUCTION |
|---|
|
|
|---|
/ß)8 barrel fold constituted by eight parallel ß-strands surrounded by eight
-helices (Chothia, 1988
/ß)8 barrel proteins have been reviewed in detail (Pujadas and Palau, 1999
The unique three-dimensional structure and stable secondary structures of globular proteins are determined by the interactions of amino acid residues among themselves along the polypeptide chain as well as with the surrounding medium. The interactions between amino acid residues have been classified into short-, medium-, and long-range, based on their distance of separation (Gromiha and Selvaraj, 1997a
). This classification has been used successfully to address the problem of protein folding and sequence recognition (Gromiha and Selvaraj, 1999
; 2001
; Miyazawa and Jernigan, 1999
). During the process of protein folding, the cooperative, noncovalent, and long-range interactions among residues provide stability to resist the local tendency for unfolding. Dosztanyi et al. (1997)
have identified clusters of such residues as stabilization centers. Further, most of these centers are found at buried positions, and have hydrophobic and aromatic side chains. Poupon and Mornon (1999)
showed the correspondence between hydrophobic positions of a given fold that constitute a folding nucleus, which is considered to play a key role in protein folding (Dobson and Karplus, 1999
), and amino acid residues responsible for specific interactions. The role of each amino acid residue toward different noncovalent interactions has been delineated in our earlier works (Ponnuswamy and Gromiha, 1994
; Gromiha and Selvaraj, 1999
).
In (
/ß)8 barrel proteins, we have identified similar and identical tertiary clusters and analyzed the importance of medium- and long-range interactions for the stabilization of (
/ß)8 barrel fold (Selvaraj and Gromiha, 1998a
,b
; Kannan et al., 2001
). It is of interest to understand the influence of long-range interactions in the mechanism underlying the common tertiary fold of this class of proteins. In the present work, we identified the hydrophobic clusters in (
/ß)8 barrel proteins and characterized the importance of medium- and long-range interactions in the formation of these centers. Further, we observed a network of long-range contacts in ß-strands belonging to (
/ß)8 barrel proteins. Based on these results, we propose a common mechanism for the folding of this class of proteins.
| MATERIALS AND METHODS |
|---|
|
|
|---|
/ß)8 barrel proteins from the information available in the literature (Farber and Petsko, 1990
/ß sandwich or roll by CATH database (Orengo et al. 1999
Computational methods
Delineation of hydrophobic domains and key residue in hydrophobic clusters
The amino acid residues in a protein molecule are represented by their
-carbon atoms and each residue is assigned with the hydrophobicity index obtained from thermodynamic transfer experiments (Nozaki and Tanford, 1971
; Jones, 1975
). The surrounding hydrophobicity is defined as the sum of hydrophobic indices of various residues that appear within 8-Å radius limit from a given residue (Manavalan and Ponnuswamy, 1977
, 1978
; Ponnuswamy, 1993
; Ponnuswamy and Gromiha, 1993
). It has been shown that the influence of each residue over the surrounding medium extends effectively only up to 8 Å (Manavalan and Ponnuswamy, 1977
) and this limit is sufficient to characterize the hydrophobic behavior of amino acid residues (Manavalan and Ponnuswamy, 1978
) and to accommodate both the local and nonlocal interactions (Gromiha and Selvaraj, 2000
; Jiang et al., 2002
). Further, 8-Å limit has been used in several studies, such as, to understand the folding rate of two-state proteins (Debe and Goddard, 1999
; Gromiha and Selvaraj, 2001
), protein stability upon mutations (Gromiha et al., 1999
), thermal stability of proteins (Gromiha, 2001
; Gromiha and Thangakani, 2001
), and to determine the transition state structures of two-state protein mutants (Gromiha and Selvaraj, 2002
).
The local regions that are comprised of a cluster of residues with high surrounding hydrophobicity (equal to or greater than twice the average value for all the residues in a protein) are considered as hydrophobic domains; the residue of the highest surrounding hydrophobicity within a domain is taken as the key residue in hydrophobic cluster (Ponnuswamy and Prabhakaran, 1980
).
Computation of medium- and long-range contacts in the hydrophobic clusters of (
/ß)8 barrel proteins
For each residue in the (
/ß)8 barrel proteins, we computed the residues coming within a sphere of 8-Å radius as described in previous section. For a given residue, the composition of surrounding residues is analyzed in terms of the location at the sequence level and the contributions from <±3 residues are treated as short-range contacts, ±3 or ±4 residues as medium-range contacts, and >±4 residues are treated as long-range contacts (Gromiha and Selvaraj, 1997a
; 2000
).
Computation of thermal unfolding character in hydrophobic clusters
We followed the method of Ponnuswamy's group to compute the thermal unfolding behavior of amino acid residues in globular proteins (Ponnuswamy et al., 1982
; Muthusamy et al., 2000
). It has been shown that the set of amino acid residues Asp, Cys, Glu, Lys, Leu, Arg, Trp, and Tyr enhance the stability and the set of residues Ala, Gly, Gln, Ser, Thr, and Val decrease the stability. Accordingly, the amount of stabilizing (X1) and destabilizing (X2) groups (in %) were obtained for each residue within the 8-Å radius volume. The thermal unfolding for a residue is calculated using the equation, Tm = 64.462 + 0.894 X1 - 0.591 X2.
| RESULTS AND DISCUSSIONS |
|---|
|
|
|---|
/ß)8 barrel proteins along with their location, surrounding hydrophobicity, and the number of medium- and long-range contacts, are presented in Table 1. We observed that, except for 1RUS, 1MUC, and 2BNH, all the (
/ß)8 barrel proteins contain hydrophobic clusters, which contribute a good packing of the protein core. There is no direct relationship between number of residues and number of hydrophobic clusters. The protein 2TMD has the highest of seven hydrophobic clusters, followed by 2TAA and 1CDG, that have five hydrophobic clusters each. Among these three proteins, all the hydrophobic clusters in Taka amylase (2TAA) fall within the (
/ß)8 barrel domain and the contact map (within a spatial distance of 8 Å) for 2TAA is shown in Fig. 1. This figure shows the presence of four hydrophobic clusters in ß-strands and one in
-helix. Interestingly, the faraway residues contribute toward the formation of hydrophobic clusters (Y12, I60, T291, and I326) in ß-strands, whereas medium-range contacts are responsible for T239 in
-helix.
|
|
/ß)8 barrel proteins.
Recently, Nagano et al. (2002)
have made an extensive analysis of the sequence, structure and function of proteins that adopt the TIM barrel fold. They classified the proteins belonging to the (
/ß)8 barrel fold into 18 families on the basis of homology. We have examined the conservation of hydrophobic clusters identified in the present work within these groups. Among the 36 considered proteins, 2TAA, 1BYB, 1CTN, 1NAR, and 1CDG belong to the Glycosidase family. Within this family, hydrophobic clusters have been identified in the first beta strands (ß1) of 2TAA, 1BYB, and 1CDG. In both 2TAA and 1CDG another hydrophobic cluster is identified in the second beta-strand (ß2). Further, both 2TAA and 1NAR have a hydrophobic cluster at the seventh beta-strand (ß7). These observations suggest that the hydrophobic clusters may be conserved within homologous family members.
Role of long-range interactions in hydrophobic clusters
We have estimated the relative contribution of medium- and long-range interactions for the key residue in the hydrophobic cluster of (
/ß)8 barrel proteins and the results are included in Table 1. The percentage of long-range contacts for each key residue is obtained by dividing the number of long-range contacts to the total number of contacts. Interestingly, we observed that the long-range interactions have the highest contribution, 38.5 to 76.5% for each protein, and the average contribution is 65.4%. This result indicates the vital role of long-range interactions to form the hydrophobic clusters and to stabilize the proteins. Further, the long-range interactions contribute an appreciably higher percentage in the hydrophobic clusters of ß-strands or turn regions near the strands, and there is no significant contribution from medium-range interactions. However, the medium-range interactions play a dominant role in the hydrophobic clusters formed by
-helices. This observation is consistent with the previous results that
-helices are influenced by medium-range interactions and ß-strands are dominated by long-range interactions (Gromiha and Selvaraj, 1998
).
Occurrence of hydrophobic clusters in secondary structures
It is generally observed that in (
/ß)8 barrel proteins, the
-helical segments connecting two ß-strands are amphipathic and the hydrophobic face of the helix interacts with the ß-strands (Pujadas and Palau, 1999
). The secondary structure content of all the 36 (
/ß)8 barrel proteins is presented in Table 2. We observed that 38% of the residues are in
-helix, 17% in ß-strand, and 45% in coil region. Our analysis on the location of key residues in hydrophobic clusters in different secondary structures reveals that 74% of the residues fall in the ß-strands of barrel fold. Similar result (71%) is observed with the data set of Kannan et al. (2001)
. Interestingly, 16% of hydrophobic clusters are identified in the coil region and only 10% are found in
-helical segments. We also computed the random probability for each key residue in the cluster to belong to these secondary structures. The values are 31% in helix, 23% in strand, and 46% in coil (Table 2). Further, we observed from Table 2 that the locations of more than 70% of the residues with high random probability are different from that in three-dimensional structure. This result demonstrates the importance of the location of key residues in different secondary structures for the folding and stability of (
/ß)8 barrel proteins.
|
As an example, in Fig. 2 we have illustrated the long-range contact network in the hydrophobic clusters of 2TAA. This protein has five hydrophobic clusters in which four of them are in ß-strands and one in
-helix of the (
/ß)8 barrel domain. For the key residue at I60 in ß-strand (S2), long-range contacts are observed with the hydrophobic residues, I11, Y12, F13, L14, M112, Y113, L114, and M115. At the N-terminal side, the farthermost long-range contact was observed at I11, belonging to the first strand (S1) and part of the key residue of first hydrophobic cluster (HC1). Further examination at I11 indicated no long-range contact in the same direction. In the C-terminal side of residue I60, we observed the farthest long-range contact at M115 (S3). Further search resulted in residues at R204, E230, L250, T291, and I326. Interestingly, all these residues are part of the ß-strands, S4, S5, S6, S7, and S8, respectively, in which the residues T291 and I326 are themselves hydrophobic clusters of 2TAA. Thus, we observed a good pattern of long-range contact network with links to all ß-strands constituting the (
/ß)8 barrel domain and also among the hydrophobic clusters. A similar pattern of long-range network was also observed in other (
/ß)8 barrel proteins.
|
-helical (H5) part of (
/ß)8 barrel. On the N-terminal side, the network was observed at K209 and E162; no further extension was possible at E162 as it has no long-range contacts. It is noteworthy that each of these residues belong to coil regions. On the C-terminal side, the network terminates at P253 that is in
-helix (H6).
Implications for (
/ß)8 barrel folding and stability
The location of hydrophobic clusters and the existence of long-range networks connecting the ß-strands of the (
/ß)8 barrel domain lead one to envisage the following probable scenario on (
/ß)8 barrel folding and stability. In our earlier works, we have shown the dominance of short- and medium-range interactions in
-helices and long-range interactions in ß-strands (Gromiha and Selvaraj, 1997a
,b
; 1998
). In (
/ß)8 barrel proteins,
-helices and ß-strands alternate each other. During the process of folding, initially short- and medium-range interactions may predominate and direct the formation of
-helices. This view has been supported by several recent theoretical and experimental studies on the importance of local and nonlocal interactions in the folding of globular proteins. In molecular dynamics simulation of a small
/ß protein, Sheinerman and Brooks (1998)
observed that
-helix forms earlier and ß-strand forms concomitantly during the process of protein folding. Further, the studies by Unger and Moult (1996)
support the idea that initial formation of local substructures is important to the folding of proteins.
The formation of
-helices in (
/ß)8 barrel domain will bring together adjacent ß-strand forming residues in spatial proximity so that these residues can interact with each other to form ß-strands, which are stabilized by both hydrogen-bonding network and hydrophobic interactions. Successive strands will also be stabilized by the enrichment of hydrophobic interactions. During the process of protein folding, Parker et al. (1996)
showed through protein engineering and relaxation kinetics experiments in an
/ß protein that sequence-local groups produced microdomains which is followed by the contacts between remote segments of secondary structures. Thus, the hydrophobic clusters and long-range contacts may steer the folding and stabilization of the (
/ß)8 barrel fold.
Thermal unfolding character of key residues in hydrophobic clusters
We have analyzed the thermal unfolding behavior of all the residues forming hydrophobic clusters as described in the Methods section. We found that most of the residues forming hydrophobic clusters have high thermal stability and withstand higher temperature. Thus these residues help to maintain the strong interactions with other residues and provide stability to the native state. It might be due to high packing effect, as packing is one of the most important characters for the enhanced stability of thermophilic proteins (Gromiha et al., 1999b
).
We have examined the thermal stability of several (
/ß)8 barrel proteins using the thermodynamic database for proteins and mutants developed by us and available in web (Gromiha et al., 2000
). We found that Taka amylase has the highest thermostability of 40.8 kcal/mol (Ooi and Oobatake, 1988
) followed by 1TIM and 1WSY. Also, we observed a direct relationship between stability (
G) and number of hydrophobic clusters.
| CONCLUSIONS |
|---|
|
|
|---|
/ß)8 barrel proteins revealed the presence of hydrophobic clusters in most of the studied proteins. The residues in hydrophobic clusters are mainly influenced by long-range interactions and these residues have high thermal stability, which may help to resist the local tendency for unfolding. Further, a network of long-range interactions is observed to link through ß-strands of (
/ß)8 barrel domain. In essence, the hydrophobic clusters and network of long-range contacts pave the way for the folding and stabilization of (
/ß)8 barrel fold. | ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on April 24, 2002; accepted for publication November 13, 2002.
| REFERENCES |
|---|
|
|
|---|
Chothia, C. 1988. Protein structure: the 14th barrel rolls out. Nature. 333:598599.[Medline]
Debe, D. A., and W. A. Goddard. 1999. First principles prediction of protein folding rates. J. Mol. Biol. 294:619625.[Medline]
Dobson, C. M., and M. Karplus. 1999. The fundamentals of protein folding: bringing together theory and experiment. Curr. Opin. Struct. Biol. 9:92101.[Medline]
Dosztanyi, Z., A. Fiser, and I. Simon. 1997. Stabilization centers in proteins: identification, characterization and predictions. J. Mol. Biol. 272:597612.[Medline]
Farber, G. K. 1993. An (
/ß) barrel full of evolutionary trouble. Curr. Opin. Str. Biol. 3:409412.
Farber, G. K., and G. A. Petsko. 1990. The evolution of
/ß barrel enzymes. Trends Biochem. Sci. 15:228233.[Medline]
Gromiha, M. M. 2001. Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins. Biophys. Chem. 91:7177.[Medline]
Gromiha, M. M., and S. Selvaraj. 2002. Important amino acid properties for determining the transition state structures of two-state protein mutants. FEBS Lett. 526:129134.[Medline]
Gromiha, M. M., and S. Selvaraj. 2001. Comparison between long-range interactions and contact order in determining the folding rates of two-state proteins: application of long-range order to folding rate prediction. J. Mol. Biol. 310:2732.[Medline]
Gromiha, M. M., and S. Selvaraj. 2000. Inter-residue interactions in the structure, folding and stability of proteins. Recent Res. Devel. Biophys. Chem. 1:114.
Gromiha, M. M., and S. Selvaraj. 1999. Importance of long range interactions in protein folding. Biophys. Chem. 77:4968.[Medline]
Gromiha, M. M., and S. Selvaraj. 1998. Protein secondary structure prediction in different structural classes. Protein Eng. 11:249251.
Gromiha, M. M., and S. Selvaraj. 1997a. Influence of medium and long range interactions in different structural classes of globular proteins. J. Biol. Phys. 23:151162.
Gromiha, M. M., and S. Selvaraj. 1997b. Influence of medium and long range interactions in (
/ß)8 barrel proteins. J. Biol. Phys. 23:209217.
Gromiha, M. M., and A. M. Thangakani. 2001. Role of medium- and long-range interactions to the stability of the mutants of T4 lysozyme. Prep. Biochem. Biotech. 31:217227.[Medline]
Gromiha, M. M., M. Oobatake, H. .H. Kono, H. Uedaira, and A. Sarai. 1999a. Role of structural and sequence information in the prediction of protein stability changes: comparison between buried and partially buried mutations. Protein Eng. 12:549555.
Gromiha, M. M., M. Oobatake, and A. Sarai. 1999b. Important amino acid properties for enhanced thermostability from mesophilic to thermophilic proteins. Biophys. Chem. 82:5167.[Medline]
Gromiha, M. M., J. An, H. Kono, M. Oobatake, H. Uedaira, P. Prabakaran, and A. Sarai. 2000. ProTherm, version 2.0: thermodynamic database for proteins and mutants. Nucleic Acids Res. 28:283285.
Holm, L., and C. Sander. 1993. Protein structure comparison by alignment of distance matrices. J. Mol. Biol. 233:123138.[Medline]
Jiang, Z., L. Zhang, J. Chen, A. Xia, and D. Zhao. 2002. Effect of amino acid on forming residue-residue contacts in proteins. Polymer. 43:60376047.
Jones, D. D. 1975. Amino acid properties and side-chain orientation in proteins: a cross correlation approach. J. Theor. Biol. 50:167183.[Medline]
Kabsch, W., and C. Sander. 1993. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 22:25772637.
Kannan, N., S. Selvaraj, M. M. Gromiha, and S. Vishveshwara. 2001. Clusters in (
/ß)8 barrel proteins: implications for protein structure, function and folding: a graph theoretical approach. Prot. Struct. Funct. Genet. 43:103112.
Manavalan, P., and P. K. Ponnuswamy. 1978. Hydrophobic character of amino acid residues in globular proteins. Nature. 275:673674.[Medline]
Manavalan, P., and P. K. Ponnuswamy. 1977. A study of the preferred environment of amino acid residues in globular proteins. Arch. Biochem. Biophys. 184:476487.[Medline]
Miyazawa, S., and R. L. Jernigan. 1999. Evaluation of short-range interactions as secondary structure energies for protein fold and sequence recognition. Proteins. 36:347356.[Medline]
Muthusamy, R., M. M. Gromiha, and P. K. Ponnuswamy. 2000. On the thermal unfolding character of globular proteins. J. Prot. Chem. 19:18.[Medline]
Nagano, N., C. A. Orengo, and J. M. Thornton. 2002. One fold with many functions: the evolutionary relationships between TIM barrel families based on their sequences, structures and functions. J. Mol. Biol. 321:741765.[Medline]
Nozaki, Y., and C. Tanford. 1971. The solubility of amino acids and two glycine peptides in aqueous ethanol and dioxane solutions. Establishment of a hydrophobicity scale. J. Biol. Chem. 246:22112217.
Ooi, T., and M. Oobatake. 1988. Intermolecular interactions between protein and other molecules including hydration effects. J. Biochem. 104:440444.
Orengo, C. A., F. M. Pearl, J. E. Bray, A. E. Todd, A. C. Martin, L. Lo Conte, and J. M. Thornton. 1999. The CATH Database provides insights into protein structure/function relationships. Nucleic Acids Res. 27:275279.
Parker, M. J., R. B. Sessions, I. G. Badcoe, and A. R. Clarke. 1996. The development of tertiary interactions during the folding of a large protein. Fold. Des. 1:145156.[Medline]
Ponnuswamy, P. K. 1993. Hydrophobic characteristics of folded proteins. Prog. Biophys. Mol. Biol. 59:57103.[Medline]
Ponnuswamy, P. K., and M. M. Gromiha. 1994. On the conformational stability of folded proteins. J. Theor. Biol. 166:6374.
Ponnuswamy, P. K., and M. M. Gromiha. 1993. Prediction of transmembrane helices from hydrophobic characteristics of proteins. Int. J. Pept. Protein Res. 42:326341.[Medline]
Ponnuswamy, P. K., and M. Prabhakaran. 1980. Properties of nucleation sites in globular proteins. Biochem. Biophys. Res. Commun. 97:15821590.[Medline]
Ponnuswamy, P. K., R. Muthusamy, and P. Manavalan. 1982. Amino acid composition and thermal stability of globular proteins. Int. J. Biol. Macromol. 4:186190.
Poupon, A., and J. P. Mornon. 1999. Predicting the protein folding nucleus from sequences. FEBS Lett. 452:283289.[Medline]
Pujadas, G. 2002. TIM barrel research in the genomic and proteomic era. In Recent Research Developments in Protein Folding, Stability and Design. M. M. Gromiha, and S. Selvaraj, editors. Research Signpost, Trivandrum, India. pp3755.
Pujadas, G., and J. Palau. 1999. TIM barrel fold: structural, functional and evolutionary characteristics in natural and designed molecules. Biologia Bratislava. 54:231254.
Selvaraj, S., and M. M. Gromiha. 1998a. An analysis of amino acid clustering pattern in (
/ß)8 barrel proteins. J. Prot. Chem. 17:407415.[Medline]
Selvaraj, S., and M. M. Gromiha. 1998b. Importance of long-range interactions in (
/ß)8 barrel fold. J. Prot. Chem. 17:691697.[Medline]
Sergeev, Y., and B. Lee. 1994. Alignment of beta-barrels in (ß/
)8 proteins using hydrogen-bonding pattern. J. Mol. Biol. 244:168182.[Medline]
Sheinerman, F. B., and C. L. Brooks 3rd. 1998. Molecular picture of folding of a small alpha/beta protein. Proc. Natl. Acad. Sci. USA. 95:15621567.
Shindyalov, I. N., and P. E. Bourne. 1998. Protein structure alignment by incremental combinatorial extension (CE) of the optimal path. Protein Eng. 11:739747.
Unger, R., and J. Moult. 1996. Local interactions dominate folding in a simple protein model. J. Mol. Biol. 259:988994.[Medline]
Wolf, Y. I., S. E. Brenner, P. A. Bash, and E. V. Koonin. 1999. Distribution of protein folds in the three superkingdoms of life. Genome Res. 9:1726.
This article has been cited by other articles:
![]() |
A. Paiardini, F. Bossa, and S. Pascarella Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes Protein Sci., November 1, 2004; 13(11): 2992 - 3005. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |