Turn scanning by site-directed mutagenesis: application to the protein folding problem using the intestinal fatty acid binding protein

We have systematically mutated residues located in turns between beta-strands of the intestinal fatty acid binding protein (IFABP), and a glycine in a half turn, to valine and have examined the stability, refolding rate constants and ligand dissociation constants for each mutant protein. IFABP is an almost all beta-sheet protein exhibiting a topology comprised of two five-stranded sheets surrounding a large cavity into which the fatty acid ligand binds. A glycine residue is located in seven of the eight turns between the antiparallel beta-strands and another in a half turn of a strand connecting the front and back sheets. Mutations in any of the three turns connecting the last four C-terminal strands slow the folding and decrease stability with the mutation between the last two strands slowing folding dramatically. These data suggest that interactions between the last four C-terminal strands are highly cooperative, perhaps triggered by an initial hydrophobic collapse. We suggest that this trigger is collapse of the highly hydrophobic cluster of amino acids in the D and E strands, a region previously shown to also affect the last stage of the folding process (Kim et al., 1997). Changing the glycine in the strand between the front and back sheets also results in a unstable, slow folding protein perhaps disrupting the D-E strand interactions. For most of the other turn mutations there was no apparent correlation between stability and refolding rate constants. In some turns, the interaction between strands, rather than the turn type, appears to be critical for folding while in others, turn formation itself appears to be a rate limiting step. Although there is no simple correlation between turn formation and folding kinetics, we propose that turn scanning by mutagenesis will be a useful tool for issues related to protein folding.     

Unique stabilizing interactions identified in the two-stranded alpha-helical coiled-coil: crystal structure of a cortexillin I/GCN4 hybrid coiled-coil peptide

We determined the 1.17 A resolution X-ray crystal structure of a hybrid peptide based on sequences from coiled-coil regions of the proteins GCN4 and cortexillin I. The peptide forms a parallel homodimeric coiled-coil, with C(alpha) backbone geometry similar to GCN4 (rmsd value 0.71 A). Three stabilizing interactions have been identified: a unique hydrogen bonding-electrostatic network not previously observed in coiled-coils, and two other hydrophobic interactions involving leucine residues at positions e and g from both g-a' and d-e' interchain interactions with the hydrophobic core. This is also the first report of the quantitative significance of these interactions. The GCN4/cortexillin hybrid surprisingly has two interchain Glu-Lys' ion pairs that form a hydrogen bonding network with the Asn residues in the core. This network, which was not observed for the reversed Lys-Glu' pair in GCN4, increases the combined stability contribution of each Glu-Lys' salt bridge across the central Asn15-Asn15' core to approximately 0.7 kcal/mole, compared to approximately 0.4 kcal mole(-1) from a Glu-Lys' salt bridge on its own. In addition to electrostatic and hydrogen bonding stabilization of the coiled-coil, individual leucine residues at positions e and g in the hybrid peptide also contribute to stability by 0.7 kcal/mole relative to alanine. These interactions are of critical importance to understanding the stability requirements for coiled-coil folding and in modulating the stability of de novo designed macromolecules containing this motif.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.     

Early formation of a beta hairpin during folding of staphylococcal nuclease H124L as detected by pulsed hydrogen exchange

Pulsed hydrogen exchange methods were used to follow the formation of structure during the refolding of acid-denatured staphylococcal nuclease containing a stabilizing Leu substitution at position 124 (H124L SNase). The protection of more than 60 backbone amide protons in uniformly (15)N-labeled H124L SNase was monitored as a function of refolding time by heteronuclear two-dimensional NMR spectroscopy. As found in previous studies of staphylococcal nuclease, partial protection was observed for a subset of amide protons even at the earliest folding time point (10 msec). Protection indicative of marginally stable hydrogen-bonded structure in an early folding intermediate was observed at over 30 amide positions located primarily in the beta-barrel and to a lesser degree in the alpha-helical domain of H124L SNase. To further characterize the folding intermediate, protection factors for individual amide sites were measured by varying the pH of the labeling pulse at a fixed refolding time of 16 msec. Protection factors >5.0 were observed only for amide positions in a beta-hairpin formed by strands 2 and 3 of the beta-barrel domain and a single site near the C-terminus. The results indicate that formation of stable hydrogen-bonded structure in a core region of the beta-sheet is among the earliest structural events in the folding of SNase and may serve as a nucleation site for further structure formation.     

Comparison of the hemocyanin beta-barrel with other Greek key beta-barrels: possible importance of the "beta-zipper" in protein structure and folding.

The Greek key beta-barrel topology is a folding motif observed in many proteins of widespread evolutionary origin. The arthropodan hemocyanins also have such a Greek key beta-barrel, which forms the core of the third domain of this protein. The hemocyanin beta-barrel was found to be structurally very similar to the beta-barrels of the immunoglobulin domains, Cu,Zn-superoxide dismutase and the chromophore carrying antitumor proteins. The structural similarity within this group of protein families is not accompanied by an evolutionary or functional relationship. It is therefore possible to study structure-sequence relations without bias from nonstructural constraints. The present study reports a conserved pattern of features in these Greek key beta-barrels that is strongly suggestive of a folding nucleation site. This proposed nucleation site, which we call a "beta-zipper," shows a pattern of well-conserved, large hydrophobic residues on two sequential beta-strands joined by a short loop. Each beta-zipper strand is near the center of one of the beta-sheets, so that the two strands face each other from opposite sides of the barrel and interact through their hydrophobic side chains, rather than forming a hydrogen-bonded beta-hairpin. Other protein families with Greek key beta-barrels that do not as strongly resemble the immunoglobulin fold--such as the azurins, plastocyanins, crystallins, and prealbumins--also contain the beta-zipper pattern, which might therefore be a universal feature of Greek key beta-barrel proteins.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.     

Salt effects on hydrophobic‐core formation in folding of a helical miniprotein studied by molecular dynamics simulations

We have investigated effects of salt ions on folding events of a helical miniprotein chicken villin headpiece subdomain HP36. Low concentrations of ions alter electrostatic interactions between charged groups of a protein and can change the populations of conformers. Here, we compare two data sets of folding simulations of HP36 in explicit water solvent with or without ions. For efficient sampling of the conformational space of HP36, the multicanonical replica‐exchange molecular dynamics method was employed. Our analyses suggest that salt alters salt‐bridging nature of the protein at later stages of folding at room temperature. Especially, more nonnative, nonlocal salt bridges are formed at near‐native conformations in pure water. Our analyses also show that such salt‐bridge formation hinders the fully native hydrophobic‐core packing at the final stages of folding. Proteins 2014; 82:933–943. © 2013 Wiley Periodicals, Inc.     

Optimization of the electrostatic interactions in proteins of different functional and folding type

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed αβ folding type are electrostatically better optimized than pure α-helical or β-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.     

Sequence and structure patterns in proteins from an analysis of the shortest helices: implications for helix nucleation

The shortest helices (three-length 3(10) and four-length alpha), most abundant among helices of different lengths, have been analyzed from a database of protein structures. A characteristic feature of three-length 3(10)-helices is the shifted backbone conformation for the C-terminal residue (phi,psi angles: -95 degrees,0 degrees ), compared to the rest of the helix (-62 degrees,-24 degrees ). The deviation can be attributed to the release of electrostatic repulsion between the carbonyl oxygen atoms at the two C-terminal residues and further stabilization (due to a more linear geometry) of an intrahelical hydrogen bond. A consequence of this non-canonical C-terminal backbone conformation can be a potential origin of helix kinks when a 3(10)-helix is sequence-contiguous at the alpha-helix N-terminal. An analysis of hydrogen bonding, as well as hydrophobic interactions in the shortest helices shows that capping interactions, some of them not observed for longer helices, dominate at the N termini. Further, consideration of the distribution of amino acid residues indicates that the shortest helices resemble the N-terminal end of alpha-helices rather than the C terminus, implying that the folding of helices may be initiated at the N-terminal end, which does not get propagated in the case of the shortest helices. Finally, pairwise comparison of beta-turns and the shortest helices, based on correlation matrices of site-specific amino acid composition, and the relative abundance of these short secondary structural elements, leads to a helix nucleation scheme that considers the formation of an isolated beta-turn (and not an alpha-turn) as the helix nucleation step, with shortest 3(10)-helices as intermediates between the shortest alpha-helix and the beta-turn. Our results ascribe an important role played by shortest 3(10)-helices in proteins with important structural and folding implications.     

Intra- versus intermolecular interactions in monellin: contribution of surface charges to protein assembly

The relative significance of weak non-covalent interactions in biological context has been much debated. Here, we have addressed the contribution of Coulombic interactions to protein stability and assembly experimentally. The sweet protein monellin, a non-covalently linked heterodimeric protein, was chosen for this study because of its ability to spontaneously reconstitute from separated fragments. The reconstitution of monellin mutants containing large surface charge perturbations was compared to the thermostability of structurally equivalent single-chain monellin containing the same sets of mutations under varying salt concentrations. The affinity between monellin fragments is found to correlate with the thermostability of single chain monellin, indicating the involvement of the same underlying Coulombic interactions. This confirms that there are no principal differences in the interactions involved in folding and binding. Based on comparison with a previous mutational study involving hydrophobic core residues, the relative contribution of Coulombic interactions to stability and affinity is modest. However, the Coulombic perturbations only affect the association rates of reconstitution in contrast to perturbations involving hydrophobic residues, which affect primarily the dissociation rates. These results indicate that Coulombic interactions are likely to be of main importance for the association of protein assembly, relevant for functions of proteins.     

The role of aromatic residues in the hydrophobic core of the villin headpiece subdomain

Small autonomously folding proteins are of interest as model systems to study protein folding, as the same molecule can be used for both experimental and computational approaches. The question remains as to how well these minimized peptide model systems represent larger native proteins. For example, is the core of a minimized protein tolerant to mutation like larger proteins are? Also, do minimized proteins use special strategies for specifying and stabilizing their folded structure? Here we examine these questions in the 35‐residue autonomously folding villin headpiece subdomain (VHP subdomain). Specifically, we focus on a cluster of three conserved phenylalanine (F) residues F47, F51, and F58, that form most of the hydrophobic core. These three residues are oriented such that they may provide stabilizing aromatic–aromatic interactions that could be critical for specifying the fold. Circular dichroism and 1D‐NMR spectroscopy show that point mutations that individually replace any of these three residues with leucine were destabilized, but retained the native VHP subdomain fold. In pair‐wise replacements, the double mutant that retains F58 can adopt the native fold, while the two double mutants that lack F58 cannot. The folding of the double mutant that retains F58 demonstrates that aromatic–aromatic interactions within the aromatic cluster are not essential for specifying the VHP subdomain fold. The ability of the VHP subdomain to tolerate mutations within its hydrophobic core indicates that the information specifying the three dimensional structure is distributed throughout the sequence, as observed in larger proteins. Thus, the VHP subdomain is a legitimate model for larger, native proteins.     

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.     

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.     

Mutational analysis of the BPTI folding pathway: I. Effects of aromatic-->leucine substitutions on the distribution of folding intermediates.

The roles of aromatic residues in determining the folding pathway of bovine pancreatic trypsin inhibitor (BPTI) were analyzed mutationally by examining the distribution of disulfide-bonded intermediates that accumulated during the refolding of protein variants in which tyrosine or phenylalanine residues were individually replaced with leucine. The eight substitutions examined all caused significant changes in the intermediate distribution. In some cases, the major effect was to decrease the accumulation of intermediates containing two of the three disulfides found in the native protein, without affecting the distribution of earlier intermediates. Other substitutions, however, led to much more random distributions of the intermediates containing only one disulfide. These results indicate that the individual residues making up the hydrophobic core of the native protein make clearly distinguishable contributions to conformation and stability early in folding: The early distribution of intermediates does not appear to be determined by a general hydrophobic collapse. The effects of the substitutions were generally consistent with the structures of the major intermediates determined by NMR studies of analogs, confirming that the distribution of disulfide-bonded species is determined by stabilizing interactions within the ordered regions of the intermediates. The plasticity of the BPTI folding pathway implied by these results can be described using conformational funnels to illustrate the degree to which conformational entropy is lost at different stages in the folding of the wild-type and mutant proteins.     

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.     

The hydrogen exchange core and protein folding.

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.     

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.     

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.     

Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.