Long-range coupling between separate docking sites in interleukin-1beta

The human cytokine interleukin-1beta (IL-1beta) interacts with the interleukin type I receptor using two large docking surfaces designated A and B. Crystallographic studies reveal that a single histidine residue (His30) in IL-1beta makes critical electrostatic interactions at the receptor/ligand interface. To study the function of this residue at site A, four mutant forms of IL-1beta (H30A, H30D, H30F and H30R) were investigated. The mutation that introduces charge repulsion at His30 destabilizes the protein, but paradoxically causes the least effect on receptor binding (H30D). Mutations that enhance hydrophobic or electrostatic interactions have little effect on protein stability yet markedly impair receptor binding (H30F, H30R). All mutations can transmit effects from site A to site B, as evidenced by changes in the binding of a single-chain antibody highly specific for site B. Dihedral scalar coupling constants for the wild-type IL-1beta and the four His mutant proteins showed changes in backbone angles in residues located around site B, some approximately 30 angstroms away from His30 in site A. A comparison of native solvent exchange in wild-type and mutated IL-1beta shows transmission of local destabilization along the hydrogen bond network of the beta-sheet. Taken together, the data indicate that a single residue in site A of IL-1beta can impact stability and function through perturbations in both local and long-range contacts.     

Identification of a key structural element for protein folding within beta-hairpin turns

Specific residues in a polypeptide may be key contributors to the stability and foldability of the unique native structure. Identification and prediction of such residues is, therefore, an important area of investigation in solving the protein folding problem. Atypical main-chain conformations can help identify strains within a folded protein, and by inference, positions where unique amino acids may have a naturally high frequency of occurrence due to favorable contributions to stability and folding. Non-Gly residues located near the left-handed alpha-helical region (L-alpha) of the Ramachandran plot are a potential indicator of structural strain. Although many investigators have studied mutations at such positions, no consistent energetic or kinetic contributions to stability or folding have been elucidated. Here we report a study of the effects of Gly, Ala and Asn substitutions found within the L-alpha region at a characteristic position in defined beta-hairpin turns within human acidic fibroblast growth factor, and demonstrate consistent effects upon stability and folding kinetics. The thermodynamic and kinetic data are compared to available data for similar mutations in other proteins, with excellent agreement. The results have identified that Gly at the i+3 position within a subset of beta-hairpin turns is a key contributor towards increasing the rate of folding to the native state of the polypeptide while leaving the rate of unfolding largely unchanged.     

Identification of kinetically hot residues in proteins.

A number of recent studies called attention to the presence of kinetically important residues underlying the formation and stabilization of folding nuclei in proteins, and to the possible existence of a correlation between conserved residues and those participating in the folding nuclei. Here, we use the Gaussian network model (GNM), which recently proved useful in describing the dynamic characteristics of proteins for identifying the kinetically hot residues in folded structures. These are the residues involved in the highest frequency fluctuations near the native state coordinates. Their high frequency is a manifestation of the steepness of the energy landscape near their native state positions. The theory is applied to a series of proteins whose kinetically important residues have been extensively explored: chymotrypsin inhibitor 2, cytochrome c, and related C2 proteins. Most of the residues previously pointed out to underlie the folding process of these proteins, and to be critically important for the stabilization of the tertiary fold, are correctly identified, indicating a correlation between the kinetic hot spots and the early forming structural elements in proteins. Additionally, a strong correlation between kinetically hot residues and loci of conserved residues is observed. Finally, residues that may be important for the stability of the tertiary structure of CheY are proposed.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.     

Conformation of amyloid fibrils of beta2-microglobulin probed by tryptophan mutagenesis

Beta2-microglobulin (beta2-m), a protein responsible for dialysis-related amyloidosis, adopts an immunoglobulin domain fold in its native state. Although beta2-m has Trp residues at positions 60 and 95, both are located near the surface of the domain. Hence, beta2-m does not have a conserved Trp common to other immunoglobulin domains, which is buried in close proximity to the disulfide bond. To study the structure of amyloid fibrils in relation to their native fold, we prepared a series of Trp mutants. Trp60 and Trp95 were both replaced with Phe, and a single Trp was introduced at various positions. Among various mutants, W39-beta2-m, in which a Trp was introduced at the position corresponding to the conserved Trp, exhibited a remarkable quenching of fluorescence in the native state, as observed for other immunoglobulin domains. An x-ray structural analysis revealed that W39-beta2-m assumes the native fold with Trp39 located in the vicinity of the disulfide bond. Comparison of the fluorescence spectra of various mutants for the native and fibrillar forms indicated that, while the Trp residues introduced in the middle of the beta2-m sequence tend to be buried in the fibrils, those located in the C-terminal region are more exposed. In addition, the fluorescence spectra of fibrils prepared at pH 2.5 and 7.0 revealed a large difference in the fluorescence intensity for W60-beta2-m, implying a major structural difference between them.     

Accommodation of single amino acid insertions by the native state of staphylococcal nuclease

Single alanine and glycine insertions were introduced at 20 randomly selected positions in staphylococcal nuclease. The resulting changes in catalytic activity and in stability to guanidine hydrochloride denaturation indicate that the native state structure is frequently able to accommodate the extra residue without great difficulty, even insertions within secondary structural elements such as alpha helices and beta sheets. On average, an inserted residue reduces the free energy of denaturation (delta GH2O) by an amount roughly comparable to an alanine or glycine substitution for one of the residues flanking the site of insertion. Several positions outside of the enzyme active site were found where insertions, but not substitutions, lead to structural changes that modify catalytic activity and the circular dichroism spectrum. Amino acid insertions represent a virtually unexplored class of genetic mutation that may prove complementary to amino acid substitutions for engineering proteins with altered functional and structural properties.     

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.     

Conserved tryptophan mutation disrupts structure and function of immunoglobulin domain revealing unusual tyrosine fluorescence

Immunoglobulin (Ig) domains are the most prevalent protein domain structure and share a highly conserved folding pattern; however, this structural family of proteins is also the most diverse in terms of biological roles and tissue expression. Ig domains vary significantly in amino acid sequence but share a highly conserved tryptophan in the hydrophobic core of this beta‐stranded protein. Palladin is an actin binding and bundling protein that has five Ig domains and plays an important role in normal cell adhesion and motility. Mutation of the core tryptophan in one Ig domain of palladin has been identified in a pancreatic cancer cell line, suggesting a crucial role for this sole tryptophan in palladin Ig domain structure, stability, and function. We found that actin binding and bundling was not completely abolished with removal of this tryptophan despite a partially unfolded structure and significantly reduced stability of the mutant Ig domain as shown by circular dichroism investigations. In addition, this mutant palladin domain displays a tryptophan‐like fluorescence attributed to an anomalous tyrosine emission at 341 nm. Our results indicate that this emission originates from a tyrosinate that may be formed in the excited ground state by proton transfer to a nearby aspartic acid residue. Furthermore, this study emphasizes the importance of tryptophan in protein structural stability and illustrates how tyrosinate emission contributions may be overlooked during the interpretation of the fluorescence properties of proteins.     

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.     

How strong are side chain interactions in the folding intermediate?

Influence of 12 nonpolar amino acids residues from the hydrophobic core of apomyoglobin on stability of its native state and folding intermediate was studied. Six of the selected residues are from the A, G and H helices; these are conserved in structure of the globin family, although nonfunctional, that is, not involved in heme binding. The rest are nonconserved hydrophobic residues that belong to the B, C, D, and E helices. Each residue was substituted by alanine, and equilibrium pH‐induced transitions in apomyoglobin and its mutants were studied by circular dichroism and fluorescent spectroscopy. The obtained results allowed estimating changes in their free energy during formation of the intermediate state. It was first shown that the strength of side chain interactions in the apomyoglobin intermediate state amounts to 15–50% of that in its native state for conserved residues, and practically to 0% for nonconserved residues. These results allow a better understanding of interactions occurring in the intermediate state and shed light on involvement of certain residues in protein folding at different stages.     

Identification of cooperative folding units in a set of native proteins.

Cooperative unfolding penalties are calculated by statistically evaluating an ensemble of denatured states derived from native structures. The ensemble of denatured states is determined by dividing the native protein into short contiguous segments and defining all possible combinations of native, i.e., interacting, and non-native, i.e., non-interacting, segments. We use a novel knowledge-based scoring function, derived from a set of non-homologous proteins in the Protein Data Bank, to describe the interactions among residues. This procedure is used for the structural identification of cooperative folding cores for four globular proteins: bovine pancreatic trypsin inhibitor, horse heart cytochrome c, French bean plastocyanin, and staphylococcal nuclease. The theoretical folding units are shown to correspond to regions that exhibit enhanced stability against denaturation as determined from experimental hydrogen exchange protection factors. Using a sequence similarity score for related sequences, we show that, in addition to residues necessary for enzymatic function, those amino acids comprising structurally important folding cores are also preferentially conserved during evolution. This implies that the identified folding cores may be part of an array of fundamental structural folding units.     

Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.     

Topological frustration and the folding of interleukin-1 beta

The cytokine, interleukin-1beta (IL-1beta), adopts a beta-trefoil fold. It is known to be much slower folding than similarly sized proteins, despite having a low contact order. Proteins are sufficiently well designed that their folding is not dominated by local energetic traps. Therefore, protein models that encode only the folded structure and are energetically unfrustrated (Gō-type), can capture the essentials of the folding routes. We investigate the folding thermodynamics of IL-1beta using such a model and molecular dynamics (MD) simulations. We develop an enhanced sampling technique (a modified multicanonical method) to overcome the sampling problem caused by the slow folding. We find that IL-1beta has a broad and high free energy barrier. In addition, the protein fold causes intermediate unfolding and refolding of some native contacts within the protein along the folding trajectory. This "backtracking" occurs around the barrier region. Complex folds like the beta-trefoil fold and functional loops like the beta-bulge of IL-1beta can make some of the configuration space unavailable to the protein and cause topological frustration.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

The native energy landscape for interleukin-1beta. Modulation of the population ensemble through native-state topology

A minimalist Go-model, with no energetic frustration in the native conformation, has been shown to describe accurately the folding pathway of the beta-trefoil protein, interleukin-1beta (IL-1beta). While it appears that these models successfully model transition states and intermediates between the unfolded and native ensembles, it is unclear how accurately they capture smaller, yet biologically relevant, structural changes within the native ensemble after energetic perturbation. Here, we address the following questions. Can a simple Go-model of interleukin-1beta, based on native topology, describe changes in structural properties of the native ensemble as the protein stability is changed? Or is it necessary to include a more explicit representation of atoms, electrostatic, hydrogen bonding, and van der Waals forces to describe these changes? The native ensemble of IL-1beta was characterized using a variety of experimental probes under native (0 M NaCl, guanidine hydrochloride (Gdn-HCl)), moderately destabilized (0 M NaCl, 0.8 M Gdn-HCl), and in moderate salt concentration (0.8 M NaCl, 0 M Gdn-HCl). Heteronuclear (1)H-(15)N nuclear Overhauser effect spectroscopy (NOESY) and heteronuclear single quantum correlation (HSQC) NMR spectra confirmed that the beta-trefoil global fold was largely intact under these three conditions. However, 25 of the 153 residues throughout the chain did demonstrate (13)C and (1)H-(15)N chemical shifts when perturbed with 0.8 M NaCl or Gdn-HCl. Despite large differences in protection factors from solvent hydrogen-deuterium exchange for all residues between stable (0 M Gdn-HCl) and destabilized (0.8 M Gdn-HCl) IL-1beta, no difference in steady-state (15)N-(1)H NOE enhancements were measured. Thus, the chemical shifts correlate with a global but limited increase in residue flexibility in the presence of Gdn-HCl. Minimalist simulations highlight the regions of greatest position shift between native and 0.8 M Gdn-HCl, which were determined experimentally. This correlation demonstrates that structural changes within the native ensemble of IL-1beta are, at least partially, governed by the principle of minimal energetic frustration.     

The role of distal S6 hydrophobic residues in the voltage-dependent gating of CaV2.3 channels

The hydrophobic locus VAVIM is conserved in the S6 transmembrane segment of domain IV (IVS6) in Ca(V)1 and Ca(V)2 families. Herein we show that glycine substitution of the VAVIM motif in Ca(V)2.3 produced whole cell currents with inactivation kinetics that were either slower (A1719G approximately V1720G), similar (V1718G), or faster (I1721G approximately M1722G) than the wild-type channel. The fast kinetics of I1721G were observed with a approximately +10 mV shift in its voltage dependence of activation (E(0.5,act)). In contrast, the slow kinetics of A1719G and V1720G were accompanied by a significant shift of approximately -20 mV in their E(0.5,act) indicating that the relative stability of the channel closed state was decreased in these mutants. Glycine scan performed with Val (349) in IS6, Ile(701) in IIS6, and Leu(1420) in IIIS6 at positions predicted to face Val(1720) in IVS6 also produced slow inactivating currents with hyperpolarizing shifts in the activation and inactivation potentials, again pointing out a decrease in the stability of the channel closed state. Mutations to other hydrophobic residues at these positions nearly restored the channel gating. Altogether these data indicate that residues at positions equivalent to 1720 exert a critical control upon the relative stability of the channel closed and open states and more specifically, that hydrophobic residues at these positions promote the channel closed state. We discuss a three-dimensional homology model of Ca(V)2.3 based upon Kv1.2 where hydrophobic residues at positions facing Val(1720) in IS6, IIS6, and IIIS6 play a critical role in stabilizing the closed state in Ca(V)2.3.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

NMR structure of an archaeal homologue of ribonuclease P protein Rpp29

A protein component of the Archaeoglobus fulgidus RNase P was expressed in Escherichia coli, purified, and structurally characterized using multidimensional NMR methods. The dominant structural feature of this 11 kDa protein is a sheet of six antiparallel beta-strands, wrapped around a core of conserved hydrophobic amino acids. Amide proton exchange and (15)N relaxation rate data provide evidence that the first 16 residues of the protein, located before the start of the first beta-strand, and the last 24 residues, located past the end of the last beta-strand, are relatively flexible; this contrasts with the relatively rigid and well-defined structure of the beta-sheet. Amino acid sequence comparisons among a diverse set of species indicate that the A. fulgidus protein is homologous to the human RNase P protein Rpp29, yeast RNase P protein Pop4, and a known archaeal RNase P protein from Methanobacter thermoautotrophicus; conserved hydrophobic residues indicate that the homologous protein in each of these species contains a similar beta-sheet structure. Conserved surface residues located in the loop connecting strands beta2 and beta3, the loop connecting strands beta4 and beta5, and in the flexible N- and C-terminal tails are most likely to have specific interactions with the RNA and other proteins of RNase P. The structural model of an RNase P protein component provided by the present work provides an essential step toward eventually understanding the overall architecture of this complex enzyme and the mechanism by which it performs its functions.     

Packing and hydrophobicity effects on protein folding and stability: effects of beta-branched amino acids, valine and isoleucine, on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers.

The aim of this study was to examine the differences between hydrophobicity and packing effects in specifying the three-dimensional structure and stability of proteins when mutating hydrophobes in the hydrophobic core. In DNA-binding proteins (leucine zippers), Leu residues are conserved at positions "d," and beta-branched amino acids, Ile and Val, often occur at positions "a" in the hydrophobic core. In order to discern what effect this selective distribution of hydrophobes has on the formation and stability of two-stranded alpha-helical coiled coils/leucine zippers, three Val or three Ile residues were simultaneously substituted for Leu at either positions "a" (9, 16, and 23) or "d" (12, 19, and 26) in both chains of a model coiled coil. The stability of the resulting coiled coils was monitored by CD in the presence of Gdn.HCl. The results of the mutations of Ile to Val at either positions "a" or "d" in the reduced or oxidized coiled coils showed a significant hydrophobic effect with the additional methylene group in Ile stabilizing the coiled coil (delta delta G values range from 0.45 to 0.88 kcal/mol/mutation). The results of mutations of Leu to Ile or Val at positions "a" in the reduced or oxidized coiled coils showed a significant packing effect in stabilizing the coiled coil (delta delta G values range from 0.59 to 1.03 kcal/mol/mutation). Our results also indicate the subtle control hydrophobic packing can have not only on protein stability but on the conformation adopted by the amphipathic alpha-helices. These structural findings correlate with the observation that in DNA-binding proteins, the conserved Leu residues at positions "d" are generally less tolerant of amino acid substitutions than the hydrophobic residues at positions "a."