Conserved residues and their role in the structure, function, and stability of acyl-coenzyme A binding protein

In the family of acyl-coenzyme A binding proteins, a subset of 26 sequence sites are identical in all eukaryotes and conserved throughout evolution of the eukaryotic kingdoms. In the context of the bovine protein, the importance of these 26 sequence positions for structure, function, stability, and folding has been analyzed using single-site mutations. A total of 28 mutant proteins were analyzed which covered 17 conserved sequence positions and three nonconserved positions. As a first step, the influence of the mutations on the protein folding reaction has been probed, revealing a folding nucleus of eight hydrophobic residues formed between the N- and C-terminal helices [Kragelund, B. B., et al. (1999) Nat. Struct. Biol. (In press)]. To fully analyze the role of the conserved residues, the function and the stability have been measured for the same set of mutant proteins. Effects on function were measured by the extent of binding of the ligand dodecanoyl-CoA using isothermal titration calorimetry, and effects on protein stability were measured with chemical denaturation followed by intrinsic tryptophan and tyrosine fluorescence. The sequence sites that have been conserved for direct functional purposes have been identified. These are Phe5, Tyr28, Tyr31, Lys32, Lys54, and Tyr73. Binding site residues are mainly polar or charged residues, and together, four of these contribute approximately 8 kcal mol-1 of the total free energy of binding of 11 kcal mol-1. The sequence sites conserved for stability of the structure have likewise been identified and are Phe5, Ala9, Val12, Leu15, Leu25, Tyr28, Lys32, Gln33, Tyr73, Val77, and Leu80. Essentially, all of the conserved residues that maintain the stability are hydrophobic residues at the interface of the helices. Only one conserved polar residue, Gln33, is involved in stability. The results indicate that conservation of residues in homologous proteins may result from a summed optimization of an effective folding reaction, a stable native protein, and a fully active binding site. This is important in protein design strategies, where optimization of one of these parameters, typically function or stability, may influence any of the others markedly.     

Hydrophobic core packing in the SH3 domain folding transition state.

How tightly packed is the hydrophobic core of a folding transition state structure? We have addressed this question by characterizing the effects on folding kinetics of > 40 substitutions of both large and small amino acids in the hydrophobic core of the Fyn SH3 domain. Our results show that residues at three positions, which we designate as the 'core folding nucleus', are tightly packed in the transition state, and substitutions at these positions cause the largest changes in the folding rate. The other six positions examined appear to be loosely packed; thus, substitutions at these positions with larger hydrophobic residues generally accelerate folding, presumably by increasing the rate of nonspecific hydrophobic collapse. Surprisingly, the folding rate can be greatly accelerated by residues that also significantly destabilize the native state structure. Furthermore, mutants with identical thermodynamic stability can differ by up to 55-fold in their folding rates. These results highlight the importance of hydrophobic core composition, as opposed to only topology, in determining the folding rate of a protein. They also provide a new explanation for the 'abnormal' phi-values observed in many protein folding kinetics studies.     

The folding of an enzyme. VI. The folding pathway of barnase: comparison with theoretical models.

The sequence of events in the refolding pathway of barnase has been analysed to search for general principles in protein folding. There appears to be a correlation between burying hydrophobic surface area and early folding events. All the regions that fold early interact extensively with the beta-sheet. These interactions involve predominantly hydrophobic interactions and the burial of very extensive hydrophobic areas in which multiple, close, hydrophobic-hydrophobic contacts are established around a central group of aliphatic residues. There is no burial of hydrophilic residues in these regions; those that are partly screened from the solvent make hydrogen bonds. All the regions or interactions that are made late in the folding pathway do not make extensive contacts with the beta-sheet. Their buried hydrophobic regions lack a central hydrophobic residue or residues around which other hydrophobic residues pack. Further, in some of these regions there is an extensive burial of hydrophilic residues. The results are consistent with one of the earlier events in protein folding being the local formation of native-like secondary structure elements driven by local hydrophobic surface burial. A possible candidate for an initiation site is a beta-hairpin between beta-strands 3 and 4 that is conserved in the microbial ribonuclease family. A comparison of structures in this family shows that those regions that can be superimposed, or have sequence homology, correspond to elements of structure that are formed and interact with each other early in the folding pathway, suggesting that some of these residues could be involved in directing the folding process. The data on barnase combined with results from other laboratories suggest the following tentative conclusions for the refolding of small monomeric proteins. (1) The refolding pathway is, at least in part, sequential and of compulsory order. (2) Secondary structure formation is driven by local hydrophobic surface burial and precedes the formation of most tertiary interactions. These elements are then stabilized and sometimes elongated by tertiary interactions. It is plausible that there are stop signals encoded in the linear sequence that prevent the elongation of isolated secondary structure elements in solution to a larger extent than is found in the folded protein. (3) Many tertiary interactions are not very constrained in the intermediate but become more and more defined as the hydrophobic cores consolidate, loop structures form and the configuration of surface residues takes place. The interactions between different elements of secondary structure are the last ones to be consolidated while the interactions within the secondary structure elements are consolidated earlier.(ABSTRACT TRUNCATED AT 400 WORDS)     

Localized, stereochemically sensitive hydrophobic packing in an early folding intermediate of dihydrofolate reductase from Escherichia coli

Mutational analysis was performed to probe the development of hydrophobic clusters during the early events in the folding of dihydrofolate reductase. Replacements were made in several hydrophobic subdomains to examine the roles of hydrophobicity and stereochemistry in the formation of kinetic intermediates. Amide protons in two of these clusters, including residues I91, I94, and I155, have been shown to be protected against solvent exchange within 13 ms of folding. Additional hydrophobic clusters were probed by substitutions at residues I2, I61, and L112; these residues are not protected from exchange until later in the folding reaction. Valine and leucine replacements at positions I91, I94, and I155 significantly diminish the formation of the burst phase kinetic intermediate, relative to the wild-type protein. In contrast, I2 and I61 are insensitive to these substitutions in the first 5 ms of the folding reaction, as is the replacement of L112 with either isoleucine or valine. These results demonstrate that the tightly packed core of dihydrofolate reductase is acquired in a non-uniform fashion, beginning in the submillisecond time frame. The progressive development of specific side-chain packing in localized hydrophobic clusters may be a common theme for complex protein folding reactions.     

Native protein sequences are designed to destabilize folding intermediates

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.     

Frequencies of amino acid strings in globular protein sequences indicate suppression of blocks of consecutive hydrophobic residues

Patterns of hydrophobic and hydrophilic residues play a major role in protein folding and function. Long, predominantly hydrophobic strings of 20-22 amino acids each are associated with transmembrane helices and have been used to identify such sequences. Much less attention has been paid to hydrophobic sequences within globular proteins. In prior work on computer simulations of the competition between on-pathway folding and off-pathway aggregate formation, we found that long sequences of consecutive hydrophobic residues promoted aggregation within the model, even controlling for overall hydrophobic content. We report here on an analysis of the frequencies of different lengths of contiguous blocks of hydrophobic residues in a database of amino acid sequences of proteins of known structure. Sequences of three or more consecutive hydrophobic residues are found to be significantly less common in actual globular proteins than would be predicted if residues were selected independently. The result may reflect selection against long blocks of hydrophobic residues within globular proteins relative to what would be expected if residue hydrophobicities were independent of those of nearby residues in the sequence.     

Visualization of the nature of protein folding by a study of a distance constraint approach in two-dimensional models

A distance constraint approach is applied to two-dimensional models of proteins in order to visualize the nature of protein folding and to examine the relative roles of different ranges of interaction. Three different native structures (I, II, and III) are considered; they have two different kinds of residues, viz., hydrophobic and hydrophilic, and different sequences of these residues. We examine how the distance constraint approach functions in the prediction of protein folding when we know the sequence of the residues, the (fixed) bond lengths, the mean distances between residues i and i + 2, and i and i + 3, and the mean distances for hydrophobic–hydrophobic, hydrophobic–hydrophilic, and hydrophilic–hydrophilic contacts between residues i and i + j, where j ≥ 4. This approach involves optimization of an object function with respect to 98 variables and is not free of the multiple-minimum problem. The optimization is always terminated if the chain is entangled and/or the segments (residues) are packed too compactly to move. In order to escape from such situations and to take the excluded-volume effect into account, a Monte Carlo method is used after the optimization is trapped in local minima. Success in the prediction of folding is found to depend on the starting conformations and on the native conformations. Fair success is obtained in predicting the helix-like structure in protein I and the overall structure of protein III, but not the β-like structures of proteins I and II. Insofar as the prediction of the structure of protein III is reasonable, it appears that some sequences of residues produce greater constraints on their conformations than others, if one considers only the hydrophobic and hydrophilic nature of the residues. These results imply that, in the folding of real proteins in three dimensions, the competition for hydrophobic (and hydrophilic) residues for inside (outside) positions in the molecule probably constitutes a necessary but not a sufficient condition to form and stabilize the native structure. The failure to predict the structure of protein II, and part of that of protein I, suggests that there are two types of long-range interactions. One (which we considered here) is nonspecific (i.e., is defined only in terms of contacts between residues of the same or different polarity) and acts at any stage of protein folding; the other (which we did not consider here) is a specific interaction between residues in pairs and contributes only when the residues in the specific pair take on the native conformation. Presumably, incorporation of such specific long-range interactions, together with the nonspecific ones, is necessary for successful protein folding, using the distance constraint approach.     

Differential stabilization of two hydrophobic cores in the transition state of the villin 14T folding reaction

We report the distribution of hydrophobic core contacts during the folding reaction transition state for villin 14T, a small 126-residue protein domain. The solution structure of villin 14T contains a central beta-sheet with two flanking hydrophobic cores; transition states for this protein topology have not been previously studied. Villin 14T has no disulfide bonds or cis-proline residues in its native state; it folds reversibly, and in an apparently two-state manner under some conditions. To map the hydrophobic core contacts in the transition state, 27 point mutations were generated at positions spread throughout the two hydrophobic cores. After each point mutation, comparison of the change in folding kinetics with the equilibrium destabilization indicates whether the site of mutation is stabilized in the transition state. The results show that the folding nucleus, or the sub-region with the strongest transition state contacts, is located in one of the two hydrophobic cores (the predominantly aliphatic core). The other hydrophobic core, which is mostly aromatic, makes much weaker contacts in the transition state. This work is the first transition state mapping for a protein with multiple major hydrophobic cores in a single folding unit; the hydrophobic cores cannot be separated into individual folding subdomains. The stabilization of only one hydrophobic core in the transition state illustrates that hydrophobic core formation is not intrinsically capable of nucleating folding, but must also involve the right specific interactions or topological factors in order to be kinetically important.     

Non-functional conserved residues in globins and their possible role as a folding nucleus.

Structure-based sequence alignment of 728 sequences of different globin subfamilies shows that in each subfamily there are two clusters of consensually conserved residues. The first is the well-known "functional" cluster which includes six heme-binding conserved residues (Phe CD1, His F8; aliphatic E11, FG5; hydrophobic F4, G5) and seven other conserved residues (Pro C2; aliphatic H19; hydrophobic B10, B13, B14, CD4, E4) that do not bind the heme but belong to its immediate neighborhood. The second cluster revealed here (aliphatic A8, G16, G12; aromatic A12; hydrophobic H8 and possibly H12) is distant from the heme. It is entirely non-polar and includes one turn (i, i+4 positions) from each of helices A, G, and H. It is known that A, G, and H helices formed at the earliest stage of apomyoglobin folding remain relatively stable in the equilibrium molten globule state, and are likely to be tightly packed with each other in this state. We have shown the existence of two similar conserved clusters in c -type cytochromes, heme-binding and distal from the heme. The second cluster in c -cytochromes includes one turn from each of the N and C-terminal alpha-helices. These N and C-terminal helices in cytochrome c are formed at the earliest stage of protein folding, remain relatively stable in the molten globule state, and are tightly packed with each other in this state, similar to the observed behavior of the globins. At least these two large protein families (c -type cytochromes and globins) have a close similarity in the existence and mutual positions of non-functional conserved residues. We assume that non-functional conserved residues are requisite for the fast and correct folding of both of these protein families into their stable 3D structures.     

Sequence evolution and the mechanism of protein folding

The impact on protein evolution of the physical laws that govern folding remains obscure. Here, by analyzing in silico-evolved sequences subjected to evolutionary pressure for fast folding, it is shown that: First, a subset of residues in the thermodynamic folding nucleus is mainly responsible for modulating the protein folding rate. Second and most important, the protein topology itself is of paramount importance in determining the location of these residues in the structure. Further stabilization of the interactions in this nucleus leads to fast folding sequences. Third, these nucleation points restrict the sequence space available to the protein during evolution. Correlated mutations between positions around these hot spots arise in a statistically significant manner, and most involve contacting residues. When a similar analysis is carried out on real proteins, qualitatively similar results are obtained.     

Single amino acid substitutions globally suppress the folding defects of temperature-sensitive folding mutants of phage P22 coat protein.

The amino acid sequence of a polypeptide defines both the folding pathway and the final three-dimensional structure of a protein. Eighteen amino acid substitutions have been identified in bacteriophage P22 coat protein that are defective in folding and cause their folding intermediates to be substrates for GroEL and GroES. These temperature-sensitive folding (tsf) substitutions identify amino acids that are critical for directing the folding of coat protein. Additional amino acid residues that are critical to the folding process of P22 coat protein were identified by isolating second site suppressors of the tsf coat proteins. Suppressor substitutions isolated from the phage carrying the tsf coat protein substitutions included global suppressors, which are substitutions capable of alleviating the folding defects of numerous tsf coat protein mutants. In addition, potential global and site-specific suppressors were isolated, as well as a group of same site amino acid substitutions that had a less severe phenotype than the tsf parent. The global suppressors were located at positions 163, 166, and 170 in the coat protein sequence and were 8-190 amino acid residues away from the tsf parent. Although the folding of coat proteins with tsf amino acid substitutions was improved by the global suppressor substitutions, GroEL remained necessary for folding. Therefore, we believe that the global suppressor sites identify a region that is critical to the folding of coat protein.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

Prediction of folding mechanisms for Ig-like beta sandwich proteins based on inter-residue average distance statistics methods

To understand the folding mechanism of a protein is one of the goals in bioinformatics study. Nowadays, it is enigmatic and difficult to extract folding information from amino acid sequence using standard bioinformatics techniques or even experimental protocols which can be time consuming. To overcome these problems, we aim to extract the initial folding unit for titin protein (Ig and fnIII domains) by means of inter-residue average distance statistics, Average Distance Map (ADM) and contact frequency analysis (F-value). TI I27 and TNfn3 domains are used to represent the Ig-domain and fnIII-domain, respectively. Beta-strands 2, 3, 5, and 6 are significant for the initial folding processes of TI I27. The central strands of TNfn3 were predicted as a primary folding segment. Known 3D structure and unknown 3D structure domains were investigated by structure or non-structure based multiple sequence alignment, respectively, to learn the conserved hydrophobic residues and predicted compact region relevant to evolution. Our results show good correspondence to experimental data, phi-value and protection factor from H-D exchange experiments. The significance of conserved hydrophobic residues near F-value peaks for structural stability using hydrophobic packing is confirmed. Our prediction methods once again could extract a folding mechanism only knowing the amino acid sequence.     

Folding behavior of chaperonin-mediated substrate protein

Chaperonin-mediated protein folding is complex. There have been diverse results on folding behavior, and the chaperonin molecules have been investigated as enhancing or retarding the folding rate. To understand the diversity of chaperonin-mediated protein folding, we report a study based on simulations using a simplified Gō-type model. By considering effects of affinity between the substrate protein and the chaperonin wall and spatial confinement of the chaperonin cavity, we study the thermodynamics and kinetics of folding of an unfrustrated substrate protein encapsulated in a chaperonin cavity. The affinity makes the hydrophobic residues of the protein bind to the chaperonin wall, and a strong (or weak) affinity results in a large (or small) effect of binding. Compared with the folding in bulk, the folding in chaperonin cavity with different strengths of affinity shows two kinds of behaviors: one with less dependence on the affinity but more reliance on the spatial confinement effect and the other relying strongly on the affinity. It is found that the enhancement or retardation of the folding rate depends on the competition between the spatial confinement and the affinity due to the chaperonin cavity, and a strong affinity produces a slow folding while a weak affinity induces a fast folding. The crossover between two kinds of folding behaviors happens in the case that the favorable effect of confinement is balanced by the unfavorable effect of the affinity, and a critical affinity strength is roughly defined. By analyzing the contacts formed between the residues of the protein and the chaperonin wall and between the residues of the protein themselves, the role of the affinity in the folding processes is studied. The binding of the residues with the chaperonin wall reduces the formation of both native contacts and nonnative contact or mis-contacts, providing a loose structure for further folding after allosteric change of the chaperonin cavity. In addition, 15 single-site-mutated mutants are simulated in order to test the validity of our model and to investigate the importance of affinity. Inspiringly, our results of the folding rates have a good correlation with those obtained from experiments. The folding rates are inversely correlated with the strength of the binding interactions, i.e., the weaker the binding, the faster the folding. We also find that the inner hydrophobic residues have larger effects on the folding kinetics than those of the exterior hydrophobic residues. We suggest that, besides the confinement effect, the affinity acts as another important factor to affect the folding of the substrate proteins in chaperonin systems, providing an understanding of the folding mechanism of the molecular chaperonin systems.     

Sequence–function relationships in folding upon binding

Folding coupled to binding is ubiquitous in biology. Nevertheless, the relationship of sequence to function for protein segments that undergo coupled binding and folding remains to be determined. Specifically, it is not known if the well-established rules that govern protein folding and stability are relevant to ligand-linked folding transitions. Upon small ligand biotinoyl-5′-AMP (bio-5′-AMP) binding the Escherichia coli protein BirA undergoes a disorder-to-order transition that results in formation of a network of packed hydrophobic side chains. Ligand binding is also allosterically coupled to protein association, with bio-5′-AMP binding enhancing the dimerization free energy by −4.0 kcal/mol. Previous studies indicated that single alanine replacements in a three residue hydrophobic cluster that contributes to the larger network disrupt cluster formation, ligand binding, and allosteric activation of protein association. In this work, combined equilibrium and kinetic measurements of BirA variants with alanine substitutions in the entire hydrophobic network reveal large functional perturbations resulting from any single substitution and highly non-additive effects of multiple substitutions. These substitutions also disrupt ligand-linked folding. The combined results suggest that, analogous to protein folding, functional disorder-to-order linked to binding requires optimal packing of the relevant hydrophobic side chains that contribute to the transition. The potential for many combinations of residues to satisfy this requirement implies that, although functionally important, segments of homologous proteins that undergo folding linked to binding can exhibit sequence divergence.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

Topology of amphipathic motifs mediating Golgi localization in ArfGAP1 and its splice isoforms

The interaction of the Arf1-directed GTPase-activating protein ArfGAP1 with the Golgi apparatus depends on motifs in its noncatalytic part that are unstructured in solution but are capable of folding into amphipathic helices in vitro upon interaction with poorly packed lipids. In previous studies a few hydrophobic residues that are critical for lipid binding and Golgi localization were identified, but the precise topology of the amphipathic motifs has not been determined. Here we present a detailed analysis of the Golgi targeting and in vitro folding features of the region encompassing the amphipathic motifs (residues 199-294). Point mutation analysis revealed that most hydrophobic residues within this region contribute to Golgi localization, whereas analysis by proline replacements and alanine insertions revealed that Golgi interaction depends on folding into two amphipathic helices with a short interrupting sequence. Analysis of splice isoforms containing 10-residue in-frame insertions within their first amphipathic motifs revealed that the insertion causes a truncation of the amphipathic helix that does not extend beyond the insertion sequence. Lastly, a lysine replacement mutant recently reported to bind to negatively charged liposomes in a curvature-independent manner showed normal cellular distribution, suggesting that Golgi targeting of Arf-GAP1 may involve factors other than sensing lipid packing.     

不同类型氨基酸网络参量与蛋白质折叠的关系

理论和实验研究表明,蛋白质天然拓扑结构对其折叠过程具有重要的影响．采用复杂网络的方法分析蛋白质天然结构的拓扑特征,并探索蛋白质结构特征与折叠速率之间的内在联系．分别构建了蛋白质氨基酸网络、疏水网、亲水网、亲水-疏水网以及相应的长程网络,研究了这些网络的匹配系数(assortativity coefficient)和聚集系数(clustering coefficient)的统计特性．结果表明,除了亲水-疏水网,上述各网络的匹配系数均为正值,并且氨基酸网和疏水网的匹配系数与折叠速率表现出明显的线性正相关,揭示了疏水残基间相互作用的协同性有助于蛋白质的快速折叠．同时,研究发现疏水网的聚集系数与折叠速率有明显的线性负相关关系,这表明疏水残基间三角结构(triangle construction)的形成不利于蛋白质快速折叠．还进一步构建了相应的长程网络,发现序列上间距较远的残基接触对的形成将使蛋白质折叠进程变慢．     

Three topologically equivalent core residues affect the transition state ensemble in a protein folding reaction

The single domain protein, interleukin-1beta, is representative of a distinct class of proteins characterized by their beta-trefoil topology. Each subdomain of this structural class is composed of a beta beta beta loop beta (betabetabetaLbeta) motif comprised of approximately 50 residues and gives the protein a pseudo- 3-fold axis of symmetry. A common feature of proteins in this topological family appears to be that they are slow folders, which reach the native state on the order of tens to 100s of seconds. Sequence analysis of interleukin-1beta indicates that three phenylalanine residues located at positions 42, 101, and 146 are well conserved, separated by approximately 50 residues in the primary sequence, located in similar positions in the pseudo-symmetric units of the trefoil, and are juxtaposed to one another in conformational space. These residues surround the hydrophobic cavity and "pin" the hairpin triplet cap to the core beta-barrel. To determine if cap-barrel interactions are involved in maintaining the structural stability and cooperativity or in controlling the slow formation of the native state, we performed a series of mutational studies. The results indicate that interleukin-1beta tolerates large increases in side-chain volume at these three topologically conserved sites with little effect on stability, while the kinetics show significant differences in both the unfolding and refolding rates. Taken together, our results indicate that these conserved core residues are essential contacts in the transition-state ensemble for folding.