Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.     

Positional dependence of the effects of negatively charged Glu side chains on the stability of two-stranded α-helical coiled-coils

The effects on protein stability of negatively charged Glu side chains at different positions along the length of the α-helix were investigated in the two-stranded α-helical coiled-coil. A native coiled-coil has been designed which consists of two identical 35 residue polypeptide chains with a heptad repeat Q_gV_aG_bA_cL_dQ_eK_f and a Cys residue at position 2 to allow the formation of an interchain 2-2′ disulphide bridge. This coiled-coil contains no intra- or interchain electrostatic interactions and served as a control for peptides in which Glu was substituted for Gln in the e or g heptad positions. The effect of the substitutions on stability was determined by urea denaturation at 20°C with the degree of unfolding monitored by circular dichroism spectroscopy. A Glu substituted for Gln near the N-terminus in each chain of the coiled-coil stabilizes the coiled-coil at pH 7, consistent with the charge–helix dipole interaction model. This stability increase is modulated by pH change and the addition of salt (KCl or guanidine hydrochloride), confirming the electrostatic nature of the effect. In contrast, Glu substitution in the middle of the helix destabilizes the coiled-coil because of the lower helical propensity and hydrophobicity of Glu compared with Gln at pH 7. Taking the intrinsic differences into account, the apparent charge–helix dipole interaction at the N-terminus is approximately 0.35 kcal/mol per Glu substitution. A Glu substitution at the C-terminus destabilizes the coiled-coil more than in the middle owing to the combined effects of intrinsic destabilization and unfavourable charge–helix dipole interaction with the negative pole of the helix dipole. The estimated destabilizing charge–helix dipole interaction of 0.08 kcal/mol is smaller than the stabilizing interaction at the N-terminus. The presence of a 2-2′disulphide bridge appears to have little influence on the magnitude of the charge–helix dipole interactions at either end of the coiled-coil. © 1997 European Peptide Society and John Wiley & Sons, Ltd.     

Stabilizing and destabilizing clusters in the hydrophobic core of long two-stranded alpha-helical coiled-coils

Detailed sequence analyses of the hydrophobic core residues of two long two-stranded alpha-helical coiled-coils that differ dramatically in sequence, function, and length were performed (tropomyosin of 284 residues and the coiled-coil domain of the myosin rod of 1086 residues). Three types of regions were present in the hydrophobic core of both proteins: stabilizing clusters and destabilizing clusters, defined as three or more consecutive core residues of either stabilizing (Leu, Ile, Val, Met, Phe, and Tyr) or destabilizing (Gly, Ala, Cys, Ser, Thr, Asn, Gln, Asp, Glu, His, Arg, Lys, and Trp) residues, and intervening regions that consist of both stabilizing and destabilizing residues in the hydrophobic core but no clusters. Subsequently, we designed a series of two-stranded coiled-coils to determine what defines a destabilizing cluster and varied the length of the destabilizing cluster from 3 to 7 residues to determine the length effect of the destabilizing cluster on protein stability. The results showed a dramatic destabilization, caused by a single Leu to Ala substitution, on formation of a 3-residue destabilizing cluster (DeltaT(m) of 17-21 degrees C) regardless of the stability of the coiled-coil. Any further substitution of Leu to Ala that increased the size of the destabilizing cluster to 5 or 7 hydrophobic core residues in length had little effect on stability (DeltaT(m) of 1.4-2.8 degrees C). These results suggested that the contribution of Leu to protein stability is context-dependent on whether the hydrophobe is in a stabilizing cluster or its proximity to neighboring destabilizing and stabilizing clusters.     

Designing heterodimeric two-stranded alpha-helical coiled-coils. Effects of hydrophobicity and alpha-helical propensity on protein folding,stability, and specificity

The E/K coil, a heterodimeric coiled-coil, has been designed as a universal peptide capture and delivery system for use in applications such as biosensors and as an expression and affinity purification tag. In this design, heterodimer formation is specified through the placement of charged residues at the e and g positions of the heptad repeat such that the E coil contains all glutamic acid residues at these positions, and the K coil contains all lysine residues at these positions. The affinity and stability of the E/K coil have been modified to allow a greater range of conditions for association and dissociation. Increasing the hydrophobicity of the coiled-coil core, by substituting isoleucine for valine, gave increases in stability of 2.81 and 3.73 kcal/mol (0.47 kcal/mol/substitution). Increasing the alpha-helical propensity of residues outside the core, by substituting alanine for serine, yielded increases in stability of 2.68 and 3.28 kcal/mol (0.41 and 0.45 kcal/mol/substitution). These sequence changes yielded a series of heterodimeric coiled-coils whose stabilities varied from 6.8 to 11.2 kcal/mol, greatly expanding their scope for use in protein engineering and biomedical applications.     

The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins

We describe here a systematic investigation into the role of position a in the hydrophobic core of a model coiled-coil protein in determining coiled-coil stability and oligomerization state. We employed a model coiled coil that allowed the formation of an extended three-stranded trimeric oligomerization state for some of the analogs; however, due to the presence of a Cys-Gly-Gly linker, unfolding occurred from the same two-stranded monomeric oligomerization state for all of the analogs. Denaturation from a two-stranded state allowed us to measure the relative contribution of 20 different amino acid side chains to coiled-coil stability from chemical denaturation profiles. In addition, the relative hydrophobicity of the substituted amino acid side chains was assessed by reversed-phase high-performance liquid chromatography and found to correlate very highly (R = 0.95) with coiled-coil stability. We also determined the effect of position a in specifying the oligomerization state using ultracentrifugation as well as high-performance size-exclusion chromatography. We found that nine of the analogs populated one oligomerization state exclusively at peptide concentrations of 50 microM under benign buffer conditions. The Leu-, Tyr-, Gln-, and His-substituted analogs were found to be exclusively three-stranded trimers, while the Asn-, Lys-, Orn-, Arg-, and Trp-substituted analogs formed exclusively two-stranded monomers. Modeling results for the Leu-substituted analog showed that a three-stranded oligomerization state is preferred due to increased side-chain burial, while a two-stranded oligomerization state was observed for the Trp analog due to unfavorable cavity formation in the three-stranded state.     

Designing heterodimeric two-stranded alpha-helical coiled-coils: the effect of chain length on protein folding, stability and specificity.

The E/K coil, a heterodimeric coiled-coil, has been designed as a universal peptide capture and delivery system for use in applications such as biosensors and affinity chromatography. In this design, heterodimer formation is specified through the placement of charged residues at the e and g positions of the heptad repeat. The affinity and stability of the E/K coil has been modified in order to allow a greater range of conditions for association and dissociation by varying the chain length to obtain three, four and five heptad coiled-coils (21, 28 and 35 residues per polypeptide chain). The effect of chain length on stability and folding was examined by circular dichroism spectroscopy, guanidine hydrochloride denaturation, and redox equilibrium experiments. It was found that increases in chain length produced increases in the stability of heterodimeric coiled-coils, but in a nonlinear fashion. The resulting disulfide-bridged heterostranded molecules and reduced heterodimers span a wide range of stabilities (deltaG=3.3-11.9 kcal/mol), greatly expanding their scope for use in protein engineering and biomedical applications.     

Protein destabilization by electrostatic repulsions in the two-stranded alpha-helical coiled-coil/leucine zipper.

The destabilizing effect of electrostatic repulsions on protein stability has been studied by using synthetic two-stranded alpha-helical coiled-coils as a model system. The native coiled-coil consists of two identical 35-residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf and a Cys residue at position 2 to allow formation of an interchain disulfide bridge. This peptide, designed to contain no intrahelical or interhelical electrostatic interactions, forms a stable coiled-coil structure at 20 degrees C in benign medium (50 mM KCl, 25 mM PO4, pH 7) with a [urea]1/2 value of 6.1 M. Four mutant coiled-coils were designed to contain one or two Glu substitutions for Gln per polypeptide chain. The resulting coiled-coils contained potential i to i' + 5 Glu-Glu interchain repulsions (denoted as peptide E2(15,20)), i to i' + 2 Glu-Glu interchain repulsions (denoted E2(20,22)), or no interchain ionic interactions (denoted E2(13,22) and E1(20)). The stabilities of the coiled-coils were determined by measuring the ellipticities at 222 nm as a function of urea or guanidine hydrochloride concentration at 20 degrees C in the presence and absence of an interchain disulfide bridge. At pH 7, in the presence of urea, the stabilities of E2(13,22) and E2(20,22) were identical suggesting that the potential i to i' + 2 interchain Glu-Glu repulsion in the E2(20,22) coiled-coil does not occur. In contrast, the mutant E2(15,20) is substantially less stable than E2(13,22) or E2(15,20) by 0.9 kcal/mol due to the presence of two i to i' + 5 interchain Glu-Glu repulsions, which destabilize the coiled-coil by 0.45 kcal/mol each. At pH 3 the coiled-coils were found to increase in stability as the number of Glu substitutions were increased. This, combined with reversed-phase HPLC results at pH 7 and pH 2, supports the conclusion that the protonated Glu side chains present at low pH are significantly more hydrophobic than Gln side chains which are in turn more hydrophobic than the ionized Glu side chains present at neutral pH. The protonated Glu residues increase the hydrophobicity of the coiled-coil interface leading to higher coiled-coil stability. The guanidine hydrochloride results at pH 7 show similar stabilities between the native and mutant coiled-coils indicating that guanidine hydrochloride masks electrostatic repulsions due to its ionic nature and that Glu and Gln in the e and g positions of the heptad repeat have very similar effects on coiled-coil stability in the presence of GdnHCl.     

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."     

Estimation of the maximum change in stability of globular proteins upon mutation of a hydrophobic residue to another of smaller size.

Although the hydrophobic effect is generally considered to be one of the most important forces in stabilizing the folded structure of a globular protein molecule, there is a lack of consensus on the precise magnitude of this effect. The magnitude of the hydrophobic effect is most directly measured by observing the change in stability of a protein molecule when an internal hydrophobic residue is mutated to another of smaller size. Results of such measurements have, however, been confusing because they vary greatly and are generally considerably larger than expected from the transfer free energies of corresponding small molecules. In this article, a thermodynamic argument is presented to show (1) that the variation is mainly due to that in the flexibility of the protein molecule at the site of mutation, (2) that the maximum destabilization occurs when the protein at the site of mutation is rigid, in which case the value of the destabilization is approximately given by the work of cavity formation in water, and (3) that the transfer free energy approximately gives the minimum of the range of variations. The best numerical agreements between the small molecule and the protein systems are obtained when the data from the small molecule system are expressed as the molarity-based standard free energies without other corrections.     

The optimal fraction of hydrophobic residues required to ensure protein collapse

The hydrophobic interaction is the main driving force for protein folding. Here, we address the question of what is the optimal fraction, f of hydrophobic (H) residues required to ensure protein collapse. For very small f (say f<0.1), the protein chain is expected to behave as a random coil, where the H residues are "wrapped" locally by polar (P) residues. However, for large enough f this local coverage cannot be achieved and the thermodynamic alternative to avoid contact with water is burying the H residues in the interior of a compact chain structure. The interior also contains P residues that are known to be clustered to optimize their electrostatic interactions. This means that the H residues are clustered as well, i.e. they effectively attract each other like the H-monomers in Dill's HP lattice model. Previously, we asked the question: assuming that the H monomers in the HP model are distributed randomly along the chain, what fraction of them is required to ensure a compact ground state? We claimed there that f approximately p(c), where p(c) is the site percolation threshold of the lattice (in a percolation experiment, each site of an initially empty lattice is visited and a particle is placed there with a probability p. The interest is in the critical (minimal) value, p(c), for which percolation occurs, i.e. a cluster connecting the opposite sides of the lattice is created). Due to the above correspondence between the HP model and real proteins (and assuming that the H residues are distributed at random) we suggest that the experimental f should lead to percolating clusters of H residues over the highly dense protein core, i.e. clusters of the core size. To check this theory, we treat a simplified model consisting of H and P residues represented by their alpha-carbon atoms only. The structure is defined by the C(alpha)-C(alpha) virtual bond lengths, angles and dihedral angles, and the X-ray structure is best-fitted onto a face-centered cubic lattice. Percolation experiments are carried out for 103 single-chain proteins using six different hydrophobic sets of residues. Indeed, on average, percolating clusters are generated, which supports our theory; however, some sets lead to a better core coverage than others. We also calculate the largest actual hydrophobic cluster of each protein and show that, on average, these clusters span the core, again in accord with our theory. We discuss the effect of protein size, deviations from the average picture, and implications of this study for defining reliable simplified models of proteins.     

Automatic recognition of hydrophobic clusters and their correlation with protein folding units.

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.     

Clusters of isoleucine,leucine, and valine side chains define cores of stability in high‐energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high‐energy states that populate their folding free‐energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high‐energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high‐energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.     

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.     

Osmolyte effects on helix formation in peptides and the stability of coiled-coils

The ability of several naturally occurring substances known as osmolytes to induce helix formation in an alanine-based peptide have been investigated. As predicted by the osmophobic effect hypothesis, the osmolytes studies here do induce helix formation. Trimethylamine-N-oxide (TMAO) is the best structure-inducing osmolytes investigated here, but it is not as effective in promoting helix formation as the common cosolvent trifluoroethanol (TFE). We also provide a semiquantitative study of the ability of TMAO to induce helix formation and urea, which acts as a helix (and protein) denaturant. We find that on a molar basis, these agents are exactly counteractive as structure inducing and unfolding agents. Finally, we extend the investigations to the effects of urea and TMAO on the stability of a dimeric coiled-coil peptide and find identical results. Together these results support the tenets of the osmophobic hypothesis and highlight the importance of the polypeptide backbone in protein folding and stability.     

Defining process design space for a hydrophobic interaction chromatography (HIC) purification step: Application of quality by design (QbD) principles

The concept of design space has been taking root under the quality by design paradigm as a foundation of in‐process control strategies for biopharmaceutical manufacturing processes. This paper outlines the development of a design space for a hydrophobic interaction chromatography (HIC) process step. The design space included the impact of raw material lot‐to‐lot variability and variations in the feed stream from cell culture. A failure modes and effects analysis was employed as the basis for the process characterization exercise. During mapping of the process design space, the multi‐dimensional combination of operational variables were studied to quantify the impact on process performance in terms of yield and product quality. Variability in resin hydrophobicity was found to have a significant influence on step yield and high‐molecular weight aggregate clearance through the HIC step. A robust operating window was identified for this process step that enabled a higher step yield while ensuring acceptable product quality. Biotechnol. Bioeng. 2010;107: 985–997. © 2010 Wiley Periodicals, Inc.     

A desolvation barrier to hydrophobic cluster formation may contribute to the rate-limiting step in protein folding.

To gain insight into the free energy changes accompanying protein hydrophobic core formation, we have used computer simulations to study the formation of small clusters of nonpolar solutes in water. A barrier to association is observed at the largest solute separation that does not allow substantial solvent penetration. The barrier reflects an effective increase in the size of the cavity occupied by the expanded but water-excluding cluster relative to both the close-packed cluster and the fully solvated separated solutes; a similar effect may contribute to the barrier to protein folding/unfolding. Importantly for the simulation of protein folding without explicit solvent, we find that the interactions between nonpolar solutes of varying size and number can be approximated by a linear function of the molecular surface, but not the solvent-accessible surface of the solutes. Comparison of the free energy of cluster formation to that of dimer formation suggests that the assumption of pair additivity implicit in current protein database derived potentials may be in error.     

The role of unstructured highly charged regions on the stability and specificity of dimerization of two-stranded alpha-helical coiled-coils: analysis of the neck-hinge region of the kinesin-like motor protein Kif3A

We investigated the folding, stability, and specificity of dimerization of the neck-hinge region (residues 356-416) of the kinesin-like protein Kif3A. We showed that the predicted coiled-coil on its own (residues 356-377) will fold autonomously in solution. We then explored the ability of oppositely charged regions to specify heterodimer formation in coiled-coils by synthesizing analogs of the neck coiled-coil region with and without various negatively and positively charged extensions to the C-terminus of the neck coiled-coil and characterizing these analogs by circular dichroism spectroscopy. The charged region alone (residues 378-416) adopted a random-coil structure and this region remained unfolded in the presence of the coiled-coil. Redox experiments demonstrated that oppositely charged regions specified the formation of a hetero-two-stranded coiled-coil. Denaturation studies with urea demonstrated a decrease in coiled-coil stability with the addition of negatively charged residues in the homostranded coiled-coil; conversely, the addition of the positively charged region (residues 403-416) of Kif3A C-terminally to the neck coiled-coil did not affect coiled-coil stability. Overall, our results suggest that electrostatic attractions drive the specificity of heterodimerization of the coiled-coil, not the removal of positive or negative charge-charge repulsions, while maintaining the stability of the heterodimer compared to that of the stablest homodimer.     

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.     

Chemically crosslinked protein dimers: stability and denaturation effects.

Nine single substitution cysteine mutants of staphylococcal nuclease (nuclease) were preferentially crosslinked at the introduced cysteine residues using three different bifunctional crosslinking reagents; 1,6-bismaleimidohexane (BMH), 1,3-dibromo-2-propanol (DBP), and the chemical warfare agent, mustard gas (bis(2-chloroethyl)sulfide; mustard). BMH and mustard gas are highly specific reagents for cysteine residues, whereas DBP is not as specific. Guanidine hydrochloride (GuHCl) denaturations of the resulting dimeric proteins exhibited biphasic unfolding behavior that did not fit the two-state model of unfolding. The monofunctional reagent, epsilon-maleimidocaproic acid (MCA), was used as a control for the effects of alkylation. Proteins modified with MCA unfolded normally, showing that this unusual unfolding behavior is due to crosslinking. The data obtained from these crosslinked dimers was fitted to a three-state thermodynamic model of two successive transitions in which the individual subunits cooperatively unfold. These two unfolding transitions were very different from the unfolding of the monomeric protein. These differences in unfolding behavior can be attributed in large part to changes in the denatured state. In addition to GuHCl titrations, the crosslinked dimers were also thermally unfolded. In contrast to the GuHCl denaturations, analysis of this data fit a two-state model well, but with greatly elevated van't Hoff enthalpies in many cases. However, clear correlations between the thermal and GuHCl denaturations exist, and the differences in thermal unfolding can be rationalized by postulating interactions of the denatured crosslinked proteins.     

Computational design of proteins stereochemically optimized in size, stability, and folding speed

Artificial proteins potentially barrier-free in the folding kinetics are approached computationally under the guidance of protein-folding theories. The smallest and fastest folding globular protein triple-helix-bundle (THB) is so modified as to minimize or eliminate its presumed barriers in folding speed. As the barriers may reside in the ordering of either secondary or tertiary structure, the elements of both secondary and tertiary structure in the protein are targeted for prenucleation with suitable stereochemically constrained amino acid residues. The required elements of topology and sequence for the THB are optimized independently; first the topology is optimized with simulated annealing in polypeptides of highly simplified alphabet; next, the sequence in side chains is optimized using the standard inverse design methods. The resultant three best-adapted THBs, variable in topology and distinctive in sequences, are assessed by comparing them with a few benchmark proteins. The results of mainly molecular dynamics (MD) comparisons, undertaken in explicit water at different temperatures, show that the designed sequences are favorably placed against the chosen benchmarks as THB proteins potentially thermostable in the native folds. Folding simulation experiments with MD establish that the designed sequences are rapid in the folding of individual helices, but not in the evolution of tertiary structure; energetic cum topological frustrations remain but could be the artifacts of the starting conformations that were chosen in the THBs in the folding simulations. Overall, a practical high-throughput approach for de novo protein design has been developed that may have fruitful application for any type of tertiary structure.